Liquefaction is a phenomenon that results in a loss of strength and stability of a saturated soil mass due to dynamic excitation such as that imposed by an earthquake. The granular nature of New Zealand soils and the location of many of our cities and towns on fluvial foundations are such that the effects of liquefaction can be very important. Research was undertaken to build on the past work undertaken at the University of Canterbury studying the effects of the 1929 Murchison earthquake, the 1968 Inangahua earthquake and the 1991 Hawks Crag earthquakes on the West Coast. Additional archival information has been gathered from newspapers and reports and from discussions with people who experienced one or all of these large earthquakes that occurred on the West Coast during the 20th Century. Further, some twenty Cone Penetrometer Tests were carried out, with varying success, in Greymouth and Karamea using the Department of Civil Engineering's Drilling Rig. These, combined with the basic site investigation information, consolidate and add to the liquefaction case history data bank at the University of Canterbury. Many of the sites have liquefied in some but not all of the three earthquakes and thus provide both upper and lower bounds for the calibration of empirical models. While a lack of knowledge of the 1929 source location reduces the value of information from that event, the data form a useful set of liquefaction case histories and will become more so as further earthquakes occur. A list of critical sites for checking of the future earthquakes is provided and recommendations are made for the installation of downhole arrays of accelerometers and pore water pressure transducers at a number of sites.
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.
As a consequence of the 2010 – 2011 Canterbury earthquake sequence, Christchurch experienced widespread liquefaction, vertical settlement and lateral spreading. These geological processes caused extensive damage to both housing and infrastructure, and increased the need for geotechnical investigation substantially. Cone Penetration Testing (CPT) has become the most common method for liquefaction assessment in Christchurch, and issues have been identified with the soil behaviour type, liquefaction potential and vertical settlement estimates, particularly in the north-western suburbs of Christchurch where soils consist mostly of silts, clayey silts and silty clays. The CPT soil behaviour type often appears to over-estimate the fines content within a soil, while the liquefaction potential and vertical settlement are often calculated higher than those measured after the Canterbury earthquake sequence. To investigate these issues, laboratory work was carried out on three adjacent CPT/borehole pairs from the Groynes Park subdivision in northern Christchurch. Boreholes were logged according to NZGS standards, separated into stratigraphic layers, and laboratory tests were conducted on representative samples. Comparison of these results with the CPT soil behaviour types provided valuable information, where 62% of soils on average were specified by the CPT at the Groynes Park subdivision as finer than what was actually present, 20% of soils on average were specified as coarser than what was actually present, and only 18% of soils on average were correctly classified by the CPT. Hence the CPT soil behaviour type is not accurately describing the stratigraphic profile at the Groynes Park subdivision, and it is understood that this is also the case in much of northwest Christchurch where similar soils are found. The computer software CLiq, by GeoLogismiki, uses assessment parameter constants which are able to be adjusted with each CPT file, in an attempt to make each more accurate. These parameter changes can in some cases substantially alter the results for liquefaction analysis. The sensitivity of the overall assessment method, raising and lowering the water table, lowering the soil behaviour type index, Ic, liquefaction cutoff value, the layer detection option, and the weighting factor option, were analysed by comparison with a set of ‘base settings’. The investigation confirmed that liquefaction analysis results can be very sensitive to the parameters selected, and demonstrated the dependency of the soil behaviour type on the soil behaviour type index, as the tested assessment parameters made very little to no changes to the soil behaviour type plots. The soil behaviour type index, Ic, developed by Robertson and Wride (1998) has been used to define a soil’s behaviour type, which is defined according to a set of numerical boundaries. In addition to this, the liquefaction cutoff point is defined as Ic > 2.6, whereby it is assumed that any soils with an Ic value above this will not liquefy due to clay-like tendencies (Robertson and Wride, 1998). The method has been identified in this thesis as being potentially unsuitable for some areas of Christchurch as it was developed for mostly sandy soils. An alternative methodology involving adjustment of the Robertson and Wride (1998) soil behaviour type boundaries is proposed as follows: Ic < 1.31 – Gravelly sand to dense sand 1.31 < Ic < 1.90 – Sands: clean sand to silty sand 1.90 < Ic < 2.50 – Sand mixtures: silty sand to sandy silt 2.50 < Ic < 3.20 – Silt mixtures: clayey silt to silty clay 3.20 < Ic < 3.60 – Clays: silty clay to clay Ic > 3.60 – Organics soils: peats. When the soil behaviour type boundary changes were applied to 15 test sites throughout Christchurch, 67% showed an improved change of soil behaviour type, while the remaining 33% remained unchanged, because they consisted almost entirely of sand. Within these boundary changes, the liquefaction cutoff point was moved from Ic > 2.6 to Ic > 2.5 and altered the liquefaction potential and vertical settlement to more realistic ii values. This confirmed that the overall soil behaviour type boundary changes appear to solve both the soil behaviour type issues and reduce the overestimation of liquefaction potential and vertical settlement. This thesis acts as a starting point towards researching the issues discussed. In particular, future work which would be useful includes investigation of the CLiq assessment parameter adjustments, and those which would be most suitable for use in clay-rich soils such as those in Christchurch. In particular consideration of how the water table can be better assessed when perched layers of water exist, with the limitation that only one elevation can be entered into CLiq. Additionally, a useful investigation would be a comparison of the known liquefaction and settlements from the Canterbury earthquake sequence with the liquefaction and settlement potentials calculated in CLiq for equivalent shaking conditions. This would enable the difference between the two to be accurately defined, and a suitable adjustment applied. Finally, inconsistencies between the Laser-Sizer and Hydrometer should be investigated, as the Laser-Sizer under-estimated the fines content by up to one third of the Hydrometer values.
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.
Surface rupture of the previously unrecognised Greendale Fault extended west-east for ~30 km across alluvial plains west of Christchurch, New Zealand, during the Mw 7.1 Darfield (Canterbury) earthquake of September 2010. Surface rupture displacement was predominantly dextral strike-slip, averaging ~2.5 m, with maxima of ~5 m. Vertical displacement was generally less than 0.75 m. The surface rupture deformation zone ranged in width from ~30 to 300 m, and comprised discrete shears, localised bulges and, primarily, horizontal dextral flexure. About a dozen buildings, mainly single-storey houses and farm sheds, were affected by surface rupture, but none collapsed, largely because most of the buildings were relatively flexible and resilient timber-framed structures and also because deformation was distributed over a relatively wide zone. There were, however, notable differences in the respective performances of the buildings. Houses with only lightly-reinforced concrete slab foundations suffered moderate to severe structural and non-structural damage. Three other buildings performed more favourably: one had a robust concrete slab foundation, another had a shallow-seated pile foundation that isolated ground deformation from the superstructure, and the third had a structural system that enabled the house to tilt and rotate as a rigid body. Roads, power lines, underground pipes, and fences were also deformed by surface fault rupture and suffered damage commensurate with the type of feature, its orientation to the fault, and the amount, sense and width of surface rupture deformation.
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.
On 4 September 2010, a magnitude Mw 7.1 earthquake struck the Canterbury region on the South Island of New Zealand. The epicentre of the earthquake was located in the Darfield area about 40 km west of the city of Christchurch. Extensive damage was inflicted to lifelines and residential houses due to widespread liquefaction and lateral spreading in areas close to major streams, rivers and wetlands throughout Christchurch and Kaiapoi. Unreinforced masonry buildings also suffered extensive damage throughout the region. Despite the severe damage to infrastructure and residential houses, fortunately, no deaths occurred and only two injuries were reported in this earthquake. From an engineering viewpoint, one may argue that the most significant aspects of the 2010 Darfield Earthquake were geotechnical in nature, with liquefaction and lateral spreading being the principal culprits for the inflicted damage. Following the earthquake, an intensive geotechnical reconnaissance was conducted to capture evidence and perishable data from this event. This paper summarizes the observations and preliminary findings from this early reconnaissance work.
Improving community resilience requires a way of thinking about the nature of a community. Two complementary aspects are proposed: the flows connecting the community with its surrounding environment and the resources the community needs for its ongoing life. The body of necessary resources is complex, with many interactions between its elements. A systems approach is required to understand the issues adequately. Community resilience is discussed in general terms together with strategies for improving it. The ideas are then illustrated and amplified by an extended case study addressing means of improving the resilience of a community on the West Coast of New Zealand to natural disasters. The case study is in two phases. The first relies on a mix of on-the-ground observations and constructed scenarios to provide recommendations for enhancing community resilience, while the second complements the first by developing a set of general lessons and issues to be addressed from observations of the Christchurch earthquakes of 2010 and 2011.
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.
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.
On 4 September 2010, a magnitude Mw 7.1 earthquake struck the Canterbury region on the South Island of New Zealand. The epicentre of the earthquake was located in the Darfield area about 40 km west of the city of Christchurch. Extensive damage occurred to unreinforced masonry buildings throughout the region during the mainshock and subsequent large aftershocks. Particularly extensive damage was inflicted to lifelines and residential houses due to widespread liquefaction and lateral spreading in areas close to major streams, rivers and wetlands throughout Christchurch and Kaiapoi. Despite the severe damage to infrastructure and residential houses, fortunately, no deaths occurred and only two injuries were reported in this earthquake. From an engineering viewpoint, one may argue that the most significant aspects of the 2010 Darfield Earthquake were geotechnical in nature, with liquefaction and lateral spreading being the principal culprits for the inflicted damage. Following the earthquake, a geotechnical reconnaissance was conducted over a period of six days (10–15 September 2010) by a team of geotechnical/earthquake engineers and geologists from New Zealand and USA (GEER team: Geo-engineering Extreme Event Reconnaissance). JGS (Japanese Geotechnical Society) members from Japan also participated in the reconnaissance team from 13 to 15 September 2010. The NZ, GEER and JGS members worked as one team and shared resources, information and logistics in order to conduct thorough and most efficient reconnaissance covering a large area over a very limited time period. This report summarises the key evidence and findings from the reconnaissance.
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.
his poster presents the ongoing development of a 3D Canterbury seismic velocity model which will be used in physics-based hybrid broadband ground motion simulation of the 2010-2011 Canterbury earthquakes. Velocity models must sufficiently represent critical aspects of the crustal structure over multiple length scales which will influence the results of the simulations. As a result, numerous sources of data are utilized in order to provide adequate resolution where necessary. Figure 2: (a) Seismic reflection line showing P-wave velocities and significant geologic horizons (Barnes et al. 2011), and (b) Shear wave profiles at 10 locations (Stokoe et al. 2013). Figure 4: Cross sections of the current version of the Canterbury velocity model to depths of 10km as shown in Figure 1: (a) at a constant latitude value of -43.6˚, and (b) at a constant longitude value of 172.64˚. 3. Ground Surface and Geologic Horizon Models Figure 3: (a) Ground surface model derived from numerous available digital elevation models, and (b) Base of the Quaternary sediments derived from structural contours and seismic reflection line elevations. The Canterbury region has a unique and complex geology which likely has a significant impact on strong ground motions, in particular the deep and loose deposits of the Canterbury basin. The Canterbury basin has several implications on seismic wave phenomena such as long period ground motion amplification and wave guide effects. Using a realistic 3D seismic velocity model in physics-based ground motion simulation will implicitly account for such effects and the resultant simulated ground motions can be studied to gain a fundamental understanding of the salient ground motion phenomena which occurred during the Canterbury earthquakes, and the potential for repeat occurrences in the Canterbury region. Figure 1 shows the current model domain as a rectangular area between Lat=[-43.2˚,-44.0˚], and Lon=[171.5˚,173.0˚]. This essentially spans the area between the foot of the Southern Alps in the North West to Banks Peninsula in the East. Currently the model extends to a depth of 50km below sea level.
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’.
The Stone Jug Fault (SJF) ruptured during the November 14th, 2016 (at 12:02 am), Mw 7.8 Kaikōura Earthquake which initiated ~40 km west-southwest of the study area, at a depth of approximately 15 km. Preliminary post-earthquake mapping indicated that the SJF connects the Conway-Charwell and Hundalee faults, which form continuous surface rupture, however, detailed study of the SJF had not been undertaken prior to this thesis due to its remote location and mountainous topography. The SJF is 19 km long, has an average strike of ~160° and generally carries approximately equal components of sinistral and reverse displacement. The primary fault trace is sigmoidal in shape with the northern and southern tips rotating in strike from NNW to NW, as the SJF approaches the Hope and Hundalee faults. It comprises several steps and bends and is associated with many (N=48) secondary faults, which are commonly near irregularities in the main fault geometry and in a distributed fault zone at the southern tip. The SJF is generally parallel to Torlesse basement bedding where it may utilise pre-existing zones of weakness. Horizontal, vertical and net displacements range up to 1.4 m, with displacement profiles along the primary trace showing two main maxima separated by a minima towards the middle and ends of the fault. Average net displacement along the primary trace is ~0.4m, with local changes in relative values of horizontal and vertical displacement at least partly controlled by fault strike. Two trenches excavated across the northern segment of the fault revealed displacement of mainly Holocene stratigraphy dated using radiocarbon (N=2) and OSL (N=4) samples. Five surface-rupturing paleoearthquakes displaying vertical displacements of <1 m occurred at: 11,000±1000, 7500±1000, 6500±1000, 3500±100 and 3 (2016 Kaikōura) years BP. These events produce an average slip rate since ~11 ka of 0.2-0.4 mm/yr and recurrence intervals of up to 5500 years with an average recurrence interval of 2750 yrs. Comparison of these results with unpublished trench data suggests that synchronous rupture of the Hundalee, Stone Jug, Conway-Charwell, and Humps faults at ~3500 yrs BP cannot be discounted and it is possible that multi-fault ruptures in north Canterbury are more common than previously thought.
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.
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 Acheron rock avalanche is located in the Red Hill valley almost 80 km west of Christchurch and is one of 42 greywacke-derived rock avalanches identified in the central Southern Alps. It overlies the Holocene active Porters Pass Fault; a component of the Porters Pass-Amberley Fault Zone which extends from the Rakaia River to beyond the Waimakariri River. The Porters Pass Fault is a dextral strike-slip fault system viewed as a series of discontinuous fault scarps. The location of the fault trace beneath the deposit suggests it may represent a possible source of seismic shaking resulting in the formation of the Acheron rock avalanche. The rock mass composition of the rock avalanche source scar is Torlesse Supergroup greywacke consisting of massive sandstone and thinly bedded mudstone sequences dipping steeply north into the centre of the source basin. A stability analysis identified potential instability along shallow north dipping planar defects, and steep south dipping toppling failure planes. The interaction of the defects with bedding is considered to have formed conditions for potential instability most likely triggered by a seismic event. The dTositional area of the rock avalanche covers 7.2 x 105 m2 with an estimated volume of 9 x 10 m3 The mobilised rock mass volume was calculated at 7.5 x 106 m3• Run out of the debris from the top of the source scar to the distal limit reached 3500m, descending over a vertical fall of almost 700m with an estimated Fahrboschung of 0.2. The run out of the rock avalanche displayed moderate to high mobility, travelling at an estimated maximum velocity of 140-160 km/hour. The rapid emplacement of the deposit is confirmed by highly fragmented internal composition and burial of forest vegetation New radiocarbon ages from buried wood retrieved from the base of Acheron rock avalanche deposit represents an emplacement age closely post-dating (Wk 12094) 1152 ± 51 years B.P. This differs significantly from a previous radiocarbon age of (NZ547) 500 ± 69 years B.P. and modal lichenometry and weathering-rind thickness ages of approximately 460 ± 10 yrs and 490 ± 50 years B.P. The new age shows no resemblance to an earthquake event around 700- 500 years B.P. on the Porters Pass-Amberley Fault Zone. The DAN run out simulation using a friction model rheology successfully replicated the long run out and velocity of the Acheron rock avalanche using a frictron angle of 27° and high earth pressure coefficients of 5.5, 5.2, and 5.9. The elevated earth pressure coefficients represent dispersive pressures derived from dynamic fragmentation of the debris within the mobile rock avalanche, supporting the hypothesis of Davies and McSaveney (2002). The DAN model has potential applications for areas prone to large-scale instability in the elevated slopes and steep waterways of the Southern Alps. A paleoseismic investigation of a newly identified scarp of the Porters Pass Fault partially buried by the rock avalanche was conducted to identify any evidence of a coseismic relationship to the Acheron rock avalanche. This identified three-four fault traces striking at 078°, and a sag pond displaying a sequence of overbank deposits containing two buried soils representing an earthquake event horizon. A 40cm vertical offset of the ponded sediment and lower buried soil horizqn was recorded, which was dated to (Wk 13112 charcoal in palosol) 653 ± 54 years B.P. and (Wk 13034 palosol) 661 ± 34 years B.P. The evidence indicates a fault rupture occurred along the Porters Pass Fault, west of Porters Pass most likely extending to the Red Lakes terraces, post-dating 700 years B.P., resulting in 40cm of vertical displacement and an unknown component of dextral strike slip movement. This event post dates the event one (1000 ± 100 years B.P) at Porters Pass previously considered to represent the most recent rupture along the fault line. This points to a probable source for resetting of the modal weathering-rind thicknesses and lichen size populations in the Red Hill valley and possibly the Red Lakes terraces. These results suggest careful consideration must be given to the geomorphic and paleoseismic history of a specific site when applying surface dating techniques and furthermore the origin of dates used in literature and their useful range should be verified. An event at 700-500 years B.P did not trigger the Acheron rock avalanche as previously assumed supporting Howard's conclusions. The lack of similar aged rupture evidence in either of the Porters Pass and Coleridge trenches supports Howard's hypothesis of segmentation of the Porters Pass Fault; where rupture occurs along one fault segment but not along another. The new rock avalanche age closely post-dating 1200-1100 years B.P. resembles the poorly constrained event one rupture age of 1700-800 years B.P for the Porters Pass Fault and the tighter constrained Round Top event of 1010 ± 50 years B.P. on the Alpine Fault. Eight other rock avalanche deposits spread across the central Southern Alps also resemble the new ages however are unable to be assigned specific earthquake events due to the large associated error bars of± 270 years. This clustering of ages does represent compelling lines of evidence for large magnitude earthquake events occurring over the central Southern Alps. The presence of a rock avalanche deposit does not signify an earthquake based on the historical evidence in the Southern Alps however clustering of ages does suggest that large Mw >7 earthquakes occurred across the Southern Alps between 1200-900 years BP.
This dissertation addresses a diverse range of applied aspects in ground motion simulation validation via the response of complex structures. In particular, the following topics are addressed: (i) the investigation of similarity between recorded and simulated ground motions using code-based 3D irregular structural response analysis, (ii) the development of a framework for ground motion simulations validation to identify the cause of differences between paired observed and simulated dataset, and (iii) the illustration of the process of using simulations for seismic performance-based assessment. The application of simulated ground motions is evaluated for utilisation in engineering practice by considering responses of 3D irregular structures. Validation is performed in a code-based context when the NZS1170.5 (NZS1170.5:2004, 2004) provisions are followed for response history analysis. Two real buildings designed by engineers and physically constructed in Christchurch before the 2010-2011 Canterbury earthquake sequence are considered. The responses are compared when the buildings are subjected to 40 scaled recorded and their subsequent simulated ground motions selected from 22 February 2011 Christchurch. The similarity of recorded and simulated responses is examined using statistical methods such as bootstrapping and hypothesis testing to determine whether the differences are statistically significant. The findings demonstrate the applicability of simulated ground motion when the code-based approach is followed in response history analysis. A conceptual framework is developed to link the differences between the structural response subjected to simulated and recorded ground motions to the differences in their corresponding intensity measures. This framework allows the variability to be partitioned into the proportion that can be “explained” by the differences in ground motion intensity measures and the remaining “unexplained” variability that can be attributed to different complexities such as dynamic phasing of multi-mode response, nonlinearity, and torsion. The application of this framework is examined through a hierarchy of structures reflecting a range of complexity from single-degree-of-freedom to 3D multi-degree-of-freedom systems with different materials, dynamic properties, and structural systems. The study results suggest the areas that ground motion simulation should focus on to improve simulations by prioritising the ground motion intensity measures that most clearly account for the discrepancies in simple to complex structural responses. Three approaches are presented to consider recorded or simulated ground motions within the seismic performance-based assessment framework. Considering the applications of ground motions in hazard and response history analyses, different pathways in utilising ground motions in both areas are explored. Recorded ground motions are drawn from a global database (i.e., NGA-West2 Ancheta et al., 2014). The NZ CyberShake dataset is used to obtain simulations. Advanced ground motion selection techniques (i.e., generalized conditional intensity measure, GCIM) are used for ground motion selection at a few intensity levels. The comparison is performed by investigating the response of an example structure (i.e., 12-storey reinforced concrete special moment frame) located in South Island, NZ. Results are compared and contrasted in terms of hazard, groundmotion selection, structural responses, demand hazard, and collapse risk, then, the probable reasons for differences are discussed. The findings from this study highlight the present opportunities and shortcomings in using simulations in risk assessment. i
In the period between September 2010 and December 2011, Christchurch (New Zealand) and its surroundings were hit by a series of strong earthquakes including six significant events, all generated by local faults in proximity to the city: 4 September 2010 (Mw=7.1), 22 February 2011 (Mw=6.2), 13 June 2011 (Mw=5.3 and Mw=6.0) and 23 December 2011 (M=5.8 and (M=5.9) earthquakes. As shown in Figure 1, the causative faults of the earthquakes were very close to or within the city boundaries thus generating very strong ground motions and causing tremendous damage throughout the city. Christchurch is shown as a lighter colour area, and its Central Business District (CBD) is marked with a white square area in the figure. Note that the sequence of earthquakes started to the west of the city and then propagated to the south, south-east and east of the city through a set of separate but apparently interacting faults. Because of their strength and proximity to the city, the earthquakes caused tremendous physical damage and impacts on the people, natural and built environments of Christchurch. The 22 February 2011 earthquake was particularly devastating. The ground motions generated by this earthquake were intense and in many parts of Christchurch substantially above the ground motions used to design the buildings in Christchurch. The earthquake caused 182 fatalities, collapse of two multi-storey reinforced concrete buildings, collapse or partial collapse of many unreinforced masonry structures including the historic Christchurch Cathedral. The Central Business District (CBD) of Christchurch, which is the central heart of the city just east of Hagley Park, was practically lost with majority of its 3,000 buildings being damaged beyond repair. Widespread liquefaction in the suburbs of Christchurch, as well as rock falls and slope/cliff instabilities in the Port Hills affected tens of thousands of residential buildings and properties, and shattered the lifelines and infrastructure over approximately one third of the city area. The total economic loss caused by the 2010-2011 Christchurch earthquakes is currently estimated to be in the range between 25 and 30 billion NZ dollars (or 15% to 18% of New Zealand’s GDP). After each major earthquake, comprehensive field investigations and inspections were conducted to document the liquefaction-induced land damage, lateral spreading displacements and their impacts on buildings and infrastructure. In addition, the ground motions produced by the earthquakes were recorded by approximately 15 strong motion stations within (close to) the city boundaries providing and impressive wealth of data, records and observations of the performance of ground and various types of structures during this unusual sequence of strong local earthquakes affecting a city. This paper discusses the liquefaction in residential areas and focuses on its impacts on dwellings (residential houses) and potable water system in the Christchurch suburbs. The ground conditions of Christchurch including the depositional history of soils, their composition, age and groundwater regime are first discussed. Detailed liquefaction maps illustrating the extent and severity of liquefaction across Christchurch triggered by the sequence of earthquakes including multiple episodes of severe re-liquefaction are next presented. Characteristic liquefaction-induced damage to residential houses is then described focussing on the performance of typical house foundations in areas affected by liquefaction. Liquefaction impacts on the potable water system of Christchurch is also briefly summarized including correlation between the damage to the system, liquefaction severity, and the performance of different pipe materials. Finally, the characteristics of Christchurch liquefaction and its impacts on built environment are discussed in relation to the liquefaction-induced damage in Japan during the 11 March 2011 Great East Japan Earthquake.
Structures of the Lowry Peaks Range - Waikari Valley district are complex. The majority comprise three members of a predominantly WSW -ENE striking major northwards-directed, leading edge imbricate thrust system, with associated angular, asymmetric fault-propagation folds. This system forms anomalously within a large NESW trending belt of structures characterising the entire east coast of north Canterbury, both onshore and offshore and terminates westwards against N-S striking, east facing fold-fault zone. The objectives of this study address the origin, geometry and kinematics of the interaction between these diversely trending systems. Stratigraphy and small-scale structures denote three periods of deformation, namely: i) Middle Cretaceous deformation of the basement rocks, ii) weak Middle Oligocene deformation associated with the inception of the plate boundary through the South Island, and iii) major Pliocene - Recent deformation that formed the majority of the above-mentioned structures. Stress tensor analyses within competent basement and limestone cover rocks suggest two sets of sub-horizontal compression, NE-SW and NW-SE, the former likely to relate to a localised earlier period of deformation, now overprinted by the latter. NW-SE oriented sub-horizontal compression correlates well with results from other parts of north Canterbury. The result of NW-SE compression on the W-E to WSW-ENE striking structures is a large component of oblique motion, which is manifest in four ways: i) movement on two, differently oriented splays rather than a single fault strand, ii) the development of a sinuous trace for a number of the major folds, whereby the ends are oriented normal to the compression direction, the centres parallel to the strike of the faults, iii) the development of a number of cross-folds, striking NNE-SSW and iv) the apparently recent development of a strike-slip component on at least one of the major thrust faults. The origin of the W-E, or WSW-ENE striking structures may be reactivation of Late Cretaceous faults, stratigraphic evidence for the existence of a "structural high" (the Hurunui High) over the majority of the area in the Late Cretaceous to Early Eocene times suggests the formation of a W-E trending horst structure, with a corresponding asymmetric graben to the south. The junction of WSW-ENE trending structures with N-S trending structures to the west centres on an alluvial-filled depression, Waikari Flat, into which the structures of the WSW-ENE trending imbricate thrust system plunge, locally curling to the SW at their ends to link with N-S trending structures to the south. Roof thrusting on two orientations, W-E and N-S, towards to SE is currently occurring above these structures. Currently the area is not highly seismically active, although a magnitude ~6.4 Ms earthquake in historic times has been recorded. The effects of tectonics on the drainage of the area does suggest that the majority of the systems, are still potentially active, albeit moving at a comparatively slow rate. The majority of the recent motion appears to be concentrated on the roof-thrusting occurring in Waikari Flat, and uplift along the Lowry Peaks Fault System. Increasing amounts of secondary movement on back-thrusts and cross fractures is also implied for western ends of the major imbricate thrust system. In contrast, the southern-most fault system appears to be largely sustaining dextral strike-slip motion, with some local folding in central portions.
The structure and geomorphology of active orogens evolves on time scales ranging from a single earthquake to millions of years of tectonic deformation. Analysis of crustal deformation using new and established remote sensing techniques, and integration of these data with field mapping, geochronology and the sedimentary record, create new opportunities to understand orogenic evolution over these timescales. Timor Leste (East Timor) lies on the northern collisional boundary between continental crust from the Australian Plate and the Banda volcanic arc. GPS studies have indicated that the island of Timor is actively shortening. Field mapping and fault kinematic analysis of an emergent Pliocene marine sequence identifies gentle folding, overprinted by a predominance of NW-SE oriented dextral-normal faults and NE-SW oriented sinistral-normal faults that collectively bound large (5-20km2) bedrock massifs throughout the island. These fault systems intersect at non-Andersonian conjugate angles of approximately 120° and accommodate an estimated 20 km of orogen-parallel extension. Folding of Pliocene rocks in Timor may represent an early episode of contraction but the overall pattern of deformation is one of lateral crustal extrusion sub-parallel to the Banda Arc. Stratigraphic relationships suggest that extrusion began prior to 5.5 Ma, during and after initial uplift of the orogen. Sedimentological, geochemical and Nd isotope data indicate that the island of Timor was emergent and shedding terrigenous sediment into carbonate basins prior to 4.5 Ma. Synorogenic tectonic and sedimentary phases initiated almost synchronously across much of Timor Leste and <2 Myr before similar events in West Timor. An increase in plate coupling along this obliquely converging boundary, due to subduction of an outlying continental plateau at the Banda Trench, is proposed as a mechanism for uplift that accounts for orogen-parallel extension and early uplift of Timor Leste. Rapid bathymetric changes around Timor are likely to have played an important role in evolution of the Indonesian Seaway. The 2010 Mw 7.1 Darfield (Canterbury) earthquake in New Zealand was complex, involving multiple faults with strike-slip, reverse and normal displacements. Multi-temporal cadastral surveying and airborne light detection and ranging (LiDAR) surveys allowed surface deformation at the junction of three faults to be analyzed in this study in unprecedented detail. A nested, localized restraining stepover with contractional bulging was identified in an area with the overall fault structure of a releasing bend, highlighting the surface complexities that may develop in fault interaction zones during a single earthquake sequence. The earthquake also caused river avulsion and flooding in this area. Geomorphic investigations of these rivers prior to the earthquake identify plausible precursory patterns, including channel migration and narrowing. Comparison of the pre and post-earthquake geomorphology of the fault rupture also suggests that a subtle scarp or groove was present along much of the trace prior to the Darfield earthquake. Hydrogeology and well logs support a hypothesis of extended slip history and suggests that that the Selwyn River fan may be infilling a graben that has accumulated late Quaternary vertical slip of <30 m. Investigating fault behavior, geomorphic and sedimentary responses over a multitude of time-scales and at different study sites provides insights into fault interactions and orogenesis during single earthquakes and over millions of years of plate boundary deformation.
This dissertation addresses a diverse range of topics in the physics-based broadband ground motion simulation, with a focus on New Zealand applications. In particular the following topics are addressed: the methodology and computational implementation of a New Zealand Velocity Model for broadband ground motion simulation; generalised parametric functions and spatial correlations for seismic velocities in the Canterbury, New Zealand region from surface-wave-based site characterisation; and ground motion simulations of Hope Fault earthquakes. The paragraphs below outline each contribution in more detail. A necessary component in physics-based ground motion simulation is a 3D model which details the seismic velocities in the region of interest. Here a velocity model construction methodology, its computational implementation, and application in the construction of a New Zealand velocity model for use in physics-based broadband ground motion simulation are presented. The methodology utilises multiple datasets spanning different length scales, which is enabled via the use of modular sub-regions, geologic surfaces, and parametric representations of crustal velocity. A number of efficiency-related workflows to decrease the overall computational construction time are employed, while maintaining the flexibility and extensibility to incorporate additional datasets and re- fined velocity parameterizations as they become available. The model comprises explicit representations of the Canterbury, Wellington, Nelson-Tasman, Kaikoura, Marlborough, Waiau, Hanmer and Cheviot sedimentary basins embedded within a regional travel-time tomography-based velocity model for the shallow crust and provides the means to conduct ground motion simulations throughout New Zealand for the first time. Recently developed deep shear-wave velocity profiles in Canterbury enabled models that better characterise the velocity structure within geologic layers of the Canterbury sedimentary basin to be developed. Here the development of depth- and Vs30-dependent para-metric velocity and spatial correlation models to characterise shear-wave velocities within the geologic layers of the Canterbury sedimentary basin are presented. The models utilise data from 22 shear-wave velocity profiles of up to 2.5km depth (derived from surface wave analysis) juxtaposed with models which detail the three-dimensional structure of the geologic formations in the Canterbury sedimentary basin. Parametric velocity equations are presented for Fine Grained Sediments, Gravels, and Tertiary layer groupings. Spatial correlations were developed and applied to generate three-dimensional stochastic velocity perturbations. Collectively, these models enable seismic velocities to be realistically represented for applications such as 3D ground motion and site response simulations. Lastly the New Zealand velocity model is applied to simulate ground motions for a Mw7.51 rupture of the Hope Fault using a physics-based simulation methodology and a 3D crustal velocity model of New Zealand. The simulation methodology was validated for use in the region through comparison with observations for a suite of historic small magnitude earthquakes located proximal to the Hope Fault. Simulations are compared with conventionally utilised empirical ground motion models, with simulated peak ground velocities being notably higher in regions with modelled sedimentary basins. A sensitivity analysis was undertaken where the source characteristics of magnitude, stress parameter, hypocentre location and kinematic slip distribution were varied and an analysis of their effect on ground motion intensities is presented. It was found that the magnitude and stress parameter strongly influenced long and short period ground motion amplitudes, respectively. Ground motion intensities for the Hope Fault scenario are compared with the 2016 Kaikoura Mw7.8 earthquake, it was found that the Kaikoura earthquake produced stronger motions along the eastern South Island, while the Hope Fault scenario resulted in stronger motions immediately West of the near-fault region. The simulated ground motions for this scenario complement prior empirically-based estimates and are informative for mitigation and emergency planning purposes.
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.
The purpose of this thesis is to conduct a detailed examination of the forward-directivity characteristics of near-fault ground motions produced in the 2010-11 Canterbury earthquakes, including evaluating the efficacy of several existing empirical models which form the basis of frameworks for considering directivity in seismic hazard assessment. A wavelet-based pulse classification algorithm developed by Baker (2007) is firstly used to identify and characterise ground motions which demonstrate evidence of forward-directivity effects from significant events in the Canterbury earthquake sequence. The algorithm fails to classify a large number of ground motions which clearly exhibit an early-arriving directivity pulse due to: (i) incorrect pulse extraction resulting from the presence of pulse-like features caused by other physical phenomena; and (ii) inadequacy of the pulse indicator score used to carry out binary pulse-like/non-pulse-like classification. An alternative ‘manual’ approach is proposed to ensure 'correct' pulse extraction and the classification process is also guided by examination of the horizontal velocity trajectory plots and source-to-site geometry. Based on the above analysis, 59 pulse-like ground motions are identified from the Canterbury earthquakes , which in the author's opinion, are caused by forward-directivity effects. The pulses are also characterised in terms of their period and amplitude. A revised version of the B07 algorithm developed by Shahi (2013) is also subsequently utilised but without observing any notable improvement in the pulse classification results. A series of three chapters are dedicated to assess the predictive capabilities of empirical models to predict the: (i) probability of pulse occurrence; (ii) response spectrum amplification caused by the directivity pulse; (iii) period and amplitude (peak ground velocity, PGV) of the directivity pulse using observations from four significant events in the Canterbury earthquakes. Based on the results of logistic regression analysis, it is found that the pulse probability model of Shahi (2013) provides the most improved predictions in comparison to its predecessors. Pulse probability contour maps are developed to scrutinise observations of pulses/non-pulses with predicted probabilities. A direct comparison of the observed and predicted directivity amplification of acceleration response spectra reveals the inadequacy of broadband directivity models, which form the basis of the near-fault factor in the New Zealand loadings standard, NZS1170.5:2004. In contrast, a recently developed narrowband model by Shahi & Baker (2011) provides significantly improved predictions by amplifying the response spectra within a small range of periods. The significant positive bias demonstrated by the residuals associated with all models at longer vibration periods (in the Mw7.1 Darfield and Mw6.2 Christchurch earthquakes) is likely due to the influence of basin-induced surface waves and non-linear soil response. Empirical models for the pulse period notably under-predict observations from the Darfield and Christchurch earthquakes, inferred as being a result of both the effect of nonlinear site response and influence of the Canterbury basin. In contrast, observed pulse periods from the smaller magnitude June (Mw6.0) and December (Mw5.9) 2011 earthquakes are in good agreement with predictions. Models for the pulse amplitude generally provide accurate estimates of the observations at source-to-site distances between 1 km and 10 km. At longer distances, observed PGVs are significantly under-predicted due to their slower apparent attenuation. Mixed-effects regression is employed to develop revised models for both parameters using the latest NGA-West2 pulse-like ground motion database. A pulse period relationship which accounts for the effect of faulting mechanism using rake angle as a continuous predictor variable is developed. The use of a larger database in model development, however does not result in improved predictions of pulse period for the Darfield and Christchurch earthquakes. In contrast, the revised model for PGV provides a more appropriate attenuation of the pulse amplitude with distance, and does not exhibit the bias associated with previous models. Finally, the effects of near-fault directivity are explicitly included in NZ-specific probabilistic seismic hazard analysis (PSHA) using the narrowband directivity model of Shahi & Baker (2011). Seismic hazard analyses are conducted with and without considering directivity for typical sites in Christchurch and Otira. The inadequacy of the near-fault factor in the NZS1170.5: 2004 is apparent based on a comparison with the directivity amplification obtained from PSHA.
The Lake Coleridge Rock Avalanche Deposits (LCRADs) are located on Ryton Station in the middle Rakaia Valley, approximately 80 km west of Christchurch. Torlesse Supergroup greywacke is the basement material and has been significantly influenced by both active tectonics and glaciation. Both glacial and post-glacial processes have produced large volumes of material which blanket the bedrock on slopes and in the valley floors. The LCRADs were part of a regional study of rock avalanches by WHITEHOUSE (1981, 1983) and WHITEHOUSE and GRIFFITHS (1983), and a single rock avalanche event was recognised with a weathering rind age of 120 years B.P. that was later modified to 150 ± 40 years B.P. The present study has refined details of both the age and the sequence of events at the site, by identifying three separate rock avalanche deposits (termed the LCRA1, LCRA2 and LCRA3 deposits), which are all sourced from near the summit of Carriage Drive. The LCRA1 deposit is lobate in shape and had an estimated original deposit volume of 12.5 x 10⁶ m³, although erosion by the Ryton River has reduced the present day debris volume to 5.1 x 10⁶ m³. An optically stimulated luminescence date taken from sandy loess immediately beneath the LCRA1 deposit provided a maximum age for the rock avalanche event of 9,720 ± 750 years B.P., which is believed to be realistic given that this is shortly after the retreat of Acheron 3 ice from this part of the valley. Emplacement of rock avalanche material into an ancestral Ryton riverbed created a natural dam with a ~17 M m³ lake upstream. The river is thought to have created a natural spillway over the dam structure at ~557 m (a.s.l), and to have existed for a number of years before any significant downcutting occurred. Although a triggering mechanism for the LCRA1 deposit was poorly constrained, it is thought that stress rebound after glacial ice removal may have initiated failure. Due to the event occurring c.10,000 years ago, there was a lack of definition for a possible earthquake trigger, though the possibility is obvious. The LCRA₂ event had an original deposit volume of 0.66 x 10⁶ m³, and was constrained to the low-lying area adjacent to the Ryton River that had been created by river erosion of the LCRA1 deposit. Further erosion by the Ryton River has reduced the deposit volume to 0.4 x 10⁶ m³. A radiocarbon date from a piece of mānuka found within the LCRA2 deposit provided an age of 668 ± 36 years B.P., and this is thought to reliably date the event. The LCRA2 event also dammed the Ryton River, and the preservation of dam-break outwash terraces downstream from the deposit provides clear evidence of rapid dam erosion and flooding after overtopping, and breaching by the Ryton River. Based on the mean annual flow of the Ryton River, the LCRA2 lake would have taken approximately two weeks to fill assuming that there were no preferred breach paths and the material was relatively impermeable. The LCRA2 event is thought to have been coseismic with a fault rupture along the western segment of the PPAFZ, which has been dated at 600 ± 100 years B.P. by SMITH (2003). The small LCRA3 event was not able to be dated, but it is believed to have failed shortly after the LCRA2 event and it may in fact be a lag deposit of the second rock avalanche event possibly triggered by an aftershock. The deposit is only visible at one locality within the cliffs that line the Ryton River, and its lack of geomorphic expression is attributed to it occurring closely after the LCRA2 event, while the Ryton River was still dammed from the second rock avalanche event. A wedge-block of some 35,000 m³ of source material for a future rock avalanche was identified at the summit of Carriage Drive. The dilation of the rock mass, combined with unfavourably oriented sub-vertical bedding in the Torlesse Supergroup bedrock, has allowed toppling-style failure on both of the main ridge lines around the source area for the LCRADs. In the event of a future rock avalanche occurring within the Ryton riverbed an emergency response plan has been developed to provide a staged response, especially in relation to the camping ground located at the mouth of the Ryton River. A long-term management plan has also been developed for mitigation measures for the Ryton riverbed and adjacent floodplain areas downstream of a future rock avalanche at the LCRAD site.
Documenting earthquake-induced ground deformation is significant to assess the characteristics of past and contemporary earthquakes and provide insight into seismic hazard. This study uses airborne light detection and ranging (LiDAR) and conducts multi-disciplinary field techniques to document the surface rupture morphology and evaluate the paleoseismicity and seismic hazard parameters of the Hurunui segment of the Hope Fault in the northern South Island of New Zealand. It also documents and evaluates seismically induced features and ground motion characteristics of the 2010 Darfield and 2011 Christchurch earthquakes in the Port Hills, south of Christchurch. These two studies are linked in that they investigate the near-field coseismic features of large (Mw ~7.1) earthquakes in New Zealand and produce data for evaluating seismic hazards of future earthquakes. In the northern South Island of New Zealand, the Australian-Pacific plate boundary is characterised by strike-slip deformation across the Marlborough Fault System (MFS). The ENE-striking Hope Fault (length: ~230 km) is the youngest and southernmost fault in the MFS, and the second fastest slipping fault in New Zealand. The Hope Fault is a major source of seismic hazard in New Zealand and has ruptured (in-part) historically in the Mw 7.1 1888 Amuri earthquake. In the west, the Hurunui segment of the Hope Fault is covered by beech forest. Hence, its seismic hazard parameters and paleoearthquake chronology were poorly constrained and it was unknown whether the 1888 earthquake ruptured this segment or not and if so, to what extent. Utilising LiDAR and field data, a 29 km-long section of the Hurunui segment of the Hope Fault is mapped. LiDAR-mapping clearly reveals the principal slip zone (PSZ) of the fault and a suite of previously unrecognised structures that form the fault deformation zone (FDZ). FDZ width measurements from 415 locations reveal a spatially-variable, active FDZ up to ~500 m wide with an average width of 200 m. Kinematic analysis of the fault structures shows that the Hurunui segment strikes between 070° and 075° and is optimally oriented for dextral strike-slip within the regional stress field. This implies that the wide FDZ observed is unlikely to result from large-scale fault mis-orientation with respect to regional stresses. The analysis of FDZ width indicates that it increases with increased hanging wall topography and increased topographic relief suggesting that along-strike topographic perturbations to fault geometry and stress states increase fault zone complexity and width. FDZ width also increases where the tips of adjacent PSZ strands locally vary in strike, and where the thickness of alluvial deposits overlying bedrock increases. LiDAR- and photogrammetrically-derived topographic mapping indicates that the boundary between the Hurunui and Hope River segments is characterised by a ~850-m-wide right stepover and a 9º-14° fault bend. Paleoseismic trenching at Hope Shelter site reveals that 6 earthquakes occurred at A.D. 1888, 1740-1840, 1479-1623, 819-1092, 439-551, and 373- 419. These rupture events have a mean recurrence interval of ~298 ± 88 yr and inter-event times ranging from 98 to 595 yrs. The variation in the inter-event times is explained by (1) coalescing rupture overlap from the adjacent Hope River segment on to the Hurunui segment at the study site, (2) temporal clustering of large earthquakes on the Hurunui segment, and/or (3) ‘missing’ rupture events. It appears that the first two options are more plausible to explain the earthquake chronologies and rupture behaviour on the Hurunui segment, given the detailed nature of the geologic and chronologic investigations. This study provides first evidence for coseismic multi-segment ruptures on the Hope Fault by identifying a rupture length of 44-70 km for the 1888 earthquake, which was not confined to the Hope River segment (primary source for the 1888 earthquake). LiDAR data is also used to identify and measure dextral displacements and scarp heights from the PSZ and structures within the FDZ along the Hurunui segment. Reconstruction of large dextrally-offset geomorphic features shows that the vertical component of slip accounts for only ~1% of the horizontal displacements and confirms that the fault is predominantly strike-slip. A strong correlation exists between the dextral displacements and elevations of geomorphic features suggesting the possibility of age correlation between the geomorphic features. A mean single event displacement (SED) of 3.6 ± 0.7 m is determined from interpretation of sets of dextral displacements of ≤ 25 m. Using the available surface age data and the cumulative dextral displacements from Matagouri Flat, McKenzie Fan, Macs Knob and Hope River sites, and the mean SED, a mean slip rate of 12.2 ± 2.4 mm/yr, and a mean recurrence interval of ~320 ± 120 yr, and a potential earthquake magnitude of Mw 7.2 are determined for the Hurunui segment. This study suggests that the fault slip rate has been constant over the last ~15000 yr. Strong ground motions from the 2010 Darfield (Canterbury) earthquake displaced boulders and caused ground damage on some ridge crests in the Port Hills. However, the 2011 Christchurch earthquake neither displaced boulders nor caused ground damage at the same ridge crests. Documentation of locations (~400 m a.s.l.), lateral displacements (8-970 cm), displacement direction (250° ± 20°) of displaced boulders, in addition to their hosting socket geometries (< 1 cm to 50 cm depth), the orientation of the ridges (000°-015°) indicate that boulders have been displaced in the direction of instrumentally recorded transient peak ground horizontal displacements nearby and that the seismic waves have been amplified at the study sites. The co-existence of displaced and non-displaced boulders at proximal sites suggests small-scale ground motion variability and/or varying boulder-ground dynamic interactions relating to shallow phenomena such as variability in soil depth, bedrock fracture density and/or microtopography on the bedrock-soil interface. Shorter shaking duration of the 2011 Christchurch event, differing frequency contents and different source characteristics were all factors that may have contributed to generating circumstances less favourable to boulder displacement in this earthquake. Investigating seismically induced features, fault behaviour, site effects on the rupture behaviour, and site response to the seismic waves provides insights into fault rupture hazards.
