The Canterbury region experienced widespread damage due to liquefaction induced by seismic shaking during the 4 September 2010 earthquake and the large aftershocks that followed, notably those that occurred on 22 February, 13 June and 23 December 2011. Following the 2010 earthquake, the Earthquake Commission directed a thorough investigation of the ground profile in Christchurch, and to date, more than 7500 cone penetration tests (CPT) have been performed in the region. This paper presents the results of analyses which use a subset of the geotechnical database to evaluate the liquefaction process as well as the re-liquefaction that occurred following some of the major events in Christchurch. First, the applicability of existing CPT-based methods for evaluating liquefaction potential of Christchurch soils was investigated using three methods currently available. Next, the results of liquefaction potential evaluation were compared with the severity of observed damage, categorised in terms of the land damage grade developed from Tonkin & Taylor property inspections as well as from observed severity of liquefaction from aerial photography. For this purpose, the Liquefaction Potential Index (LPI) was used to represent the damage potential at each site. In addition, a comparison of the CPT-based strength profiles obtained before each of the major aftershocks was performed. The results suggest that the analysis of spatial and temporal variations of strength profiles gives a clear indication of the resulting liquefaction and re-liquefaction observed in Christchurch. The comparison of a limited number of CPT strength profiles before and after the earthquakes seems to indicate that no noticeable strengthening has occurred in Christchurch, making the area vulnerable to liquefaction induced land damage in future large-scale earthquakes.
The North Canterbury and Marlborough regions of Aotearoa | New Zealand were severely impacted by almost 30,000 landslides triggered during the 2016 Kaikōura Earthquake. Of these landslides approximately 200 dammed rivers. In the study area near Waiau, rupture of The Humps and Leader faults (and associated ground motions) initiated at least 42 co-seismic landslides. The Leader Landslide is the largest of these landslides, with an area of approximately 600,000 m2 and a volume of 6-8 million m3. The landslide buried approximately 980 m of active Leader River bed length and dammed the river. The dam produced four lakes, with two remaining today and two having been breached by partial landslide collapse and knickpoint migration in the year following the earthquake. As of 2025, the landslide dam has not been completely breached and Lake Rebekah remains. The Leader Landslide dam presents a unique opportunity to chart the evolution of the active riverbed pre- and post-earthquake, for up to 2 km downstream of Lake Rebekah. The river’s evolutionary timeline was observed using LiDAR, satellite aerial imagery, and drone surveys from 2001 to 2024 to develop maps and topographic difference models. Key timeframes for riverbed change events were also constrained with information and dated photography gathered from previous communications with the landowners at Woodchester Station, where the landslide is located. Finally, Schmidt Hammer testing of the Pliocene-Miocene Greta Siltstone Formation was conducted to investigate the role of bedrock strength on the rate of riverbed erosion. I present the history of evolution of the Leader River, pre- and post-earthquake, and consider factors impacting riverbed morphology changes. Despite the stability of Lake Rebekah, these data show that the position and morphology of the Leader River has changed significantly in response to the landslide, with the formation of two knickpoint waterfalls up to 14 m-high, four waterbodies, and diversion of the river around the landslide toe. Evolution of the river is characterised by longer periods of stasis (e.g., months to years) punctuated by rapid changes in riverbed morphology (e.g., hours to weeks) associated with incision and aggradation. In particular, the knickpoints migrated upstream at variable spatial and temporal rates. Factors controlling the rates of processes include; rain-storm events, partial lake outburst flooding, spatial changes in Pliocene-Miocene siltstone bed induration and landowner intervention to stabilise the landslide dam. An overarching conclusion of this thesis is that landforms can develop rapidly (i.e., hours to weeks) and in the absence of historical accounts, could be interpreted to have formed over hundreds to thousands of years.
Spatial variations in river facies exerted a strong influence on the distribution of liquefaction features observed in Christchurch during the 2010-11 Canterbury Earthquake Sequence (CES). Liquefaction and liquefaction-induced ground deformation was primarily concentrated near modern waterways and areas underlain by Holocene fluvial deposits with shallow water tables (< 1 to 2 m). In southern Christchurch, spatial variations of liquefaction and subsidence were documented in the suburbs within inner meander loops of the Heathcote River. Newly acquired geospatial data, geotechnical reports and eye-witness discussions are compiled to provide a detailed account of the surficial effects of CES liquefaction and ground deformation adjacent to the Heathcote River. LiDAR data and aerial photography are used to produce a new series of original figures which reveal the locations of recurrent liquefaction and subsidence. To investigate why variable liquefaction patterns occurred, the distribution of surface ejecta and associated ground damage is compared with near-surface sedimentologic, topographic, and geomorphic variability to seek relationships between the near-surface properties and observed ground damages. The most severe liquefaction was concentrated within a topographic low in the suburb of St Martins, an inner meander loop of the Heathcote River, with liquefaction only minor or absent in the surrounding areas. Subsurface investigations at two sites in St Martins enable documentation of fluvial stratigraphy, the expressions of liquefaction, and identification of pre-CES liquefaction features. Excavation to water table depths (~1.5 m below the surface) across sand boils reveals multiple generations of CES liquefaction dikes and sills that cross-cut Holocene fluvial and anthropogenic stratigraphy. Based on in situ geotechnical tests (CPT) indicating sediment with a factor of safety < 1, the majority of surface ejecta was sourced from well-sorted fine to medium sand at < 5 m depth, with the most damaging liquefaction corresponding with the location of a low-lying sandy paleochannel, a remnant river channel from the Holocene migration of the meander in St Martins. In the adjacent suburb of Beckenham, where migration of the Heathcote River has been laterally confined by topography associated with the volcanic lithologies of Banks Peninsula, severe liquefaction was absent with only minor sand boils occurring closest to the modern river channel. Auger sampling across the suburb revealed thick (>1 m) clay-rich overbank and back swamp sediments that produced a stratigraphy which likely confined the units susceptible to liquefaction and prevented widespread ejection of liquefied material. This analysis suggests river migration promotes the formation and preservation of fluvial deposits prone to liquefaction. Trenching revealed the strongest CES earthquakes with large vertical accelerations favoured sill formation and severe subsidence at highly susceptible locations corresponding with an abandoned channel. Less vulnerable sites containing deeper and thinner sand bodies only liquefied in the strongest and most proximal earthquakes forming minor localised liquefaction features. Liquefaction was less prominent and severe subsidence was absent where lateral confinement of a Heathcote meander has promoted the formation of fluvial stratum resistant to liquefaction. Correlating CES liquefaction with geomorphic interpretations of Christchurch’s Heathcote River highlights methods in which the performance of liquefaction susceptibility models can be improved. These include developing a reliable proxy for estimating soil conditions in meandering fluvial systems by interpreting the geology and geomorphology, derived from LiDAR data and modern river morphology, to improve the methods of accounting for the susceptibility of an area. Combining geomorphic interpretations with geotechnical data can be applied elsewhere to identify regional liquefaction susceptibilities, improve existing liquefaction susceptibility datasets, and predict future earthquake damage.
