DEM-Based Stream and Watershed Delineation

1. DEM-Based Stream and Watershed Delineation Francisco Olivera, Ph.D., P.E. Texas A&M University Department of Civil Engineering

2. Watershed Delineation • Watershed delineation is the process of identifying the drainage area of a point or set of points. • For many years, Civil Engineers have been using USGS paper maps for watershed delineation. • The maps are in Universal Transverse Mercator (UTM) coordinates and in 1:24,000 scale.

3. Watershed Delineation • Water flows in the direction of the terrain steepest downhill slope. • Streamlines are orthogonal to the elevation contour lines.

4. Watershed Delineation • Streamlines do not flow towards drainage divides, and do not intersect them. • Drainage divides are found along the highest points of the terrain.

5. DEMs • Digital Elevation Models (DEMs) are grids of elevation. • DEMs store the same type of information contour lines do, but with a different data structure. • Watershed delineation can be based on DEMs rather than contour lines.

6. 30-Meter DEMs • For the United States, 30-meter DEMs have been developed by USGS. • 30-meter DEMs are consistent with 1:24,000 paper maps.

7. 720 720 740 720 700 680 720 700 740 680 30-Meter DEMs

8. 30-Meter and 3" DEMs ~ 90m 30m

9. 32 64 128 16 1 8 4 2 Hydrologic Functions • The eight-direction pour point algorithm (D-8) assigns a flow direction code to each cell, based on the steepest downhill slope as defined by the DEM. • The flow direction code indicates the cell – out of the eight neighbor cells – towards which the water flows. • Flow direction cannot be defined for cells within a terrain depression. DEMs have to be filled before determining flow directions. Flow direction codes

10. 67 56 49 53 44 37 58 55 22 Hydrologic Functions 1 1 67 56 49 53 44 37 58 55 22 Slope:

11. 78 72 69 71 58 2 2 2 4 4 74 67 56 49 46 2 2 2 4 4 69 53 44 37 38 1 1 2 4 8 64 58 55 22 31 128 128 1 2 4 68 61 47 21 16 128 128 1 1 4 Hydrologic Functions Function: Flow direction Argument: DEM Flow network (not a grid) Flow direction grid DEM

12. Flow Direction • Water flows to one of its neighbor cells according to the direction of the steepest descent. • Flow direction takes one out eight possible values.

13. 78 72 69 71 58 74 67 56 49 46 69 53 44 37 38 64 58 55 22 31 68 61 47 21 16 0 0 0 0 0 0 1 1 2 1 0 3 8 5 2 0 1 1 20 0 0 0 0 1 24 Hydrologic Functions Function: Flow accumulation Argument: Flow direction Flow network (not a grid) Flow accumulation grid DEM

14. Flow Accumulation • Flow accumulation is a measure of the drainage area in units of grid cells. • The cell itself is not included.

15. 78 72 69 71 58 74 67 56 49 46 69 53 44 37 38 64 58 55 22 31 68 61 47 21 16 Hydrologic Functions Process: Stream delineation Argument: Flow accumulation 0 0 0 0 0 0 1 1 2 1 0 3 8 5 2 0 1 1 20 0 0 0 0 1 24 Flow accumulation (> 2 cells) DEM

16. Stream Network • All grid cells draining more than a user-defined threshold value (blue streams) are part of the stream network. • All grid cells located downstream of user-defined cells (red streams) are also part of the stream network.

17. 78 72 69 71 58 74 67 56 49 46 69 53 44 37 38 64 58 55 22 31 68 61 47 21 16 Hydrologic Functions Function: Stream links Argument: Delineated streams DEM Streams links Delineated streams • Stream links are assigned an identification number at random.

18. Stream links • Stream links are the segments of a stream channel connecting two successive junctions, a junction and an outlet, or a headwater and a junction.

19. 78 72 69 71 58 74 67 56 49 46 69 53 44 37 38 64 58 55 22 31 68 61 47 21 16 Hydrologic Functions Process: Outlets Argument: Stream links DEM Outlets Stream links • Outlets are assigned the identification number of their link.

20. Watershed Outlets • The most downstream cells of the stream segments (yellow/brown cells) are watershed outlets. • User-defined points (red dots) are also watershed outlets.

21. Modified Stream Links • User-defined outlets modify the stream links by subdividing the segment in which they are located into two segments.

22. 78 72 69 71 58 74 67 56 49 46 69 53 44 37 38 64 58 55 22 31 68 61 47 21 16 Hydrologic Functions Function: Watershed Argument: Flow direction and Outlets DEM Watersheds Outlets • Watersheds are assigned the identification number of their outlet.

23. Watershed Delineation • The drainage area of each watershed outlet is delineated.

24. Raster to Vector Conversion • The grids of stream links and watersheds are vectorized to obtain stream line and watershed polygon datasets. • Both datasets include a field Grid-Code in their attribute tables, which stores the stream link and watershed identification number. • Since the identification number of the streams and their corresponding watershed is the same, the attribute tables can be related.

25. Dissolving Spurious Polygons • Cells connected to the main watershed polygon through a corner are defined as a separate polygon (spurious polygon). • These polygons are dissolved into the main polygon.

26. Merging Watersheds • Adjacent watershed polygons can be merged into a single polygon, if they share their outlet or one flows into the other.

27. Merging Watersheds • The attribute tables are updated so that the watershed code (WshCode) of the streams is modified and the area of the new polygon is calculated.