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.
Object Overview of 'Greendale Fault: investigation of surface rupture characteristics for fault avoidance zonation (Villamor et al, 2011).'
Oblique convergence of the Pacific and Australian Plates is accommodated in the northern South Island by the Marlborough Fault System. The Hope Fault is the southern of four major dextral strike-slip faults of this system. Hanmer Basin is a probable segment boundary between the Hope River and Conway segments of the Hope Fault. The Conway segment is transpressional and shows increasing structural complexity near the segment boundary at Hanmer Basin, with multiple Late Quaternary traces, and fault-parallel folding in response to across-fault shortening. Between Hossack Station and Hanmer Basin a crush zone in excess of one kilometre wide is exposed in incised streams and rivers. The crush zone has an asymmetrical geometry about the active trace of the Hope Fault, being only 100-300 metres wide south of the fault, and more than 500 metres wide north of the fault. The most intense deformation of Torlesse bedrock occurs at the south side of the fault zone, indicating that strain is accommodated against the fault footwall. North of the fault deformation is less intense, but occurs over a wider area. The wide fault zone at Hossack Station may reflect divergence of the Hanmer Fault, a major splay of the Hope Fault. At Hossack Station, the Hope Fault has accommodated at least 260 metres of dextral displacement during the Holocene. Dating of abandoned stream channels, offset by the Hope Fault, indicated a Late Holocene dextral slip-rate of 18±8 mm-¹ for the west end of the Conway segment. Using empirical formulae and inferred fault parameters, the expected magnitude of an earthquake generated by the Conway segment is M6.9 to M7.4; for an exceedence probability of 10%, the magnitude is M7.7 to M7.9. Effects associated with coseismic rupture of the Conway segment include shaking of up to MMIX along the ruptured fault and at Hanmer Basin. Uplift at the east end of Hanmer Basin, in conjunction with subsidence at the southwest margin of the basin, is resulting in the development of onlapping stratigraphy. Seismic reflection profiles support this theory. Possible along-fault migration of the basin is inferred to be a consequence of non-parallelism of the master faults.
This study contains an evaluation of the seismic hazard associated with the Springbank Fault, a blind structure discovered in 1998 close to Christchurch. The assessment of the seismic hazard is approached as a deterministic process in which it is necessary to establish: 1) fault characteristics; 2) the maximum earthquake that the fault is capable of producing and 3) ground motions estimations. Due to the blind nature of the fault, conventional techniques used to establish the basic fault characteristics for seismic hazard assessments could not be applied. Alternative methods are used including global positioning system (GPS) surveys, morphometric analyses along rivers, shallow seismic reflection surveys and computer modelling. These were supplemented by using multiple empirical equations relating fault attributes to earthquake magnitude, and attenuation relationships to estimate ground motions in the near-fault zone. The analyses indicated that the Springbank Fault is a reverse structure located approximately 30 km to the northwest of Christchurch, along a strike length of approximately 16 km between the Eyre and Ashley River. The fault does not reach the surface, buy it is associated with a broad anticline whose maximum topographic expression offers close to the mid-length of the fault. Two other reverse faults, the Eyrewell and Sefton Faults, are inferred in the study area. These faults, together with the Springbank and Hororata Faults and interpreted as part of a sys of trust/reverse faults propagating from a decollement located at mid-crustal depths of approximately 14 km beneath the Canterbury Plains Within this fault system, the Springbank Fault is considered to behave in a seismically independent way, with a fault slip rate of ~0.2 mm/yr, and the capacity of producing a reverse-slip earthquake of moment magnitude ~6.4, with an earthquake recurrence of 3,000 years. An earthquake of the above characteristics represents a significant seismic hazard for various urban centres in the near-fault zone including Christchurch, Rangiora, Oxford, Amberley, Kaiapoi, Darfield, Rollestion and Cust. Estimated peak ground accelerations for these towns range between 0.14 g to 0.5 g.
Recent surface-rupturing earthquakes in New Zealand have highlighted significant exposure and vulnerability of the road network to fault displacement. Understanding fault displacement hazard and its impact on roads is crucial for mitigating risks and enhancing resilience. There is a need for regional-scale assessments of fault displacement to identify vulnerable areas within the road network for the purposes of planning and prioritising site-specific investigations. This thesis employs updated analysis of data from three historical surface-rupturing earthquakes (Edgecumbe 1987, Darfield 2010, and Kaikoūra 2016) to develop an empirical model that addresses the gap in regional fault displacement hazard analysis. The findings contribute to understanding of • How to use seismic hazard model inputs for regional fault displacement hazard analysis • How faulting type and sediment cover affects the magnitude and spatial distribution of fault displacement • How the distribution of displacement and regional fault displacement hazard is impacted by secondary faulting • The inherent uncertainties and limitations associated with employing an empirical approach at a regional scale • Which sections of New Zealand’s roading network are most susceptible to fault displacement hazard and warrant site-specific investigations • Which regions should prioritise updating emergency management plans to account for post-event disruptions to roading. I used displacement data from the aforementioned historical ruptures to generate displacement versus distance-to-fault curves for various displacement components, fault types, and geological characteristics. Using those relationships and established relationships for along-strike displacement, displacement contours were generated surrounding active faults within the NZ Community Fault Model. Next, I calculated a new measure of 1D strain along roads as well as relative hazard, which integrated 1D strain and normalised slip rate data. Summing these values at the regional level identified areas of heightened relative hazard across New Zealand, and permits an assessment of the susceptibility of road networks using geomorphon land classes as proxies for vulnerability. The results reveal that fault-parallel displacements tend to localise near the fault plane, while vertical and fault-perpendicular displacements sustain over extended distances. Notably, no significant disparities were observed in off-fault displacement between the hanging wall and footwall sides of the fault, or among different surface geology types, potentially attributed to dataset biases. The presence of secondary faulting in the dataset contributes to increased levels of tectonic displacement farther from the fault, highlighting its significance in hazard assessments. Furthermore, fault displacement contours delineate broader zones around dip-slip faults compared to strike-slip faults, with correlations identified between fault length and displacement width. Road ‘strain’ values are higher around dip-slip faults, with notable examples observed in the Westland and Buller Districts. As expected, relative hazard analysis revealed elevated values along faults with high slip rates, notably along the Alpine Fault. A regional-scale analysis of hazard and exposure reveals heightened relative hazard in specific regions, including Wellington, Southern Hawke’s Bay, Central Bay of Plenty, Central West Coast, inland Canterbury, and the Wairau Valley of Marlborough. Notably, the Central West Coast exhibits the highest summed relative hazard value, attributed to the fast-slipping Alpine Fault. The South Island generally experiences greater relative hazard due to larger and faster-slipping faults compared to the North Island, despite having fewer roads. Central regions of New Zealand face heightened risk compared to Southern or Northern regions. Critical road links intersecting high-slipping faults, such as State Highways 6, 73, 1, and 2, necessitate prioritisation for site-specific assessments, emergency management planning and targeted mitigation strategies. Roads intersecting with the Alpine Fault are prone to large parallel displacements, requiring post-quake repair efforts. Mitigation strategies include future road avoidance of nearby faults, modification of road fill and surface material, and acknowledgement of inherent risk, leading to prioritised repair efforts of critical roads post-quake. Implementing these strategies enhances emergency response efforts by improving accessibility to isolated regions following a major surface-rupturing event, facilitating faster supply delivery and evacuation assistance. This thesis contributes to the advancement of understanding fault displacement hazard by introducing a novel regional, empirical approach. The methods and findings highlight the importance of further developing such analyses and extending them to other critical infrastructure types exposed to fault displacement hazard in New Zealand. Enhancing our comprehension of the risks associated with fault displacement hazard offers valuable insights into various mitigation strategies for roading infrastructure and informs emergency response planning, thereby enhancing both national and global infrastructure resilience against geological hazards.
This report provided information on the location and character of the Ostler Fault Zone near Twizel. The fault traces, and associated recommended fault avoidance zones, were mapped in detail for inclusion in a District Plan Change for the Twizel area. The Ostler Fault Zone was mapped in detail because of the higher likelihood of movement on that fault than others in the district, and the potential for future development across the fault zone because of its proximity to Twizel. See Object Overview for background and usage information. The report recommended that the information be incorporated into the District Plan Change and that site-specific investigations be undertaken before development is allowed within the fault avoidance zones. These recommendations were taken up by Mackenzie District Council.
The previously unknown Greendale Fault ruptured to the ground surface, causing up to 5 metres horizontal and 1 metre vertical permanent offset of the ground, during the September 2010 Darfield (Canterbury) earthquake. Environment Canterbury commissioned GNS Science, with help from the University of Canterbury, to define a fault avoidance zone and to estimate the fault recurrence interval. There is little evidence for past movement on the fault in the past 16,000 years. However, because of the uncertainties involved, a conservative approach was taken and the fault has been categorised as a Recurrence Interval Class IV fault (a recurrence interval of between 5,000 and 10,000 years). A PhD study by a University of Canterbury student will work towards refining the Recurrence Interval Class over the next three years. Taking a risk-based approach, the Ministry for the Environment Active Fault Guidelines recommend that normal residential development be allowed within the fault avoidance zone for faults of this Recurrence Interval Class, but recommends restrictions for larger community buildings or facilities with post-disaster functions. The report is assisting Selwyn District Council in granting consents for rebuilding houses on or near the Greendale Fault that were damaged by permanent distortion of the ground due to the fault rupture in the September 2010 earthquake. The report provides specific recommendations for building on or close to the Greendale Fault, which are being implemented by Selwyn District Council. See Object Overview for background and usage information.
The Eastern Humps and Leader faults, situated in the Mount Stewart Range in North Canterbury, are two of the ≥17 faults which ruptured during the 2016 MW7.8 Kaikōura Earthquake. The earthquake produced complex, intersecting ground ruptures of these faults and the co-seismic uplift of the Mount Stewart Range. This thesis aims to determine how these two faults accommodated deformation during the 2016 earthquake and how they interact with each other and with pre-existing geological structures. In addition, it aims to establish the most likely subsurface geometry of the fault complex across the Mount Stewart Range, and to investigate the paleoseismic history of the Leader Fault. The Eastern Humps Fault strikes ~240° and dips 80° to 60° to the northwest and accommodated right- lateral – reverse-slip, with up to 4 m horizontal and 2 m vertical displacement in the 2016 earthquake. The strike of the Leader Fault varies from ~155 to ~300°, and dips ~30 to ~80° to the west/northwest, and mainly accommodated left-lateral – reverse-slip of up to 3.5 m horizontal and 3.5 m vertical slip in the 2016 earthquake. On both the Eastern Humps and Leader faults the slip is variable along strike, with areas of low total displacement and areas where horizontal and vertical displacement are negatively correlated. Fault traces with low total displacement reflect the presence of off-fault (distributed) displacement which is not being captured with field measurements. The negative correlation of horizontal and vertical displacement likely indicates a degree of slip partitioning during the 2016 earthquake on both the Eastern Humps and Leader faults. The Eastern Humps and Leader faults have a complex, interdependent relationship with the local bedrock geology. The Humps Fault appears to be a primary driver of ongoing folding and deformation of the local Mendip Syncline and folding of the Mount Stewart Range, which probably began prior to, or synchronous with, initial rupture of The Humps Fault. The Leader Fault appears to use existing lithological weaknesses in the Cretaceous-Cenozoic bedrock stratigraphy to rupture to the surface. This largely accounts for the strong variability on the strike and dip of the Leader Fault, as the geometry of the surface ruptures tend to reflect the strike and dip of the geological strata which it is rupturing through. The Leader Fault may also accommodate some degree of flexural slip in the Cenozoic cover sequence of the Mendip Syncline, contributing to the ongoing growth of the fold. The similarity between topography and uplift profiles from the 2016 earthquake suggest that growth of the Mount Stewart Range has been primarily driven by multiple (>500) discrete earthquakes that rupture The Humps and Leader faults. The spatial distribution of surface displacements across the Mount Stewart Range is more symmetrical than would be expected if uplift is driven primarily by The Humps and Leader faults alone. Elastic dislocation forward models were used to model potential sub-surface geometries and the resulting patterns of deformation compared to photogrammetry-derived surface displacements. Results show a slight preference for models with a steeply southeast-dipping blind fault, coincident with a zone of seismicity at depth, as a ‘backthrust’ to The Humps and Leader faults. This inferred Mount Stewart Fault accommodated contractional strain during the 2016 earthquake and contributes to the ongoing uplift of the Mount Stewart Range with a component of folding. Right-lateral and reverse shear stress change on the Hope Fault was also modelled using Coulomb 3.3 software to examine whether slip on The Humps and Leader faults could transfer enough stress onto the Hope Fault to trigger through-going rupture. Results indicate that during the 2016 earthquake right-lateral shear and reverse stress only increased on the Hope Fault in small areas to the west of the Leader Fault, and similar ruptures would be unlikely to trigger eastward propagating rupture unless the Hope Fault was close to failure prior to the earthquake. Paleoseismic trenches were excavated on the Leader Fault at four locations from 2018 to 2020, revealing near surface (< 4m depth) contractional deformation of Holocene stratigraphy. Three of the trench locations uncovered clear evidence for rupture of the Leader Fault prior to 2016, with fault displacement of near surface stratigraphy being greater than displacement recorded during the 2016 earthquake. Radiocarbon dating of in-situ organic material from two trenches indicate a date of the penultimate earthquake on the Leader Fault within the past 1000 years. This date is consistent with The Humps and Leader faults having ruptured simultaneously in the past, and with multi-fault ruptures involving The Humps, Leader, Hundalee and Stone Jug faults having occurred prior to the 2016 Kaikōura earthquake. Overall, the results contribute to an improved understanding of the Kaikōura earthquake and highlight the importance of detailed structural and paleoseismic investigations in determining controls on earthquake ‘complexity’.
This study investigates evidence for linkages and fault interactions centred on the Cust Anticline in Northwest Canterbury between Starvation Hill to the southwest and the Ashley and Loburn faults to the northeast. An integrated programme of geologic, geomorphic, paleo-seismic and geophysical analyses was undertaken owing to a lack of surface exposures and difficulty in distinguishing active tectonic features from fluvial and/or aeolian features across the low-relief Canterbury Plains. LiDAR analysis identified surface expression of several previously unrecognised active fault traces across the low-relief aggradation surfaces of the Canterbury Plains. Their presence is consistent with predictions of a fault relay exploiting the structural mesh across the region. This is characterised by interactions of northeast-striking contractional faults and a series of re-activating inherited Late Cretaceous normal faults, the latter now functioning as E–W-striking dextral transpressive faults. LiDAR also allowed for detailed analysis of the surface expression of individual faults and folds across the Cust Anticline contractional restraining bend, which is evolving as a pop-up structure within the newly established dextral shear system that is exploiting the inherited, now re-activated, basement fault zone. Paleo-seismic trenches were located on the crest of the western arm of the Cust Anticline and across a previously unrecognised E–W-striking fault trace, immediately southwest of the steeply plunging Cust Anticline termination. These studies confirmed the location and structural style of north-northeast-striking faults and an E–W-striking fault associated with the development of this structural culmination. A review of available industry seismic reflection lines emphasised the presence of a series of common structural styles having the same underlying structural drivers but with varying degrees of development and expression, both in the seismic profiles and in surface elevations across the study area. Based on LiDAR surface mapping and preliminary re-analysis of industry seismic reflection data, four fault zones are identified across the restraining bend structural culminations, which together form the proposed Oxford–Cust–Ashley Fault System. The 2010–2012 Canterbury Earthquake Sequence showed many similarities to the structural pattern established across the Oxford–Cust–Ashley Fault System, emphasising the importance of identification and characterization of presently hidden fault sources, and the understanding of fault network linkages, in order to improve constraints on earthquake source potential. Improved understanding of potentially-interactive fault sources in Northwest Canterbury, with the potential for combined initial fault rupture and spatial and temporal rupture propagation across this fault system, can be used in probabilistic seismic hazard analysis for the region, which is essential for the suitability and sustainability of future social and economic development.
New research shows that if a major earthquake ruptures the Alpine Fault Christchurch is likely to suffer more intense damage than some areas closer to the fault.
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.
This report provides information on the locations and character of active geological faults and folds in Ashburton District. The faults are mapped at a district scale and the information is intended to highlight areas where there is a risk of permanent fault movement at the ground surface, and where more detailed investigations should be done if development is proposed in that area (depending on the potential activity of the fault and the type of development proposed). See Object Overview for background and usage information. Most of the faults and folds identified at the ground surface in Ashburton District are in rural or very sparsely populated areas. In addition, most of the faults have relatively long recurrence intervals (long-term average time between fault movements) in the order of several thousand years. Following the Ministry for the Environment Active Fault Guidelines, normal residential development would be allowed on or near faults with recurrence intervals this long. There are no recommendations associated with this report. The information in the report will be reviewed as required, after the remaining district reports are completed in the region.
This report provides information on the locations and character of active geological faults and folds in Mackenzie District. The faults are mapped at a district scale and the information is intended to highlight areas where there is a risk of fault movement, and where more detailed investigations should be done if development is proposed in that area(depending on the potential activity of the fault and the type of development proposed). Most of the faults and folds identified at the ground surface in Mackenzie District are in rural or very sparsely populated areas. In addition, most of the faults have relatively long recurrence intervals (long-term average time between fault movements) in the order of several thousand years. Following the Ministry for the Environment Active Fault Guidelines, normal residential development would be allowed on or near faults with recurrence intervals this long. There are no recommendations associated with this report. The information in the report will be reviewed as required, after the remaining district reports are completed in the region. See Object Overview for background and usage information.
Surface rupture and slip from the Mw 7.8 2016 Kaikōura Earthquake have been mapped in the region between the Leader and Charwell rivers using field mapping and LiDAR data. The eastern Humps, north Leader and Conway-Charwell faults ruptured the ground surface in the study area. The E-NE striking ‘The Humps’ Fault runs along the base of the Mt Stewart range front, appears to dip steeply NW and intersects the NNW-NNE Leader Fault which itself terminates northwards at the NE striking Conway-Charwell Fault. The eastern Humps Fault is up to the NW and accommodates oblique slip with reverse and right lateral displacement. Net slip on ‘The Humps’ Fault is ≤4 m and produced ≤4 m uplift of the Mt Stewart range during the earthquake. The Leader Fault strikes NNW-NNE with dips ranging from ~10° west to 80° east and accommodated ≤4 m net slip comprising left-lateral and up-to-the-west vertical displacement. Like the Humps west of the study area, surface-rupture of the Leader Fault occurred on multiple strands. The complexity of rupture on the Leader Fault is in part due to the occurrence of bedding-parallel slip within the Cretaceous-Cenozoic sequence. Although the Mt Stewart range front is bounded by ‘The Humps’ Fault, in the study area neither this fault nor the Leader Fault were known to have been active before the earthquake. Fieldwork and trenching investigations are ongoing to characterise the geometry, kinematics and paleoseismic history of the mapped active faults.
Oblique-convergent plate collision between the Pacific and Australian plates across the South Island has resulted in shallow, upper crustal earthquake activity and ground surface deformation. In particular the Porters Pass - Amberley Fault Zone displays a complex hybrid zone of anastomosing dextral strike-slip and thrust/reverse faulting which includes the thrust/reverse Lees Valley Fault Zone and associated basin deformation. There is a knowledge gap with respect to the paleoseismicity of many of the faults in this region including the Lees Valley Fault Zone. This study aimed to investigate the earthquake history of the fault at a selected location and the structural and geomorphic development of the Lees Valley Fault Zone and eastern rangefront. This was investigated through extensive structural and geomorphic mapping, GPS field surveying, vertical aerial photo interpretation, analysis of Digital Elevation Models, paleoseismic trenching and optically stimulated luminescence dating. This thesis used a published model for tectonic geomorphology development of mountain rangefronts to understand the development of Lees Valley. Rangefront geomorphology is investigated through analysis of features such as rangefront sinuosity and faceted spurs and indicates the recently active and episodic nature of the uplifted rangefront. Analysis of fault discontinuity, fault splays, distribution of displacement, fault deformation zone and limited exposure of bedrock provided insight into the complex structure of the fault zone. These observations revealed preserved, earlier rangefronts, abandoned and uplifted within the eastern ranges, indicating a basinward shift in focus of faulting and an imbricate thrust wedge development propagating into the footwall of the fault zone and along the eastern ranges of Lees Valley. Fault scarp deformation analysis indicated multiple events have produced the deformation present preserved by the active fault trace in the northern valley. Vertical deformation along this scarp varied with a maximum of 11.5 m and an average of 5 m. Field mapping revealed fan surfaces of various ages have been offset and deformed, likely during the Holocene, based on expected relative surface ages. Geomorphic and structural mapping highlighted the effect of cross-cutting and inherited structures on the Lees Valley Fault, resulting in a step-over development in the centre of the eastern range-bounding trace. Paleoseismic trenching provided evidence of at least two earthquakes, which were constrained to post 21.6 ± 2.3 ka by optically stimulated luminescence dating. Single event displacements (1.48 ± 0.08 m), surface rupture earthquake magnitudes (Mw 6.7 ± 0.1, with potential to produce ≥ 7.0), and a minimum recurrence interval (3.6 ± 0.3 ka) indicated the Lees Valley Fault is an active structure capable of producing significant earthquake events. Results from this study indicate that the Lees Valley Fault Zone accommodates an important component of the Porters Pass - Amberley Fault Zone deformation and confirms the fault as a source of potentially damaging, peak ground accelerations in the Canterbury region. Remnants of previous rangefronts indicate a thrust wedge development of the Lees Valley Fault Zone and associated ranges that can potentially be used as a model of development for other thrust-fault bounded basins.
The Mw 7.1 Darfield earthquake generated a ~30 km long surface rupture on the Greendale Fault and significant surface deformation related to related blind faults on a previously unrecognized fault system beneath the Canterbury Plains. This earthquake provided the opportunity for research into the patterns and mechanisms of co-seismic and post-seismic crustal deformation. In this thesis I use multiple across-fault EDM surveys, logic trees, surface investigations and deformation feature mapping, seismic reflection surveying, and survey mark (cadastral) re-occupation using GPS to quantify surface displacements at a variety of temporal and spatial scales. My field mapping investigations identified shaking and crustal displacement-induced surface deformation features south and southwest of Christchurch and in the vicinity of the projected surface traces of the Hororata Blind and Charing Cross Faults. The data are consistent with the high peak ground accelerations and broad surface warping due to underlying reverse faulting on the Hororata Blind Fault and Charing Cross Fault. I measured varying amounts of post-seismic displacement at four of five locations that crossed the Greendale Fault. None of the data showed evidence for localized dextral creep on the Greendale Fault surface trace, consistent with other studies showing only minimal regional post-seismic deformation. Instead, the post-seismic deformation field suggests an apparent westward translation of northern parts of the across-fault surveys relative to the southern parts of the surveys that I attribute to post-mainshock creep on blind thrusts and/or other unidentified structures. The seismic surveys identified a deformation zone in the gravels that we attribute to the Hororata Blind Fault but the Charing Cross fault was not able to be identified on the survey. Cadastral re-surveys indicate a deformation field consistent with previously published geodetic data. We use this deformation with regional strain rates to estimate earthquake recurrence intervals of ~7000 to > 14,000 yrs on the Hororata Blind and Charing Cross Faults.
Object Overview of 'Fault mapping studies.'
New research says Canterbury would suffer severe damage in an Alpine Fault earthquake.
"There are major faults in low-seismicity areas that could produce a major event."
The timing of large Holocene prehistoric earthquakes is determined by dated surface ruptures and landslides at the edge of the Australia-Pacific plate boundary zone in North Canterbury, New Zealand. Collectively, these data indicate two large (M > 7) earthquakes during the last circa 2500 years, within a newly formed zone of hybrid strike-slip and thrust faulting herein described as the Porter's Pass-to-Amberley Fault Zone (PPAFZ). Two earlier events during the Holocene are also recognized, but the data prior to 2500 years are presumed to be incomplete. A return period of 1300–2000 years between large earthquakes in the PPAFZ is consistent with a late Holocene slip rate of 3–4 mm/yr if each displacement is in the range 4–8 m. Historical seismicity in the PPAFZ is characterized by frequent small and moderate magnitude earthquakes and a seismicity rate that is identical to a region surrounding the structurally mature Hope fault of the Marlborough Fault System farther north. This is despite an order-of-magnitude difference in slip rate between the respective fault zones and considerable differences in the recurrence rate of large earthquakes. The magnitude-frequency distribution in the Hope fault region is in accord with the characteristic earthquake model, whereas the rate of large earthquakes in the PPAFZ is approximated (but over predicted) by the Gutenberg-Richter model. The comparison of these two fault zones demonstrates the importance of the structural maturity of the fault zone in relation to seismicity rates inferred from recent, historical, and paleoseismic data.
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.
The Leader Fault was one of at least 17 faults that ruptured the ground surface across the northeastern South Island of New Zealand during the Mw 7.8 2016 Kaikōura Earthquake. The southern ~6 km of the Leader Fault, here referred to as the South Leader Fault (SLF), ruptured the North Canterbury (tectonic) Domain and is the primary focus of this study. The main objective of the thesis is to understand the key factors that contributed to the geometry and kinematics of the 2016 SLF rupture and its intersection with The Humps Fault (HF). This thesis employs a combination of techniques to achieve the primary objective, including detailed mapping of the bedrock geology, geomorphology and 2016 rupture, measurement of 2016 ground surface displacements, kinematic analysis of slip vectors from the earthquake, and logging of a single natural exposure across a 2016 rupture that was treated as a paleoseismic trench. The resulting datasets were collected in the field, from terrestrial LiDAR and InSAR imagery, and from historical (pre-earthquake) aerial photographs for a ~11 km2 study area. Surface ruptures in the study area are a miniature version of the entire rupture from the earthquake; they are geometrically and kinematically complex, with many individual and discontinuous segments of varying orientations and slip senses which are distributed across a zone up to ~3.5 km wide. Despite this variability, three main groups of ruptures have been identified. These are: 1) NE-SW striking, shallow to moderate dipping (25-45°W) faults that are approximately parallel to Cenozoic bedding with mainly reverse dip-slip and, and for the purposes of this thesis, are considered to be part of the SLF. 2) N-S striking, steeply dipping (~85°E) oblique sinistral faults that are up to the west and part of the SLF. 3) E-NE striking, moderate to steeply dipping (45-68°N) dextral reverse faults which are part of the HF. Bedding-parallel faults are interpreted to be flexural slip structures formed during folding of the near-surface Cenozoic strata, while the steeply dipping SLF ruptured a pre-existing bedrock fault which has little topographic expression. Groups 1 and 2 faults were both locally used for gravitational failure during the earthquake. Despite this non-tectonic fault movement, the slip vectors for faults that ruptured during the earthquake are broadly consistent with NCD tectonics and the regional ~100-120° trend of the principal horizontal stress/strain axes. Previous earthquake activity on the SLF is required by its displacement of Cenozoic formations but Late Quaternary slip on the fault prior to 2016 is neither supported by pre-existing fault scarps nor by changes in topography across the fault. By contrast, at least two earthquakes (including 2016) appear to have ruptured the HF from the mid Holocene, consistent with recurrence intervals of no more than ~7 kyr, and with preliminary observations from trenches on the fault farther to the west. The disparity in paleoearthquake records of the two faults suggests that they typically do not rupture together, thus it is concluded that the HF-SLF rupture pattern observed in the Kaikōura Earthquake rarely occurs in a single earthquake.
The objective of this study is to examine the influence of near-fault motions on liquefaction triggering in Christchurch and neighboring towns during the 2010-2011 Canterbury earthquake sequence (CES). The CES began with the 4 September 2010, Mw7.1 Darfield earthquake and included up to ten events that triggered liquefaction. However, most notably, widespread liquefaction was induced by the Darfield earthquake and the Mw6.2, 22 February 2011 Christchurch earthquake. Of particular relevance to this study is the forward directivity effects that were prevalent in the motions recorded during the Darfield earthquake, and to a much lesser extent, during the Christchurch earthquake. A 2D variant of the Richart-Newmark fatigue theory was used to compute the equivalent number of cycles (neq) for the ground motions, where volumetric strain was used as the damage metric. This study is unique because it considers the contribution and phasing of both the fault-normal and fault-parallel components of motion on neq and the magnitude scaling factor (MSF). It was found that when the fault-normal and fault-parallel motions were treated individually, the former yielded a lower neq than the latter. Additionally, when the combined effects of fault-normal and fault-parallel components were considered, it was found that the MSF were higher than those commonly used. This implies that motions containing near-fault effects are less demanding on the soil than motions that do not. This may be one of several factors that resulted in less severe liquefaction occurring during the Darfield earthquake than the Christchurch earthquake.
The Porters Pass fault (PPF) is a prominent element of the Porters Pass-Amberley Fault Zone (PPAFZ) which forms a broad zone of active earth deformation ca 100 km long, 60-90 km west and north of Christchurch. For a distance of ca 40 km the PPF is defined by a series of discontinuous Holocene active traces between the Rakaia and Waimakariri Rivers. The amount of slip/event and the timing of paleoearthquakes are crucial components needed to estimate the earthquake potential of a fault. Movement was assumed to be, coseismic and was quantified by measuring displaced geomorphic features using either tape measure or surveying equipment. Clustering of offset data suggests that four to five earthquakes occurred on the PPF during the Holocene and these range between ca 5-7 m/event. Timing information was obtained from four trenches excavated across the fault and an auger adjacent to the fault. Organic samples from these sites were radiocarbon dated and used in conjunction with data from previous studies to identify the occurrence of at least four earthquakes at 8500 ± 200, 5300 ± 700, 2500 ± 200 and 1000 ± 100 years B.P. Evidence suggests that an additional event is also possible at 6200 ± 500 years B.P. The ~1000, 5300 and 6200 years B.P. paleoearthquakes were previously unrecognised, while the 500 year event previously inferred from rock-avalanche data has been discarded. The present data set produces recurrence intervals of ~2000-2500 years for the Holocene. The identification of only one Holocene PPF rupture to the west of Red Lakes indicates the presence of a segment boundary that prevents the propagation of rupture beyond this point. This is consistent with displacement data and results in slip rates of 0.5-0.7 mm/yr and 2.5-3.4 mm/yr to the west and east of Red Lakes respectively. It is possible that the nearby extensional Red Hill Fault influences PPF rupture propagation. The combination of geometric, slip rate and timing data has enabled the magnitude of prehistoric earthquakes on the PPF to be estimated. These magnitudes range from an average of between 6.9 for a fault rupture from Waimakariri River to Red Lakes, to a maximum of 7.4 that ruptures the entire length of the PPAFZ, including the full length of the PPF. These estimates are approximately consistent with previous magnitude estimates along the full length of the PPAFZ of between 7.0 and 7.5.
A map of earthquake events in Canterbury.
An incomplete infographic showing a proposed fault rupture avoidance zone where the Alpine Fault runs under Franz Joseph township.
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.
The previously unknown Greendale Fault was buried beneath the Canterbury Plains and ruptured in the September 4th 2010 moment magnitude (Mw) 7.1 Darfield Earthquake. The Darfield Earthquake and subsequent Mw 6 or greater events that caused damage to Christchurch highlight the importance of unmapped faults near urban areas. This thesis examines the morphology, age and origin of the Canterbury Plains together with the paleoseismology and surface-rupture displacement distributions of the Greendale Fault. It offers new insights into the surface-rupture characteristics, paleoseismology and recurrence interval of the Greendale Fault and related structures involved in the 2010 Darfield Earthquake. To help constrain the timing of the penultimate event on the Greendale Fault the origin and age of the faulted glacial outwash deposits have been examined using sedimentological analysis of gravels and optically stimulated luminescence (OSL) dating combined with analysis of GPS and LiDAR survey data. OSL ages from this and other studies, and the analysis of surface paleochannel morphology and subsurface gravel deposits indicate distinct episodes of glacial outwash activity across the Canterbury Plains, at ~20 to 24 and ~28 to 33 kyr separated by a hiatus in sedimentation possibly indicating an interstadial period. These data suggest multiple glacial periods between ~18 and 35 kyr which may have occurred throughout the Canterbury region and wider New Zealand. A new model for the Waimakariri Fan is proposed where aggradation is mainly achieved during episodic sheet flooding with the primary river channel location remaining approximately fixed. The timing, recurrence interval and displacements of the penultimate surface-rupturing earthquake on the Greendale Fault have been constrained by trenching the scarp produced in 2010 at two locations. These excavations reveal a doubling of the magnitude of surface displacement at depths of 2-4 m. Aided by OSL ages of sand lenses in the gravel deposits, this factor-of-two increase is interpreted to indicate that in the central section of the Greendale Fault the penultimate surface-rupturing event occurred between ca. 20 and 30 kyr ago. The Greendale Fault remained undetected prior to the Darfield earthquake because the penultimate fault scarp was eroded and buried during Late Pleistocene alluvial activity. The Darfield earthquake rupture terminated against the Hororata Anticline Fault (HAF) in the west and resulted in up to 400 mm of uplift on the Hororata Anticline immediately above the HAF. Folding in 2010 is compared to Quaternary and younger deformation across the anticline recorded by a seismic reflection line, GPS-measured topographic profiles along fluvial surfaces, and river channel sinuosity and morphology. It is concluded that the HAF can rupture during earthquakes dissimilar to the 2010 event that may not be triggered by slip on the Greendale Fault. Like the Greendale Fault geomorphic analyses provide no evidence for rupture of the HAF in the last 18 kyr, with the average recurrence interval for the late Quaternary inferred to be at least ~10 kyr. Surface rupture of the Greendale Fault during the Darfield Earthquake produced one of the most accessible and best documented active fault displacement and geometry datasets in the world. Surface rupture fracture patterns and displacements along the fault were measured with high precision using real time kinematic (RTK) GPS, tape and compass, airborne light detection and ranging (LiDAR), and aerial photos. This allowed for detailed analysis of the cumulative strike-slip displacement across the fault zone, displacement gradient (ground shear strain) and the type of displacement (i.e. faulting or folding). These strain profiles confirm that the rupture zone is generally wide (~30 to ~300 metres) with >50% of displacement (often 70-80%) accommodated by ground flexure rather than discrete fault slip and ground cracking. The greatest fault-zone widths and highest proportions of folding are observed at fault stepovers.
A map showing the location of a study dating historic earthquakes on the Alpine Fault.
Object Overview of 'Hunter Hills Fault Zone study – Earthquake hazard assessment (Yetton, 2008).'