Mention the words "earthquake" in the same brief as "remediation" and it's enough to strike fear in the hearts of all New Zealanders, particularly those in Christchurch and other earthquake prone areas of the country. Now we find the chances of the ground shaking more violently in a quake is much higher than previously thought for large parts of the country. In some places it has doubled or even trebled. What are the ramifications of this new found knowledge? Joining the show to discuss is Michelle Grant, President of the Structural Engineering Society New Zealand, and Matt Gerstenberger, Principal Scientist and Seismologist at GNS Science
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’.
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.
