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.
Rock mass defect controlled deep-seated landslides are widespread within the deeply incised landscapes formed in Tertiary soft rock terrain in New Zealand. The basal failure surfaces of deep-seated slope failures are defined by thin, comparatively weak and laterally continuous bedding parallel layers termed critical stratigraphic horizons. These horizons have a sedimentary origin and have typically experienced some prior tectonically induced shear displacement at the time of slope failure. The key controls on the occurrence and form of deep-seated landslides are considered in terms of rock mass defect properties and tectonic and climatic forcing. The selection of two representative catchments (in southern Hawke's Bay and North Canterbury) affected by tectonic and climatic forcing has shown that the spatial and temporal initiation of deep-seated bedrock landslides in New Zealand Tertiary soft rock terrain is a predictable rather than a stochastic process; and that deep-seated landslides as a mass wasting process have a controlling role in landscape evolution in many catchments formed in Tertiary soft rock terrain. The Ella Landslide in North Canterbury is a deep-seated (~85 m) translational block slide that has failed on a 5 - 10 mm thick, kaolinite-rich, pre-sheared critical stratigraphic horizon. The residual strength of this sedimentary horizon, (C'R 2.6 - 2.7 kPa, and Ѳ'R = 16 - 21°), compared to the peak strength of the dominant lithology (C' = 176 kPa, and Ѳ' = 37°) defines a high strength contrast in the succession, and therefore a critical location for the basal failure surface of deep-seated slope failures. The (early to mid Holocene) Ella Landslide debris formed a large landslide dam in the Kate Stream catchment and this has significantly retarded rates of mass wasting in the middle catchment. Numerical stability analysis shows that this slope failure would have most likely required the influence of earthquake induced strong ground motion and the event is tentatively correlated to a Holocene event on the Omihi Fault. The influence of this slope failure is likely to affect the geomorphic development of the catchment on a scale of 10⁴ - 10⁵ years. In deeply incised catchments at the southeastern margin of the Maraetotara Plateau, southern Hawke's Bay, numerous widespread deep-seated landslides have basal failure surfaces defined by critical stratigraphic horizons in the form of thin « 20 mm) tuffaceous beds in the Makara Formation flysch (alternating sandstone and mudstone units). The geometry of deep-seated slope failures is controlled by these regularly spaced (~70 m), very weak critical stratigraphic horizons (C'R 3.8 - 14.2 kPa, and Ѳ'R = 2 - 5°), and regularly spaced (~45 m) and steeply dipping (-50°) critical conjugate joint/fault sets, which act as slide block release surfaces. Numerical stability analysis and historical precedent show that the temporal initiation of deep-seated landslides is directly controlled by short term tectonic forcing in the form of periodic large magnitude earthquakes. Published seismic hazard data shows the recurrence interval of earthquakes producing strong ground motions of 0.35g at the study site is every 150 yrs, however, if subduction thrust events are considered the level of strong ground motion may be much higher. Multiple occurrences of deep-seated slope failure are correlated to failure on the same critical stratigraphic horizon, in some cases in three adjacent catchments. Failure on multiple critical stratigraphic horizons leads to the development of a "stepped" landscape morphology. This slope form will be maintained during successive accelerated stream incision events (controlled by long term tectonic and climatic forcing) for as long as catchments are developing in this specific succession. Rock mass defect controlled deep seated landslides are controlling catchment head progression, landscape evolution and hillslope morphology in the Hawke's Bay study area and this has significant implications for the development of numerical landscape evolution models of landscapes formed in similar strata. Whereas the only known numerical model to consider deep seated landslides as an erosion process (ZSCAPE) considers them as stochastic in time and space, this study shows that this could not be applied to a landscape where the widespread spatial occurrence of deep-seated landslides is controlled by rock mass defects. In both of the study areas for this project, and by implication in many catchments in Tertiary soft rock terrain, deep-seated landslides controlled by rock mass defect strength, spacing and orientation, and tectonic and climatic forcing have an underlying control on landscape evolution. This study quantifies parameters for the development of numerical landscape evolution models that would assess the role of specific parameters, such as uplift rates, incision rates and earthquake recurrence in catchment evolution in Tertiary soft rock terrain.
The Canterbury Region is susceptible to a variety of natural hazards, including earthquakes, landslides and climate hazards. Increasing population and tourism within the region is driving development pressures and as more and more development occurs, the risk from natural hazards increases. In order to avoid development occurring in unacceptably vulnerable locations, natural hazard assessments are required. This study is a reconnaissance natural hazard assessment of Lakes Lyndon, Coleridge and Tekapo. There is restricted potential for development at Lake Lyndon, because the land surrounding the lake is owned by the Crown and has a number of development restrictions. However, there is the potential for conservation or recreation-linked development to occur. There is more potential for development at Lake Coleridge. Most of the land surrounding the lake is privately owned and has less development restrictions. The majority of land surrounding Lake Tekapo is divided into Crown-owned pastoral leases, which are protected from development, such as subdivision. However, there are substantial areas around the lake, which are privately owned and, therefore, have potential for development. Earthquake, landslide and climate hazards are the main natural hazards threatening Lakes Lyndon, Coleridge and Tekapo. The lakes are situated in a zone of active earth deformation in which large and relatively frequent earthquakes are produced. A large number of active faults lie within 15 km of each lake, which are capable of producing M7 or larger earthquakes. Ground shaking, liquefaction, landslides, tsunami and seiches are among the consequences of earthquakes, all of which have the potential to cause severe damage to lives, lifelines and infrastructure. Landslides are also common in the landscape surrounding the lakes. The majority of slopes surrounding the lakes are at significant risk from earthquake-induced failure under moderate to strong earthquake shaking. This level of shaking is expected to occur in any 50 year period around Lakes Lyndon and Coleridge, and in any 150 year period around Lake Tekapo. Injuries, fatalities and property damage can occur directly from landslide impact or from indirect effects such as flooding from landslide-generated tsunami or from landslide dam outbreaks. Lakes Lyndon, Coleridge and Tekapo are also susceptible to climate hazards, such as high winds, drought, heavy snowfall and heavy rainfall, which can lead to landslides and flooding. Future climate change due to global warming is most likely going to affect patterns of frequency and magnitudes of extreme weather events, leading to an increase in climate hazards. Before development is permitted around the lakes, it is essential that each of these hazards is considered so that unacceptably vulnerable areas can be avoided.
Case study analysis of the 2010-2011 Canterbury Earthquake Sequence (CES), which particularly impacted Christchurch City, New Zealand, has highlighted the value of practical, standardised and coordinated post-earthquake geotechnical response guidelines for earthquake-induced landslides in urban areas. The 22nd February 2011 earthquake, the second largest magnitude event in the CES, initiated a series of rockfall, cliff collapse and loess failures around the Port Hills which severely impacted the south-eastern part of Christchurch. The extensive slope failure induced by the 22nd February 200 earthquake was unprecedented; and ground motions experienced significantly exceeded the probabilistic seismic hazard model for Canterbury. Earthquake-induced landslides initiated by the 22nd February 2011 earthquake posed risk to life safety, and caused widespread damage to dwellings and critical infrastructure. In the immediate aftermath of the 22nd February 2011 earthquake, the geotechnical community responded by deploying into the Port Hills to conduct assessment of slope failure hazards and life safety risk. Coordination within the voluntary geotechnical response group evolved rapidly within the first week post-earthquake. The lack of pre-event planning to guide coordinated geotechnical response hindered the execution of timely and transparent management of life safety risk from coseismic landslides in the initial week after the earthquake. Semi-structured interviews were conducted with municipal, management and operational organisations involved in the geotechnical response during the CES. Analysis of interview dialogue highlighted the temporal evolution of priorities and tasks during emergency response to coseismic slope failure, which was further developed into a phased conceptual model to inform future geotechnical response. Review of geotechnical responses to selected historical earthquakes (Northridge, 1994; Chi-Chi, 1999; Wenchuan, 2008) has enabled comparison between international practice and local response strategies, and has emphasised the value of pre-earthquake preparation, indicating the importance of integration of geotechnical response within national emergency management plans. Furthermore, analysis of the CES and international earthquakes has informed pragmatic recommendations for future response to coseismic slope failure. Recommendations for future response to earthquake-induced landslides presented in this thesis include: the integration of post-earthquake geotechnical response with national Civil Defence and Emergency Management; pre-earthquake development of an adaptive management structure and standard slope assessment format for geotechnical response; and emergency management training for geotechnical professionals. Post-earthquake response recommendations include the development of geographic sectors within the area impacted by coseismic slope failure, and the development of a GIS database for analysis and management of data collected during ground reconnaissance. Recommendations provided in this thesis aim to inform development of national guidelines for geotechnical response to earthquake-induced landslides in New Zealand, and prompt debate concerning international best practice.
The 22 February 2011, Mw6.2 Christchurch earthquake is the most costly earthquake to affect New Zealand, causing an estimated 181 fatalities and severely damaging thousands of residential and commercial buildings. This paper presents a summary of some of the observations made by the NSF-sponsored GEER Team regarding the geotechnical/geologic aspects of this earthquake. The Team focused on documenting the occurrence and severity of liquefaction and lateral spreading, performance of building and bridge foundations, buried pipelines and levees, and significant rockfalls and landslides. Liquefaction was pervasive and caused extensive damage to residential properties, water and wastewater networks, high-rise buildings, and bridges. Entire neighborhoods subsided, resulting in flooding that caused further damage. Additionally, liquefaction and lateral spreading resulted in damage to bridges and to stretches of levees along the Waimakariri and Kaiapoi Rivers. Rockfalls and landslides in the Port Hills damaged several homes and caused several fatalities.
One of the less understood geotechnical responses to the cyclic loading from the MW6.2 Christchurch Earthquake, on the 22nd of February 2011, is the fissuring in the loessial soil-mantled, footslope positions of the north-facing valleys of the Port Hills. The fissures are characterized by mostly horizontal offset (≤500mm), with minor vertical displacement (≤300mm), and they extend along both sides of valleys for several hundred metres in an approximately contour-parallel orientation. The fissure traces correspond to extensional features mapped in other studies. Previous studies have suggested that the fissures are the headscarps of incipient landslides, but the surface and subsurface features are not typical of landslide movement. Whilst there are some features that correlate with landslide movement, there are many features that contradict the landslide movement hypothesis. Of critical importance to this investigation was the fact that there are no landslide flanks, there has been no basal shear surface found, there is little deformation in the so-called ‘landslide body’, and there have been no recorded zones of low shear strength in the soil deposit that are indicative of a basal shear surface. This thesis is a detailed geotechnical study on the fissures along part of Ramahana Road in the Hillsborough Valley, Christchurch. Shallow and deep investigation methods found that the predominant soil is loess-colluvium, to depths of ~20m, and this soil has variable geotechnical characteristics depending on the layer sampled. The factor that has the most influence on shear strength was found to be the moisture content. Direct shear-box testing of disturbed, recompacted loess-colluvium found that the soil had a cohesion of 35-65kPa and a friction angle of 38-43° when the soil moisture content was at 8-10%. However when the moisture content was at 19-20% the soil’s cohesion decreased to 3-5kPa and its friction angle decreased to 33-38°, this moisture content is at or slightly above the plastic limit. An electrical resistivity geophysical survey was conducted perpendicular to multiple fissure traces and through the compressional zone at 17 Ramahana Road. The electrical resistivity line found that there was an area of high resistivity at the toe of the slope, and an area of high conductivity downslope of this and at greater depths. This area correlated to the compressional zone recorded by previous studies. Moisture content testing of the soil in these locations showed that the soil in the resistive area was relatively dry (9%) compared to the surrounding soil (13%), whilst the soil in the conductive area was relatively wet (22%)compared to the surrounding soil (19%). Density tests of the soil in the compressional zone recorded that the resistive area had a higher dry density than the surrounding soil (~1790 kg/m3 compared to ~1650 kg/m3). New springs arose downslope of the compressional zone contemporaneously with the fissures, and it is interpreted that these have arisen from increased hydraulic head in the Banks Peninsula bedrock aquifer system, and earthquake induced-bedrock fracturing. A test pit was dug across an infilled fissure trace at 17 Ramahana Road to a depth of 3m. The fissure trace had an aperture of 450-470mm at the ground surface, but it gradually lost aperture with depth until 2.0-2.1m where it became a segmented fissure trace with 1-2mm aperture. A mixed-colluvium layer was intercepted by the fissure trace at 2.4m depth, and there was no observable vertical offset of this layer. The fissure trace was at an angle of 78° at the ground surface, but it also flattened with depth, which gave it a slightly curved appearance. The fissure trace was at an assumed angle of 40-50° near the base of the test pit. Rotational slide, translational slide and lateral spread landslide movement types were compared and contrasted as possibilities for landslide movement types, whilst an alternative hypothesis was offered that the fissures are tensile failures with a quasi-toppling motion involving a cohesive block of loessial soil moving outwards from the slope, with an accommodating compressional strain in the lower less cohesive soil. The mechanisms behind this movement are suggested to be the horizontal earthquake inertia forces from the Christchurch Earthquake, the static shear stress of the slope, and bedrock uplift of the Port Hills in relation to the subsidence of the Christchurch city flatlands. Extremely high PGA is considered to be a prerequisite to the fissure trace development, and these can only be induced in the Hillsborough Valley from a Port Hills Fault rupture, which has a recurrence interval of ~10,000 years. The current understanding of how the loess-colluvium soil would behave under cyclic loading is limited, and the mechanisms behind the suggested movement type are not completely understood. Further research is needed to confirm the proposed mechanism of the fissure traces. Laboratory tests such as the cyclic triaxial and cyclic shear test would be beneficial in future research to quantitatively test how the soil behaves under cyclic loading at various moisture contents and clay contents, and centrifuge experiments would be of great use to qualitatively test the suggested mode of movement in the loessial soil.
1. INTRODUCTION. Earthquakes and geohazards, such as liquefaction, landslides and rock falls, constitute a major risk for New Zealand communities and can have devastating impacts as the Canterbury 2010/2011 experience shows. Development patterns expose communities to an array of natural hazards, including tsunamis, floods, droughts, and sea level rise amongst others. Fostering community resilience is therefore vitally important. As the rhetoric of resilience is mainstreamed into the statutory framework, a major challenge emerges: how can New Zealand operationalize this complex and sometimes contested concept and build ‘community capitals’? This research seeks to provide insights to this question by critically evaluating how community capitals are conceptualized and how they can contribute to community resilience in the context of the Waimakariri District earthquake recovery and regeneration process.
Following the 22 February 2011, MW 6.2 earthquake located on a fault beneath the Port Hills of Christchurch, fissuring of up to several hundred metres in length was observed in the loess and loess-colluvium of foot-slope positions in north-facing valleys of the Port Hills. The fissuring was observed in all major valleys, occurred at similar low altitudes, showing a contour-parallel orientation and often accompanied by both lateral compression/extension features and spring formation in the valley floor below. Fissuring locations studied in depth included Bowenvale Valley, Hillsborough Valley, Huntlywood Terrace–Lucas Lane, Bridle Path Road, and Maffeys Road–La Costa Lane. Investigations into loess soil, its properties and mannerisms, as well as international examples of its failure were undertaken, including study of the Loess Plateau of China, the Teton Dam, and palaeo-fissuring on Banks Peninsula. These investigations lead to the conclusion that loess has the propensity to fail, often due to the infiltration of water, the presence of which can lead to its instantaneous disaggregation. Literature study and laboratory analysis of Port Hills loess concluded that is has the ability to be stable in steep, sub-vertical escarpments, and often has a sub-vertically jointed internal structure and has a peak shear strength when dry. Values for cohesion, c (kPa) and the internal friction angle, ϕ (degrees) of Port Hills loess were established. The c values for the 40 Rapaki Road, 3 Glenview Terrace loess samples were 13.4 kPa and 19.7 kPa, respectively. The corresponding ϕ values were thought unusually high, at 42.0° and 43.4°.The analysed loess behaved very plastically, with little or no peak strength visible in the plots as the test went almost directly to residual strength. A geophysics resistivity survey showed an area of low resistivity which likely corresponds to a zone of saturated clayey loess/loess colluvium, indicating a high water table in the area. This is consistent with the appearances of local springs which are located towards the northern end of each distinct section of fissure trace and chemical analysis shows that they are sourced from the Port Hills volcanics. Port Hills fissuring may be sub-divided into three categories, Category A, Category B, and Category C, each characterised by distinctive features of the fissures. Category A includes fissures which display evidence of, spring formation, tunnel-gullying, and lateral spreading-like behaviour or quasi-toppling. These fissures are several metres down-slope of the loess-bedrock interface, and are in valleys containing a loess-colluvium fill. Category B fissures are in wider valleys than those in Category A, and the valleys contain estuarine silty sediments which liquefied during the earthquake. Category C fissures occurred at higher elevations than the fissures in the preceding categories, being almost coincident with bedrock outcropping. It is believed that the mechanism responsible for causing the fissuring is a complex combination of three mechanisms: the trampoline effect, bedrock fracturing, and lateral spreading. These three mechanisms can be applied in varying degrees to each of the fissuring sites in categories A, B, and C, in order to provide explanation for the observations made at each. Toppling failure can describe the soil movement as a consequence of the a three causative mechanisms, and provides insight into the movement of the loess. Intra-loess water coursing and tunnel gullying is thought to have encouraged and exacerbated the fissuring, while not being the driving force per se. Incipient landsliding is considered to be the least likely of the possible fissuring interpretations.
The M7.8 Kaikoura Earthquake in 2016 presented a number of challenges to science agencies and institutions throughout New Zealand. The earthquake was complex, with 21 faults rupturing throughout the North Canterbury and Marlborough landscape, generating a localised seven metre tsunami and triggering thousands of landslides. With many areas isolated as a result, it presented science teams with logistical challenges as well as the need to coordinate efforts across institutional and disciplinary boundaries. Many research disciplines, from engineering and geophysics to social science, were heavily involved in the response. Coordinating these disciplines and institutions required significant effort to assist New Zealand during its most complex earthquake yet recorded. This paper explores that effort and acknowledges the successes and lessons learned by the teams involved.
On 14 November 2016, the Mw 7.8 Kaikōura earthquake caused widespread damage along the east coast of the South Island, New Zealand. Kaikōura town itself was isolated from the rest of the country by landslides blocking off major roads. While impacts from the Kaikōura earthquake on large, urban population centres have been generally well documented, this thesis aims to fill gaps in academic knowledge regarding small rural towns. This thesis investigates what, where and when critical infrastructure and lifeline service disruption occurred following the 2016 Kaikōura earthquake in a selection of small towns, and how the communities in these areas adapted to disruption. Following a robust review of literature and news media, four small rural towns were selected from North Canterbury (Culverden & Waiau) and Marlborough (Seddon & Ward) in the South Island, New Zealand. Semi-structured interview sessions with a special focus on these towns were held with infrastructure managers, emergency response and recovery officials, and organisation leaders with experience or expertise in the 2016 Kaikōura earthquake. Findings were supplemented with emergency management situation reports to produce hazard maps and infrastructure exposure maps. A more detailed analysis was conducted for Waiau involving interdependence analyses and a level of service timeline for select lifeline services. The earthquake impacted roads by blocking them with landslides, debris and surface rupture. Bridges where shaken off their abutments, breaking infrastructure links such as fibre landlines as they went. Water supplies and other forms of infrastructure relied heavily on the level of service of roads, as rough rural terrain left few alternatives. Adapting to an artificial loss of road service, some Waiau locals created their own detour around a road cordon in order to get home to family and farms. Performance of dwellings was tied to socioeconomic factors as much as proximity to the epicentre. Farmers who lost water access pulled out fences to allow stock to drink from rivers. Socioeconomic differences between farmland and township residents also contributed to resilience variations between the towns assessed in this study. Understanding how small rural towns respond and adapt to disaster allows emergency management officials and policy to be well informed and flexible with planning for multiple size classes of towns.
On 14 November 2016 a magnitude Mw 7.8 earthquake struck the upper South Island of New Zealand with effects also being observed in the capital city, Wellington. The affected area has low population density but is the largest wine production region in New Zealand and also hosts the main national highway and railway routes connecting the country’s three largest cities of Auckland, Wellington and Christchurch, with Marlborough Port in Picton providing connection between the South and North Islands. These transport facilities sustained substantial earthquake related damage, causing major disruptions. Thousands of landslides and multiple new faults were counted in the area. The winery facilities and a large number of commercial buildings and building components (including brick masonry veneers, historic masonry construction, and chimneys), sustained damage due to the strong vertical and horizontal acceleration. Presented herein are field observations undertaken the day immediately after the earthquake, with the aim to document earthquake damage and assess access to the affected area.
<strong>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.</strong>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.
This project was initiated by ENGEO Limited and KiwiRail Holdings Limited to assess the stability of Slovens Creek Viaduct (specifically its western abutment) and a 3km section of rail corridor between Slovens Creek Viaduct and Avoca on the Midland Line (MDL). Commonly known as the scenic TranzAlpine rail journey (through Arthurs Pass National Park) the MDL connects Greymouth to Christchurch via Rolleston, where the MDL meets the Main South Line into Christchurch. The project area is approximately 40km southeast of Arthurs Pass Township, in the eastern extension of the Castle Hill Basin which is part of the Waimakariri Catchment and Canterbury Foothills. The field area is underlain by Rakaia Terrane, which is part of the Torlesse Composite Terrane forming the basement rock unit for the field area. Cretaceous-Tertiary rocks of the Castle Hill Basin overlie the basement strata and record a transgression-regression sequence, as well as mid-Oligocene submarine volcanism. The stratigraphic sequence in the Castle Hill Basin, and its eastern extension to Avoca, comprises two formations of the Eyre group, the older Broken River Formation and the younger Iron Creek Formation. Deep marine Porter Group limestones, marls, and tuffs of Oligocene age succeed the Iron Creek Formation of the Eyre Group, and probably records the maximum of the transgression. The Enys Formation lies disconformably on the Porter Group and is overlain unconformably by Late Pleistocene glacifluvial and glacial deposits. The Tertiary strata in the Slovens-Avoca rail corridor are weak, and the clay-rich tuff derived from mid-Oligocene volcanism is particularly prone to slaking. Extensive mapping carried out for this project has identified that some 90 percent of the surface along the length of the Slovens-Avoca corridor has been subject to mass movement. The landslides of the Slovens-Avoca rail corridor are clearly younger than the Last Glaciation, and Slovens Creek has been downcutting, with associated faulting and uplift, to form the present day geomorphology of the rail corridor. Deep-seated landslides in the rail corridor extend to Slovens Creek, locally deflecting the stream course, and a generic ground failure model for the rail corridor has been developed. Exploratory geotechnical investigations, including core drilling, installation of an inclinometer and a piezometer, enabled the construction of a simple ground model and cross section for the Slovens Creek Viaduct western abutment. Limit-equilibrium and pseudo-static slope stability analyses using both circular and block critical slip surface search methods were applied to the ground model for the western abutment of Slovens Creek Viaduct. Piezometric and strength data obtained during laboratory testing of core material have been used to constrain the western abutment stability assessment for one representative section line (C-C’). Prior to pseudo-static sensitivity analyses peak ground acceleration (PGA) for various Ultimate Limit State (ULS) design return periods, defined by an equation given in NZS1170.5:2004, were calculated and have been used as a calibration technique to find and compare specific PGA values for pseudo-static analyses in the Slovens Creek Viaduct area. The main purpose has been to provide an indication of how railway infrastructure could be affected by seismic events of various return periods defined by ULS design standards for the area. Limit equilibrium circular slip surface search methods, both grid search and auto refine search, indicated the slope is stable with a FoS greater than 1.0 returned from each, although one particular surface returned the lowest FoS in each. This surface is in the lower portion of the slope, adjacent to Slovens Stream and northeast of the MDL. As expected, pseudo-static analyses returned a lower FoS overall when compared to limit equilibrium analyses. The PGA analyses suggest that partial ground failure at the Slovens Creek Viaduct western abutment could occur in a 1 in 25-year return period event within materials on the slower slope beyond the immediate rail corridor. A ULS (1 in 500-year) event in the Slovens Creek Viaduct area would likely produce a PGA of ~0.9g, and the effects on the western abutment and rail infrastructure would most likely be catastrophic. Observed ground conditions for the western abutment of the Slovens Creek Viaduct suggest there is no movement within the landslide at depth within the monitoring timeframe of this project (22 May 2015 – 4 August 2015). Slope stability monitoring is recommended to be continued in two parts: (1) the inclinometer in BH1 is to be monitored on a six monthly basis for one year following completion of this thesis, and then annually unless ground movements become evident; and (2) surface movement monitoring should be installed using a fixed datum on the stable eastern abutment. Long-term stability management strategies for the Slovens Creek Viaduct western abutment are dependent upon future observed changes and ongoing monitoring. Hazard and risk assessment using the KiwiRail Qualitative Risk Assessment Framework (QRA) is recommended, and if slope stability becomes problematic for operation of the Midland Line consideration should be given to deep slope drainage. In the event of a large magnitude or high PGA earthquake all monitoring should be reviewed.
At 00:02 on 14th November 2016, a Mw 7.8 earthquake occurred in and offshore of the northeast of the South Island of New Zealand. Fault rupture, ground shaking, liquefaction, and co-seismic landslides caused severe damage to distributed infrastructure, and particularly transportation networks; large segments of the country’s main highway, State Highway 1 (SH1), and the Main North Line (MNL) railway line, were damaged between Picton and Christchurch. The damage caused direct local impacts, including isolation of communities, and wider regional impacts, including disruption of supply chains. Adaptive measures have ensured immediate continued regional transport of goods and people. Air and sea transport increased quickly, both for emergency response and to ensure routine transport of goods. Road diversions have also allowed critical connections to remain operable. This effective response to regional transport challenges allowed Civil Defence Emergency Management to quickly prioritise access to isolated settlements, all of which had road access 23 days after the earthquake. However, 100 days after the earthquake, critical segments of SH1 and the MNL remain closed and their ongoing repairs are a serious national strategic, as well as local, concern. This paper presents the impacts on South Island transport infrastructure, and subsequent management through the emergency response and early recovery phases, during the first 100 days following the initial earthquake, and highlights lessons for transportation system resilience.
At 00:02 on 14th November 2016, a Mw 7.8 earthquake occurred in and offshore of the northeast of the South Island of New Zealand. Fault rupture, ground shaking, liquefaction, and co-seismic landslides caused severe damage to distributed infrastructure, and particularly transportation networks; large segments of the country’s main highway, State Highway 1 (SH1), and the Main North Line (MNL) railway line, were damaged between Picton and Christchurch. The damage caused direct local impacts, including isolation of communities, and wider regional impacts, including disruption of supply chains. Adaptive measures have ensured immediate continued regional transport of goods and people. Air and sea transport increased quickly, both for emergency response and to ensure routine transport of goods. Road diversions have also allowed critical connections to remain operable. This effective response to regional transport challenges allowed Civil Defence Emergency Management to quickly prioritise access to isolated settlements, all of which had road access 23 days after the earthquake. However, 100 days after the earthquake, critical segments of SH1 and the MNL remain closed and their ongoing repairs are a serious national strategic, as well as local, concern. This paper presents the impacts on South Island transport infrastructure, and subsequent management through the emergency response and early recovery phases, during the first 100 days following the initial earthquake, and highlights lessons for transportation system resilience.
Nowadays the telecommunication systems’ performance has a substantial impact on our lifestyle. Their operationality becomes even more substantial in a post-disaster scenario when these services are used in civil protection and emergency plans, as well as for the restoration of all the other critical infrastructure. Despite the relevance of loss of functionality of telecommunication networks on seismic resilience, studies on their performance assessment are few in the literature. The telecommunication system is a distributed network made up of several components (i.e. ducts, utility holes, cabinets, major and local exchanges). Given that these networks cover a large geographical area, they can be easily subjected to the effects of a seismic event, either the ground shaking itself, or co-seismic events such as liquefaction and landslides. In this paper, an analysis of the data collected after the 2010-2011 Canterbury Earthquake Sequence (CES) and the 2016 Kaikoura Earthquake in New Zealand is conducted. Analysing these data, information gaps are critically identified regarding physical and functional failures of the telecommunication components, the timeline of repair/reconstruction activities and service recovery, geotechnical tests and land planning maps. Indeed, if these missing data were presented, they could aid the assessment of the seismic resilience. Thus, practical improvements in the post-disaster collection from both a network and organisational viewpoints are proposed through consultation of national and international researchers and highly experienced asset managers from Chorus. Finally, an outline of future studies which could guide towards a more resilient seismic performance of the telecommunication network is presented.
The greater Wellington region, New Zealand, is highly vulnerable to large earthquakes. While attention has been paid to the consequences of earthquake damage to road, electricity and water supply networks, the consequences of wastewater network damage for public health, environmental health and habitability of homes remain largely unknown for Wellington City. The Canterbury and Kaikōura earthquakes have highlighted the vulnerability of sewerage systems to disruption during a disaster. Management of human waste is one of the critical components of disaster planning to reduce faecal-oral transmission of disease and exposure to disease-bearing vectors. In Canterbury and Kaikōura, emergency sanitation involved a combination of Port-a-loos, chemical toilets and backyard long-drops. While many lessons may be learned from experiences in Canterbury earthquakes, it is important to note that isolation is likely to be a much greater factor for Wellington households, compared to Christchurch, due to the potential for widespread landslides in hill suburbs affecting road access. This in turn implies that human waste may have to be managed onsite, as options such as chemical toilets and Port-a-loos rely completely on road access for delivering chemicals and collecting waste. While some progress has been made on options such as emergency composting toilets, significant knowledge gaps remain on how to safely manage waste onsite. In order to bridge these gaps, laboratory tests will be conducted through the second half of 2019 to assess the pathogen die-off rates in the composting toilet system with variables being the type of carbon bulking material and the addition of a Bokashi composting activator.
Farming and urban regions are impacted by earthquake disasters in different ways, and feature a range of often different recovery requirements. In New Zealand, and elsewhere, most earthquake impact and recovery research is urban focused. This creates a research deficit that can lead to the application of well-researched urban recovery strategies in rural areas to suboptimal effect. To begin to reduce this deficit, in-depth case studies of the earthquake impacts and recovery of three New Zealand farms severely impacted by the 14th November 2016, M7.8 Hurunui-Kaikōura earthquake were conducted. The initial earthquake, its aftershocks and coseismic hazards (e.g., landslides, liquefaction, surface rupture) affected much of North Canterbury, Marlborough and the Wellington area. The three case study farms were chosen to broadly represent the main types of farming and topography in the Hurunui District in North Canterbury. The farms were directly and indirectly impacted by earthquakes and related hazards. On-farm infrastructure (e.g., woolsheds, homesteads) and essential services (e.g., water, power), frequently sourced from distributed networks, were severely impacted. The earthquake occurred after two years of regional drought had already stressed farm systems and farmers to restructuring or breaking point. Cascading interlinked hazards stemming from the earthquakes and coseismic hazards continued to disrupt earthquake recovery over a year after the initial earthquake. Semi-structured interviews with the farmers were conducted nine and fourteen months after the initial earthquake to capture the timeline of on-going impacts and recovery. Analysis of both geological hazard data and interview data resulted in the identification of key factors influencing farm level earthquake impact and recovery. These include pre-existing conditions (e.g., drought); farm-specific variations in recovery timelines; and resilience strategies for farm recovery resources. The earthquake recovery process presented all three farms with opportunities to change their business plans and adapt to mitigate on-going and future risk.
The 22 February 2011, Mw6.2-6.3 Christchurch earthquake is the most costly earthquake to affect New Zealand, causing 181 fatalities and severely damaging thousands of residential and commercial buildings, and most of the city lifelines and infrastructure. This manuscript presents an overview of observed geotechnical aspects of this earthquake as well as some of the completed and on-going research investigations. A unique aspect, which is particularly emphasized, is the severity and spatial extent of liquefaction occurring in native soils. Overall, both the spatial extent and severity of liquefaction in the city was greater than in the preceding 4th September 2010 Darfield earthquake, including numerous areas that liquefied in both events. Liquefaction and lateral spreading, variable over both large and short spatial scales, affected commercial structures in the Central Business District (CBD) in a variety of ways including: total and differential settlements and tilting; punching settlements of structures with shallow foundations; differential movements of components of complex structures; and interaction of adjacent structures via common foundation soils. Liquefaction was most severe in residential areas located to the east of the CBD as a result of stronger ground shaking due to the proximity to the causative fault, a high water table approximately 1m from the surface, and soils with composition and states of high susceptibility and potential for liquefaction. Total and differential settlements, and lateral movements, due to liquefaction and lateral spreading is estimated to have severely compromised 15,000 residential structures, the majority of which otherwise sustained only minor to moderate damage directly due to inertial loading from ground shaking. Liquefaction also had a profound effect on lifelines and other infrastructure, particularly bridge structures, and underground services. Minor damage was also observed at flood stop banks to the north of the city, which were more severely impacted in the 4th September 2010 Darfield earthquake. Due to the large high-frequency ground motion in the Port hills numerous rock falls and landslides also occurred, resulting in several fatalities and rendering some residential areas uninhabitable.
The Amuri Earthquake of September 1, 1888 (magnitude M = 6.5 to 6.8) occurred on the Hope River Segment of the Hope Fault west of Hanmer Plains. The earthquake was felt strongly in North Canterbury and North Westland and caused considerable property damage and landsliding in the Lower Hope Valley. However, damage reports and the spatial distribution of felt intensities emphasize extreme variations in seismic effects over short distances, probably due to topographic focusing and local ground conditions. Significant variations in lateral fault displacement occurred at secondary fault segment boundaries (side-steps and bends in the fault trace) during the 1888 earthquake. This historical spatial variation in lateral slip is matched by the Late Quaternary geomorphic distribution of slip on the Hope River Segment of the Hope Fault. Trenching studies at two sites on the Hope Fault have also identified evidence for five pre-historic earthquakes of similar magnitude to the 1888 earthquake and an average recurrence interval of 134 ± 27 years between events. Magnitude estimates for the 1888 earthquake are combined with a. strong ground motion attenuation expression to provide an estimate of potential ground accelerations in Amuri District during-future earthquakes on the Hope River Segment of the Hope Fault. The predicted acceleration response on bedrock sites within 20 km of the epicentral region is between 0.23 g and 0.34 g. The close match between the historic, inferred pre-historic and geomorphic distribution of lateral slip indicates that secondary fault segmentation exerts a strong structural control on rupture propagation and the expression of fault displacement at the surface. In basement rocks at depth the spatial variations in slip are inferred to be distributed within zones of pervasive cataclastic shear, on either side of the fault segment boundaries. The large variations in surface displacement across fault segment boundaries means that one must know the geometry of the fault in order to evaluate slip-rates calculated from individual locations. The average Late Quaternary slip-rate on the Hope Fault at Glynn Wye Station is between 15.5 mm/yr and 18.25 mm/yr and the rate on the subsidiary Kakapo Fault is between 5.0 mm/yr and 7.5 mm/yr. These rates have been determined from sites which are relatively free of structural complication.
Following the 22nd February 2011, Mw 6.2 earthquake located along a previously unknown fault beneath the Port Hills of Christchurch, surface cracking was identified in contour parallel locations within fill material at Quarry Road on the lower slopes of Mount Pleasant. GNS Science, in the role of advisor to the Christchurch City Council, concluded that these cracks were a part of a potential rotational mass movement (named zone 11A) within the fill and airfall loess material present. However, a lack of field evidence for slope instability and an absence of laboratory geotechnical data on which slope stability analysis was based, suggested this conclusion is potentially incorrect. It was hypothesised that ground cracking was in fact due to earthquake shaking, and not mass movement within the slope, thus forming the basis of this study. Three soil units were identified during surface and subsurface investigations at Quarry Road: fill derived from quarry operations in the adjacent St. Andrews Quarry (between 1893 and 1913), a buried topsoil, and underlying in-situ airfall loess. The fill material was identified by the presence of organic-rich topsoil “clods” that were irregular in both size (∼10 – 200 mm) and shape, with variable thicknesses of 1 – 10 m. Maximum thickness, as indicated by drill holes and geophysical survey lines, was identified below 6 Quarry Road and 7 The Brae where it is thought to infill a pre-existing gully formed in the underlying airfall loess. Bearing strength of the fill consistently exceeded 300 kPa ultimate below ∼500 mm depth. The buried topsoil was 200 – 300 mm thick, and normally displayed a lower bearing strength when encountered, but not below 300 kPa ultimate (3 – 11 blows per 100mm or ≥100 kPa allowable). In-situ airfall loess stood vertically in outcrop due to its characteristic high dry strength and also showed Scala penetrometer values of 6 – 20+ blows per 100 mm (450 – ≥1000 kPa ultimate). All soils were described as being moist to dry during subsurface investigations, with no groundwater table identified during any investigation into volcanic bedrock. In-situ moisture contents were established using bulk disturbed samples from hand augers and test pitting. Average moisture contents were low at 9% within the fill, 11 % within the buried topsoil, and 8% within the airfall loess: all were below the associated average plastic limit of 17, 15, and 16, respectively, determined during Atterberg limit analysis. Particle size distributions, identified using the sieve and pipette method, were similar between the three soil units with 11 – 20 % clay, 62 – 78 % silt, and 11 – 20 % fine sand. Using these results and the NZGS soil classification, the loess derived fill and in-situ airfall loess are termed SILT with some clay and sand, and the buried topsoil is SILT with minor clay and sand. Dispersivity of the units was found using the Emerson crumb test, which established that the fill can be non- to completely dispersive (score 0 – 4). The buried topsoil was always non-dispersive (score 0), and airfall loess completely dispersive (score 4). Values for cohesion (c) and internal friction angle (φ) of the three soil units were established using the direct shear box at field moisture contents. Results showed all soil units had high shear strengths at the moisture contents tested (c = 18 – 24 kPa and φ = 42 – 50°), with samples behaving in a brittle fashion. Moisture content was artificially increased to 16% within the buried topsoil, which reduced the shear strength (c = 10 kPa, φ = 18°) and allowed it to behave plastically. Observational information indicating stability at Quarry Road included: shallow, discontinuous, cracks that do not display vertical offset; no scarp features or compressional zones typical of landsliding; no tilted or deformed structures; no movement in inclinometers; no basal shear zone identified in logged core to 20 m depth; low field moisture contents; no groundwater table; and high soil strength using Scala penetrometers. Limit equilibrium analysis of the slope was conducted using Rocscience software Slide 5.0 to verify the slope stability identified by observational methods. Friction, cohesion, and density values determined during laboratory were input into the two slope models investigated. Results gave minimum static factor of safety values for translational (along buried topsoil) and rotational (in the fill) slides of 2.4 – 4.2. Sensitivity of the slope to reduced shear strength parameters was analysed using c = 10 kPa and φ = 18° for the translational buried topsoil plane, and a cohesion of 0 kPa within the fill for the rotational plane. The only situation that gave a factor of safety <1.0 was in nonengineered fill at 0.5 m depth. Pseudostatic analysis based on previous peak ground acceleration (PGA) values for the Canterbury Earthquake Sequence, and predicted PGAs for future Alpine Fault and Hope Fault earthquakes established minimum factor of safety values between 1.2 and 3.3. Yield acceleration PGAs were computed to be between 0.8g and 1.6g. Based on all information gathered, the cracking at Quarry Road is considered to be shallow deformation in response to earthquake shaking, and not due to deep-seated landsliding. It is recommended that the currently bare site be managed by smoothing the land, installing contour drainage, and bioremediation of the surface soils to reduce surface water infiltration and runoff. Extensive earthworks, including removal of the fill, are considered unnecessary. Any future replacement of housing would be subject to site-specific investigations, and careful foundation design based on those results.
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.
The Porter's Pass-Amberley Fault Zone (PPAFZ) is a complex zone of anastomosing faults and folds bounding the south-eastern edge of the transition from subducting Pacific Plate to continental collision on the Australia Plate boundary. This study combines mapping of a 2000 km2 zone from the Southern Alps northeast to the coast near Amberley, 40 km north of metropolitan Christchurch, with an analysis of seismicity and a revision of regional seismic hazard. Three structural styles: 1) a western strike-slip, and 2) a more easterly thrust and reverse domain, pass into 3) a northwest verging fold belt on the northern Canterbury Plains, reflecting the structural levels exposed and the evolving west to east propagation. Basal remnants of a Late Cretaceous-Cenozoic, largely marine sedimentary cover sequence are preserved as outliers that unconformably overlie Mesozoic basement (greywacke and argillite of the Torlesse terrain) in the mountains of the PPAFZ and are underlain by a deeply leached zone which is widely preserved. Structure contouring of the unconformity surface indicates maximum, differential uplift of c.2600 m in the southwest, decreasing to c.1200 m in the coastal fold belt to the northeast. Much lower rates (or reversal) of uplift are evident a few kilometres southeast of the PPAFZ range-front escarpment. The youngest elements of the cover sequence are basement-derived conglomerates of Plio-Pleistocene age preserved on the SE margin. The source is more distant than the intervening mountains of the PPAFZ, probably from the Southern Alps, to the west and northwest. The absence of another regional unconformity on Mesozoic basement, older than Pleistocene, indicates that this uplift is post-Pliocene. Late Pleistocene(<100 kyr) differential uplift rates of c.0.5-2.7 m/kyr from uplifted marine terraces at the east coast, and rates of 2.5-3.3 m/kyr for tectonically-induced river-down cutting further west, suggest that uplift commenced locally during the last 1 Ma, and possibly within the last 0.5 Ma, if average rates are assumed to be uniform over time. Analysis of seismicity, recorded during a 10 week regional survey of micro earthquakes in 1990, identified two seismic zones beneath North Canterbury: 1) a sub-horizontal zone of activity restricted to the upper crust (≤12 km); and 2) a seismic zone in the lower crust (below a ceiling of ≤17 km), that broadens vertically to the north and northwest to a depth of c.40 km, with a bottom edge which dips 10°N and 15°NW, respectively. No events were recorded at depths between 12 km and 17 km, which is interpreted as a relatively aseismic, mid-crustal ductile layer. Marked differences (up to 60°) in the trend of strain axes for events above and below the inferred ductile layer are observed only north of the PPAFZ. A fundamental, north-to-south increase in the Wave-length of major geological structures occurs across the PPAFZ, and is interpreted as evidence that the upper crust beneath the Canterbury Plains is coupled to the lower crust, whereas the upper crust further north is not. Most of the recorded micro earthquakes <12 km deep beneath the PPAFZ have strike-slip mechanisms. It is probable that faults splay upward into the thrusts and folds at the surface as an evolving transpression zone in response to deep shear in basement. There have been no historic surface ruptures of the PPAFZ, but the zone has been characterised historically by frequent small earthquakes. Paleoseismic data (dated landslides and surface ruptures) compiled in this study, indicate a return period of 1500-1900 years between the last two M>7-7.5 earthquakes, and 500-700 years have elapsed since the last. The magnitudes of these events are estimated at c.M7.5, which represents a probable maximum magnitude for the PPAFZ. There are insufficient data to determine whether or not the frequency of large earthquakes conforms to a recognised model of behaviour, but comparison of the paleoseismic data with the historic record of smaller earthquakes, suggests that the magnitudes of the largest earthquakes in this zone are not exponentially distributed. A seismicity model for the PPAFZ (Elder et al., 1991) is reviewed, and a b-value of 1.0 is found to be consistent with the newly acquired paleoseismic data. This b-value reduces the predicted frequency of large earthquakes (M≥7.0) in this zone by a factor of 3.5, while retaining a conservative margin that allows for temporal variations in the frequency of large events and the possibility that the geological database is incomplete, suggesting grounds for revising the hazard model for Christchurch.
