Identification and Mitigation of Large Landslide Risks in Europe Advances in Risk Assessment CH
Identification and Mitigation of Large Landslide Risks in Europe Advances in Risk Assessment CH
IMIRILAND PROJECT
Page ii
European Commission
IMIRILAND Project
ACKNOWLEDGEMENTS
Although all care is taken to ensure the integrity and quality of this
publication and the
information herein, no responsibility is assumed by the publishers
nor the author for any damage
to property or persons as a result of operation or use of this
publication and/or the information
contained herein.
POLITECNICO DI TORINO/TECHNICAL
UNIVERSITY OF TURIN (Italy)
Arpa Piemonte—(ARPA)
10134 Torino—Italy
Tel: +39–011–3169335
Fax: +39–011–3169340
e-mail: [email protected]
URL: www.arpa.piemonte.it
10129 Torino—Italy
Tel: +39–011–5644823
Fax: +39–011–5644899
e-mail: [email protected]
URL: http://www.polito.it/ricerca/dipartimenti/distr/
10123 Torino—Italy
Tel: +39–011–530652
Fax: +39–011–530652
e-mail: [email protected]
URL: http://www.igg.cnr.it/
Boulevard Lefebvre 58
75732 Paris—France
Tel: +33–01–40435246
Fax: +33–01–40436516
URL: www.lcpc.fr/
Tel: +33–04–72143215
Fax: +33–04–72143035
e-mail: [email protected]
URL: http://www.equipement.gouv.fr/
08034 Barcelona—Spain
Tel: +34–93–4016511
Fax: +34–93–4017251
e-mail: [email protected]
URL: http://www.etcg.upc.es/
1040 Wien—Austria
Tel: +43–1–5880120–319
Fax: +43–1–5880120–399
e-mail: [email protected]
URL: http://www.tuwien.ac.at/
1015 Lausanne—Switzerland
Tel: +41–21–6932312
Fax: +41–21–6934153
e-mail: [email protected]
URL: http://lmswww.epfl.ch/
Contents
Preface xiii
1. Introduction 1
M.Castelli, & F.Forlati
3.1 Introduction 13
3.2 Definitions 15
4.1 Introduction 45
5.1 Introduction 89
References 307
Preface
The disastrous events of the last few years have also clearly shown
that attainable security in natural hazard management remains
limited despite effective preventive measures. This is largely due to
the vulnerability and the sensitivity of our infrastructures (Fig. 2), as
well as our society and the environment which we live in. In addition,
it is more and more questionable to distinguish between natural risk
and that induced by human activity.
The editors:
Christophe Bonnard, Ferruccio Forlati, Claudio Scavia.
Page xviii
1
Introduction
M.Castelli1, F.Forlati2
1Politecnico di Torino; 2Arpa Piemonte
Table 1.1. Some major landslides which have occurred over the past
centuries in the Alps.
Since then, many technical and scientific papers have been written
on this topic (i.e. Brabb, 1984; Brand, 1988; Bunce et al., 1997;
Carrara, 1992; Carrara et al., 1992; Cherubini et al., 1993; Cruden &
Fell, 1997; Cruden & Varnes, 1994; Einstein, 1988; Fell, 1994; Fell et
al., 2000; Finlay & Fell, 1997; Hungr, 1981; Hungr et al., 1993;
Hutchinson, 1992; Leroi, 1997; Mulder, 1991).
Finally, vulnerability and risk analyses have been developed for the
same landslides considering short and long-term perspectives, direct
and indirect consequences, as well as technical and social impacts
(Fig. 1.2). The combination of the results obtained in relation to
hazard analysis and through the risk approach will allow the
development of a new practical and quantified risk assessment
program which will be applied to several sites.
A general scheme for data description was created and compiled for
each site (Table 1.2). Each partner was then able to access data and
to manage them in a homogeneous manner. The objective was also
to check the reliability of the available data and to choose which
were to be considered as significant for the risk analysis.
Hydroelectric power
plant
Valley downstream
International Torino-
Modane railway
Valley downstream
Prequartera and
Campioli villages
Communities
downstream
Valley downstream
Town of Lanquart
Valley downstream
Note: The volumes indicated above refer to the modelled active landslide
masses, whereas the areas indicated refer to the globally involved instability
phenomena.
1. Introduction
2
The meaning of risk assessment related to large
landslides
Ch.Bonnard
Ecole Polytechnique Fédérale de Lausanne
is often limited to a localised zone near the scarp or at the toe. In the
case of a marked reactivation, the movement of the mass may
evolve into a rockfall phenomenon or a debris slide near the toe,
which, according to the size of the blocks, might be able to dam the
valley and cause flooding first upstream, and then downstream if a
sudden release of stored water occurs following a dam break event.
In this case, the impact of the phenomenon is much larger than its
own surface extension and may affect hundreds of km2 (see the
Sechilienne landslide commented on in Chapter 9). Such a case
occurred in 1993 in southern Ecuador at La Josefina as presented in
Chap. 10.2 below.
3
A key approach: the IMIRILAND project method
3.1 INTRODUCTION
– Where?
– How often?
Action
5. Estimation of vulnerabilities
– Calculation of the risk from the hazard, elements at risk and the
vulnerability of the elements at risk.
3.2 DEFINITIONS
Total Physical risk is the sum of the respective risks for all the
scenarios, whenever relevant.
Total Economic risk is the sum of the respective risks for all the
scenarios, whenever relevant.
Total Social risk is the sum of the respective risks for all the
scenarios, whenever relevant.
Total Environmental risk is the sum of the respective risks for all the
scenarios, whenever relevant.
Page 16
1. Hazard analysis
2. Consequence analysis
3. Risk computation.
– Danger characterisation,
– Process intensity,
– Occurrence probability.
Page 17
—Comparison of
different
measurement
typologies.
– geomorphological constraints;
– geostructural constraints;
Morphostructural model
In recent years, investigations on large landslides have more
frequently concerned the whole slope area or even the entire valley
on which they are located, in order to discover common (ubiquitous)
elements that could have driven the gravitational instability.
When these conditions are satisfied, it can be assumed that the “on
site” and “around site” interpretations are consistent. In this case,
there is the possibility of extrapolating the “on site” interpretations to
broader areas or vice-versa.
Ancient cadastral maps, defining the limits of the plots of land which
indicate the location of old houses, also help establish past rates of
movement, as in the case of the La Frasse landslide in Switzerland,
where ancient maps date back to 1768 (Bonnard, 1984).
• Rockfall: 1940, Oct. 1971, Apr. 1977, Oct. 2000, Jun. 2002
• Movement of the entire mass (or of a large part of it): 312, 843,
1816, 2000
Numerical triggering models can indeed represent the rock slope as:
The most widely used numerical methods are the limit equilibrium
method (LEM), the finite element method (FEM), the boundary
element method (BEM) and the finite difference method (FDM).
– Empirical methods
– Analytical methods
– Path profile
– Deposition area
– Mass balance
– Granular methods
• The Particle Flow Code (PFC) from the ITASCA (1999a and 1999b)
Consulting Group in cooperation with FLAC (Itasca) as the Ball Wall
model and stand alone as the All Ball model;
In contrast to the Ball Wall model, in the All Ball model, the relatively
stationary bedrock, the failure mechanism of the slope and the
detachment mechanism are modeled by balls. This allows the all in
one calculation of failure mechanisms, detachment and run-out. For
calculating the failure mechanism, however, quite a demanding
calibration of materials is necessary; this is particularly due to the
fact that the modeling of the bedrock by balls determines a rather
rough surface. When comparing the volume of the detached rock
with FLAC, both methods correspond closely.
One limitation of PFC is due to the fact that it is not able to model the
influence of water (e.g. pore pressure). On the other hand, as the
balls are supposed to model simple rock blocks, and as the number
of balls is limited by the computation capacity, a certain restriction
exists in the modeling of large landslides of dozens of millions of m3.
Finally, the All Ball model allows an all in one calculation but the
computational effort required obliges the modeling of the detached
mass with less and bigger balls than those used in the Ball Wall
model.
The main data used as inputs for the DAN code are the slope profile
and the top profile of the initial mass. In particular, to perform a
three-dimensional analysis, the width of the path must be determined
at each path profile point by assuming that the entire envelope
considered is enveloped by the moving mass.
The final outputs are: run-out distance and run-up on the opposite
slope, landslide velocity and expiring time, limits of propagation and
depth of the displaced mass once it stops.
The main limitation of the Hungr model is due to the fact that it is
approximate, as it involves the reduction of a complex and
heterogeneous three-dimensional problem into an extremely simple
formulation. The simplicity of the model is, however, an advantage in
making possible an immediate and rapid numerical simulation of real
cases.
– Rockfall: 1940, Oct. 1971, Apr. 1977, Oct. 2000, Jun. 2002
– Movement of the entire mass (or of a large part of it): 312, 843,
1816, 2000.
Apr. 1977
Oct. 2000
Jun. 2002
312
1×106 m3 30–200
Geological— • • • •••
structural
Geomorphological•• •• • •••
Monitoring •• — • •
investigation
Geomechanical •• • — ••
modeling:
triggering
3.6.1 Introduction
Once the hazard scenarios have been identified, each one being
characterised by an area in which the distribution of the process
intensity is known, as well as the corresponding probability of
Page 35
• vulnerability evaluation,
Dead persons 1
W=[Rm(Mm−Em)]*Nab+Ned*Ced+Cstr+Cmorf
In the following tables (Tables 3.6 and 3.7), relative values related to
some elements at risk categories are shown. For each element at
risk, a “relative value index” for assets, economic activity interruption
and environment has been indicated. The human lives relative
“value” has been applied on the basis of the persons involved (the
ranges can vary according to the area studied).
Page 37
Residential areas 4 1 1
Extra-municipal infrastructures 4 4 2
and plants
Tourist accommodation— 3 3 1
buildings
Strategic lifelines 3 4 1
Tourist accommodation—camp 2 2 1
grounds
Local services 2 2 1
Secondary roads 1 2 1
Secondary lifelines 1 2 1
hospitals 3 4 1
schools 2 4 1
garbage dumps 2 3 4
2÷9 2
10÷19 3
>20 4
“Although the state of the art for identifying the elements at risk and their
characteristics is relatively well developed, the state of the art for the
assessment of vulnerability is in general relatively primitive.
Thus, the main criteria for determining the value of the environmental
vulnerability coefficients are:
V=f (T, I)
In the first case (1), 5 classes of percentage loss (0, 0.25, 0.5, 0.75,
1) have been applied to each vulnerability category. Through
detailed studies it could be possible, however, to attribute various
weights at various percentage classes based on “utility curves” [e.g.:
in the road interruption case—economic vulnerability—some studies
have evidenced a logarithmic type curve (Fig. 3.7a), as the
vulnerability diminishes with the passing of time for the opening of
the alternative roads (e.g.: variable indices–0, 0.5, 0.8, 0.95, 1–
according to elapsed time)].
C=VE×V
R=H×C
R=P×C
4
The Cassas landslide
4.1 INTRODUCTION
Between 12th and 14th June 1957 in the Susa Valley, a huge flood
caused great damage to the entire valley. On that occasion, the
western area of the Cassas landslide had a paroxystic phase:
millions of m3 of rock mass slid and caused a debris flow that
reached the valley floor. Re-activations of the landslide were
recorded until the middle Sixties (Peretti, 1967, 1969; Brovero et al.,
1996a) (Fig. 4.2).
Figure 4.2. General overview of the Cassas landslide in a photo from
1965, at the conclusion of the paroxystic phase of the period 1950–
1965 (Brovero et al., 1996a).
Page 47
Figure 4.4. Location of the weather stations in the upper part of the
Susa Valley used for climatic analysis.
4.2.2 Hydrological data
The main historical floods of the Dora Riparia have been recorded
since the 18th century; they are sometimes associated with
disastrous events (Brovero et al., 1996a; Fratianni & Motta, 2002)
listed as follows:
In the lower part of landslide, two springs run (estimated flow on July
1997: from 2×10−4 to 8×10−4 m3/s). There are also some springs
along the scarp that cut across the debris cover between the
altitudes 1’350 and 1’550 m a.s.l. (Brovero et al., 1996a; Paro,
1997).
Page 50
Figure 4.8. View of the Susa Valley right slope, from Salbertrand to
Serre la Voute, with indications of the main morphological elements:
1. Cassas landslide main scarp, June 1957 event; 2. main scraps of
gravitational phenomena; 3. secondary scraps of gravitational origin;
4. scraps by fluvial-torrential erosion; 5. landslide accumulations; 6.
main rock outcrops; 7. actual fluvial-torrential deposits; 8. alluvial
fans; 9. fluvial torrential terraced deposits; 10. talus; 11. landslide
scar and transition area; 12. area characterised by morphological
features related to deep-seated gravitational deformations.
In the valley floor, some alluvial fans are observed, such as the large
fan formed by the Gran Plenei stream (right tributary of the Dora
Riparia river). Some of these, of smaller dimensions, are located at
the feet of landslide accumulations. Among those, the fan placed at
the foot of the “1957 Cassas”, was involved in the paroxystic event of
June 1957. Successively, this fan was modified by the construction
of the service station for the A32 Motorway (C in Fig. 4.15), which
also rests on the alluvial deposits of the Dora Riparia river.
Figure 4.10. 1934 map at scale of 1:25’000 (IGM sheet n.54, II SE,
Oulx). The shaded zones indicate the main instability features.
– during the event of June 1957, the upper part of the Cassas slope
(L in Fig. 4.11) underwent a total displacement of about 250 m,
mostly by sliding. A sub-vertical scarp, tens of meters high, was
formed on the right flank from which a part of the existing debris
flowed (S in Fig. 4.12).
– prevalent calcschists;
– carbonatic massive calcschists with interlayered, cm-mm thick
micaceous schists;
The available mesostructural data of the upper and lower part of the
landslide permitted the observation of six main fracture (s.l.) systems
(Fig. 4.16a):
– the main schistosity (S1) dips to the S-SE with an inclination of 40–
50°. In some places, the schistosity gently dips to the N; in this case
it is labelled S2;
Class V IV III II I
SMR—RMR+(F1·F2·F3)+F4
The RMR is the Rock Mass Rating (Bieniawski, 1974); F1, F2, F3
and F4 are corrective factors that define the type of slope failure
resulting from geometric relationships between discontinuity systems
and the slope (Romana, 1990).
Page 61
The collected data are summarised in Table 4.3 and in the relative
histogram (Fig. 4.18), where the RMR, SMR and RQD (rock quality
designation) are reported.
The comparison between the SMR value and the Romana table put
the rock mass of the Cassas landslide between classes of quality V
and III.
13 Calcschists 82 48 20.5
14 Carbonatic Calcschists 75 39 45
28 Calcschists 32 44 50
Figure 4.19. Histogram of the geomechanical data reported in Table
3.
The amount of available regional data was high; those used here
have been critically revised, since they were originally produced for
very different purposes from those of the IMIRILAND project.
the presence of drag folds or tilting of rock blocks. Since the s1 has a
relatively constant orientation, it can be easily extrapolated along
strike and at depth for at least several hundreds of meters.
The frequencies for the classes of lineament length are more widely
distributed and permit the identification of four lineament systems
(see Fig. 4.24):
Five drilling surveys have been carried out until the present:
The rock wall at the left flank of the landslide is monitored by 12 joint
gauges and by 1 geophone.
To the West of the Cassas landslide a stable area may be seen (G)
where superficial landslides are located; the monitoring data of the
sub-vertical rock wall above the motorway service station do not
indicate movement.
– N-S minor faults (F3) widespread all around the landslide; they are
sub-parallel to the trenches and scarps in the middle and upper parts
of the landslide. The most extensive scarp among those observed in
the landslide body separates Sector B from Sector C.
Page 70
The continued movements measured over the last ten years in the
upper part of the slope emphasised an acceleration starting from
June 2001 (approximately 20 to 80 mm/y), caused by the heavy
rainfalls of October 2000 and by the winter/spring snow melting.
Page 74
– the water quantity would condition the typology and the velocity of
the debris movement in the middle-lower area and subsequently the
ways of accumulation in the valley floor;
– blocks or rock masses falling from the main scarp of the Cassas
landslide (velocity of some m/s) with an accumulation zone at the
scarp foot or at the slope. Each block has a maximum volume of
some tens of m3; the rockfall has a maximum volume of 105 m3;
– collapse of parts of the slope with debris slide and flow (variable
velocity from some m/min to some tens of m/month) with
accumulation partly on the slope, mostly on the valley floor.
Maximum volume inferior to 10×106 m3 (indicative volume useful for
the delimitation of the area affected). It is excluded, at present, that
the phenomenon could evolve as a debris avalanche;
– for the complex debris slide and flow phenomena (<107 m3):
morphologically based hypothesis for volumes of about 1÷2×106 m3;
– for the complex debris slide and flow phenomena (>107 m3):
morphologically based hypothesis for volumes of about 1.5×107 m3.
Geophones
Extensometers
Table 4.6. Debris slide & flow scenario (< 107 m3): element at risk
identification.
Element at risk Name Notes
Forests (private and Bosco Chapel Area affected: 0.23 km2 Persons
public properties) (Piccolo bosco) affected: 2 (occasional presence of
monitoring operators)
Table 4.7. Debris slide & flow scenario (1÷1.5×10 m3): element at
risk identification.
Persons affected: 10
Geophones
Extensometers
In the debris slide and flow scenarios, a maximum intensity has been
given to the recognised processes, as the involved volumes are
huge (millions or tens of millions m3). Considering the elements at
risk involved in the gravitational phenomenon (buildings,
infrastructures, roads and railway), a vulnerability value of 1 (100%)
was considered: in the area affected (failure, depletion and
accumulation areas) the destruction is considered total. Also, at the
valley floor, when the velocity tends to zero, the accumulation can be
so big that all the assets can be completely destroyed, heavily
damaged or buried and, therefore, permanently useless. The velocity
and thickness distribution in the accumulation area depend on the
initial velocity and on viscosity (water quantity in the material in
movement).
As shown in Figures 4.35, 4.36, 4.37 and 4.38, for the debris flow
and slide scenarios, the major risks are connected to the landslide
body on the slope: this situation is possible because the landslide on
the slope is active and, therefore, with high frequency (≥1).For the
correct considerations of the real risk, only that which exists in the
accumulation area on the valley floor must be evaluated.
Page 81
Figure 4.35. Map obtained by G.I.S. techniques concerning the
physical risk for the 3 scenarios analysed.
In the two debris slide and flow scenarios, the major values of
physical risk are, first, those of the service station “Gran Bosco” and,
second, those of the motorway and the areas classified as
productive areas in the land planning.
Page 82
Figure 4.36. Map obtained by G.I.S. techniques concerning the
economic risk for the 3 scenarios analysed.
As it can be observed in the images, the areas with major risk values
are the A32 motorway (Turin-Frejus tunnel) and the service station.
Particularly important is the Dora Riparia river, as well as the intrinsic
value of the river, its waters are captured downstream of the Cassas
landslide by a power plant.
Page 83
Figure 4.37. Map obtained by G.I.S. techniques concerning the
environmental risk for the 3 scenarios analysed.
The main areas with major risk values are, obviously, the slope
areas (in the Natural Park “Gran Bosco di Salbertrand”) and the
valley floor areas where the Dora R. river runs.
Page 84
Figure 4.38. Map obtained by G.I.S. techniques concerning the
social risk for the 3 scenarios analysed.
Concerning the rockfall scenario, the risk values are not zero as
periodically 2 monitoring system operators are at the landslide.
For the debris slide and flow scenarios, the areas with elevated
social risk are the Turin-Frejus motorway tunnel and the related
service station. A detailed analysis shows how traffic variation due to
time and seasons leads to different considerations about risk
evaluation.
Page 85
As already detailed for the risk analysis, the high risk index obtained
for the slope sectors involved in the active phenomenon must not be
taken into account. This derives from the fact that for the Cassas
landslide, a high occurrence probability value has been considered
and this leads to a high risk value, even if linked to a low value of the
affected assets. Therefore, for a correct risk analysis, it is necessary
to consider the affected areas at the valley floor whose occurrence
probability values are 10 times inferior to those on the slope and that
reach high risk values for the important consequences evaluated.
Concerning the assets evaluation, the analyses are not detailed, but
the relative values obtained by the methodologies have been
applied. Probably, to obtain a correct framework of the risk situations
related specifically to the Cassas area, careful social and economic
evaluation should be done on the basis of the considerations of the
next paragraph (§ “Indirect risk and regional impacts”).
– How much material will deposit on the slope and how much will
reach the valley floor?
This will depend on the water quantity and on the related viscosity
and velocity of the mass in movement.
The motorway and the railway play an important role in the economic
development of these mountain areas, as they are linked to tourist
resorts. The Susa Valley and the Chisone Valley are two important
ski resorts connected to other ski stations in France: international ski
competitions take place in this area every year and at this time new
infrastructures are developing for the winter Olympic games of 2006.
The motorway and the service station “Gran Bosco” at the foot of the
Cassas landslide raise some problems related to risk mitigation,
problems never taken into account either in relation to the
Salbertrand village or in relation to other infrastructures in this area.
The S.I.T.A.F. (the society that manages the motorway) has been
studying the best protection both for the motorway and the service
station since the end of the ’80s of the last century.
5
The Rosone landslide
5.1 INTRODUCTION
The entire southern side of the ridge bounded by the Orco and
Piantonetto Rivers, in the Province of Turin, has been undergoing a
slow process of deep-seated gravitational slope deformation
(DSGSD). This deformation process, a vast landslide movement, is
historically known as the “Rosone landslide” because it frequently
involved the small village of Rosone located at the toe of the slope.
At present, this phenomenon also involves the 99 MWh hydroelectric
power plant of the Electricity Agency (AEM) of the city of Turin (Fig.
5.1). Coming from the Ceresole Reale dam, the water reaches the
AEM facilities through a 17 km long pipe after spanning the entire
length of the gravitational deformation, where it falls towards the
power plant via a penstock with a drop of 813 m. At the toe of the
slope, the only National Road (No. 460) from Turin to Ceresole
Reale is located.
Between 1957 and the early 1960’s, the movements accelerated and
phenomena similar to those recorded in 1953 were observed, though
less severe. Further phenomena were recorded in the fall of 1963
and in the spring of 1964 and 1969, 1988, 1993.
May: the highest monthly rainfall value (160 mm), the greatest
average number of rainy days (12 days) and the highest cumulative
rainfall values are reached in three or four days (up to nearly 195
mm).
It is of interest to point out that two perpetual springs are found near
Perebella and Bertodasco and some temporary springs at the toe of
the slope and near the right flank of the Bertodasco sector. The flow
of water in the slope occurs through the main discontinuities, which
are generally open and persistent. Considering that the upper part of
the slope acts as a reservoir, the largest part of the groundwater is
cumulated when rainfalls coincide with the melting of the snow cover.
The Orco Valley is located in the central part of the Gran Paradiso
Massif (Fig. 5.5). The geological unit belongs to the Pennine domain
and it consists of a composite crystalline basement and a Permo-
Liassic cover, locally preserved in the peripheral areas. The Gran
Paradiso Massif consists of three different complexes (Compagnoni
et al., 1974):
All of these units are locally covered by glacial deposits and they are
cut by two main joint sets, mutually orthogonal, corresponding to E-
W and N-S striking sub-vertical faults.
In this area, the deep seated slope deformations have reached the
maximum evolution level. The gravitational processes, having lasted
here for a long time, have erased any evidence of glacial
morphology and strongly modified the profile of the slope. In ancient
times, the slope failed and rockslides occurred.
Figure 5.6. Aerial view of the high part of the Ronchi sector. The
boundary main scarp, on the right, as well as the intense disruption
of rock mass, is depicted.
Page 96
Figure 5.8. Up-hill facing scarps in the upper part of the Perebella
sector.
Page 97
Figure 5.9. Aerial view of the medium part of the Perebella sector:
cross-cutting relations of the main steeply dipping joint systems are
well exposed.
– Zone A includes the upper part of the Bertodasco sector, where the
movements of the slope (as a whole) and related morphological
evidences, are poorly defined. Moreover, although disrupted, this
zone still preserves its original structural features. This suggests the
presence of minor translational movements, as also confirmed by the
AEM monitoring system (Forlati et al., 1993).
The eastern lower part of Zone A shows more prominent
movements, as indicated by opened fractures, and is affected by
undulations and local bulges. The rocks here are split into blocks,
although still aligned along the main joint systems.
The movements in the upper part of this zone have caused, during
the last fifty years, the break-down of part of the Bertodasco village,
which is now completely abandoned.
Figure 5.10. View of the Bertodasco sector where the origin and
deposit zones of debris flows are well exposed.
Two main joint systems have been recognised (Fig. 5.11). They are
the K1 and K2 steeply dipping (75°–85°) systems, which strike N10
and N90–N100 respectively (Fig. 5.9).
Two other minor joint systems are also present. They are the K3
system, that dips 70° to N240 and the K4 system, that is steeply
dipping and strikes N40 (Geoengineering, 1984).
Page 99
Figure 5.12. Structural sketch map of the upper Orco valley from
1:10.000 field mapping and 3D digital photo-interpretation.
Since 1985, attention has been focused on the whole slope, with
special emphasis on the Bertodasco sector: a systematic campaign
has been conducted on the entire rock mass involved in the sliding
activities.
Page 101
3.55–31.3 m: rock
with sampling
>80%
Failure surfaces
intercepted at 3.1 m
depth
Water: 26.2 m
22.9–39.2 m: rock
with sampling
>80%
Failure surfaces
intercepted at 22.6
m depth
57.7–81 m: rock
with sampling
>80%
21.5–30.4 m: rock
with sampling
>80%
Failure surfaces
intercepted at 20.1
m depth
Water: 25.4 m
Lugeon
permeability test:
depth [m]/UL
extrapolated to 10
bars
62.95–64.45/68–78
72.70–74.70/10–15
Page 105
A3 12/12/19915/10/1993 71.31Closed — —
5.5.4.4 Scenarios
On the basis of geomorphological observation and field investigation,
the rock slope foot, over deepened by glacial erosion and weakened
by severe stress concentration, does not seem of particularly
efficient significance in contrasting a global movement of the slope.
As the deformation affecting the middle part of the slope will spread
througt the basal part, a quite defined failure surface might be
formed, leading eventually to a huge landslide and possibly severe
rockfall and/or rock avalanche phenomena.
The rock volume involved by the inferred failure surface coupling its
position with the topographic surface, obtained from DEM, has been
also calculated. The side boundaries of the moving mass have been
obtained from the geomorphological evidence (§ 5.5.2.2): hence,
vertical joints have been introduced to reproduce the boundaries of
Zones C and B (Figure 2). The estimated volume is about
11’000’000 m3, in good agreement with the volume evaluated on the
basis of the geomorphological studies (§ 5.5.4.4). In Figure 5.21, the
3D representation of the unstable volume obtained as described
above is shown.
Page 110
The following step was the assessment of the value of the creep
coefficient C. Activating the Map3D creep option, a simulation was
carried out with an estimated value of C and the results were
compared with the measured displacements. If these values were
not comparable, more simulations were carried out in order to find a
suitable C. The final set of mechanical parameters used for the
numerical analyses is shown in Table 5.4.
The obtained results are very far from giving the probability of failure
needed for the risk analysis, which is to be obtained through a
probabilistic procedure. For this purpose however the variation of
pore pressure in time and the decay of rock strength parameters in
time should be known. These data are seldom known and are not
known in the case of the Rosone landslide.
Page 111
Equivalent continuous
Failure surface
Friction angle, φ 25 o
– Scenario C.
– Scenario C+B.
– cenario C+B+A.
C 1’864’000 2’423’000
ISMES methodology
– the dynamic friction coefficient (µ) and the ratio of involved mass
(M) to aerodynamic strength coefficient (D) are considered constants
and their value is defined through the back-analyses of historical
cases.
The Perla model was also applied to estimate the maximum lateral
travel distance (Rl), this aspect supposes that, once the bottom of
the valley is reached, the river direction is immediately followed by
the moving mass (Table 5.7).
c 1776 2245
Hungr method
Rheology Frictional
(unit weight)
The analyses carried out underline that in the triggering case the
displaced mass would determine the river damming with disastrous
consequences.
698.0 4.70
723.2 12.50
743.4 19.71
759.4 25.27
775.0 28.27
788.9 30.69
801.2 31.82
812.4 31.77
823.0 30.65
833.0 28.49
842.6 25.58
851.4 22.38
859.5 18.20
867.4 12.86
874.5 6.70
880.9 1.68
Table 5.11. DAN output data.
Final=2428449.50
Final=103757.18
From the historical data analysis (see Fig. 5.4), a frequency inferior
to 10 years for both rockfall/debris flow phenomena and the
Bertodasco sector instability is found. The latest signs were seen
about 40 years ago. This indicates a relationship among these
phenomena: the movements in Zone C cause the blockfall from the
scarp over the National road 460. This situation is particularly
apparent during heavy rainfalls.
phases; the last ones could cause sudden collapses of the slope
with the development of rock avalanche phenomena.
The risk analyses of the scenarios are based on land planning of the
Locana Municipality and on the landslide run-out areas obtained
using the ISMES method (Friz & Pinelli, 1993).
Persons affected: 0
(R.D.445/1908).
32 buildings of
new Rosone
community
Persons affected:
88
Aghettini Buildings
affected: 8
Persons affected:
12
(R.D.445/1908). 62
buildings of new
Rosone community.
Persons affected:
186
Persons affected: 12
Persons affected: 18
The results of the risk assessment are shown in the following figures
(Fig. 5.26, 5.27, 5.28 and 5.29) obtained by GIS layout and
described in the captions.
Page 121
Page 122
Figure 5.26. Maps by GIS of physical risk for the three rock
avalanche scenarios. The maps underline the areas with highest risk
value, i.e. the National road no. 450 and part of the new Rosone
community. Particularly relevant is the run-out area of the A+B+C
scenario: part of the accumulation involves the hydroelectric power
plant, part of the failure area involves the penstocks.
Page 123
Page 124
Page 125
Figure 5.27. Maps by GIS of economic risk for the three rock
avalanche scenarios. The maps underline the areas with highest risk
value, i.e. the National road no. 450 and the lifelines (electric power
line “Valle Locana” and secondary electric power line for the funicular
service). From an economic point of view the electric power plant of
Rosone and the related infrastructures are very important, as they
supply energy to a large part of Turin.
Page 126
Page 127
Page 128
Figure 5.28. Maps by GIS of environmental risk for the three rock
avalanche scenarios. The maps underline the areas with highest risk
value, i.e. the wooded areas and the valley floor. The Orco river is
completely buried by rock avalanche accumulation, therefore, the
risk value of the affected zones has no meaning.
Page 129
Page 130
Page 131
Figure 5.29. Maps by GIS of social risk for the three rock avalanche
scenarios. The maps underline the areas with highest risk value, i.e.
the urban areas of the new Rosone and Casetti communities. Also
the National road no. 450 involved in the rock avalanche
phenomenon has high risk values, being characterised by intense
traffic during the tourist season.
Figure 5.31. Maps by GIS of total risk. The maps underline the areas
with highest risk value for each category physical, economic,
environmental and social.
Page 137
6
The Oselitzenbach landslide
6.1 INTRODUCTION
The area on the northern slopes of the so-called “Carnic Alps” in the
southern part of Austria’s Carinthia country is affected by numerous
mass movements.
The risks concern the destruction of the Naßfeld Road (National road
B 90, crossing the Alps to Italy and connecting to a large ski area),
as well as the possible spreading out of landslide material by the
Oselitzenbach torrent (which occurred in part after heavy rainfalls
and a consecutive flood in 1983).
The main scarp of the whole Reppwand landslide (Fig. 6.2 and 6.3)
shows a typical profile of the Southern Alpine sedimentary
sequences from Devonian sandy—marly calcareous sediments up to
lower Triassic carbonates.
Figure 6.1. Geological sketch of Austria with the location of the site.
Figure 6.3. Panoramic view of the area (Kahler & Prey, 1963).
Especially the toe zone of the Reppwand landslide is still quite active
and is responsible for debris-generating slope failures and the
destruction of the Naßfeld Road. After sustained regional rainstorms
in Sept. 1983 with numerous embankment failures, an extensive
project of construction work and research was initiated along the
Oselitzenbach torrent and the sagging slope above.
Thus, the area of the eastern toe of the Reppwand landslide has
been the object of numerical investigations as the displacements of
this area are the largest and are, therefore, intensively monitored.
During World War I, a military road was constructed. Since then, the
movements along the road have been known, but have not been
investigated in detail.
In the seventies, the skiing area at the Naßfeld Pass was opened
and then enlarged and the road was upgraded to a primary road,
connecting Austria and Italy.
There are three levels of springs: The highest level lies between an
elevation of 1150 m and 1170 m. At this level, there are two lakes,
the “Great Bodensee” and the “Small Bodensee”, also there are
springs west of the lake “Great Bodensee” and water outlets in the
water logging zone east of the lake “Small Bodensee”.
The area of the second level lies 20 m deeper (1130–1150 m), in the
north of the lake “Great Bodensee”.
At the third level, there are many selective water outlets below the
Naßfeld Road (at an elevation of 700–800 m).
In the main valley of the Gailtal, near the site, the continental tectonic
divide of the “Periadriatic lineament” passes through (dividing the
Northern from the Southern Alps from the Piemont to Slovenia). This
leads to parallel faults of the periadriatic lineament (North West-
South East), to some Riedel shears (diagonal) and extension joints
(A-C- Joints, North-south) (see Fig. 6.7).
Above the main scarp of the mass, there are outcrops of massive
limestones with steep slopes. Within the deformed mass, there are
many small trenches (Fig. 6.8). Within the sagging mass, there are
many secondary scarps, forming different elements of the deformed
mass. As shown in § 6.3.3, there are different movement rates in the
various sections.
In the Reppwand (main scarp), the general bedding has a dip of 30–
50° with a dip direction of 200° to the south. The bedding planes are
usually flat, partly polished and may contain mylonitic layers. The
thickness of the bedding is in the range of cm-dm.
The main joint system dips 80° to the SW (dip direction 230°), is
medium rough and has an average trace length of 5–10 m. Some
clay seams may reach a thickness up to several cm.
A second joint system dips 75° to the NW (dip direction 310°), has
usually rough joint surfaces and can be covered by mm thick seams
of clay.
Page 147
The main fault zones are (sub-) parallel to the Periadriatic lineament
(carving the main valley of the Gail River), striking WNW.
Figure 6.13. Joint diagrams in the wider area (Läufer et al., 1997);
the location of the landslide is indicated by a red ellipse.
Page 148
– Naßfeld series: the rocks are heavily disturbed, therefore, they can
be defined as a cohesive soil
– 1 inclinometer
Figure 6.14. Location of the monitoring instruments.
Page 149
– The most active zone in both periods is the upper area of the
Quellenbach landslide, which, however, shows a stabilisation as well
in Period B from 10–15 cm/y down to 7–15 cm/y. Figure. 15 shows
the map of the movements, which were measured between 1991
and 2000.
– The lowest layer does not show any displacements at the moment
although it is influenced by the Reppwand landslide.
2. Bodensee landslide
Figure 6.18. Sliding planes of the possible scenarios at the toe zone
of the landslide.
Scenario 2: The run out of the Bodensee slide will bury the
Oselitzenbach ravine by 3’375’000 m3 of debris flowing over the
Oselitzenbach debris cone, burying 200’000 m2 (Fig. 6.20).
DRAC model
The values selected for the material properties are listed in Table
6.2. As mentioned above, the Naßfeld schists were divided into three
layers due to their different degree of fracturation.
Naßfeld schists
FLAC3D model
Model 1
The blocks of the Naßfeld schists and of the Hochwipfel schists are
homogeneous. The mesh of the model (Fig. 6.29) is based on the
digital surface model of the BEV (Fig. 6.24). The grid distance is 25
meters. A Mohr-Coulomb constitutive model was investigated using
the following material properties:
2,5 0,2 20 14
9 0,1 40 1e3
The density of the material is 2500 kg/m3. Figure 6.29 shows the
distribution of the shear strain rate and the displacement vectors. In
general, the displacements decrease continuously with depth
resulting from creep (slope sagging). The distribution of the shear
strain rate indicates a zone of maximum shear strain rate in a certain
depth. Below this zone, displacements are zero, above they have a
value increasing to the surface. Thus, the zone of maximum shear
strain rate is the “sliding” zone.
Model 2
The main objective of the analyses with DRAC and FLAC3D was the
comparison of the different computer codes for the back-analysis of
the displacements measured during 1991 and 2000. DRAC and
FLAC are not directly comparable, as DRAC doesn’t provide any
plasticity models (e.g. Mohr-Coulomb). An elastic analysis will not
present any failure mechanism, but it will provide some insight into
the landslide behaviour, as areas with large displacements or tensile
stresses may be potential zones of failure.
Page 162
Generally PFC is not yet able to model the influence of water (e.g.
pore pressure) on the runout, whereas DAN needs an estimate of
the detached rock mass and of the run-out direction (Hungr, 1995).
An estimation of the detached rock mass is also needed for the PFC
Ball Wall model, which was provided by a FLAC3D investigation (see
§ 6.3.5.1). Thus, the methods described should be used in
combination, which makes a comprehensive assessment of the
runout as close to reality as possible.
In the PFC Ball Wall model, the bedrock is simulated by linear (2D)
and planar (3D) elements. In contrast to the All Ball model, where
the relatively stationary bedrock is modelled by balls as well in order
to model the failure mechanism of the slope and the detachment
mechanism also, in the Ball Wall model only the detached rock mass
is modelled by balls. Therefore, in the Ball Wall model, an estimate
of the failure mechanism of the slope and of the detachment
mechanism is needed as an input parameter. Consequently, in the
Ball Wall model, the detached mass can be modelled with the help of
more and smaller balls with the same computational effort. One goal
of the investigations of WP4, therefore, was to compare the two
different approaches.
The PFC Ball Wall model offers the possibility of making use of
know-how related to run-out relevant resistances (factors of
restitution, absorption, friction, etc.) applied in rock fall programs
and, consequently, makes a realistic calculation of the run-out
possible.
Figure 6.34. Deposit in the valley after the rock mass fall.
D1 346
D2 397
D3 573
D4 556
W2 377
The PFC All Ball model allows the all in one calculation of failure
mechanisms, detachment and run-out. For calculating the failure
mechanism, however, quite a demanding calibration of materials is
necessary. When comparing the volume of the detached rock with
FLAC both methods correspond closely, which verifies the All Ball
model. The surface of the model is rather rough due to the modelling
of the bedrock by balls, which has to be considered at the calibration
of the run-out parameters. This aspect needs further investigation.
Figure 6.38. Initial state of the 3D model.
Page 166
Figure 6.39 and Figure 6.40 show the final state of the All Ball
model: the run-out is described with the help of four horizontal
measurement lines (D1–D4), the path width by the lines W1 and W2
(see Table 6.8).
D1 330
D2 389
D3 373
D4 378
W1 204
W2 364
Page 167
Figure 6.40. Run-out distance (D1–D4) & path width (W1, W2).
The area of the maximum displacement rates (light blue) in the All
Ball model corresponds closely to the direction and width of the run-
out in the Ball Wall model. The Ball Wall model, however, indicates a
far greater travel distance (Fig. 6.37). This is due to the collision of
moving particles with stationary ones (bedrock) in the All Ball model,
with the moving particles losing energy additionally. This aspect also
needs further investigation.
Page 168
Since in the DAN method the run-out direction and the path width
are assigned a priori, in the investigations of WP4 they have been
determined with the aid of a 3D rockfall program called ROTOMAP.
The comparison of the results of the PFC2D Ball Wall model, which
also needs an estimate of the runout direction, and the results of the
DAN model shows that two different run-out directions have been
chosen and that these profile directions have an enormous influence
on the results.
The run-out directions determined by Rotomap correspond very well
to those of the 3D All Ball model in the West, whereas the run-out in
the 3D All Ball model in the East indicates a smaller width (Fig.
6.44). The same applies to the 3D Ball Wall model (Fig. 6.43).
Page 169
– the FLAC3D simulation as the basis of the DAN and of the 3D Ball
Wall model give approximately the same detached rock volume as
the All Ball model,
(244’300–571’300 depending
on the shape factor assumed)
470–833
Maximum 29 4 21
travel velocity
[m/s]
Page 171
– the run-out widths and the affected areas obtained by the models
correspond more or less and,
These data mean that the Nassfeld road will be destroyed or buried
with a maximum height of 7 m by the run-out in case of a slope
failure over a length of 400 m. The Oselitzenbach torrent will be
dammed up as well over a length of 460 m with a maximum height of
14 m, thus endangering the villages of Tröpolach and Watsching,
including the adjoining agricultural and forest areas by debris flows.
Occurrence probability
by mudflow: 0,896
Persons affected: 1
Persons affected: 3
Pylons affected: 4
by
mudflow:
0,137
Persons
affected: 15
by
landslide:
0,536
by
mudflow:
1,652
Persons
affected: 1
by
mudflow:
0,204
Persons
affected: 3
Pylons
affected: 9
by
mudflow:
0,500
Persons
affected:
129
Tröpolach Area
affected
[km2]:
by
mudflow:
0,120
by
inundation:
0,248
Persons
affected:
180
RattendorfArea
affected
[km2]:
by
inundation:
0,271
Persons
affected:
50
by
landslide:
2,911
by
mudflow:
1,804
by
inundation:
0,2
Persons
affected: 1
(contd)
Page 173
by mudflow: 0,204
Heavy traffic or strategic roads National road B90 Length affected [m]:
by landslide: 3450
Persons affected: 3
by mudflow: 2422
by inundation: 3440
by mudflow: 1935
by inundation: 2628
Pylons affected:
by mudflow: 15
by inundation: 11
Tourist accommodation—buildings 3 3 1 1
Lifeline 1 2 1 0
Vulnerability [V]
Residential areas 4 1 1 3
Lifeline 1 2 1 0
Vulnerability [V]
Residential areas 4 1 1 4
Secondary railway 1 2 1 0
Tourist accommodation—buildings 3 3 1 2
Lifeline 1 2 1 0
Page 175
Vulnerability [V]
6.5.1 Viability
In order to prevent serious debris flow in the Oselitzenbach torrent
and in the receiving stream “Gail”, the following measures are
possible:
6.5.2 Facilities
An extensive project of construction work and research was started
in 1988 to prevent debrisgenerating slope failures (especially at the
toe zone of the Reppwand slide) and the destruction of the “Naßfeld
Road”. As a part of the construction measures, a 400 m long channel
was excavated in massive Hochwipfel formations (comp. 1.2) and a
landfill was made at the toe of the excavated material as well (about
170’000 m3).
7
The Ceppo Morelli rockslide
7.1 INTRODUCTION
The Ceppo Morelli rockslide extends on the left side of the middle
Anzasca valley in the Pennine Alps, northern Piedmont (Italy), a few
kilometres from the border with Switzerland (Fig. 7.1). The national
road N° 549 runs on the valley bottom for about 30km, from the
Piedimulera village (Toce River plain) to the small village of Pecetto
di Macugnaga (at the foot of the Monte Rosa Massif). The road is the
only link to Macugnaga, a popular tourist resort close to Monte Rosa.
The Anza River drains the valley.
Major mass falls may reach two small villages dating back to the XVI
century: Prequartera and Campioli located at the bottom of the
valley.
Figure 7.3. Rock fall triggered during the October 2000 heavy
meteorological event. The boulder (right red arrow) damaged the
national road (red line) and stopped near a house.
The landslide has been subdivided into three areas (see Fig. 7.5):
A. Detachment area
This fan-shaped area lies below the supposed landslide toe and
consists of debris featuring big rock blocks. On the basis of the
geomorphological features, it has been subdivided into two sectors:
Mesostructural analyses
These analyses pointed out that the available data (Fig. 7.6),
although collected only in the upper part of the slope, are significant
for the whole landslide body and that they can consequently be
assumed as statistically representative.
Regional schistosity
The schistosity regularly dips to the SSW, both within and outside
the landslide area. Nevertheless, in the SW part of the landslide, it is
characterised by a steeper inclination (50–60°); in this sector it has
consequently been labelled as “S2”. On the other hand, the more
widely diffused low to medium angle schistosity has been indicated
as S1. The schistosity also dips to the NE in some places within the
landslide, probably due to the presence of open folds or the local
tilting of isolated rock blocks.
One of the most prominent is the K1 (Fig. 7.8 and 7.10a) system that
is considered here as a comprehensive discontinuity system whose
strike ranges from NNE to NE, and which has been synthetically
represented in Figures 7.6 and 7.8 as a single average plane. K1 is
indeed characterised by different sub-systems:
• faults, well developed in the lower part of the slope, namely where
the slope shows a greater inclination (at about 1250 m, see Fig. 7.5
Sector B3). These last surfaces strike from NE to ENE and consist of
very steep planes, tens of meters long and dipping both SW and NE.
It is assumed here that these faults belong to a fault zone; this
interpretation is mainly based on the geometric array of F elements
(scattered and anastomosed discontinuities) at both the regional and
local scales.
The western boundary fault of the landslide (Fb) also strikes sub-
parallel to the K1 system. This fault is characterised by a cataclastic
belt tens of meters long; minor slip planes within this belt indicate S-
ward right-lateral movements.
Other metric NNE reverse faults (Fi) have been also locally observed
(Fig. 7.10a). The NNE striking steeply dipping Fb and Fi displace the
less inclined ZT shear zones.
There are evident similarities between the aerial photos and the
satellite lineament data. The former shows a distribution range from
N20 to N60 with a maximum at N40. These lineaments are mainly
concentrated in the NE. The greatest lineament length was found
along the N80–90 direction.
The satellite data display the same distribution for both frequency
and length diagrams with three average maxima at N290, N350–20
and N50. The lineaments can also be distinguished in terms of their
length to determine the relations between frequency and length.
Frequencies were divided into three length classes according to
statistical distribution.
In summary, in both the aerial photo and the satellite data sets, four
main lineament systems can be recognised: the NE-SW (Ln1CM)
and WNW-ESE (Ln2CM) striking systems display higher azimuthal
frequency and longer lineaments; the NW-SE (Ln3CM) and N-S
(Ln4CM) striking systems show, instead, lower azimuthal frequency
and lengths.
• 1 rain gauge.
the highest total displacement (262 mm), and Points 8, 14, 21 and
23, located on the scarps of the detaching niche, which seem to be
still;
Sector B2: suffers minor translations and it is less damaged than the
other sectors.
Sector B3: is more damaged than the others, since it is cut by the
closely spaced K1 and F systems that subdivide it into small rock
blocks. The greater inclination of the slope induces relatively strong
gravity instabilities, as shown by toppling of blocks.
Ft joints are closely spaced and intensively affect the rocks, mainly
within the landslide. These joints, probably originated by unloading
processes (release joints, Engelder, 1985) show inclinations ranging
from few to 30 degrees.
K1 joints and maybe faults which are more closely spaced than the
other systems in Sector B3. Their mutual array, at different scales,
suggests that they could be segments of a NE striking fault zone.
In the lower part of the landslide (Sector B3), both sliding and
toppling have been found. Sliding could occur toward SSW (Fig.
7.15), while toppling could be originated by the partial reactivation of
pre-existing vertical fractures (K1 and K2) and cracking of the rock
bridges along the Ft system.
Page 200
7.3.4.5 Scenarios
On the basis of the previous considerations, three different hazard
scenarios are thus possible, with decreasing occurrence probability
and increasing impact on land planning:
ancient times. The sites from which rock blocks are released are
located all along the lower boundary of the unstable Sector B3.
Rockfalls range from small cobbles to large boulders hundreds of
cubic metres in size.
The 2D analyses were carried out using the Finite Element Method
(code DRAC, Prat et al., 1993). The approach applied during the
numerical modelling consisted of two stages: a first stage during
which the calculation of the initial stress conditions of the slope was
carried out using topographical information (the actual morphology,
where possible) and incorporating the effect of the glacier weight; a
second stage in which the glacier was removed and the stress field
computed. The objectives of the simulations were twofold. First, the
general stress field of the slope was calculated in order to obtain
information on the failure mechanism. The interaction of the upper
part (B1) that consists of a rather well-structured rock mass and the
lower part (B3) that is characterised by a disjointed rock body, was
an additional aim. Second, the influence of different material
properties and mechanical models for the cracks (interfaces in the
discontinuous model) were analysed.
Gneiss 28 50 0.3
Figure 7.19. Opening of the interfaces (red lines) at the final stage of
the Ceppo Morelli model. Width of the red lines is proportional to the
opening displacement of interfaces.
Figure 7.20. The unstable volume considered in the 3D numerical
analyses (Sector B3).
Page 204
The following step was the assessment of the value of the creep
coefficient C. Activating the Map3D creep option, a simulation was
made with an assumed value of C and the results were compared
with the measured displacements. If these values were not close to
each other, more simulations were carried out in order to find a
proper C. The final set of mechanical parameters used for the
numerical analyses is shown in Table 7.3.
Equivalent continuum
Failure surface
The obtained results are very far from giving the probability of failure
needed for the risk analysis, which is to be obtained through a
probabilistic procedure. For this purpose however the variation of
pore pressure in time and the decay of rock strength parameters in
time should be known. These data are seldom known and are not
known in the case of the Ceppo Morelli landslide.
Scenario 1
The aim of the analysis is to reproduce the rock falls occurred in the
past and to forecast new rock falls affecting a larger area.
To get realistic results the assumed range of velocity for the released
blocks changes if they originate from Line 1 (Fig. 7.21) or from the
Line 2 (Fig. 7.22) of the considered unstable area as indicated in
Table 7.4.
Page 205
1 1.0 1.8
2 0.6 1.0
In addition to rock fall trajectories (Fig. 7.21 and 7.22), the analyses
allow the distribution of the kinetic energy of the blocks (Fig. 7.23
and 7.24): the slope can then be subdivided in homogeneous
sectors, each associated to a specific value of the medium energy.
The distribution obtained in terms of energy represents an important
tool to be taken into account for risk and vulnerability evaluations.
Due to the fact that ROTOMAP is a “lumped mass” code, the energy
values are formulated as E=v2/2 [J/kg]. Hence, in order to link the
obtained results to a possible real situation, an indicative mass
should be hypothesised for the boulders. The maximum mass of the
boulders observed at the bottom of the valley (around 2.7×106 kg,
corresponding to a volume of 1000 m3) seems a reasonable value.
Figure 7.21. Scenario 1: rock fall trajectories from Line 1 (in red:
blocks in flight).
Figure 7.22. Scenario 1: rock fall trajectories from Line 2 (in red:
blocks in flight).
Page 206
Scenario 2
Together with the run-out area (Fig. 7.26), the value of the avalanche
front velocity is defined along the section taken into account (Fig.
7.27); the front velocity distribution is then supposed to be the same
for all the points located at the same height.
As the energy distribution for the rock fall case, the homogeneous
areas obtained in terms of velocity distribution are an important
aspect, that should contribute to the risk and vulnerability evaluation.
Scenario 3
Conclusions
distribution per unit mass, and the avalanche analyses give both the
definition of the run-out areas and the front velocity distribution.
Oct. 1971
Apr. 1977
Oct. 2000
Jun. 2002
843
1816
2000
Number of
buildings
affected: 5
Persons
affected: 12
Area
affected:
6300 m2
Prequartera Number of
buildings
affected: 18
Persons
affected: 42
Annual
traffic:
100000
vehicles
Annual
traffic: 1000
vehicles
Lifelines Power line Length
Valle affected:
Anzasca 15 1025 m
kV
Pylons
affected: 4
Number of
buildings
affected: 3
Persons
affected: 3
Annual
traffic:
100000
vehicles
Annual
traffic: 1000
vehicles
15 kV Pylons
affected: 1
Number of
buildings
affected: 3
Persons
affected: 3
Area
affected:
6300 m2
PrequarteraNumber of
buildings
affected: 18
Persons
affected: 42
Annual
traffic:
100000
vehicles
Annual
traffic: 1000
vehicles
15 kV Pylons
affected: 3
In the quantitative risk analysis for the Ceppo Morelli rock slide,
detailed economic evaluations were not considered for elements at
risk and arbitrary relative value indices have been temporarily
applied for each category (Table 7.9).
Secondary roads 1 2 1 2
Secondary lifelines 1 2 1 0
Figure 7.45. Flood effect in alpine area due to landslide dam (Randa,
CH, 9th May 1991). Effects also on buildings 1 km far from landslide
site. (source: www.crealp.ch).
Figure 7.46. Areas that might be affected by a flash flood due to
rockslide dam breaking (blue colored). The zones in orange squares
are detailed in Figure 7.47.
Page 223
7.5.1 Viability
After the many rockfalls of October 2000, the risk conditions of the
National Road were considered too high to be accepted. Since the
road is the only possible connection for Macugnaga, a well-known
ski resort, a provisional by-pass road was constructed along the right
banks of the Anza River, in order to avoid the most dangerous zone.
The road is about 1 km long and construction took about 30 days;
the total cost was on the order of 1.25 M€. The river is crossed twice
by means of Bayley-type bridges provided by army engineers.
Before construction, the Regional Geologists made a running
evaluation of the risk condition of the provisional road. Results
showed that the road is protected from rockfalls but may be affected
by major mass falls (e.g. a rock avalanche). A rock avalanche would
obviously destroy the provisional road.
Opening of the road was thus made conditional on the results of the
monitoring system. If some defined threshold values (both for rainfall
and displacements) are exceeded, the road is blocked by local
authorities. Up until now this has happened only once, during the
heavy rains of November 2002.
The wall is about 120 m and 5 m high; average base width is about 8
m. The wall was designed in order to withstand impact energies up
to 8 MJ; the total cost was about 0.25 M€.
As for Campioli, the situation is much more critical, for it is just on the
edge of the zone of possible invasion in the case of failure of the
lower part of the landslide body. The evacuated part of the small
village consists of seven houses. One of them, built in the seventies,
may also be affected in the case of common rockfalls; in fact, a 70
m3 boulder stopped a few meters from the building in October 2000.
The house will be demolished and relocated elsewhere.
The remaining six buildings pose more problems, in that they form a
historic nucleus built around 1650, a marvellous and perfectly
preserved sample of typical Walser architecture, protected by the
Sovraintendendenza ai Beni Architettonici (National agency for the
protection of the Architectural heritage).
8
The Sedrun landslide
F.Noverraz
Consulting geologist, Lausanne, Switzerland
8.1 INTRODUCTION
The unstable area, covering 1.5 km2, can be quite well determined
on the basis of its morphology, the presence of major cracks more or
less opened by the movements in the upper part and due to
numerous monitoring data (see § 8.3.3). As no boreholes have been
carried out, its volume can only be roughly assessed at some 100
million m3. The important movements recorded cause the formation
of typical uphill facing scarps in the upper part of Cuolm da Vi (see
Fig. 1.3).
This landslide site has been studied particularly between 1993 and
1997 as a part of a large research project funded by the Swiss
National Science Foundation, within the PNR 31 national project on
climate change and natural disasters (VERSINCLIM project—
Noverraz et al., 1998; Bonnard et al, 1996).
8.1.2 Historical background
The temporal origin of the movements observed in the Cuolm da Vi
unstable zone is not easy to determine. Due to the late installation of
the Swiss triangulation network in this area in 1942, and to the lack
of a renewed determination of the position of the landmarks until
1973, past information concerning the evolution of movements is
quite limited. No damage or special catastrophic event has taken
place in Cuolm da Vi itself, as no building or infrastructure exists in
that zone. The population was aware that some movements were
occurring in this zone but did not pay much attention and was not
overly concerned. The occasional snow avalanches originating from
Val Strem sparked much greater concern.
The extreme rainfall values for a return period of 100 years are 70
mm for a rain duration of 1 hour and 175 mm for a rain duration of 24
hours. The maximum annual precipitation amounts recorded at the
Sedrun station were 1’844 mm in 1935, 1’809 mm in 1999 and 1’836
mm in 2002, this last value being mostly influenced by the extremely
high precipitation in November (544 mm), but which fell mainly as
snow in the area of Cuolm da Vi. Since the beginning of the
measurements, no significant long-term tendency of variation has
been observed, except a long dry period between 1959 and 1964
and a wet period between 1999 and 2002 (133% of the average long
term rainfall over a period of four years). Therefore, it can be inferred
that rainfall is not directly the cause of the increasing movements
over the last 20 years.
Most of the extreme rainfall events in the area are produced by the
same type of meteorological process: warm and wet air masses
coming from the Mediterranean encountering the Alps are forced to
reach high altitudes, leading to sudden condensation and discharge
of their water content.
In the vicinity of the Cuolm da Vi area, there are two main streams
determining the limiting groundwater conditions. The first one is the
Val Strem stream, which runs along the western side of the unstable
zone of Cuolm da Vi, but at a much lower level, producing a deep
drainage effect for the landslide mass. To the east, the Drun Tobel
creek, which is less important in terms of flow, but presents a very
steep slope of 60% over 1’000 m, limits the unstable zone (Fig. 8.2).
The Drun Tobel creek displays a temporary flow, being especially
important in the snow melting season and
Figure 8.2. Drun Tobel creek seen from the village of Sedrun. The
unstable area extends up to the top of the mountain at the back.
Page 230
On the right side of the Drun Tobel, as well as along the east slope of
Val Strem, a series of temporary water outlets can be inferred by the
presence of several gullies located between elevations 1’900 m and
2’000 m. In the Cuolm da Vi area itself, absolutely no drainage
pattern is visible, which means that the net precipitation (rainfall and
snow less evapotranspiration) infiltrates completely in the moving
mass, in particular during snowmelt. However, no groundwater level
data are available.
The toe of the slope is incised by a deep erosion creek (Drun Tobel)
resulting from a recent torrential erosion. It displays kakiritic and
mylonitic rocks and favours the destabilization of the slope (Fig. 8.2).
The villages of Sedmn and Camischolas are located on a vast
alluvial fan, formed both by the Val Strem stream and the Drun Tobel
creek. In most of the area, the forest cover is absent.
2’458 m a.s.l.) and the bottom of the Drun Tobel outlet, with its floor
at 1’550 m a.s.l., leads to an average slope of 60% over a distance
of 1’500 m. The exposure of the unstable zone is oriented towards
the south (Fig. 8.3).
8.3.1.2 Rockslides
To the west, along Val Strem, there are some instances of rock
sliding, most of all at the top of the slope. Unfortunately there are no
monitoring points in this area, but secondary scarps clearly appear
and the rock mass is deeply cut by numerous faults and cracks.
Limited phenomena have already occurred as the toe of the slope in
Val Strem is covered with scree (Fig. 8.3).
To the east side of Cuolm da Vi, along the northern slopes of Drun
Tobel, this phenomenon is even larger, in terms of spatial extension
and activity, with several rockslides and rockfalls of limited depth
being reactivated by seepage and rainfall, as this zone shows the
presence of formations of mylonites and kakirites. A horizontal limit
below which an activity is not perceptible is clearly visible in Drun
Tobel (Fig. 8.2).
Page 233
8.3.1.4 Rockfalls
Although no specific evidence of large rockfall is visible in the
Sedrun Landslide area, except for some screes of larger blocks at
the toe of the left bank of Val Strem and for minor events visible in
Figure 8.5, in this zone, it is clear that main rockfall events are in a
preparatory phase to the west of Cuolm Parlets Dadens at an
elevation of approximately 2’400 m a.s.l., in the site called Gion
Giachen. The very high drop down to the bottom of Val Strem
(approx. 600 m), the very steep slope (some 80%) and the intensely
cracked rock mass have to be interpreted as major features for a
potential rockfall. However, it is difficult to forecast if these rockfalls
could occur as a single main event or as a succession of smaller
events.
The rock in the northern part of the site displays a rough schistosity
and a high strength (50–100 MPa). In the southern part, the rock
displays a more developed schistosity and high tectonisation. This
tectonisation leads to the formation of deep series of mylonites and
kakirites (superficial tectonical crushing inducing the formation of a
pulverulent rock without recrystallisation).
The S1 and S2 faults are present in the same zones and are
apparently not related to a local evolution of the same system.
Finally, some NS isolated faults can be observed.
These various fault families defining the local structural setting,
detected by field survey in 1994 and on the basis of aerial
photographs (Fig. 8.5), can be briefly described as follows:
Family S1
Family S2
Family S3
Family S4
The first geodetic measurements took place in 1983, when the ETHZ
began to reconfirm the coordinates of the 4th order triangulation net
in the area. The purpose of this campaign was to analyse the
movements induced by probable slide processes and not those
related to active faults already described in the zone for a research
project for the National Science Foundation, by Eckardt, Funk and
Labhart (1980 and 1983). This was the reason why this study
focused on the zone of Cuolm da Vi, where a series of some 20
complementary points were installed. During these measurements,
the movements of the triangulation points of the 4th order, namely
Nos. 514, 516 and 615 were verified. Additional measurements were
carried out again in 1988, then partially in 1989 and 1990 (only
points 6 and 26 for this last year, by the survey office of the canton of
Graubünden—see Fig. 8.6) [Bovay-Huguenin surveyors, 1994].
Figure 8.8. Bend of Drun Toble Creek where the debris flow
considered in Scenario 1 will flow directly towards the village of
Sedrun (the first visible houses are only barns).
alluvial fan of Val Strem, reaching the course of this river to the west
of Sedrun; it will also extend along the Drun Tobel valley where it
crosses the village of Sedrun down to the Rhine river. The narrow
zone of impact in this last stretch may be widened by erosion
phenomena which would destabilise the creek banks (Fig. 8.9).
The affected zone has been determined on the basis of the local
topographical conditions of the slope in the lower part of Val Strem,
referring to other similar situations which have occurred in
Switzerland, in particular in August 1987 in Poschiavo and in
September 1993 in Brig. For such a case, and considering the large
uncertainties affecting the volume as well as the debris flow
mechanism, mathematical modelling would not provide more reliable
results to determine the expected zone of impact. Indeed these two
first scenarios will be used in a comparison of quantitative risks with
respect to phenomena of high and medium probabilities of
occurrence (100 to 300 year return period for Scenario 2).
seepage develops when the lake level reaches the lowest point of
the deposited rockfall mass, it would probably not prevent an
overflow which would cause a fast erosion and thus create a
massive flood downstream. Considering the velocity of such flow
induced by the steep inclination of the stream bed in Val Strem
(approx. 10%), which would exceed 10 m/s, the flood would dash
down the river bed and pass over its banks, devastating a part of
Sedrun (Community of Camischolas) down to the Rhine river. This
catastrophic scenario, assessed to have a return period of some 300
to 1’000 years, may have serious consequences on the Rhine valley
downstream, in particular for a hydro-electric plant as well as one of
the building sites of the Gotthard base tunnel and for its respective
construction facilities. The present evolution of the affected zone
near Cuolm Parlets Dadens may even lead to the assumption that
the probability of occurrence might increase in the future.
For the risk computation, the following return periods have thus been
assessed in a conservative way, especially for the last scenario.
Restaurant
Garage
Barn, farm
Hydro-electric plant
Bridge (value/m)
Forest (value/m2)
River (value/m)
Vehicles, movables
Housing (value/day)
Unemployment (value/day)
Tourism (value/day)
NumberValue VulnerabilityCost
People 25 ?
Vehicles,
movables
Housing
(value/day)
Unemployment
(value/day)
NumberValue VulnerabilityCost
People 150 ?
Vehicles,
movables
Housing
(value/day)
Unemployment
(value/day)
Figure 8.14. Bridge for the cantonal road crossing the Val Strem
creek that would be damaged by the debris flow considered in
Scenario 2.
NumberValue VulnerabilityCost
People 120 ?
1.344 mio €). The damage to buildings will represent the higher part
of the expected costs (approx. 3.43 mio CHF or 2.2 mio €). As far as
socio-economic activities are concerned, a limited amount has been
considered for closed roads, because the main access road to
Sedrun along the Rhine valley to the east should not be affected and
as it will be possible to re-establish a by-pass road through
Camischolas for the traffic coming from the Oberalp Pass. No
amount has been computed for the loss of tourist activity as the
major facilities will not be concerned; but the reputation of the resort
might be affected.
Scenario considered 1 2 3
As the three scenarios are independent of one another and will not
necessarily occur within a short period during which damage will not
yet have been repaired, the total risk implied corresponds to the sum
of the three quoted figures, namely a monetary risk of some 205’000
€ (yearly amount to spare in order to cover future damage) and a risk
to life of 4.1 victims per year. Preventive measures are thus definitely
required.
• River protection works along the lower part of Val Strem, in the
area of Camischolas-Zarcuns, down to the hydro-electric plant near
the Rhine.
8.6 CONCLUSIONS
9
The Séchilienne landslide
Philippe Marchesini
Direction Départementale de l’Équipement de l’Isère
9.1 INTRODUCTION
The landslide can be divided into three principal areas (Fig. 9.3): the
frontal mass, the most active (several decimetres per year; source of
many rockfalls) and disrupted part, the volume of which can be
estimated at approximately 3 hm3; an intermediate zone with
medium activity; the upper and north-western part, corresponding to
an elliptic (probably) post-glacial sagging, with low velocities.
Page 254
The south-facing slope of Les Ruines has an angle of 45° in its lower
part and about 20° in its upper part and near the crest. The main
feature of the slope is the sagging of Mont Sec, bounded
Page 257
These faults divide the moving zone into four parts (Fig. 9.7). The
geodetic measurements carried out over the past 20 years are fully
coherent with these structural features.
– a large slowly moving zone (the ancient sagging): 0.02 to 0.04 m/y.
Figure 9.14 shows another result for the spreading of the debris
using a digital elevation model and a mechanical process with
N=1000 pieces of rock sliding one after the other, given a law of
energy dissipation.
Figure 9.14. Debris cone resulting from a 1 hm3 volume i.e., 1.15
hm3 of debris simulation (courtesy of J.-F.Serratrice). It can be seen
that the debris cone does not reach the earth protection dam in this
case.
– Group 1 (short term i.e. <10 years): toppling and falls originating
from the frontal zone; many rockfalls; failure of the whole frontal
zone (about 3 hm3), probably in several steps. There are two main
possibilities: continuation of the current behaviour (rockfalls of a few
cubic meters to hundreds or thousands of cubic meters) or a
significant rockslide involving the whole fastmoving zone: the volume
could be several million m3 (2 to 3 in the more recent evaluation).
From the knowledge of the slope and the results of these studies,
various scenarios were defined. These scenarios make it possible to
plan the actions of the authorities, either preventive ones (i.e., land
use control) or curative ones (emergency plans). For each scenario,
the authorities take into account a reasonably elevated evaluation of
the risk, to guarantee a high level of safety, while not spending
financial resources in a useless or premature way.
– In the case of high damming and rising water behind the natural
dam formed by the debris: flooding of the major part of the village of
Séchilienne (670 inhabitants),
– In the case of overtopping and rapid erosion of the dam: flooding
downstream in the Romanche valley affecting the town of Vizille and
its surroundings (10’000 inhabitants) and chemical industries near
Grenoble, etc.
Romanche Cone scree reaches the river, likely locally diverting the flow
River towards the earth dam, inducing a beginning of diversion by
the new bed. Pollution (suspended matter) of the Romanche
River.
Inhabited Dust, etc. Possible evacuation during the crisis for the closest
zones areas. No significant consequence.
upstream
Inhabited Dust, etc. Evacuation during the crisis (several weeks) for the
zones closest sectors. Possibility of flood for some of the houses
upstream located at the lowest level.
Inhabited Dust, etc. Evacuation during the crisis (several weeks) for the
zones closest areas: 5’000 to 15’000 inhabitants affected. Damage to
downstream numerous buildings (several hundred). Damage to roads,
(2–5 km) bridges (destruction) within a radius of 20 km (local or
national level).
Element ScenarioScenarioScenarioScenario
1 2 3 4
National road RN 91 and other roads, 0,05 0,2 1 5
railway
One should also add that the social and economic impact already
exists: the development of the Séchilienne village has been stopped
and it is also probable that the chemical industries of Claix and Jarrie
downstream in the Romanche valley are no longer investing at their
sites. Moreover, large expenses have been devoted over the past
twenty years to investigation and monitoring of the landslide.
Figure 9.16. Ile-Falcon viewed from the frontal zone of Les Ruines.
The first emergency plan was worked out in 1993. A second version
was prepared in 1999, corresponding to Scenario 3. Because of the
recent analyses of the slope stability, and of the need for having an
action plan adapted to the short term risks, a new version has
recently been prepared. It corresponds to Scenarios 1 or 2.
9.6 CONCLUSIONS
The history of Les Ruines de Séchilienne has not yet come to its
end. The scientific experts have sometimes been mistaken, but the
observations and data collected from the site allow continual
understanding of the phenomenon and our ability to define the most
probable scenarios, and therefore to prepare the mitigation decisions
which are the responsibility of the political authorities.
Page 270
10
Specific situations in other contexts
Oldrich Hungr
University of British Columbia, Vancouver, Canada
10.1.1 In troduction
Some 1.2 million km2 of Canadian territory, an area equal to more
than half of Western Europe, contains mountainous terrain with
sustained local relief in excess of 1000 m. Glacial erosion during the
Holocene and Early Pleistocene left deep valleys and oversteepened
mountain slopes. Large rockslides are relatively common in this
area. Surprisingly, damage due to large rockslides is not a common
occurrence. Evans (2003)—see references at the end of the book—
reports that 238 lives have been lost to rockslides and rockfall in all
of Canadian history, more than half of which occurred in only two
events (the 1903 Frank Slide in Alberta and the 1915 slide at Jane
camp, near Vancouver). The probable reason for this rather limited
mortality is the remoteness and low population density of the
mountain lands. The Pandemonium Creek rock avalanche (Evans et
al., 1989) is an example that illustrates the point. The rock avalanche
probably occurred in 1958 in a mountainous area approximately 300
km north of Vancouver. Several million of m3 of rock failed from a
cirque headwall, passed over a glacier and traveled 11 km along a
mountain valley, devastating some 8 km2 of area (Fig. 10.1.1). Not
only were there no fatalities, but the landslide was not even noticed
by anyone! The extensive damage to the natural landscape was first
identified on air photos only in the 1970’s.
10.1.2 Examples
10.1.2.1 Large rockslides on hydroelectric dam reservoirs
B.C.Hydro, a crown corporation of the British Columbia Provincial
Government, owns 60 hydro-electric dams. The reservoirs behind
these dams form thousands of km of bank slopes in mountainous
Page 272
terrain. Ever since the Vaiont disaster of 1963, there has been
awareness of the danger posed to dams by large landslides
displacing reservoir water. In the case of the dams on the Columbia
River, for example, a landslide-generated displacement wave could
potentially destroy a large earth or rockfill dam, generating a
cascade of dam failures along a river corridor 1000 km long,
reaching into the territory of the United States. Consequently,
B.C.Hydro operates an intensive program of investigation and
monitoring of reservoir slopes. Over the last 30 years, more than
2000 km of slopes have been subjected to systematic landslide
hazard evaluations by B.C.Hydro and its consultants (Imrie and
Moore, 1997).
Some very large specific rockslides have been studied in depth and
subjected to remedial construction. The first of these, the Downie
Slide, situated on the Columbia River 70 km upstream of Revelstoke
Dam, is one of the largest recognized landslides in Canada. It
involves a mass of 1.4×109 m3 of metamorphic rock, situated above
a fault dipping at about 20° towards the reservoir and more than 200
m beneath the slope surface (Imrie et al., 1992). The rockslide
underwent a displacement of about 100 m in prehistoric time, but
was in a dormant state until disturbed by water impoundment. During
filling of the reservoir, a part of the slide mass began to move at
rates of up to 100 mm/year. Although such slow movement posed no
direct danger to the dam, there was concern about a Vaiont-like
sudden displacement. B.C.Hydro carried out a detailed investigation
with drill holes, testing and detailed suface mapping and then
designed a dewatering system with a tunnel and fans of horizontal
and inclined drainage holes. The total cost of the drainage system
was 30 million $CDN. The drainage decreased the slope movement
rates to about 10 mm/year. Design changes were implemented at
the dam, in order to accommodate a displacement wave, predicted
by physical modeling. Full time remote monitoring of displacement
rates and piezometric levels is continuing.
Figure 10.1.3. A zone of cracks and scarps on the south ridge of the
Cheekye valley. The valley slope is visible in the upper left corner of
the nearly-vertical photograph.
The study comprised two parts. The first was a detailed surface and
geophysical investigation of the two ridges above the Cheekye
Valley. Stability analyses based on the data collected were used to
predict the likely maximum volume of detachment corresponding to
the extent of the cracked zones. The second part was a detailed
stratigraphic study of the Cheekye Fan deposits to establish the
magnitude and frequency of past debris flows.
The stability investigation in the headwaters showed that the
cracking was caused by an incipient translational failure of the
volcanic breccia deposits over a zone of heavily altered basement
rock, just below the gently sloping contact between the two units. It
was concluded that a sufficiently large volume could detach during a
future failure retrogression, with the possibility of sending several
million m3 of rock into the valley. Damming of the Cheekye stream
and subsequent erosion of alluvial and colluvial deposits from the
valley floor and sides could generate a large debris flow in the river.
However, it was not possible to predict the probability of occurrence
of such an event, especially because the tension features on the
ridge were found to be currently inactive.
was likely to extend into the future. Using the curve, a prediction of
the probability of occurrence of three magnitude classes of debris
flows was made. With the guidance of the identified deposits, a run-
out analysis of debris flows between the fan apex and the hazard
area was conducted. Based on this, a hazard intensity map was
constructed as shown in Figure 10.1.4. Each of the hazard zones on
the map is associated with estimated hazard intensity parameters
(flow velocity, thickness of deposits) for each of the three classes of
debris flow magnitude. This map in itself can be used by the
planning authorities to determine whether a certain level of
development is suitable or not.
The right flank of the slide contains a large mass of rock which
moved several meters during the slide, opening large tension cracks,
but did not accelerate and remained on the slope (Mathews and
McTaggart, 1978). The British Columbia Ministry of Highways carried
out a detailed investigation of this remnant mass. Run-out studies
were conducted using a fragmental rockfall model. The highway was
relocated to the opposite side of the valley based on the results.
Figure 10.1.5. The source area of the 1964 Hope Slide. The
disturbed rock mass is visible in the upper left corner of the image.
Figure 10.1.6. An overall picture of the 1903 Frank Slide in Southern
Alberta. A disturbed rock mass exists near the peak of the Turtle
Mountain, in the upper left part of the picture (Photo courtesy Dr.
D.M.Cruden, University of Alberta).
Page 277
10.1.3 Conclusions
Over the last 30 years, a body of experience with the management
of large rockslide hazards has been gradually developing in Western
Canada. It is, however, a difficult field and no general
Page 278
P.Basabe
United Nations International Strategy for Disaster Reduction,
Geneva, Switzerland
10.2.1 Introduction
The location of Ecuador, a small country in South America located
along a convergence zone of tectonic plates, has played a principal
role in producing not only its beautiful geography, but also a high
susceptibility to natural disasters such as earthquakes, tsunami and
volcanic eruptions. Moreover, the influence of ocean currents
exposes this country to intense hydrometeorological phenomena
which cause landslides, debris flows and floods. These natural
phenomena turn into hazards as population and activities increase
without appropriate planning actions, causing a growing occupation
of land in inappropriate zones, so that the vulnerability and risk
factors continue to increase.
The economic impact of “natural disasters” in Ecuador is very
important, as shown in Figure 10.2.1 for the period of 1980 to 2001.
The gross domestic product (GDP) increase in percentage per year
normally varies between 2.0 and 4.3; nevertheless, most of the
GDP’s decrease coincides with the occurrence of the main disasters
having natural geneses in the same period. In 1982–83, the year of
the El Niño phenomenon, the GDP decreased and reached a
variation of 1.2% and −2.8%, respectively. In 1997–98 the impact of
the “El Niño” phenomenon which caused floods, landslides and
related damage was responsible for direct losses of 2’882 million
dollars (CAF-CEPAL, 2000), i.e. 15% of the GNP, which is the
highest value recorded in a South American country (Basabe, 1998).
The GDP variation decreased to less than 2.0% in 1998 and even
−7.0% in 1999.
But the major aspect of the project included the development and
application of non-structural disaster risk reduction measures. A
large inter-institutional cooperation was developed with national,
regional and local partners to create a multidisciplinary approach.
Above all, primary importance was given to the professional training,
in order to ensure the continuation of monitoring, analysis and
implementation of disaster reduction measures at the end of the
project. The measures also included Civil Defence strengthening, as
well as the training of local and regional professionals operating in
the application of hazard assessment and mapping within land use
planning actions, as well as the promotion of regulations and laws as
part of sustainable development and related programmes.
The main triggering causes are the intense and long-lasting rainfall,
as well as the anthropic actions such as deforestation, inappropriate
cuts and fills for roads, quarry exploitation and undue use of soil.
Qualitative Quantitative
A High 0.66<x<1.0
M Medium 0.33<x<0.66
B Low 0.00<x<0.33
The first two years of the project were focussed in mapping and
hazard assessment of landslides by covering 384 km2 in 1994 and
512 km2 in 1995 (Fig. 10.2.3). The communities located in main
hazard prone areas were informed through a joint effort with Civil
Defence and the Ministry of Education, regional branch. Specific
booklets and comics (40’000 copies each) were created to inform the
population.
The task of getting involved and not creating panic in the population
is not easy. The project took advantage of a specific programme
between Civil Defence (CD) and the Ministry of Education in which
bachelor students need to develop a thesis or a social work on Civil
Defence during one year. Every year some 800 hundred students
were informed on landslide assessment and trained on CD and risk
disaster reduction management. Every student had to apply her/his
knowledge for the community. The work consisted to take care of 3
or 4 families living in the landslide prone areas to inform them about
hazards and risk globally, as well as on landslide hazards, with
emphasis on disaster risk reduction and CD measures. At the end of
the social work the student delivered a report to the Ministry of
Education and CD office who took care of the follow-up.
Page 286
In the case of Paccha, the biggest landslide prone area (13 km2)
discovered (see § 10.2.3.3), the vulnerability assessment was based
on the compilation of information through the census of population,
housing, infrastructures, economic activities and social conditions
(Fig. 10.2.7). The census was developed with volunteers of Civil
Defence (CD) and the Ministry of Education, in conjunction with the
social work, to inform the population on hazards, vulnerability,
landslides and simple disaster risk reduction measures.
At the end of the project, the region of Paccha, which had been
designated within the expansion urban zone of the city of Cuenca
before the project, was redesignated as agricultural and recreational
forest area. The same designation was done for flood prone areas
along the streams; in this case the ordinance indicates longitudinal
parks.
• National level
The studies carried out, as well as the maps, the GIS and the
monitoring systems were handed over to the national institutions
participating in the project, namely the National Direction of Civil
Defence, the National Institute of Electrification and the National
Meteorological and Hydrological Institute (INAMHI). These
institutions have continued to use the results obtained and, in the
case of INAMHI, to be part of the maintenance, data treatment and
management of the hydrometeorological monitoring networks and
forecast.
• Regional level
The studies carried out by PRECUPA project allowed the
development of other projects in the region of Cuenca, in particular
the one, financed by the European Union, which was set forth by the
institution created for the rehabilitation of the zone affected by the La
Josefina Landslide and which used the results obtained by the
PRECUPA Project. Thus, the destroyed road was rebuilt,
Page 287
• Provincial level
The provincial council of Azuay, namely the provincial government,
received constant support by the PRECUPA project in order to
assess the situation of landslide-prone areas in rural surroundings,
even outside the zone of study. Such advice was directly useful to
the population for opening new roads, building protection works or
providing support to local communities. The results and maps of the
project were also handed to the Provincial Council of Cañar, north of
Azuay Province.
• University level
Several activities of the project were carried out with the participation
of the regional universities, which also contributed to the training of
students and assistants. For instance, the maintenance and
operation of the seismic monitoring network, data base treatment,
hazard and vulnerability assessment and further investigations are
managed by the Cuenca University.
• Local level
10.2.5 Conclusions
The participation of national and local institutions in the project, the
knowledge development, the training of operational professionals
and future decision makers, as well as the community and civil
defence work towards a large dissemination of the results obtained,
have guaranteed the continuity of activities and the implementation
of specific risk reduction measures (Basabe et al., 1996). Thus much
of the data and results could be incorporated in municipal ordinances
and landuse planning. The conscience of the importance of natural
hazard and risk management has clearly increased, allowing a better
use of the land. This is expressed by a widespread popular slogan:
“The knowledge and respect of Nature is fundamental to protect
ourselves”.
Page 289
11
Suggestions, guidelines and perspectives of
development
are not ready to accept living with such risk (unless they have good
reasons to do so, or do not believe in the threat).
In the case of the Sedrun landslide, for which the hazard level is not
considered as severe by the inhabitants and, consequently, the risks
do not affect the market prices in this tourist resort, the value
considered for the buildings has been assessed on the basis of
average construction prices in the area, which is quite appropriate as
a large majority of the houses are recent and fairly standard (3
storied chalets or 5 storied buildings with 8 apartments). But, of
course, the risk is assessed without considering the variation of
prices over time. As far as the cost related to traffic interruption, it
can be computed using the tables established for that purpose by
the railway companies and road administrations.
The possible attitudes and actions when facing such a situation are
very different when a master plan or a local plan is considered. At
the national or regional level, landslide risk situations, once identified
and recognised, are considered as local constraints which will be
duly taken into account in the planning of new infrastructures
(motorways, railway lines), as well as in specific development
programs, for instance for flood protection. But such constraints do
not substantially modify the development of projects, as they can be
solved technically and, socially speaking, do not represent a major
hindrance.
As far as the authorities are concerned, the situation varies from one
country to another, depending on the respective competencies of the
local and regional authorities for the execution of safety measures
(see § 11.4). It is clear that the municipal authorities, who have to be
re-elected every four to six years, are not induced to take drastic
safety measures limiting the use of the land, as this will displease the
threatened owners and induce a feeling of injustice between the
affected and non-affected citizens. Otherwise, they will take partial
measures and then declare that the landslide risks are “under
control”, whereas the level of safety is only slightly or locally
increased.
The fact that the adaptation of local plans in order to include new
hazard zones or modify their requirements must be approved by
municipal councils, makes it difficult to proceed to swift changes
when the nature of the phenomenon evolves. Of course, the mayor
generally has full power for the enacting of safety or evacuation
measures when a critical situation occurs, but such actions are of
limited duration. Therefore, it is not easy to reach an approved
consensus or the rapid application of restrictive measures in order to
reduce the landslide risks. But when the local authorities, and in
particular the mayor, have not taken the required safety measures to
avoid a disaster, they are legally considered to be responsible for the
damage, in particular if human victims are reported.
11.4.1 Introduction
In the face of a natural hazard, risk management can be divided in
several stages:
The public authorities have to choose which kind of policy for risk
prevention and mitigation they will apply. One can define three
degrees of increasing involvement by the authorities:
From a general point of view, some final remarks can be drawn from
the comparison of different risk management policies in Europe:
11.4.2 Italy
Before the 1990’s, the policy in Italy was characterised by two
principal behaviours:
Figure 11.1. Large debris flow which affected the town of Sarno.
11.4.3 France
Since the 1980’s, several laws related to natural hazard prevention
have been passed. Citizens are more and more concerned with their
own safety and the idea of “fate” is no longer accepted. The
prevention of natural or technological hazards has become a
politically important topic.
– the “Plans d’Exposition aux Risques” which regulate and plan the
land use at the municipal level; they were changed into “Plans de
Prevention des Risques” (PPR) in 1995; more than 10’000
communities are involved (5000 PPR have been produced up to
now). The PPR are the main tool for the Government to encourage
better land use in relation to natural hazard. In the most dangerous
zones, building is forbidden; in zones of medium hazard, preventive
measures have to be taken.
11.4.4 Switzerland
The 1979 Federal Law on Land Planning (LAT) requires that the
areas seriously endangered by natural hazards be designated
specifically in cantonal master plans. Their management, as well as
the type of qualification and representation, then depends on the
cantons. But a federal recommendation has been produced in 1997
specifying how to take landslides into account in land planning
(OFAT, OFEE, OFEFP, 1997).
(macroscopic) properties of the rock mass that could account for the
observed deformation and morphology. This clearly relies upon a
detailed structural field survey and should involve geologists in
collaboration with mechanical and numerical modelling specialists,
within multidisciplinary teams.
The quantitative data related to the exposed objects may exist but
are difficult to find and to collect. In the near future, data bases and
geographical information systems will probably facilitate the
obtaining of information related to the spatial distribution, the
functionality and the cost of the main objects. As an example, the
Swisstopo database includes a GIS layer containing all of the
buildings in Switzerland. The evaluation of indirect costs, because of
the closure of a road for instance, has yet to be formalised;
economists should be involved in such analyses in conjunction with
engineers.
To take into account the large landslide risk means to change the
proper assessment scale, i.e. to extend it from a local perspective
due to the directly exposed elements at risk to a more extensive
regional one, often not easy to identify, but usually of relevant
impact. These risks are related to the indirect effects, e.g. chain
effects, implying further impacts on social and economic assets
etc…To assess the global risk thus means to take into account more
vulnerability aspects, which often interact with one another, often
implying different zones of influence both in terms of dimensions and
spatial identity. These interconnected aspects also include the
quality and availability of regional socio-economic, environmental
and cultural data, as well as the capacity for understanding and
managing this topic correctly.
While hazard analysis and physical aspects of vulnerability
evaluation have been substantially facilitated and improved due to
the use of GIS techniques, the inclusion of social, economic and
environmental variables into GIS conceptual models remains a major
methodological challenge. The need to assign a quantitative value to
the variables analysed in the spatial model used by GIS is not
always feasible for some social/economic dimensions of
vulnerability. For instance, how can the social impact of a destroyed
international railroad or the impact on cultural aspects be quantified?
The spatial scale of the examined instability problem and its related
effects, the role exerted by climatic and meteorological triggering
factors and the time-dependent rock mass behaviour influence the
temporal scale of the analysis. What is the likely future behaviour of
the slope (acceleration leading to a catastrophic failure, deceleration
then stabilisation, change of speed with seasonal fluctuations, etc.)?
Even though this question is raised systematically, at the present
time it is rarely solved. Risk analysis requires the assessment of
occurrence probability in time; however, the time element is of great
uncertainty in landslide hazard/risk analysis. Occurrence probability
which can be defined as the chance or probability that a landslide
hazard will occur can be expressed in relative (qualitative) or
probabilistic (quantitative) terms. A rigorous, quantitative, procedure
could include the application of formal mechanical-probabilistic
methods, taking into account the uncertainty in all the geometrical
and mechanical parameters, including the variability of the boundary
conditions (for instance pore pressure distribution in the slope, initial
state of stresses etc.).
The risk management quality depends not only on the quality of each
phase but also on the communication among the experts of the
consecutive phases. Practical difficulty exists since each phase is
performed by experts of different disciplines and the integration of
the works may suffer from communication problems.
From the experience gained in the IMIRILAND project and also from
the previous experience of the groups involved, it is possible to
improve the integration of knowledge in order to guarantee
straightforward and efficient communication among the disciplines.
REFERENCES
Basabe, P., Neumann, A., Almeida, E., Herrera, B., García, E. &
Ontaneda, R, 1998. Informe Final del proyecto PRECUPA de
cooperación entre el CSS, DNDC, INECEL, INAMHI, Munic. de
Cuenca/ETAPA, U.Cuenca y Consejo de Programación por “La
Josefina”, para la prevención de desastres naturales en la cuenca
del río Paute. Temas: Topografía/geodesia, Geología/geot-ecnia,
Hidrometeorología, Sismología, Limnología, Capacitación/difusión y
apoyo a Defensa Civil, Cuerpo Suizo de Socorro, 369 p. más
anexos, Cuenca-Ecuador.
Bogge, A., 1975. L’alluvione del 1728 in Val di Susa. Centro Studi
Piemontesi 4(2): 379–396.
Bonnard, Ch., Noverraz, F., Lateltin, O. & Raetzo, H., 1995. Large
Landslides and Possibilities of Sudden Reactivation. Proc. 44th
Geomechanics Colloquy, Salzburg. Felsbau No 6/95, pp. 401−407.
Carol, I., Prat, P.C. & López, C.M., 1997. A normal/shear cracking
model. Application to discrete crack analysis. Journal of Engineering
Mechanics (123–8): 765–773.
Carrara, A., 1992. Landslide hazard assessment. Proc, Ist Symp. Int.
Sensores Remotos y Sistema de Inform. Geogr. para el Estudio de
Riesgos Natur., March 10–12, 1992, Bogotá, 329–355.
Castelli, M., Forlati, F. & Scavia, C., 2001. Landslides: the questions
of the decision-maker and the answer of the modelisation, Proc. 1st
Int. Conf. Albert Caquot, Paris, 3–5 October 2001.
Page 309
Compagnoni, R. & Lombardo, B., 1974. The Alpine age of the Gran
Paradiso eclogites. Rend. Soc. It. Min. Petr., 30 (1): 223–237.
Corominas, J., Copons, R., Vilaplana, J.M., Altimir, J. & Amigò, J.,
2002. Integrated landslide susceptibility analysis and hazard
assessment in the principality of Andorra. Natural Hazard, Kluwer
Academics Publisher: 1–15.
Cruden, D.M. & Varnes, D.J., 1996. Landslide types and processes.
In: “Landslides: Investigation and Mitigation”. Transportation
Research Board. National Academy of Sciences, pp. 36–75.
Davis, G.H. & Reynolds, S.J., 1996. Structural geology of rocks and
regions. 776 p. New York: Wiley & Sons.
Dela Pierre, F, Lozar, F. & Polino, R., 1997. L’utilizzo della tettono-
stratigrafia per la rappresentazione cartografica delle successioni
meta-sedimentarie nelle aree di catena. Mem. Sci. Geol. (49): 195–
206.
Evans, S.G., Clague, J.J., Woodsworth, G.J. & Hungr, O., 1989. The
Pandemonium Creek Rock Avalanche, British Columbia. Canadian
Geotechnical Journal, 26:427–446.
Fell, R., 1994. Landslide risk assessment and acceptable risk. Can.
Geot. J., (31): 261–272.
Fell, R., Hungr, O., Leroueil, S. & Riemer, W., 2000. Keynote
Lecture015—Geotechnical engineering of the stability of natural
slopes, and cuts and fills in soil. Proc. GeoEng 2000 Conf,
Melbourne, Vol. l, pp. 21–120.
Finlay, P.J. & Fell, R., 1997. Landslides: risk perception and
acceptance. Can. Geot. J., (34): 169–188.
Forlati, F, Ramasco, M., Susella, G., Barla, G., Marino, M. & Mortara,
G., 1993. La deformazione gravitativi di Rosone. Un approccio
conoscitivo per la definizione di una metodologia di studio. Studi
Trentini di Scienze Naturali, Acta Geologica 68:71–108.
Fratianni, S. & Motta, L., 2002. Andamento climatico in Alta Val Susa
negli anni 1990–1999. Regione Piemonte, Università degli Studi di
Torino. Studi climatologici in Piemonte (4).
Geo & Soft International, 1999. ROTOMAP for Windows. Turin, Italy.
Geo & Soft International, 2003. ISOMAP & ROTOMAP for Windows
(3D Surface Modelling & Rockfall Analysis). User’s Guide.
Hungr, O., 1995. A model for the runout analysis of rapid flow slides,
debris flows, and avalanches. Canadian Geotechnical Journal,
32:610–623.
Imrie, A.S., Moore, D.P. & Enegren, E.G., 1992. Performance and
maintenance of the drainage system at Downie Slide. Proc. 6th. Int.
Symp. on Landslides, Christchurch, New Zealand, Vol. 1:751–757.
Rotterdam: Balkema.
Läufer, A.L., Frisch, W., Steinitz, G. & Loeschke, J., 1997. Exhumed
fault-bounded Alpine blocks along the Periadriatic lineament: the
Eder unit (Carnic Alps, Austria). In: Geol. Rundschau, Bd. 86, S.
612–616.
Leone, F., Asté, J.-P. & Leroi, E., 1996. Vulnerability assessment of
elements exposed to massmovement: working toward a better risk
reception. Proc. VIIth Int. Symp. on Landslides, Trondheim, Norway,
Vol. 1, p. 263–270. Rotterdam: Balkema.
Luino, F., Ramasco, M. & Susella, G., 1993. Atlante dei centri abitati
instabili piemontesi. Gruppo per la difesa delle catastrofi
idrogeologiche. Prog. Spec. Studio centri abitati instabili, pubb. n°
964, CNR-IRPI, Regione Piemonte, Settore Prevenzione Rischio
Geol. Meteorol. e Sismico. p.245.
Mathews, W.H. & McTaggart, K.C., 1978. The Hope rock slides,
British Columbia, Canada. In Rockslides and Avalanches, Voight, B.
Ed., 1:197–258. Amsterdam: Elsevier.
Mayoraz, F., Cornu, Th., Djukic, D. & Vulliet, L., 1997. Neural
networks: A tool for prediction of slope movements. Proc. 14th Int.
Conf. on Soil Mechanics and Foundation Engineering, Hamburg. Vol.
1, pp. 703–706.
McConnell, R.G. & Brock, R.W., 1904. The great landslide at Frank,
Alberta. Canadian Parliament Sessional Paper, No.25, Department
of the Interior, Government of Canada.
Moore, D.P., Imrie, A.S. & Baker, D.G., 1991. Rockslide risk
reduction using monitoring. Proc. Canadian Dam Safety Conference.
Vancouver: BiTech Publishers.
Nichol, S. & Hungr, O., 2002. Brittle and ductile toppling of large rock
slopes. Canadian Geotechnical Journal, 39:1–16.
Perello, P, Delle Piane, L., Piana, F., Stella, F. & Damiano, A., 2004.
Brittle post-metamorphic tectonics in the Gran Paradiso massif
(north-eastern Italian Alps). Geodinamica Acta, in press.
Prat, P.C., Gens, A., Carol, L, Ledesma, A. & Gili, J.A., 1993. DRAC:
A computer software for the analysis of rock mechanics problems. H.
Liu (Ed.), Application of computer methods in rock mechanics. Xian
(China): Shaanxi Science and Technology Press. 1361–1368.
Puma, F., Ramasco, M., Stoppa, T. & Susella, G., 1984. Carta dei
movimenti di massa nelle alte valli del Chisone e di Susa.—in
Compagnoni et al., Geotraversa della Zona Piemontese nella Valle
di Susa, escursione pre-congresso (10÷11/9/1984) 72° Congresso
Soc. Geol. It., Serv. Geol. Reg. Piemonte.
Puma, F., Ramasco, M., Stoppa, T. & Susella, G., 1989. Movimenti di
massa nelle alte valli di Susa e Chisone. Boll Soc. Geol It. (108):
391–399.
Zettler, A.H., Poisel, R., Roth, W. & Preh, A., 1999. Slope stability
analysis based on the shear reduction technique in 3D. Proc. of the
International Flac Symposium on Numerical Modeling in
Geomechanics, Minnesota, Minneapolis, pp. 11–16.
Page 316