Beach ridge stratigraphy can provide an important record of both sustained coastal progradation and responses to events such as extreme storms, as well as evidence of earthquake induced sediment pulses. This study is a stratigraphic investigation of the late Holocene mixed sand gravel (MSG) beach ridge plain on the Canterbury coast, New Zealand. The subsurface was imaged along a 370 m shore-normal transect using 100 and 200 MHz ground penetrating radar (GPR) antennae, and cored to sample sediment textures. Results show that, seaward of a back-barrier lagoon, the Pegasus Bay beach ridge plain prograded almost uniformly, under conditions of relatively stable sea level. Nearshore sediment supply appears to have created a sustained sediment surplus, perhaps as a result of post-seismic sediment pulses, resulting in a flat, morphologically featureless beach ridge plain. Evidence of a high magnitude storm provides an exception, with an estimated event return period in excess of 100 years. Evidence from the GPR sequence combined with modern process observations from MSG beaches indicates that a paleo storm initially created a washover fan into the back-barrier lagoon, with a large amount of sediment simultaneously moved off the beach face into the nearshore. This erosion event resulted in a topographic depression still evident today. In the subsequent recovery period, sediment was reworked by swash onto the beach as a sequence of berm deposit laminations, creating an elevated beach ridge that also has a modern-day topographic signature. As sediment supply returned to normal, and under conditions of falling sea level, a beach ridge progradation sequence accumulated seaward of the storm feature out to the modern-day beach as a large flat, uniform progradation plain. This study highlights the importance of extreme storm events and earthquake pulses on MSG coastlines in triggering high volume beach ridge formation during the subsequent recovery period.
Liquefaction affects late Holocene, loose packed and water saturated sediment subjected to cyclical shear stress. Liquefaction features in the geological record are important off-fault markers that inform about the occurrence of moderate to large earthquakes (> 5 Mw). The study of contemporary liquefaction features provides a better understanding of where to find past (paleo) liquefaction features, which, if identified and dated, can provide information on the occurrence, magnitude and timing of past earthquakes. This is particularly important in areas with blind active faults.
The extensive liquefaction caused by the 2010-2011 Canterbury Earthquake Sequence (CES) gave the geoscience community the opportunity to study the liquefaction process in different settings (alluvial, coastal and estuarine), investigating different aspects (e.g. geospatial correlation with landforms, thresholds for peak ground acceleration, resilience of infrastructures), and to collect a wealth geospatial dataset in the broad region of the Canterbury Plains.
The research presented in this dissertation examines the sedimentary architecture of two environments, the alluvial and coastal settings, affected by liquefaction during the CES. The novel aim of this study is to investigate how landform and subsurface sedimentary architecture influence liquefaction and its surface manifestation, to provide knowledge for locating studies of paleoliquefaction in future.
Two study cases documented in the alluvial setting showed that liquefaction features affected a crevasse splay and point bar ridges. However, the liquefaction source layer was linked to paleochannel floor deposits below the crevasse splay in the first case, and to the point bar deposits themselves in the second case.
This research documents liquefaction features in the coastal dune system of the Canterbury Plains in detail for the first time. In the coastal dune setting the liquefiable layer is near the surface. The pore water pressure is vented easily because the coastal dune soil profile is entirely composed of non-cohesive, very well sorted sandy sediment that weakly resists disturbance from fluidised sediment under pressure. As a consequence, the liquefied flow does not need to find a specific crack through which the sediment is vented at the surface; instead, the liquefied sand finds many closely spaced conduits to vent its excess of pore water pressure. Therefore, in the coastal dune setting it is rare to observe discrete dikes (as they are defined in the alluvial setting), instead A horizon delamination (splitting) and blistering (near surface sills) are more common. The differences in styles of surface venting lead to contrasts in patterns of ejecta in the two environments. Whereas the alluvial environment is characterised by coalesced sand blows forming lineations, the coastal dune environment hosts apparently randomly distributed isolated sand blows often associated with collapse features.
Amongst the techniques tested for the first time to investigate liquefaction features are: 3D GPR, which improved the accuracy of the trenching even six years after the liquefaction events; thin section analysis to investigate sediment fabric, which helped to discriminate liquefied sediment from its host sediment, and modern from paleoliquefaction features; a Random Forest classification based on the CES liquefaction map, which was used to test relationships between surface manifestation of liquefaction and topographic parameters. The results from this research will be used to target new study sites for future paleoliquefaction research and thus will improve the earthquake hazard assessment across New Zealand.
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.
