skip to main content


The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 11:00 PM ET on Friday, September 29 until 11:59 PM ET on Saturday, September 30 due to maintenance. We apologize for the inconvenience.

Title: Landslides caused by the Mw7.8 Kaikōura earthquake and the immediate response
Tens of thousands of landslides were generated over 10,000 km2 of North Canterbury and Marlborough as a consequence of the 14 November 2016, Mw7.8 Kaikōura Earthquake. The most intense landslide damage was concentrated in 3500 km2 around the areas of fault rupture. Given the sparsely populated area affected by landslides, only a few homes were impacted and there were no recorded deaths due to landslides. Landslides caused major disruption with all road and rail links with Kaikōura being severed. The landslides affecting State Highway 1 (the main road link in the South Island of New Zealand) and the South Island main trunk railway extended from Ward in Marlborough all the way to the south of Oaro in North Canterbury. The majority of landslides occurred in two geological and geotechnically distinct materials reflective of the dominant rock types in the affected area. In the Neogene sedimentary rocks (sandstones, limestones and siltstones) of the Hurunui District, North Canterbury and around Cape Campbell in Marlborough, first-time and reactivated rock-slides and rock-block slides were the dominant landslide type. These rocks also tend to have rock material strength values in the range of 5-20 MPa. In the Torlesse ‘basement’ rocks (greywacke sandstones and argillite) of the Kaikōura Ranges, first-time rock and debris avalanches were the dominant landslide type. These rocks tend to have material strength values in the range of 20-50 MPa. A feature of this earthquake is the large number (more than 200) of valley blocking landslides it generated. This was partly due to the steep and confined slopes in the area and the widely distributed strong ground shaking. The largest landslide dam has an approximate volume of 12(±2) M m3 and the debris from this travelled about 2.7 km2 downslope where it formed a dam blocking the Hapuku River. The long-term stability of cracked slopes and landslide dams from future strong earthquakes and large rainstorms are an ongoing concern to central and local government agencies responsible for rebuilding homes and infrastructure. A particular concern is the potential for debris floods to affect downstream assets and infrastructure should some of the landslide dams breach catastrophically. At least twenty-one faults ruptured to the ground surface or sea floor, with these surface ruptures extending from the Emu Plain in North Canterbury to offshore of Cape Campbell in Marlborough. The mapped landslide distribution reflects the complexity of the earthquake rupture. Landslides are distributed across a broad area of intense ground shaking reflective of the elongate area affected by fault rupture, and are not clustered around the earthquake epicentre. The largest landslides triggered by the earthquake are located either on or adjacent to faults that ruptured to the ground surface. Surface faults may provide a plane of weakness or hydrological discontinuity and adversely oriented surface faults may be indicative of the location of future large landslides. Their location appears to have a strong structural geological control. Initial results from our landslide investigations suggest predictive models relying only on ground-shaking estimates underestimate the number and size of the largest landslides that occurred.  more » « less
Award ID(s):
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more » ; ; ; ; ; ; ; ; ; ; ; ; « less
Date Published:
Journal Name:
Bulletin of the New Zealand Society for Earthquake Engineering
Page Range / eLocation ID:
106 to 116
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Bedrock landsliding, including the formation of landslide dams, is a predominant geomorphic process in steep landscapes. Clarifying the importance of hydrologic and seismic mechanisms for triggering deep‐seated landslides remains an ongoing effort, and formulation of geomorphic metrics that predict dam preservation is crucial for quantifying secondary landslide hazards. Here, we identify >200 landslide‐dammed lakes in western Oregon and utilize dendrochronology and enhanced14C dating (“wiggle matching”) of “ghost forests” to establish slope failure timing at 20 sites. Our dated landslide dataset reveals bedrock landsliding has been common since the last Cascadia Subduction Zone earthquake in January 1700 AD. Our study does not reveal landslides that date to 1700 AD. Rather, we observe temporal clustering ofat leastfour landslides in the winter of 1889/1890 AD, coincident with a series of atmospheric rivers that generated one of the largest regionally recorded floods. We use topographic and field analyses to assess the relation between dam preservation and topographic characteristics of the impounded valleys. In contrast to previous studies, we do not observe systematic scaling between dam size and upstream drainage area, though dam stability indices for our sites correspond with “stable” dams elsewhere. Notably, we observe that dams are preferentially preserved at drainage areas of ∼1.5 to 13 km2and valley widths of ∼25 to 80 m, which may reflect the reduced downstream influence of debris flows and the accumulation of mature conifer trees upstream from landslide‐dammed lake outlets. We suggest that wood accumulation upstream of landslide dams tempers large stream discharges, thus inhibiting dam incision.

    more » « less
  2. The magnitude Mw7.8 ‘Kaikōura’ earthquake occurred shortly after midnight on 14 November 2016. This paper presents an overview of the geotechnical impacts on the South Island of New Zealand recorded during the post-event reconnaissance. Despite the large moment magnitude of this earthquake, relatively little liquefaction was observed across the South Island, with the only severe manifestation occurring in the young, loose alluvial deposits in the floodplains of the Wairau and Opaoa Rivers near Blenheim. The spatial extent and volume of liquefaction ejecta across South Island is significantly less than that observed in Christchurch during the 2010-2011 Canterbury Earthquake Sequence, and the impact of its occurrence to the built environment was largely negligible on account of the severe manifestations occurring away from the areas of major development. Large localised lateral displacements occurred in Kaikōura around Lyell Creek. The soft fine-grained material in the upper portions of the soil profile and the free face at the creek channel were responsible for the accumulation of displacement during the ground shaking. These movements had severely impacted the houses which were built close (within the zone of large displacement) to Lyell Creek. The wastewater treatment facility located just north of Kaikōura also suffered tears in the liners of the oxidation ponds and distortions in the aeration system due to ground movements. Ground failures on the Amuri and Emu Plains (within the Waiau Valley) were small considering the large peak accelerations (in excess of 1g) experienced in the area. Minor to moderate lateral spreading and ejecta was observed at some bridge crossings in the area. However, most of the structural damage sustained by the bridges was a result of the inertial loading, and the damage resulting from geotechnical issues were secondary. 
    more » « less
  3. Following the November 14 2016 Mw7.8 Kaikoura earthquake, field expeditions were undertaken using Unmanned Aerial Vehicles (UAVs) to map 25 sites of scientific interest with a plan area of 7.2 km2. A total of 23,172 images collected by the UAVs were used as input in Structure-from-Motion (SfM) to create 3D models of the target areas with a focus on landslides and fault rupture. Two sites are presented in more detail as examples of the data generated; a section of the Kekerengu fault that ruptured during the earthquake, and the Limestone Hills landslide. The sites were mapped at high resolution with ground sampling distance that varied from 0.5 to 7.0 cm/pixel. The developed SfM models were compared to 1-m aerial LiDAR data and the results were found to be comparable. However, the higher resolution of the SfM digital surface model (DSM), paired with the imagery facilitated more detailed interpretations, highlighting the usefulness of the UAV-enabled SfM as a mobile and effective technique for documenting perishable post-earthquake reconnaissance data. 
    more » « less
  4. Abstract

    The effects of strong ground shaking on hillslope stability can persist for many years after a large earthquake, leading to an increase in the rates of post earthquake land sliding. The factors that control the rate of post‐earthquake land sliding are poorly constrained, hindering our ability to reliably forecast how landscapes and landslide hazards and risk evolve. To address this, we use a unique data set comprising high‐resolution terrestrial laser scans and airborne lidar captured during and after the 2010–2011 Canterbury Earthquake Sequence, New Zealand. This earthquake sequence triggered thousands of rock falls, and rock and debris avalanches (collectively referred to as “rockfall”), resulting in loss‐of‐life and damage to residential dwellings, commercial buildings and other infrastructure in the Port Hills of Christchurch, New Zealand. This unique data set spans 5 years and includes five significant earthquakes. We used these data to (a) quantify the regional‐scale “rockfall” rates in response to these earthquakes and the postearthquake decay in rockfall rates with time; and (b) investigate the site‐specific factors controlling the location of seismic and nonseismic rockfalls using frequency ratios and logistic regression techniques. We found that rockfall rates increased significantly in response to the initial earthquake that generated the strongest shaking in the sequence—The MW6.2 22 February 2011 event—Compared to the long‐term background rates derived from the dating of pre‐2010 talus piles at the toe of the slopes. Non seismic rockfall rates also increased immediately after the 22 February 2011 earthquake and decayed with time following a power‐law trend. About 50% of the decay back to the pre‐earthquake rockfall rates occurred within 1–5 years after the 22 February 2011 earthquake. Our results show that the short‐term decay in rockfall rates over time, after the initial earthquake, was attributed to the subsequent erosion of seismically damaged rock mass materials caused by environmental processes such as rain. For earthquake‐induced rockfall at the regional‐scale, the peak ground accelerations is the most significant variable in forecasting rockfall volume, followed by the relative height above the base of the slope. For both earthquake and non‐seismic conditions at the site‐specific scale, the probability of rockfall increases when the adjacent areas have failed previously, indicating that accrued damage preconditions localized areas of the slope for subsequent failure. Such preconditioning is a crucial factor driving subsequent rockfalls; that is, future rockfalls are likely to cluster near areas that failed in the past.

    more » « less
  5. Abstract

    We generated dense, high‐resolution 3‐D ground displacement maps for the 2016 MW7.8 Kaikōura, New Zealand earthquake—the most geometrically and kinematically complex rupture yet recorded—from stereo WorldView optical satellite imagery using a new methodology that combines subpixel image correlation with a ray‐tracing approach. Our analysis reveals fundamental new details of near‐field displacement patterns, which cannot easily be obtained through other methods. From our detailed correlation maps, we measured fault slip in 3‐D along 19 faults at 500‐m spacing. Minimum resolvable horizontal slip is ~0.1 m, and vertical is ~0.5 m. Net slip measurements range from <1 to ~12 m. System‐level kinematic analysis shows that slip on faults north of the Hope fault was oriented primarily subparallel to the Pacific‐Australian plate motion direction. In contrast, slip on faults to the south was primarily at high angle to the plate motion and secondarily parallel to plate motion. Fault kinematics are in some locations consistent with long‐term uplift patterns, but inconsistent in others. Deformation within the Seaward Kaikōura Range may indicate an attempt by the plate boundary fault system to geometrically simplify. Comparison of published field measurements along the Kekerengu fault with our correlation‐derived measurements reveals that ~36% of surface displacement was accommodated as distributed off‐fault deformation when considering only field measurements of discrete slip. Comparatively, field measurements that project previously linear features (e.g., fence lines) into the fault over apertures >5–100 m capture nearly all (~90%) of the surface deformation.

    more » « less