ITCILWISApplications Guide6. Deterministic landslide hazard zonation


6. Deterministic landslide hazard zonation

In this exercise, a simple slope stability model (the infinite slope model) is used to calculate safety factor maps for different conditions. The effect of groundwater depth and seismic acceleration is evaluated using input maps of these factors for different return periods of rainfall (related to the groundwater level) and earthquakes (related to the seismic acceleration).

In ILWIS, the model is represented by a user-defined function. Different scenarios are calculated by changing the variables of this function. The model is applied on a data set of the city of Manizales, in central Colombia.

Introduction and data preparation

Two of the input maps are a slope map and a map of ash thickness:

Map SLOPE, slopes in degrees Map ASHT, thickness of volcanic ashes

The degree of slope hazard can be expressed by the Safety Factor (F) which is the ratio of the forces that make a slope fail and those that prevent a slope from failing.

To simplify matters, the infinite slope model is used. This is a 2-dimensional model describing the stability of slopes with an infinitively large failure plane.

The depth of the possible failure plane is taken at the contact of the volcanic ashes and the underlying material. The consequence of that is that safety factors will not be calculated for the entire area, but only for the areas where there is volcanic ash at the surface.

The formula to calculate the Safety Factor (Brunsden and Prior, 1979) reads:

Safety Factor (Brunsden and Prior, 1979)


Known average values (laboratory) and variables:

c’ = effective cohesion (Pa= N/m2).


g = unit weight of soil (N/m3).

variable Gamma

m = zw/z (dimensionless).

variable m

gw = unit weight of water (N/m3).


z = depth of failure surface below the surface (m).

map ASHT

zw = height of watertable above failure surface (m).

b = slope surface inclination (°).


j' = effective angle of shearing resistance (°).

30°, thus tan j' = 0.58

With some Mapcalc calculations in which you will use the inbuilt ILWIS functions DEGRAD( ), SIN( ), COS( ), and SQ( ), first some maps are prepared that will be frequently used in this application:

Creating a user-defined function for the infinite slope model, and calculating Safety Factors for various groundwater scenarios

The Safety Factor formula as presented above will be transformed into a user-defined function FS. This function already contains the known parameters (maps ASHT, SI, CO, CO2 and the known constants) but it also contains the variables Gamma and m.

The function can then be easily applied for various heights of the watertable (zw), by filling out the variables Gamma and m. This will result in a number of Safety Factor maps for specific watertable heights.

Function FS reads:

(10000+((Gamma-m*10000)*ASHT*CO2*0.58)) / (Gamma*ASHT*SI*CO)Understanding relations between m, zw and z

The maps are (re)classified according to:

Dry condition Saturated condition

Subsequently, function FS is applied repeatedly on the command line and Safety Factor maps are calculated for a number of watertable height scenarios. For each scenario, a value for Gamma is provided and the values for m are found in a number of existing groundwater maps (M016, M1, M20, M50).

Safety Factors are calculated for groundwater peaks that occur:

during two months a year

(rather dry):



once a year:



once in every 20 years:



once in every 50 years

(very wet):



Bar graph showing stability classes for various groundwater peaksOf all maps, histograms are calculated and the area percentages of the stability classes are joined into one table. From the table, a bar graph is presented. The graph displays the relative size of the areas that are found to be unstable, critical and unstable slopes for different groundwater level scenarios.

Evaluating the effect of seismic acceleration

Besides groundwater fluctuations, also earthquake acceleration will be taken into account.

The formula for the infinite slope model changes into:

Calculating Safety Factor including earthquake acceleration

Compared to the original formula (see above), the new parameters are:

zw can be replaced with m*z, as m=zw/z.

This formula thus contains the known parameters (maps ASHT, SI, CO, CO2 and the known constants) but it also contains the variables gamma, m, rho and alpha.

Another user-defined function has to be created (FSNEW), which may read:

(10000+(ASHT*Gamma*CO2-ASHT*Rho*Alpha*CO*SI-10000*m*ASHT*CO2)*0.58) / (ASHT*Gamma*SI*CO+ASHT*Rho*Alpha*CO2)

With user-defined function FSNEW and some additional information, Safety Factors are then calculated for the following scenarios:

  1. Dry condition, without earthquake.
  2. Ground water return period of 20 years, without earthquake.
  3. Completely saturated situation, without earthquake.
  4. Dry condition, and earthquake of 7 Magnitude.
  5. Groundwater with return period of 0.16 years, and earthquake of 7 Magnitude.
  6. Groundwater with return period of 20 years, and earthquake of 7 Magnitude.

In order to compare the output maps, all maps are (re)classified according to:

Results of scenarios 5 and 6 are displayed below:

Slope stability classes for scenario 5 Slope stability classes for scenario 6

For more information on this case study, contact:

C.J. van Westen

Department of Earth Systems Analysis,

International Institute for Geo-Information Science and Earth Observation (ITC),

P.O. Box 6, 7500 AA Enschede, The Netherlands.

Tel: +31 53 4874263, Fax: +31 53 4874336, e-mail: