Mcginty Methods
Mcginty Methods
ENVS421 lab 3
2/18/19
Hydrology Toolset
Preprocessing
Fill
A sink is a cell which has no cells adjacent to it that are lower in elevation. This is a
problem for hydrology modeling because if all the cells surrounding a cell are higher in
elevation, the model cannot determine which direction the water will flow once this sink fills up.
Therefore, I used the fill tool to bring the elevation of every sink up to the elevation of its lowest
neighbor before running any other hydrology tool. A z limit can be set to exclude sinks of a
certain depth from being filled, but I did not set one for this case. The input and output are both
DEMs.
The model begins with this tool.
Flow direction
The flow direction tool starts with a DEM and calculates the direction that water will
move from each cell. There are three methods this tool can use; I chose D8. The D8 method
models flow from each cell to its lowest neighboring cell. The output flow direction raster shows
flow in one of eight directions for each cell. The output drop raster shows the steepness of the
flow from cell to cell. I used my filled DEM and produced both outputs.
Figure 3. The resulting output flow direction raster from the flow direction tool and a DEM of the
south fork Nooksack River. Depicts the direction in which water will flow from each cell based on
elevation.
Flow accumulation
The flow accumulation tool calculates, for each cell, the number of cells that flow into it,
from the flow direction raster. I chose output data type “integer”. I did not use an input weight
raster and accepted the default D8 flow direction type.
Figure 4. The resulting output flow accumulation raster from the flow accumulation tool and the
flow direction raster. Depicts the number of 100 m2 cells that drain into each cell, converted to
mi2.
Raster calculator
Raster calculator is not a hydrology tool. However, I used it to define streams from the
flow accumulation output. In order to quantify this delineation, I changed the symbology type of
my flow accumulation raster to classify using two classes and overlaid my streams feature
class. In this scenario, the upper class of values represented streams and the lower class
represented not streams. By changing the upper value of the lower class and comparing to the
streams feature class, I was able to determine the value which was the break between streams
and not streams in the flow accumulation raster. In raster calculator I created the expression
“SetNull(“%SFNR_acc%” < 700, 1”, where SFNR_acc was the name of my flow accumulation
raster and 700 was the threshold value I chose to designate streams from non stream land. The
output resulted in a raster with a single value of 1 equating to every cell in the flow accumulation
raster whose value had been greater than or equal to 700. All values that had been less than
700 were set to null in the output streams raster. From here, I converted my streams raster to a
streams line feature class.
Note: if you would like to use a different threshold value to delineate your streams, you
must do so through the editing pane of the model.
Figure 5. The flow accumulation raster symbolized with two classes with a break point of 700
cells to delineate stream cells versus non stream cells, as described above.
Stream link
Stream order
Flow length
I used the flow length tool to calculate the distance of the longest continuous flow path to
the outlet of the watershed from each cell with the flow direction raster using the downstream
option. The output is a continuous raster visualizing how far water has to travel to reach the
outlet of the watershed.
Figure 8. The resulting output flow length raster from the flow length tool and the flow direction
raster. Depicts the number of cells water will travel through to reach the final drainage point of
the south fork Nooksack River watershed, from each particular cell.
Basin
The basin tool uses the flow direction raster to calculate drainage basins, or the entire
land area that flows to a single outlet. The output raster contained the entire south fork
Nooksack River watershed. I converted the raster to a polygon feature class. Some of the land
area in this feature class did not drain to the south fork Nooksack River, so I deleted these small
polygons from the dataset. This was easy to do by viewing the attribute table. In the attribute
table, one polygon had a much larger area than the rest. This polygon was the south fork
Nooksack River basin, so all other polygons were area that did not drain to the south fork
Nooksack River. I used this edited basin feature class to clip out streams that were not
incorporated in the south fork Nooksack River system from the streams feature class from the
Raster Calculator step.
Use of the model ends here. The final steps were run by hand.
Watershed
The watershed tool determines the area of cells that contribute to a particular cell. With
the pour point raster and the flow direction raster as inputs, I used the watershed tool to
delineate the subwatersheds within the area of interest. Then, with flow direction and stream link
rasters as inputs, I used the watershed tool again to delineate catchments within the basin.
Figure 9. The resulting output watershed raster from the watershed tool, the pour point raster,
and the flow direction raster. Depicts the subwatersheds of the south fork Nooksack River
watershed. Additionally, for context, the contributing streams and the pour points of each
subwatershed are displayed.
Figure 10. The resulting output watershed raster from the watershed tool, the stream link raster,
and the flow direction raster. Depicts the catchments of the south fork Nooksack River
watershed that contribute to each individual stream. Additionally, for context, the contributing
streams are displayed.
Map Credits
Data sources: USGS, WWU
Projection: NAD 1983 StatePlane Washington North FIPS 4601