Rockfalls and rockslides during the past 12 000 years have deposited
bouldery debris cones on the seafloor beneath massive rock slopes throughout
the inner part of Knight Inlet. The 885 m high rock slope, located across
from a former First Nations village destroyed in the late 1500s by a
slide-induced wave, exposes the contact between a Late Cretaceous dioritic
pluton and metamorphic rocks of the Upper Triassic Karmutsen Formation. The
pluton margin is strongly foliated parallel to primary and secondary fabrics
in the metamorphic rocks, resulting in highly persistent brittle structures.
Other important structures include a set of sheeting joints and highly
persistent mafic dykes and faults. Stability analysis indicates that planar
and wedge rock slope failures up to about 500 000
Inlets on the British Columbia (BC) coast are likely to see future development in the form of fish farms, power generation, transmission lines, port facilities, pipelines, roads, and other coastal infrastructure, yet the hazards to this development remain underappreciated because the region, until recently, has been remote and sparsely populated. Historical observations in similar mountainous coastal and lacustrine environments have shown that tsunamis generated by subaerial landslides are common and can be highly destructive. For instance, Huber (1982) reported over 500 fatalities associated with about 50 subaerial landslides that generated displacement waves in lakes and reservoirs in Switzerland over the past 600 years. Slide-induced waves have claimed more than 200 lives in the past 320 years in Norway (Blikra et al., 2005) and 10 lives in April 2007 in southern Chile (Sepúlveda et al., 2010). Damaging waves induced by subaerial landslides have also occurred on the coasts of Alaska (Miller, 1960a, b) and Greenland (Dahl-Jensen et al., 2004).
Knight Inlet, which extends to within 40 km of the highest glacier-clad peaks in the Coast Mountains (Fig. 1), is one of the longest and deepest fjords on the BC coast. Its precipitous, high walls are susceptible to rockfalls and rockslides, and several landslide deposits have been identified on the seafloor beneath these slopes (Fig. 2; van Zeyl, 2009). This evidence suggests that the inlet has been prone to landslide-generated tsunamis over the past 12 000 years. A recent study integrating archaeological and geological observations with oceanographic analysis suggests that the former village of Kwalate, located near the shoreline of Knight Inlet, was destroyed by a slide-induced wave in the late 1500s, with the possible loss of as many as 100 lives (Bornhold et al., 2007). In addition to this event, a slide from Three Finger Peak in November 1999 generated a wave that jostled log booms about 20 km to the north at the head of Knight Inlet.
Location of Knight Inlet (K) on the southern British Columbia coast. Glaciers coincide with the highest parts of the Coast Mountains. V: Vancouver.
Terrestrial and submarine topography of the middle part of Knight Inlet, showing many colluvial deposits beneath steep rock slopes and side-entry deltas at the mouths of tributary streams. Many of the seafloor deposits mapped in this figure are hybrid colluvial–alluvial landforms; the alluvial deposit beneath Three Finger Creek, for example, contains material from a rock slope failure at Three Finger Peak that produced a wave in November 1999 that was observed 20 km to the north at the head of Knight Inlet. SONAR data used to generate DEM were provided by the Canadian Hydrographic Service (2008).
Expanding on the work of van Zeyl (2009), this paper documents the geology and stability of the rock slope (Fig. 3) that is the likely source for the late-1500s tsunami and is across Knight Inlet from Kwalate village (Fig. 2). The study provides a glimpse into the geological and morphological conditions encountered in an area of the BC coast suspected to be particularly prone to tsunamis generated by subaerial landslides. Apart from boat traffic, for example – associated with ecotourism, logging operations, and seasonal wild fish harvesting – there are few elements in Knight Inlet that are at risk from tsunamis, but this situation could change in the future with developments such as fish farming operations or port facilities. Currently there are no fish farms in the inner part of Knight Inlet, but in recent years the number of such operations in the inner parts of other inlets on the BC coast has increased.
The Coast Mountains form a belt of crystalline and metamorphic rocks 1600 km long and up to 150 km wide, bordering the Pacific Ocean from Alaska to the Fraser Lowland near Vancouver. They rise abruptly from the sea and towards the axis of the range are characterized by rugged peaks and saw-toothed ridges. Numerous fjords, some more than 600 m deep, dissect the western margin of the Coast Mountains (Fig. 1).
The steep-walled fjords of the BC coast are former Tertiary river valleys that were greatly deepened and widened by glaciers during the Pleistocene (Clague, 1991; Mathews, 1991). The last (Late Wisconsinan) glaciers of the Cordilleran ice sheet retreated eastward from the coast about 13 000–14 000 years ago, and Knight Inlet was probably ice-free by about 11 000 years ago (Clague, 1981; Dyke, 2004).
With a length of 125 km, Knight Inlet is one of the longest inlets on the BC coast (Fig. 1). It has a maximum depth of about 525 m and typical widths in the range of 2–4 km. The inlet extends 70 km eastward from the open coast and then bends to the north and winds its way through terrain of considerably higher relief to the mouths of the Klinaklini and Franklin rivers.
Van Zeyl (2009) identified postglacial submarine rockfall and rockslide deposits on the seafloor within the inner part of Knight Inlet (Fig. 2). The blocky conical deposits are present beneath steep mountainsides extending as high as 1200 m a.s.l. within one lateral kilometre of the shoreline. The deposits must have formed after deglaciation and thus during the past 12 000 years, but nothing more is known about their age or history. One possibility is that rockfall and rockslide activity that produced many of these deposits was higher during the early Holocene as an immediate response to glacial debuttressing (cf. Evans and Clague, 1994). Another possibility is that activity increased in the second half of the Holocene after thousands of years of progressive rock slope deformation and fatigue (Bjerrum and Jørstad, 1968) before the fjord walls could collapse. Knight Inlet lies within one of the most seismically active regions in Canada (Adams and Atkinson, 2005) – earthquakes are concentrated along or near the North American plate boundary, associated with the Cascadia subduction zone (Fig. 1) 200 km west of Knight Inlet, and the Queen Charlotte transform fault 400 km to the northwest.
The maritime temperate climate of the southern mainland BC coast is
characterized by average annual temperatures and monthly precipitation totals
at sea level of about 8–10
We conducted fieldwork at the rock slope at Adeane Point during two site visits. We mapped major structures exposed along the crest of the slope, performed four scanlines in the central gully, and made additional structural observations at outcrops in the central and west gullies and at the base of the east wall (Figs. 3, 4). The scanlines involved measuring and describing each discontinuity of at least 25 cm in length that intersected a tape measure extending across an outcrop. We also made structural observations during traverses on the west slope, extending from the west gully to the crest of the rock slope. Georeferenced air photographs, digital elevation models (DEMs), stereoscopic aerial photographs from the 1950s through the 1990s, and scaled oblique photographs provided additional structural information.
We assessed rock mass conditions using the rock quality designation (RQD; Deere et al., 1967) and the geological strength index (GSI; Marinos et al., 2005). RQD is the percentage of intact core pieces larger than 10 cm to the length of core run or scanline. The GSI is used to derive field-scale rock mass properties for use in continuum rock mechanics modelling (Hoek et al., 2002). We estimated intact rock uniaxial compressive strength in the field using a rock hammer (ISRM, 1978).
Oblique photograph of the rock slope at Adeane Point. WS: west slope; c: central gully; w: west gully.
Vertical image showing the subaerial and submarine parts of the rock slope at Adeane Point. The offshore area is represented by a slope map generated from swath bathymetry. WS: west slope; W: west wall; w: west gully; c: central gully; E: east wall; ES: east slope.
Data on the morphology of the fjord were derived from digital elevation models. We used Canadian Digital Elevation Data (CDED), 20 m resolution elevation models (GeoBase, 2008) to analyze the terrestrial parts of the fjord. A 2 m resolution elevation model generated from multibeam echo-sounding data provided by the Canadian Hydrographic Service (2008) allowed us to analyze the submarine parts of Knight Inlet.
We performed a kinematic analysis to identify possible modes of rock slope
failure. Geologic structure and slope morphology were displayed in
stereographic projections to evaluate the potential for planar, wedge, and
toppling failure (Hoek and Bray, 1981; Wyllie and Mah, 2003). We then
analyzed selected planar and wedge failure modes identified from kinematic
analysis using the program Swedge (Rocscience, 2008) to compute the factor of
safety of rigid blocks bounded by intersecting discontinuities. This analysis
allowed us to rank the key failure types in terms of relative stability and
to help estimate possible landslide volumes. The procedure involved computing
the factor of safety for each failure type over the full range of slopes and
aspects observed at the site. The analysis assumed no tension crack, an upper
slope of 15
Slope and aspect maps. The stereonet included with the aspect map shows the 3-D shape of the slope.
The north part of the unnamed mountain at Adeane Point is a bowl-shaped
erosional basin that has a distinct U shape in plan, with a radius of
curvature at the crest of the rock slope of about 400 m (Figs. 4, 5). The
surface area of the terrestrial portion of this erosional basin is about
1.5
A north–south line extending through the centre of the rock slope can be used to divide the basin into four main geometric components (Figs. 3–5): a NNE-dipping forested slope (“west slope”); a north-dipping rock wall (“west wall”); a WNW-dipping rock wall (“east wall”); and a WNW-dipping forested slope (“east slope”). The west slope and its submarine extension are controlled by primary and secondary fabrics in the underlying metamorphic rocks, whereas the orientation of the east wall is controlled by a set of sheeting joints in the underlying diorite, and the west wall is controlled mainly by a complex of sub-vertical structures, which include mafic dykes and faults.
The maximum elevation at Adeane Point is 885 m at the east peak. Based on
the DEM, relief from seafloor to this summit is about 1400 m, with an
average slope of about 45
The submarine debris cone at the base of the slope dips 20–24
We examined slope deposits in the central and west gullies (Figs. 3 and 4). The west gully lies largely within metamorphic rocks, and its source area is the west wall. Grass, lichens, and seedlings are present on the deposit, suggesting minimal recent rockfall activity from its source area. In contrast, vegetation is generally absent in the central gully below the east wall; no lichens and mosses were observed on the rock blocks, and the blocks are less weathered than those in the west gully, suggesting greater rockfall activity from the east wall. Oblique photographs also indicate a larger number of fresh rockfall scars on the east wall than on the west wall.
Profile of the rock slope at Adeane Point, including a schematic representation of submarine debris cone and fjord infill sediments.
Left: oblique view of the rock slope. Right: map showing major structures exposed in the slope.
Oblique photographs showing major structures at the crest of the rock slope. See text for explanation of discontinuity sets.
Rocks at Adeane Point consist of greenstone of the Upper Triassic Karmutsen Formation and Late Cretaceous diorite (Roddick and Woodsworth, 2006). The greenstone is massive to foliated and contains intercalations of marble, sandstone, argillite, phyllite, biotite-chlorite schist, and amphibolite (Bancroft, 1913; Roddick, 1977). The dioritic pluton exposed in the rock slope contains screens of amphibolites and schist and is in sharp contact with the greenstone (structure 1a, Fig. 7); both the diorite and greenstone are foliated near the contact.
Portions of some of the major structures exposed in the rock slope (Figs. 7, 8) were accessed from the top of the slope, and additional information on major structures was obtained from oblique photographs taken from a helicopter and from ground locations.
The most apparent and persistent structures, which we term set 1, are
parallel to the foliation. At the crest of the slope, the four largest and
most evident of these structures – labelled 1a, 1b, 1c, and 1d – are spaced
about 50 m apart and are at least 100 m long (Figs. 7, 8). Field
measurements of three of them indicate an average attitude of
54/028
A singular structure, labelled 2a, is exposed in a rockfall scar near the slope crest. It has roughly the same orientation as a distinct set of sheeting joints (set 2) and extends out of the scar and persists for at least 100 m, forming a potential sliding plane beneath a large rock mass encompassing the west peak (Figs. 7, 8). The lower part of this structure, beneath the west peak, is noticeably curvi-planar and does not appear to intersect structure 1a.
Foliation joints (set 1) in the intrusive and metamorphic rocks are some of
the most persistent structures exposed in the rock slope. The average
orientation of the foliation is 55/027
The attitudes of the subaerial portion and submarine extension of the west
slope, as measured from the digital elevation models (Figs. 5, 9b), are
similar to outcrop foliation measurements (Fig. 9f). For instance, the SONAR
data show that the average orientation of a large portion of the submarine
extension of the west slope is 47/023
The foliation at Adeane Point is typical of foliation elsewhere in the southern Coast Belt, where structures in plutons and pendants are dominated by northwest-trending foliation that commonly dips steeply northeast (Woodsworth et al., 1991) and most plutonic rocks exhibit some degree of foliation, which is best developed near pluton margins (Roddick and Hutchison, 1974).
Bedding in marble in the west gully and on the west slope generally dips
steeply to the north and east (Fig. 9g). The variability in dip is due to
ductile deformation of marble around black argillite lenses in the west
gully. The average bedding attitude of 48/044
Synthesis of field and remotely sensed structural data.
Given that bedding is parallel to foliation and that weak rocks such as argillite and phyllite have been observed in the metamorphic belt exposed at Adeane Point (Roddick, 1977), weak layers parallel to the foliation are probably present within the rock slope, and faulting may have occurred along these layers during emplacement of the pluton. Indeed, we observed fault-like structures beneath extensive foliation-parallel beds of massive greenstone in the central gully, and many more may exist than were observed in the field. Regional mappers have pointed out that it is almost impossible to determine whether faults exist in areas where the stratified rocks occur in thin beds (Bancroft, 1913).
A single structure (2b) associated with a distinct set of sheeting joints
(set 2) exposed near the crest has an attitude of 55/273
Terzaghi (1962) referred to sheeting joints as “valley joints”, because they generally conform to the shape of the valley in which they occur. Bjerrum and Jørstad (1968) identified them as an important control on rockslides and rockfalls in the fjords of Norway, and they are one of the main causes of rockfalls and rockslides along BC Highway 99 between Vancouver and Whistler (Hoek and Bray, 1981; Gilbert, 2008).
Mafic dykes cut both metamorphic and intrusive rocks. The dykes observed in
the field range in thickness from about 15 cm to 2 m and consist of very
strong to extremely strong, greenish-grey basalt. The dykes and the
structures along which they were emplaced display great persistence; some of
them were traced through the entire height of the west wall. The high degree
of fracturing in the dykes compared to the adjacent host rock likely renders
them as hydraulic conduits, which is supported by the presence of springs
emanating from dykes at the apex of the central gully. Some of the dykes are
parallel to the foliation (44/036
Northwest- and northeast-trending dykes and brittle faults of Tertiary age are common along this part of the coast (Woodsworth et al., 1991). According to Bancroft (1913, p. 116), the dykes in this area “are frequently irregular in their strike, but the greater number of them either assume a direction which is parallel or transverse to the trend of the Coast Range”, and “where the walls of the fiord are devoid of vegetation, these dykes often appear as dark, ribbon-like bands of a remarkably constant width, extending from the shoreline to the summit of mountains several thousand feet high”.
To facilitate stability analysis, we developed a simple representation of the main discontinuity sets at this site based on the structures described above and on some additional structural details provided below (Table 1, Fig. 9i).
Foliation, bedding, and associated faults form S1, and sheeting joints form
S2. Scanlines and spot measurements provide evidence of a north-trending
sub-vertical set (S3) of joints and quartz veins, as well as a set of joints
and quartz veins (S4) dipping about 45
Because the orientations of the dykes and linear topographic features at Adeane Point differ substantially, we defined two sub-vertical discontinuity sets to incorporate this variability into the stability analysis (S5 and S6). We defined S5 based on ENE-trending, sub-vertical joints and dykes intersected in scanlines; ENE-trending linear cliffs 100–150 m high on the submarine sidewall (Fig. 9a); and distinctive lineaments extending across the rock slope at 300–600 m elevation (Fig. 9). An important set of northeast-trending structures (S6) extends through the crest of the rock slope south of the east peak (Fig. 9); one of these structures intersects structure 1c, forming a wedge-shaped cavity (Figs. 7, 8).
Summary of interpreted discontinuity sets.
Most rocks at this site – including diorite, greenstone, and schist – are very
strong (100–250 MPa) to extremely strong (
Average discontinuity spacing computed from scanline data are in the range of 15–50 cm for all sets, similar to the average block size observed in the submarine debris cone (10–40 cm). These values represent close to moderate spacing (ISRM, 1978), which qualify the rock mass as very blocky to blocky according to Cai et al. (2004) and Marinos et al. (2005). Typical discontinuity surfaces are rough, slightly rough, or smooth; we did not observe clay-filled joints in our surveys. The typical rock mass at this site has a GSI ranging from 40 to 90, reflecting favourable rock mass conditions. Similarly, scanlines indicate RQD values in the range of 75–82 %, reflecting “fair” to “good” rock mass conditions (Deere et al., 1967).
The observed outcrop-scale discontinuity spacing values in both dioritic and metamorphic rocks are typical of those in similar rocks described elsewhere in the Coast Mountains. Roddick (1965), for example, noted that hornblende-rich plutonic rocks, like the diorite at Adeane Point, are more closely jointed and sheeted than biotite-rich plutonic rocks. Of the hornblende-rich rocks, medium- to coarse-grained rocks have larger joint spacing (60–120 cm) than finer-grained varieties (15–30 cm), and the smaller spacing in the finer-grained intrusive rocks is similar to that in adjacent metamorphic belts.
Although average outcrop-scale joint persistence for all sets is in the range of 2–4 m, reflecting low to medium persistence (ISRM, 1978), oblique photographs of the rock slope suggest much greater persistence. Persistence “censoring” was an unavoidable reality at the outcrop scale, although in some cases, joint terminations were observed either in intact rock or against other discontinuities. Some of the more persistent structures observed in oblique photographs may represent closely spaced discontinuities separated by intact rock bridges or step-path joints.
The strong to extremely strong intact rock strength and generally favourable rock mass conditions at this site, together with the high frequency of very persistent structures, suggest that structurally controlled failure is important. Although progressive stress-induced fracturing of intact rock bridges may be required to permit larger structurally controlled failures, we would expect the dominant shape of failure surfaces to follow structures rather than involve a pseudo-circular failure surface. As a result, stability analyses based on structurally controlled failure mechanisms are considered more appropriate for this site.
We performed a kinematic analysis to identify possible failure modes using the discontinuity sets summarized in Table 1. Although at the outcrop scale only three or four of the six sets might be present, all sets were considered in order to evaluate the true range of potential failure modes within different parts of the slope.
We compared the average orientation of the east and west walls to the interpreted rock mass structure to test for feasible modes of failure (Fig. 10). The most probable failure modes are planar sliding on S2 in the east wall (Fig. 10a) and wedge failure on the intersections of S1–S6 and S1–S3 in the west wall (Fig. 10d). However, if we account for the variability in discontinuity orientations and slope attitudes, wedge failure on the above intersections becomes feasible in the east wall and additional wedge intersections are feasible in both walls. Toppling on S6 in the east wall (Fig. 10a) and on S4 and S5 in the west wall (Fig. 10c) are also possible if variations in discontinuity orientations are considered. Based on evidence from oblique photographs, block toppling could be an important mechanism in localized over-steepened slope sections.
Given that wedge intersections in addition to those described above daylight on relatively steep slopes, the very steep slopes in the upper part of the rock slope would appear to be the most susceptible to structurally controlled failures. Confining pressures are also lowest on this part of the slope. In this context, it is interesting to note that rockfall scars are most common near the crest of the slope.
We analyzed selected failure types identified in the kinematic analysis using
Swedge. The orientations of S1, S2, S3, and S6, as summarized in Table 1,
were used as input in the stability analysis. The results of the analysis are
expressed as the stability (
As is evident in Table 2, planar failure is feasible over a small range of
aspects on the west slope and east wall. We call attention to the fact that
the dip of S2 used in this analysis (Table 1) is based on an average for the
site, and joints from this set display steeper dips near the top of the
slope. Analysis with a dip angle of 55
Equal area lower hemisphere stereonets showing
Wedge failure is possible throughout the centre of the rock slope and on the
periphery, covering a much wider range of aspects than planar failure. In
this context, S1 is particularly important because it is involved in all
three wedge combinations. S3 and S6 intersect with S1 to form steeply
plunging intersections (about 50
Results of stability analysis.
The sketch shown in Fig. 11 illustrates two different scales of possible
wedge failure at the slope crest. Considering a slope 200 m wide and 75 m
high, with an orientation intermediate between that of the east and west
walls (85/330
As shown in Figs. 8 and 11, a large wedge, perhaps about 0.5
An important factor that we did not consider in our kinematic or limit equilibrium analyses is the possibility that failures could be triggered by earthquakes or intense precipitation events. These events could also contribute to progressive failure through the destruction of rock bridges along partially formed failure surfaces. Another possibility is that earthquakes or intense precipitation events might trigger numerous small failures over a wide area of the rock slope. In this way, multiple small failures could produce a larger cumulative volume of debris than a single large, structurally controlled failure. Considering these factors, it is possible that the cumulative size of a slope failure could be significantly larger than that of one simple, structurally controlled wedge or planar failure.
Sketch showing two scales of wedge-shaped blocks at the crest of the rock slope at Adeane Point. Wedges formed by S1–S6 have dimensions of 50–75 m. The large wedge formed by S1–S2 has dimensions of about 75 and 250 m. Azimuths are the dip directions of the rock slope at the respective locations.
Several different analytical methods for estimating the heights of tsunamis triggered by subaerial landslides have been described in the literature (e.g. Noda, 1970; Slingerland and Voight, 1979; Huber and Hager, 1997; Fritz et al., 2004). Typically, these methods are used to estimate maximum wave heights within a short distance of the landslide impact site. They commonly produce conservative results that are sensitive to the input parameter values, such as the thickness of the sliding mass and the characteristic water depth beyond the impact site. Many of these values cannot be accurately specified because of the irregular geometry of the slide mass and variable submarine topography. Different methods may produce such a wide range of estimated wave heights that the analysis is unhelpful for hazard evaluation.
Bornhold et al. (2007) hindcasted wave amplitudes from a past landslide at
Adeane Point, assuming that the 3.5
Our study focuses on the potential for future landslides, and our analysis of
the structure, morphology, and stability of the rock slope at Adeane Point
suggests the potential for failures of sizes ranging up to about
0.5
An example of a potential development that would be vulnerable to landslide-induced waves in Knight Inlet is fish farms. The tsunami caused by the earthquake-triggered landslide in southern Chile in April 2007 (Sepúlveda et al., 2010) illustrates the type of consequences that such events can have on fish farm operations, including loss of life, damage to facilities, and the release of a large number of farmed fish into a natural aquatic ecosystem.
The November 1999 Three Finger Peak landslide-induced wave further indicates the potential for tsunamis in Knight Inlet.
The abundance of steep rock slopes towering above rockfall and rockslide deposits on the seafloor in the inner part of Knight Inlet suggests a history of tsunami-generating landslides over the past 12 000 years. A recent study integrating archaeological and geological observations with oceanographic analysis suggests that the former village of Kwalate in Knight Inlet was destroyed in the late 1500s by a landslide-induced wave. In light of this disaster, we studied the structure and stability of the rock slope across Knight Inlet from Kwalate village to estimate likely rock slope failure mechanisms and the size of possible future slope failures.
The contact between a Late Cretaceous dioritic pluton and a metamorphic belt consisting of greenstone of the Upper Triassic Karmutsen Formation is exposed in the rock slope at Adeane Point. The greenstone and adjacent diorite are strongly foliated, resulting in highly persistent brittle structures. Other important structures include a set of sheeting joints and highly persistent mafic dykes and faults. Stability analysis showed the potential for planar and wedge failures and highlighted the potential for wedge failures near the crest of the slope.
We performed kinematic and limit equilibrium analyses to identify and rank
the importance of different failure mechanisms and to estimate possible
failure volumes. Based on the observed structure and morphology of the rock
mass at Adeane Point, we see the potential for failures up to about
0.5
Landslides of the nature and size described above would generate waves that might damage coastal infrastructure or vessels underway in the fjord. Given that there are several similar steep rock slopes with large submarine debris cones in the inner part of Knight Inlet, this is an area where tsunami hazards should be considered in coastal infrastructure development and land-use planning.
Financial assistance for the project was provided by a Natural Sciences and Engineering Research Council of Canada grant to D. Stead and a graduate fellowship from Simon Fraser University to D. P. van Zeyl. The authors are grateful to Steve Israel, Dan Gibson, and Bert Struik for discussions on the structural geology of the rock slope at Adeane Point. We thank journal reviewers V. Kveldsvik, C. B. Harbitz, and F. Zaniboni for their helpful comments and suggestions, which have improved the manuscript considerably. Edited by: I. Didenkulova Reviewed by: C. B. Harbitz and F. Zaniboni