ENCN205 — Lecture 08: Hydrological Analysis

Interactive Demo: From DEM to Catchment — sinks, D8 flow direction, flow accumulation, streams & pour points
~10 min activity
Computing flow routing…
Click a stream to place a pour point · scroll to zoom · drag to pan
▲
N
0500 m

Charteris Bay & Mt Herbert

Real LINZ 8 m DEM · 500 × 500 cells · every step computed live

View

Hover any cell to see its 3 × 3 neighbourhood and where D8 sends its water (below).

1Condition the DEM

Vertical precision the DEM is stored at:
Sinks in raw DEM (no downhill neighbour):–
Cells raised to drain:–
Deepest fill:–
Sinks are filled: every cell drains, step by step, to the sea or the map edge.

2Extract streams

Stream threshold: –
10 cells (dense)50,000 cells (main stems)
Stream network length:–
Drainage density:–

3Delineate a catchment

Click on the map to drop a pour point (outlet). All cells that drain to it are traced upstream.

Cell inspector (D8)

Hover over the map…
The workflow you are running:
DEM→Fill sinks→D8 flow direction→Flow accumulation→Threshold→Streams→Pour point→Catchment
Each step uses the output of the one before, so an error early on (an unfilled sink, a pour point off the channel) carries through to the catchment boundary.

Guided Questions

Q1 — Accumulation

Switch to Flow accumulation. Hover over a ridge crest, then a valley floor. What values do you get? Now look at Flow direction in the inspector: which neighbour does D8 pick, and why that one?

Q2 — Threshold

Turn on the OpenStreetMap streams. Find the threshold where the extracted network best matches them. Write down the threshold and the network length. Where does the DEM put a stream that the map doesn't have, and which one is wrong?

Q3 — Sinks

Untick Fill sinks and note the number of sinks. Now switch the precision to 0.1 m, then 1 m. It is the same terrain, so where did the extra sinks come from? Place a pour point near the coast on the largest stream and compare its area with filling on and off. Are most sinks real depressions or DEM artefacts?

Q4 — Pour points

Drop a pour point on a stream. Then untick Snap and click one or two cells beside the channel. Why does the catchment collapse? How does this explain a delineated catchment boundary that runs straight through a gully?

Q5 — Culvert design

Imagine a coastal road crossing the main stream. Delineate its catchment and read the rational-method peak flow. Move the culvert 500 m upstream: how do the area and Q change? Which input (C, i or A) would you trust least, and why?