Three words that are not synonyms

Most flawed risk maps go wrong at the vocabulary stage. Susceptibility is the spatial likelihood that a hazard could occur somewhere, given the static terrain conditions — it answers where. Hazard adds magnitude and temporal probability — where, how big, how often. Risk multiplies hazard by exposure (the people and assets present) and vulnerability (how badly they would be affected) — what is actually at stake. A susceptibility map labelled "risk map" has skipped two whole dimensions. Terrain risk mapping basics start here: know which of the three you are producing, and say so. This article walks through building a susceptibility surface from terrain and then extending toward genuine risk.

The conditioning factors

Terrain susceptibility (for landslides, erosion, or flooding-related processes) is modelled from conditioning factors — the stable terrain attributes that predispose ground to a hazard. The usual set, mostly derived from a DEM, is:

  • Slope — the dominant control for mass movement. Derive in degrees: gdaldem slope dem_utm.tif slope.tif -alg Horn.
  • Aspect — proxies for moisture, vegetation, and freeze-thaw exposure.
  • Curvature — concave hollows concentrate water; convex spurs shed it.
  • Distance to drainage / streams — toe erosion and saturation.
  • Geology and soils — material strength and permeability.
  • Land cover — vegetation root reinforcement versus bare or cleared ground.
  • Topographic wetness index — where water accumulates.

All raster factors must be reprojected to a single metric CRS (e.g. EPSG:32633) and resampled to a common cell size before they can be combined.

Weighted overlay: the workhorse method

When you do not have a large inventory of past events to train a statistical model, the transparent default is a weighted linear overlay:

  1. Reclassify each factor onto a common ordinal scale (say 1–5, low to high contribution). For slope you might map 0–5 degrees → 1, up to >35 degrees → 5; for distance to drainage, near → high.
  2. Assign weights that sum to 1, reflecting each factor's relative influence (slope usually highest). The Analytic Hierarchy Process is a structured way to derive weights from pairwise expert comparison, and it produces a consistency ratio so you can check the weights are internally coherent.
  3. Combine: susceptibility = Σ (weight × reclassified factor). Use the QGIS Raster calculator, the SAGA weighted overlay, or ArcGIS Weighted Overlay.
# QGIS raster calculator, illustrative weights:
0.45*"slope_rc" + 0.20*"geology_rc" + 0.15*"wetness_rc" +
0.10*"landcover_rc" + 0.10*"dist_drain_rc"

Where you do have an event inventory, statistical methods (frequency ratio, weights-of-evidence, logistic regression) derive the weights from data rather than expert judgement and let you quantify performance.

From susceptibility to risk

To go beyond where a hazard is likely:

  • Hazard: combine susceptibility with magnitude/frequency information — for example, return-period flood depths from a hydraulic model, or landslide recurrence from a dated inventory.
  • Exposure: overlay assets — buildings, roads, lifelines, population grids.
  • Vulnerability: assign damage functions or vulnerability classes to those assets.

Risk is then, conceptually, hazard × exposure × vulnerability, evaluated cell by cell or per asset. Keep the components as separate layers so each can be inspected and revised; a single fused "risk" raster with no traceable inputs cannot be defended.

Validation against an inventory

A susceptibility map is only credible if it predicts known events. With an inventory of past landslides (or flood extents), the standard check is the ROC / AUC approach: rank cells by predicted susceptibility, and measure what fraction of known events fall in the highest-susceptibility classes. An AUC near 0.5 is no better than chance; values approaching 0.8 and above indicate a useful model. Always hold back some events for validation rather than fitting and testing on the same data.

Worked example: a regional landslide susceptibility screen

  1. Reproject GLO-30 to EPSG:32633; derive slope, curvature, and topographic wetness index.
  2. Add geology and land cover; compute distance to drainage with the QGIS Proximity tool.
  3. Reclassify each factor to 1–5.
  4. Derive weights by AHP (slope dominant), check the consistency ratio.
  5. Combine in the raster calculator into a 1–5 susceptibility surface.
  6. Overlay the landslide inventory; compute the success-rate curve and AUC.
  7. Symbolise in five classes, masking out flat valley floors, and document every weight and breakpoint.

Common pitfalls and why they happen

  • Calling susceptibility "risk" — it omits exposure and vulnerability and overstates what the map says.
  • Slope in the wrong units or on a geographic CRS — distorts the single most important factor.
  • No validation — without an inventory check, weights are unaccountable opinion.
  • Mismatched grids — combining factors at different cell sizes or CRSs silently misaligns them.
  • Over-smoothing or over-resolving — a 90 m DEM blurs the slopes that drive shallow failures; a noisy 1 m DEM injects spurious steepness. Match resolution to the process.

QA and validation checklist

  • All factor layers share CRS, extent, and cell size.
  • Reclassification breakpoints and weights are documented and justified.
  • Susceptibility validated against a held-out event inventory with a reported AUC.
  • The legend states clearly that the product is susceptibility, hazard, or risk — not loosely "risk".
  • Source, date, and scale recorded for every input.

Bathyl perspective

We keep hazard, susceptibility, and risk as separate, validated layers, and we publish the weights and the inventory check alongside the map. A terrain risk screen earns trust by being honest about which question it answers and by showing it predicts the events that already happened.

Related reading

Sources