Brardinoni Et Al 2002
Brardinoni Et Al 2002
Brardinoni Et Al 2002
Received 1 August 2002; received in revised form 17 November 2002; accepted 29 November 2002
Abstract
Landslide inventories are routinely compiled by means of aerial photo interpretation (API). When examining photo pairs, the
forest canopy (notably in old-growth forest) hides a population of ‘‘not visible’’ landslides. In the present study, we attempt to
estimate how important is the contribution of landslides not detectable from aerial photographs to the global mass of sediment
production from mass failures on forested terrain of the Capilano basin, coastal British Columbia. API was coupled with
intensive fieldwork for identification and measurement of all landslides. A 30-year framework was adopted. We show that ‘‘not
visible’’ landslides can represent up to 85% of the total number of failures and account for 30% of the volume of debris
mobilised. Such percentages display high sub-basin variability with rates of sediment production varying by one order of
magnitude between two sub-basins of the study area. This is explained qualitatively by GIS-based analysis of slope frequency
distributions, drainage density, and spatial distribution of surficial materials. Such observations find further support in the
definitions of transport-limited and supply-limited basins. As a practical consideration to land managers, we envisage that
supplementary fieldwork for landslide identification is mandatory in transport-limited systems only. Fieldwork has
demonstrated that gully-related failures have a greater importance than one could expect from API.
D 2002 Elsevier Science B.V. All rights reserved.
Keywords: Landslide identification; Air-photo interpretation; Interbasin variability; Forested terrain; Coastal British Columbia
0169-555X/02/$ - see front matter D 2002 Elsevier Science B.V. All rights reserved.
doi: 10.1016/S0169-555X(02)00355-0
ARTICLE IN PRESS
2 F. Brardinoni et al. / Geomorphology 1321 (2002) 1–18
Fig. 3. Aerial photographs taken in 1996 showing part of (A) Sisters Creek and (B) East Cap Creek basins. Note the higher degree of landscape
dissection in Sisters Creek when compared to East Cap Creek. Examples of types of forest cover are marked with capital letters: A = recent cut-
blocks, B = old-growth forest, C = old extensive clearcut. Examples of visible mass movement scars are marked with numbers: 1 = debris slides
in recent cut-blocks, 2 = debris flows in old-growth forest reaching Strachan Creek, 3 = debris slide on old clearcut crossing a logging road and
reaching Sisters Creek.
ARTICLE IN PRESS
6 F. Brardinoni et al. / Geomorphology 1321 (2002) 1–18
Table 1
British Columbia Terrain Stability Classificationa
Terrain stability class Subjective rating system Interpretation
I Flood plains and level to undulating coastal plain areas No significant stability problems exist
Most terrain with slopes < 20%
II Most gently sloping (20 – 40%), poorly to well drained, Very low likelihood of slides following
lower slope landforms logging or road construction
Moderately sloping (40 – 60%), well to rapidly drained Minor slumping is expected along road cuts,
surficial deposits especially for 1 or 2 years following
construction
III Moderately sloping (40 – 60%), imperfectly to poorly Minor stability problems can develop
drained surficial deposits that are not marine or Timber harvesting should not significantly
lacustrine reduce terrain stability; there is a low
Level to gently sloping (0 – 40%), imperfectly likelihood of landslide initiation following
to poorly drained, deep marine clays and lacustrine logging
deposits Minor slumping is expected along road cuts
Moderately sloping, deeply gullied surficial Low likelihood of landslide initiation
deposits that are not of lacustrine or marine origin following road building
IV Steeply sloping (>60%), well drained, deeply Expected to contain areas with a moderate
gullied surficial deposits likelihood of slide initiation following
Steeply sloping, poorly drained surficial deposits logging or road construction
Moderately sloping, deeply gullied, or imperfectly
to poorly drained lacustrine or marine deposits
V Any areas where natural landslide scars are visible Expected to contain areas with high
on air photographs or in the field likelihood of slide initiation following
Very steeply sloping (>70%), imperfectly to logging or road construction
poorly drained, deeply gullied surficial deposits
a
Modified from British Columbia Ministry of Forests, 1995.
guidelines (Table 1). A polygon was considered as an were minimised by comparing polygons of equal land
experimental unit, and each polygon that received the use and stability (according to the British Columbia
same classification (‘‘treatment,’’ i.e., survey type, Terrain Stability Classification).
land use, stability class) was regarded as a replication. A two-step statistical analysis was performed.
Two different types of survey to calculate landslide First, one-way ANOVA-like tests were conducted to
density and the associated amount of mobilised vol- evaluate how the survey method affected landslide
ume of debris were conducted: (i) exclusively by API; frequency per unit area (LS/ha) and amount of mobi-
and (ii) by API coupled with intensive field survey. lised debris (m3/ha). Second, the experimental design
Control on experimental error was achieved by select- was made more complex (two-way ANOVA-like
ing study sites with homogeneous climate, geological tests) so that the location effect (Sisters vs. East
setting, and vegetation cover. Confounding factors Cap) on survey capability could be incorporated.
We took advantage of replications. The greater the
number of replications, the larger the number of
Table 2
Land use areas in East Cap and Sisters Creeks
degrees of freedom in the error term and the smaller
the experimental error, which represents the unex-
Land use East Cap Creek Sisters Creek
(ha) (ha) plained variability of the dependent variable. Two-
way layouts possess a greater explanatory power than
Old growth (undisturbed) 409.3 (61.3%) 271 (44.1%)
Old growth (including fire) 578.8 (86.6%) 345.1 (56.1%) one-way analyses as they incorporate both factors
Old harvesting 0 (0) 261.3 (42.5%) (location and survey type) over the entire database.
Recent harvesting 89.2 (13.4%) 8.6 (1.4%) In this way, treatment interactions can be evaluated,
Fire 169.5 (25.3%) 74.1 (12%) and unexplained variability (experimental error) is
Total surveyed 668 615
minimised.
ARTICLE IN PRESS
F. Brardinoni et al. / Geomorphology 1321 (2002) 1–18 7
Table 3
Analysis of variance and nonparametric tests of landslide density and volume
(a) One-way ANOVA-like of the entire database (survey)
Factor Dependent variable One-way ANOVA Kruskal – Wallis ANOVA Median test
Significance
Survey Volume (m3/ha) ns 0.0003 0.0001
Number (LS/ha) 0.002 0.0000 0.0000
The study of the spatial distribution of surficial englacial or supraglacial conditions. It is therefore still
materials across slope stability classes focused partic- adjusting to a subaerial, nonglacial environment. Con-
ularly on glacial deposits (or till) and post-glacial versely, colluvium deposits, products of subaerial
deposits (or colluvium). The former is considered to mass movement, are arranged in more stable config-
be less stable, having been deposited under subglacial, urations. Typically, in each terrain polygon, the nature
Table 4
Sediment yields from landsliding and landslide densities during an f 30-year window (1968 – 2000) as obtained from both air-photo
interpretation coupled and not coupled with detailed field survey
Field Number of LS Mobilised volume (m3) LS density Annual LS Denudationa Yieldb
surveyed APIc Field A&Fd API Field A&F (#LS/ha) density (m3/ha) (m3/ha/year)
area (ha) API A&F (#LS/ha/year) API A&F API A&F
API A&F
East Cap 668.4 14 20 34 7869 412 8281 0.021 0.051 0.0007 0.0017 11.7 12.4 .39 .41
Creek (41.2%) (58.8%) (95%) (5%)
Sisters 615 20 117 137 52,480 21,999 74,479 0.033 0.223 0.0011 0.0074 85.3 121.1 2.84 4.04
Creek (14.6%) (85.4%) (70.5%) (29.5%)
Sum 1283.4 34 137 171 60,349 22,411 82,760 0.027 0.133 0.0009 0.0044 46.9 64.5 1.56 2.15
a
Denudation = mobilised volume/surveyed area.
b
Yield=(mobilised volume/surveyed area)/time window.
c
API = Air-photo interpretation.
d
A&F = Air-photo coupled with fieldwork.
ARTICLE IN PRESS
F. Brardinoni et al. / Geomorphology 1321 (2002) 1–18 9
4.1. Landslide visibility and implications landslide scar in recent openings is about 150 m2
(Fig. 5A), but note the gap between the largest ‘‘not
Fig. 5A and B reports areas of visible (from air visible’’ scar (150 m2) and the smallest visible scar
photo) and ‘‘not visible’’ landslides that were meas- (275 m2). This is due to the small number of obser-
ured during field survey in, respectively, recent cut- vations (n = 15), and the critical area for visibility lies
blocks (logged < 15 years ago) and forests older than somewhere between these two values.
50 years. The maximum area of a ‘‘not visible’’ In forests older than 50 years, the largest landslide
scar not detected during API has an area of about 650
m2 (Fig. 5B). One landslide smaller than this (350 m2)
was identified on an air photo; this exemplifies how
the issue of landslide visibility involves various fac-
tors, such as topographic location and landscape
dissection, other than land use. As a first practical
application, then, if one wishes to perform an
unbiased (although incomplete) air-photo-based com-
parison in this area between managed and wild forest,
one should consider only failures that are> 650 m2.
Inevitably, landslides missed during API produce a
supplementary amount of sediment. The plotting of
the cumulative percent volumes against landslide
volume (Fig. 6) shows that additional amounts of
debris inventoried by means of fieldwork in East
Cap and Sisters Creeks range from 10% to 18%
among volume classes. This is clear evidence that
‘‘not-visible’’ failures do contribute significant vol-
umes of debris to the amount of mobilised material
already detected on air photos.
An important supplementary question is whether
missed events contributed significantly to the sedi-
ment loading of montane stream channels. This aspect
Fig. 5. Landslide visibility in relation to landslide area (initiation
has clear management implications in terms of water
plus transportation zone) in (A) recently harvested cut-blocks and quality (e.g., in the Capilano Reservoir) and fish
(B) forests older than 50 years. habitat; and it is also important for estimating the
ARTICLE IN PRESS
10 F. Brardinoni et al. / Geomorphology 1321 (2002) 1–18
error associated with remotely sensed sediment budg- and from an applied perspective, it shows that the rate
ets. Landslides have been classified as stream-con- of sediment in-filling of the Capilano Reservoir has
nected, gully-connected, and not connected, depend- been underestimated.
ing on whether debris was delivered to a permanent The combined effect of fieldwork and air-photo-
stream, to a gully, or remained in an unchannelled inventory was to increase the percentage of gully-
portion of the landscape. According to the field data, connected failures and the associated mobilised vol-
about 58% of the mobilised volume (Fig. 7B) enters ume by about four times (from 8.8% to 32.2% and
the channel network directly; and 35% is delivered to from 3.2% to 11.8%, Fig. 7A and B). The proportion
gullies and most probably will be evacuated to the of stream-connected failures decreased by nearly one-
stream network via debris torrent in a decennial time third, and the mobilised volume decreased only 5.5%.
scale (from 10 to 40 years, cf. Thurber Engineering, Percentage of unconnected slides showed minor
1996). However, API indicated that only 3.2% of the changes. Integrating field with air-photo measure-
mobilised volumes were delivered to gullies, and here ments, 72.6% of mobilised volumes (Fig. 7B) were
lies the highest discrepancy between air-photo and delivered to the stream channels. In other words, 49%
field data. The implication is that most of the small of ‘‘not visible’’ events were connected to stream
failures are connected to the drainage network. From a channels and 38% to gullies (Fig. 7A). The gully
process perspective, this finding elucidates the real network appears to have been very stable in the last 30
degree of geomorphic coupling of montane streams; years. Only the formation of one new gully (scoured
Fig. 7. Sediment delivery to streams as obtained by air-photo interpretation (API), fieldwork (FW), and field-coupled air-photo interpretation
(API and FW). Lumped data for East Cap and Sisters Creeks are reported in percentages (A and B). Data for East Cap Creek (C and D) and for
Sisters Creek (E and F) are reported in numbers.
ARTICLE IN PRESS
F. Brardinoni et al. / Geomorphology 1321 (2002) 1–18 11
block-by-treatment interaction (i.e., location-by-sur- As for the volumes of debris per hectare (m3/ha),
vey), while such interaction is of primary interest they show no significant survey-by-location interac-
when blocking is used to broaden the scope of tion; hence, one can comment on survey and location
inference (León and Mee, 2000). separately. The former is not significant; the latter is
As for the full factorial experiment, the advantage (Table 3b) with Sisters greater than East Cap. The
of having a factorial arrangement of treatments (or nonparametric equivalent of a two-way ANOVA, the
factors) over a block design lies in that one can Friedman test, requires the same assumptions as the
analyse interactions and, where interactions are not fixed block design (i.e., null survey – location inter-
statistically significant, one can reduce the number of action); and, as such, it can be performed only on
treatment means in the multiple comparison test. mobilised volumes of debris (Table 3b). The test
According to the full factorial ANOVA (Table 3b), indicates that at least one of the survey – location
the interaction is significant for landslide density. combinations is different from the others; this is the
Therefore, no comment can be made on locations or field-coupled API conducted in Sisters Creek.
survey separately; instead, having detected the sig- We performed the Tukey’s Honest Significant
nificant interaction, its meaning has to be graphically Difference (HSD, a generalisation of Tukey’s test
evaluated by inspecting the plot of the means (Fig. 8A for unequal sample sizes) multiple comparison test,
and B). The plots show that the coupling of API with which is considered to be of intermediate conserva-
intensive fieldwork has sensibly increased average tism. The test, in line with the Friedman’s results,
values of landslide density and volumes in Sisters reports that the coupling of API with intensive field-
Creek. In East Cap Creek, such increase appears to be work in Sisters Creek gives a significantly greater
negligible. This difference between the two locations amount of mobilised debris per unit area (m3/ha) than
is particularly large in terms of landslide density does any other survey – location combination. This
(significant location-by-survey interaction). result stresses the importance of conducting fieldwork
in the Sisters Creek area.
The database was then split in two location-wise;
hence a series of one-way ANOVA-like tests was
performed. Regardless of the dependent variable
examined, in East Cap, survey types were not sig-
nificantly different (K – W and Median tests); con-
versely, in Sisters Creek field-coupled survey was
significantly greater than the air-photo-based survey
(Table 3c).
Table 6
Management effects on rates of landsliding and denudation: comparison between air-photo-based and field-coupled landslide inventories
Land use Annual LS density Yield Management effects Yield
(#LS/ha 1/year 1) (m3/ha 1/year 1) Annual LS density (m3/ha 1/year 1)
API A&F API A&F (#LS/ha 1/year 1) API A&F
API A&F
Sisters Creek Undisturbed forest 0.0023 0.0093 1.92 3.49 1.0 1.0 1.0 1.0
Old logging 0.0005 0.0052 0.66 2.59 0.2 0.6 0.3 0.7
East Cap Creek Undisturbed forest 0.0005 0.0016 0.33 0.34 1.0 1.0 1.0 1.0
Recent logging 0.0026 0.0034 1.46 1.53 5.2 2.1 4.4 4.5
factors remain fairly constant (Table 6). From the land 40j in forested terrain, and landslide initiation was
manager’s perspective, fieldwork application does not found very unusual below such values (Brardinoni et
change the situation; on the other hand, fieldwork al., 2002). Percentages of slope steeper than 35j are
provides the geomorphologist with better sediment significantly different (Z-test), with Sisters Creek hav-
source data to construct a more realistic sediment ing 34% of DEM cells with slope gradient steeper than
budget. 35j and East Cap Creek only 19.5%.
The influence of slope aspect on landsliding
4.5. GIS as a tool for explaining interbasin variability deserves particular attention; several studies have
of landslide density and sediment yield examined this issue in the study area. Approximately
75% of the landslides inventoried by O’Loughlin
In seeking a qualitative explanation for the one (1972) in Howe Sound, Capilano watershed, and
order of magnitude difference in landslide denudation Seymour watershed exhibited a southerly aspect. He
rates between the two study areas, a GIS-based topo- justified such landslide preferential orientation by not-
graphic analysis was performed. Elevation, slope ing that ‘‘north-facing slopes are rocky and broken, a
gradient, slope aspect, drainage density, and spatial condition which discourages landslide formation,
partition of surficial materials across terrain stability while south-facing slopes are relatively uniform and
classes are the variables that were considered. underlain by an extensive unweathered till substratum’’
East Cap Creek is located at a significantly higher (p. 31). This spatial arrangement of surficial materials
elevation (832 m F 2) than Sisters Creek (664 m F 2). is likely to derive from the north-to-south flow that
Accordingly, maximum frequency was reached at the characterised the movement of Pleistocene ice sheets in
600 – 800-m category in Sisters, whereas elevation the Pacific Ranges (Armstrong and Brown, 1954).
frequency peaked at the 1000 –1200-m category in Relations between surficial materials and aspect have
East Cap Creek (Fig. 9A). Higher elevation might been investigated more systematically for the entire
have slightly contributed to a faster evacuation of Capilano watershed (Greater Vancouver Regional Dis-
glacial till deposits from East Cap upper- and mid- trict, 1999). They confirmed that rock dominates north-
slope locations, as higher elevation induces higher erly aspects; and colluvium and till are most extensive
precipitation and greater sediment transport. on west-, south-, and east-facing aspects.
Average slope in Sisters Creek is significantly Heaviest rains in the Pacific Ranges are brought by
higher (31.1j F 0.1) than in East Cap Creek southwesterly air flows; south- and west-facing slopes
(28.1j F 0.1). Although it is difficult to generalise are therefore more directly exposed; north exposures
about a critical slope gradient for slope failure, it is are likely to remain more uniformly damp; and alter-
reasonable to assert that terrain with slopes over 35j nation of wet/dry cycles is enhanced on south-facing
can be regarded as shallow, mass movement prone, slopes. In light of these climatic controls, south- and
when soil is present (Sidle et al., 1985). Data from west-oriented slopes would be expected to be more
Howe Sound and Chapman Creek, located in the same favourable for landsliding. Fig. 9C shows that south-
physiographic region, indicate a critical gradient for facing slopes are of comparable extent in East Cap
sliding of about 30j to 35j in cut-blocks and of 35j to and Sisters Creeks but that west-facing slopes are
ARTICLE IN PRESS
14 F. Brardinoni et al. / Geomorphology 1321 (2002) 1–18
(iii) higher drainage density: 3.6 km/km2 in Sisters Analysis of landslide density and denudation rates
Creek compared to 2.2 km/km2 in East Cap Creek. was tackled by considering two main factors affecting
mass failure detection: type of survey and location.
The former is external to the analysed system, and the
5. Discussion latter is an intrinsic property and can be regarded as
the spatial heterogeneity of the system’s propensity to
Estimation of the maximum size of undetectable fail. Intensive ground checking of the study areas has
landslide scars on air photos is of fundamental impor- shown that ‘‘not visible’’ landslides accounted for
tance for both scientists and forest managers. Never- about one-third of the total volume of debris mobi-
theless, little research has addressed the topic (e.g., lised via mass movements during the last 30 years.
Robison et al., 1999); this could be due to under- Even more interestingly, the two subwatersheds
estimation of the relative importance of small failures exhibited very different rates of sediment production
and to the physical- and cost-prohibitive requirements from landsliding (about one order of magnitude dis-
of fieldwork. crepancy), with Sisters Creek being the more active
In Fig. 11, we summarised the factors that have (Table 4). Although these two Capilano River tribu-
been identified to affect landslide visibility during taries possess many similar biophysical characteris-
API. Factors that are manifestly related to the quality tics, they behave in strikingly different ways.
of photography itself [such as nominal scale, sensor The overall survey-by-location picture can be sum-
type (colour, black and white, etc.), weather condi- marised by saying that the same treatment (survey
tions (snow cover, clouds)] are not included. Rela- type) performed differently in different locations
tively large failures that frequently were not detected (blocks). Field-coupled surveys exhibited significantly
in antecedent API were generally located on the lower greater landslide densities and denudation in Sisters
portion of very steep (>40j) old-growth forested Creek. The same variables in East Cap Creek were
slopes. Due to the progressive downhill increase of shown to be rather indifferent to the conduct of field-
pore water pressure, clusters of debris slides are often work as a supplement to API, thus confirming the
found in these locations. On recently harvested cut- simple counts. Friedman’s test showed that field-
blocks, the problem in detecting slides lies in the coupled API in Sisters Creek gives significantly
shallow nature of the events, which facilitates a very greater mobilised volumes than any other survey –
fast process of revegetation. Furthermore, gully- location combination did. This finding emphasises
related failures connected to permanent streams are the importance of intensive ground checking. Solely
more difficult to detect. This is because stream-con- from air photos, landslide sediment production in East
nected failures have their deposition zone (e.g., flow Cap and Sisters Creeks would have not been signifi-
fans) readily washed out by the stream they have cantly different. Apparently, some types of forested
impacted. terrain hide important numbers of ‘‘not visible’’ land-
slides, while others do not. In other words, the issue of One of the main issues that puts air-photo-based
‘‘not visible’’ events can be complex even within a landslide studies under suspicion is the relative bias
‘‘medium-size’’ drainage basin. This finding has introduced (because of different conditions of visibil-
important practical implications: in East Cap-like ity) when one compares mass wasting activity in
areas, fieldwork is practically unnecessary and sedi- logged and undisturbed terrain. This case study has
ment budget evaluation requires virtually no correction shown that the discrepancies between air-photo-based
factor to account for not visible events. Conversely, and field-coupled surveys (Table 6) have relatively
Sisters-like areas require intensive fieldwork in order slight management implications. Fieldwork has the
to evaluate the ‘‘invisible’’ volume of mobilised debris. effect of shifting the old-logging – old-growth ratio
One aspect to consider in comparing fieldwork from negative accelerations to near constancy. Man-
landslide detection efficiency is the relative abun- agement effects of recent logging on denudation rates
dance of different land use cover—in our case, the remained constant survey-wise, while landsliding
relative extent of old-growth, old-logged, and recently ratios exhibited a three-time decrease as a result of
logged areas. In Sisters, the presence of a large portion fieldwork coupling. This is due to the higher number of
of old harvested terrain and only a small area of ‘‘hidden’’ slides in old-growth forest than in recent cut-
recently harvested cut-blocks (Table 2) does not blocks, whereas canopies of old-logged and old-growth
favour landslide detection from air photos. Con- forests have a similar effect on landslide visibility.
versely, the absence of the old-logging land use In order to explain the large discrepancy in denu-
category and the large, recently developed clearcut dation rate that was found between the two creeks, we
areas aid remote recognition of mass failures in East performed a GIS analysis. The outcome depicted
Cap (Table 2). We consider historic fires to have no Sisters and East Cap basins as two contrasting geo-
effect on landslide visibility; canopy height in fire- morphic environments (Fig. 3A and B). Several
affected forests was comparable to that in old-logged physiographic characteristics make Sisters basin a
and old-growth forested terrain. highly landslide-prone environment. This includes
A last benefit of conducting intensive fieldwork is steep slopes (one-third >35j), and till-mantled slopes
the recognition that gully-related events have a greater (f 40% class V) are a large source of potential non-
importance than one could expect from API. They coherent material. Also, high drainage density allows
constitute more than one-third (both in terms of efficient evacuation of the sediment load brought into
number of events and of mobilised volumes) of the the channel network via sidewall debris slides and
‘‘missed’’ events, while from air photos alone, they windthrow. Conversely, in East Cap, slopes are gen-
accounted for just below 10% of the total number of erally gentler; till is mainly located on less steep
failures and volumes. From aerial photos and field polygons (class IV); and drainage density is signifi-
traversing of Sisters and East Cap Creeks, the gully cantly lower, thus imparting a lower connectivity to
network appears to have been very stable in the last 30 the system.
years. The stream channel network functions as a These observations are in accordance with the
preferential transportation pathway for debris flows/ framework that subdivides systems into weathering-
torrents through the landscape. Unconnected failures and transport-limited systems (also termed supply-
account for relatively low percentages ( < 16%) of the limited and -unlimited) (Carson and Kirkby, 1972).
total number and volume of landslides (see Fig. 7). In In the first case, the limiting factor is sediment
this mountain environment, the drainage network is production; in the other, sediment mobilisation. In
extremely efficient in evacuating the sediment load. the mountainous, forested environment of coastal
This efficiency is achieved through alternation of a British Columbia, transport-limited basins typically
‘‘normal fluvial regime,’’ where supply-limited have a high density of headwater channels incised
amounts of small material systematically are eroded into thick glacial drift or closely jointed bedrock. This
from the stream channel, and of a ‘‘catastrophic ensures virtually unlimited debris supply, in addition
regime,’’ which entails removal of large quantities to many unstable trigger points for debris slides and
of sediment and woody debris via debris flows/tor- flows. Supply-limited basins denote slower recharge
rents (Nistor, 1996). rates and fewer zones of instability. This is usually
ARTICLE IN PRESS
F. Brardinoni et al. / Geomorphology 1321 (2002) 1–18 17
due to more massive bedrock or a thinner cover of photo-based and field-coupled surveys, vary very
glacial drift (Bovis and Jakob, 1999). In this sense, little. The impact of recent logging on denuda-
East Cap and Sisters Creeks seem to be good exam- tion rates remained constant survey-wise: land-
ples of supply-limited and transport-limited basins, sliding ratios exhibited a three-time decrease as a
respectively (Fig. 3A and B). result of fieldwork coupling because of the
higher number of ‘‘invisible’’ slides in old-
growth forest. Canopies of forest logged more
6. Conclusions than 50 years ago and those of undisturbed forest
have a similar effect on landslide visibility.
(i) This study provided for the first time a system- (v) Large differences between Sisters and East Cap
atic examination of the assumption that the Creeks in landslide density and denudation rates
relation between the numbers of landslides were qualitatively explained via GIS-based topo-
detectable from air-photo inventories and ground graphic analysis. The outcome depicted the two
surveys is the same as that for landslides that are sub-basins as contrasting geomorphic environ-
‘‘non-visible’’ on air photos. In the Capilano ments. Accordingly, the different landslide
River basin, the number of slope failures that activity between the study areas was clarified
remained undetected from air photos and the with reference to the distinction between supply-
relative volume of debris mobilised showed limited (i.e., East Cap Creek) and transport-
great sub-basin variability. In basins like Sisters limited (i.e., Sisters Creek) basins.
Creek, where below threshold scale failures are
disproportionately extensive, the forest canopy
Acknowledgements
does hide an important population of ‘‘not
visible’’ landslides; in others, such as the East
The research was supported by the U.B.C. Chair of
Cap Creek basin, it does not. It follows that, if
South – North Studies (Professor Olav Slaymaker) and
one wants to obtain complete information for
by scholarships awarded to Francesco Brardinoni by
sediment budget evaluation and terrain stability
the Faculty of Science of the Università Cà Foscari di
assessment, supplementary fieldwork is manda-
Venezia and by the Government of Canada (GOCA).
tory in Sisters Creek-like areas. In East Cap
The Greater Vancouver Regional District, through the
Creek-like areas, the supplementary fieldwork
offices of Mr. Derek Bonin and Mr. David Dunkley,
would not be justifiable.
kindly allowed and facilitated access to the restricted
(ii) Maximum area of ‘‘not visible’’ failures in forest
water supply area. Russell White, Liz-Anne Strik, and
older than 50 years was 650 m2, larger than what
Stephanie Sork assisted with fieldwork. June Ryder
has been assumed in previous studies. In recently
and Michael Church provided many insightful sug-
harvested cut-blocks, the value dropped to about
gestions and comments that greatly improved the
150 m2. Factors that have proven to affect
paper. Chris Ayles and Erik Schiefer commented on an
landslide visibility were land use, gully relation
early draft of the paper. We thank David Alexander,
of failure, slope gradient, valley width, slope
Richard Marston, and an anonymous reviewer for their
position, and stream connection. Lower portions
constructive comments on this paper.
of steep, old-growth forested slopes located in
narrow valleys are most likely to hide a relevant
number of slope failures. References
(iii) Fieldwork has demonstrated that gully-related
events have a greater importance than one could Armstrong, J.E., Brown, W.L., 1954. Late Wisconsin marine drift
detect from API. They constitute more than one- and associated sediments of the Lower Fraser Valley, British
third of the missed events, while from air photos, Columbia, Canada. Bulletin of the Geological Society of Amer-
ica 65, 349 – 364.
they accounted for < 10%. Bovis, M.J., Jakob, M., 1999. The role of debris supply conditions
(iv) This case study has demonstrated that forest in predicting debris flow activity. Earth Surface Processes and
management effects, as perceived between air- Landforms 24, 1039 – 1054.
ARTICLE IN PRESS
18 F. Brardinoni et al. / Geomorphology 1321 (2002) 1–18
Brardinoni, F., Hassan, M.A., Slaymaker, O., 2002. Complex mass Roddick, J.A., 1965. Vancouver North, Coquitlam and Pitt Lake
wasting response of drainage basins to forest management in Map Areas, British Columbia. Memoir - Geological Survey of
coastal British Columbia. Geomorphology 49, 109 – 124. Canada, vol. 335. Ottawa, Canada.
British Columbia Ministry of Environment, 2002. Web site http:// Rollerson, T., Millard, T., Jones, C., Trainor, K., Thomson, B.,
srmwww. gov.bc.ca/risc/PUBS/TEECOLO/TerClass. 2001. Predicting Post-Logging Landslide Activity Using Ter-
British Columbia Ministry of Forests, 1995. British Columbia rain Attributes: Coast Mountains, British Columbia. Forest
Forest Practices Code. British Columbia Government, Victoria. Service British Columbia, Technical Report, vol. 011. Nanaimo,
216 pp. 20 pp.
Carrara, A., Guzzetti, F., Cardinali, M., Reichenbach, P., 1999. Use Rood, K.M., 1984. An aerial photograph inventory of the frequency
of GIS technology in the prediction and monitoring of landslide and yield of mass wasting on the Queen Charlotte Islands, Brit-
hazard. Natural Hazards 20 (2 – 3), 117 – 135. ish Columbia. Land Management Report, vol. 34. B.C. Ministry
Carson, M.A., Kirkby, M.J., 1972. Hillslope Form and Processes. of Forests, Victoria. 55 pp.
Cambridge Univ. Press, Cambridge, UK. 475 pp. Schwab, J.W., 1988. Mass wasting impacts to forest lands: forest
Donati, L., Turrini, M.C., 2002. An objective method to rank the management implications, Queen Charlotte Timber Supply
importance of the factors predisposing to landslides with the Area. In: Lousier, J.D., Still, G.W. (Eds.), Degradation of Forest
GIS methodology: application to an area of the Apennines Land: Forest Soil at Risk. Proc. 10th B.C. Soil Science Work-
(Valnerina; Perugia, Italy). Engineering Geology 63 (3 – 4), shop, Feb. 1986. Land Management Report, vol. 56. B.C. Min-
277 – 289. istry of Forests, Victoria, pp. 104 – 115.
Erskine, C.F., 1973. Landslides in the Vicinity of the Fort Randall Sidle, R.C., Pearce, A.J., O’Loughlin, C.L., 1985. Hillslope Stabil-
Reservoir, S. Dakota. Prof. Pap. - Geol. Surv. (U.S.), vol. 675, ity and Land Use. Water Resources Monograph Series, vol. 11.
p. 64. Reston, VA. Am. Geophysical Union, Washington, DC. 140 pp.
Flageollet, J.-C., 1996. The time dimension in the study of mass Smith, R.B., Commandeur, P.R., Ryan, M.W., 1983. Natural reve-
movements. Geomorphology 15, 185 – 190. getation, soil development, and forest growth on the Queen
Greater Vancouver Regional District, 1999. Annex CD to GVRD Charlotte Islands. Progress report. Working Paper 7/83. B.C.
Analysis Report Watershed Management Plan #5. Vancouver, Ministry of Forests, Fish/Forestry Interaction Program, Victoria.
BC. 44 pp.
Holland, S.S., 1976. Landforms of British Columbia, a Physio- Smith, R.B., Commandeur, P.R., Ryan, M.W., 1986. Soils, vegeta-
graphic Outline. Bulletin No. 48. BC Ministry of Mines and tion, and forest growth on landslides and surrounding logged
Petroleum Resources, Victoria. 138 pp. and old growth areas on the Queen Charlotte Islands. Land
León, R.V., Mee, R.W., 2000. Blocking multiple sources of error in Management Report, vol. 41. B.C. Ministry of Forests, Victoria.
small analytic studies. Quality Engineering 12 (4), 497 – 502. 95 pp.
Nistor, C.J., 1996. Temporal patterns in the normal-regime fine- Terlien, M.T.J., van Westen, C.J., van Asch, T.W.J., 1995. Deter-
sediment cascade in Russell Creek Basin, Vancouver Island. ministic modelling in GIS-based landslide hazard assessment.
MS thesis, The University of British Columbia, Vancouver. In: Carrara, A., Guzzetti, F. (Eds.), Geographical Information
236 pp. Systems in Assessing Natural Hazards. Kluwer Academic Pub-
O’Loughlin, C.L., 1972. An investigation of the slope stability of lishing, Dordrecht, The Netherlands, pp. 57 – 77.
the steepland forest soils in the Coast Mountains, southwest Thurber Engineering, 1996. Capilano Reservoir Investigation of
British Columbia. PhD dissertation, The University of British Shoreline Instability. Report to Greater Vancouver Regional Dis-
Columbia, Vancouver. 147 pp. trict, Vancouver, BC. 31 pp.
Pyles, M.R., Froehlich, H.A., 1987. Discussion of ‘‘Rates of land- Wieczorek, G.F., 1984. Preparing a detailed landslide-inventory
sliding as impacted by timber management activities in north- map for hazard evaluation and reduction. Bulletin of the Asso-
western California’’ by M. Wolfe and J. Williams. Bulletin of ciation of Engineering Geologists 21 (3), 337 – 342.
the Association of Engineering Geologists 24 (3), 425 – 431. Yanai, S., Usui, G., 1987. The Evaluation of Slope Stability by
Reid, L.M., Dunne, T., 1996. Rapid Evaluation of Sediment Budg- Tephrochronology in Central Hokkaido, Japan. ANZSLIDE
ets. GeoEcology Paperback. Catena Verlag, Reiskirchen, Ger- 87, 5th International Conference and Field Workshop on Land-
many. 164 pp. slides, Christchurch, New Zealand, Conference Proceedings,
Robison, E.G., Mills, K., Paul, J., Dent, L., Skaugset, A., 1999. pp. 165 – 169.
Storm Impacts and Landslides of 1996: Final Report. Forest Zaruba, Q., Mencl, V., 1969. Landslides and Their Control. Aca-
Practices Technical Report, vol. 4. Oregon Department of For- demia and Elsevier, Prague, Czechoslovakia. 205 pp.
estry, Corvallis. 145 pp.