NetMap's Technical Help Guide

NetMap Synthetic Stream Layer Derivation

NetMap Synthetic Stream Layer Derivation
"Synthetic" refers to a stream network derived directly from digital elevation models (DEMs) as described below. Advantages include that derived river networks are more complete, more consistent (drainage density from one watershed to another) and are designed numerically specifically with applications in mind (resource management, hazard assessments, risk mitigation, restoration and conservation).
 
NetMap uses several criteria to determine the headward extent of channels within a basin including: 1) critical drainage area (drainage area per unit contour length) or simply drainage area, as a threshold, 2) plan curvature and 3) minimum flow length. NetMap’s River Builder tool allows for delineation of channel networks and the calibration of the channel network extent or channel heads.
 
Critical drainage area is drainage area per unit contour length and various thresholds can be used to identify provisional locations of channel heads or channel extent (Figure 1).  Critical area is equal to area multiplied by slope gradient-squared (AS2). Threshold values vary with steep and less steep topography, for example, less than 25%, greater than 25% etc.  Plan curvature is the curvature of the surface on a cell-by-cell basis, as fitted through that cell and its eight surrounding neighbors. Curvature is the second derivative of the surface, or the slope-of-the-slope. The plan curvature is perpendicular to the direction of the maximum slope (Figure 2). Algorithms for flow direction and channel delineation are described by Clarke et al. 2008 but see latest Technical Help covering the River Builder.
 
 
 
Figure 1. NetMap’s channel delineation method requires calculating the drainage area per unit contour length, distinct for both steep and low gradient areas. Example is from northern California.
 
 
Figure 2. The plan curvature is shown for a basin defined as the second derivative of the surface, or the slope-of-the-slope. The plan curvature is perpendicular to the direction of the maximum slope. The plan curvature is combined with the specific drainage area (Figure 1) and a minimum channel length threshold to define channels and the upward extent of channels.  Algorithms for flow direction and channel delineation are described by Clarke et al. 2008 but see latest Technical Help covering the River Builder.
 
 
Two flow routing algorithms are used in NetMap. First D-infinity (Figure 3) is used to simulate flow on hillsides (and its convergence) and once a channel head is delineated, the D-8 algorithm is used (Figure 3).
 
Figure 3. Flow routing on hillsides (unchanneled areas) employs the D-infinity algorithm (top) and flow routing in channels uses the D-8 algorithm (from Mitasova, Geospatial Analysis and Modeling MEA592).
 
Using both specific drainage area per unit contour length and pan curvature thresholds, a user can manually calibrate the upward extent of channel heads, that is, define the channel density. Too few channels may be created or too many channels. When threshold values are incorrect, there will be “channel feathering” which is manifest as numerous, closely spaced channels extending up on planar hillslopes (Figure 3). The analyst must “calibrate” the thresholds to eliminate or reduce channel feathering (using visual clues of channel locations in landscapes) while providing for a high density of channels. This is done by defining the inflection point between specific drainage area and channel density, given that the threshold for plan curvature is met (Figure 4). The pixel cells that are classified as being a “channel” must extend contiguously for a certain minimum distance, often ranging between 30 to 50 meters, set by the analyst.
 
 
 
Figure 4. Synthetic stream networks derived from DEMs can have too high a drainage density. Here is an example of a synthetic channel network with “channel feathering”, defined as closely spaced channels on planar hillslopes.
 
 
Figure 5. To determine the headward extent of the channel network in NetMap requires calculating the relationship between specific drainage area and channel density. There is an inflection point that forms that indicates the transition to an environment of channel feathering (e.g., Figure 3).
 
 
Once the appropriate thresholds are set for plan curvature, specific drainage area (Figure 5) and critical channel length, the channel network is defined within NetMap (Figure 6).
 
Figure 6. An appropriate channel network for a basin in northern California is created by adjusting thresholds for plan curvature and specific drainage area (Figure 5) so that no channel feathering has occurred. Contrast this network with the inaccurate network shown in Figure 4. The NHD is also shown for comparison, and the NetMap network has a higher network density, and thus more headwater channels.
 
The overall strategy in NetMap is to create high density drainage networks so as not to eliminate headwater streams, since information about network extent (or channel heads) is generally absent for most landscapes. With a liberal network within NetMap tools, a user can remove headwater channels that do not exist, and the network is renumbered. The ‘Channel Head’ tool is located in the Fluvial Morphology Module. In this way, a user can customize NetMap’s stream layer to better match specific landscape locations. Another way to tailor NetMap’s stream layer, in terms of drainage density (channel heads) and other features such as segment breaks and segment lengths scales is to use the River Builder tool. These tools allow for more flexibility and tailoring of a stream layer to specific landscapes than non-synthetic networks, like the NHD (described in more detail below).
 
DEM derived stream networks depend on topographic enforcement of the derived channel network, and thus synthetic networks are more accurate in higher relief terrain with topographic controls on channel position.  In areas of low relief such as in wide low gradient meandering or braided rivers, DEM derived networks may be inaccurate, since it is difficult for a computer model to determine exact locations of channels in wide valleys of little relief (e.g., Figure 6). Thus, NetMap’s stream layer can be mirrored onto other stream layers, such as the cartographic NHD (derived from imagery by mapping personnel).  The analyst must define what gradient threshold to use to “enforce” NetMap’s stream layer using another stream layer, e.g., the NHD. However, because of the dynamic shifting behavior of river channels across wide floodplains, even the NHD become spatially inaccurate over time (Figure 7).
 
 
 
Figure 7. (Top) NetMap’s stream layer compared to actual stream location in the Hoh River basin, Washington. (Bottom) The NHD can be inaccurate also in low gradient and meandering or braided river channels because of the dynamic shifting of the river over time.
 
When developing a stream layer in NetMap (using the River Builder), a user can select the length scale of stream segments and the criteria to delineate segments. For options, see Technical Help for the River Builder.  Often, segment lengths of 50 to 200 meters are created, and their lengths vary, determined by a user selected option that allows for channel segments to be delineated based on homogeneity in physical properties such as channel gradient, drainage area and width. Thus, many segments are broken at small to large tributary junctions (Figure 8). The network is routed, e.g., each segment “knows” its position in the network, with respect to both upstream and downstream positions. This routed network, allows for information at the channel reach scale to be summarized both upstream and downstream; information can include channel attributes (such as fish habitat potential or road density), hillslope attributes such as erosion potential, vegetation age or fire risk indices), or hillslope information captured by varying buffer widths (see Cumulative Upstream Aggregation and Cumulative Downstream Aggregation for additional information.
 
Figure 8. NetMap’s synthetic stream network contains stream segments (variable length), tributary junction points and drainage wings. Drainage wings identify the local contributing area to the reach on each side of the valley (right and left side are delineated).
 
Based on the flow direction grid, “drainage wings” are constructed on both sides of the channel segment thereby delineating the local cumulative drainage area on both sides of the channel (Figure 8). This allows information that is contained in the terrestrial landscape, such as erosion potential, vegetation age, fire risk indices, burn severity rankings (all raster or grid information), and roads to be reported to individual channel segments (Figure 9). This NetMap feature allows unique spatial overlap analyses, such as identifying where the best aquatic habitats are spatially overlapping with the highest risk hillslopes (see Risk Analysis-Overlap tool). Information reported to streams is routed downstream revealing patterns of terrestrial attributes (or channel attributes) at any spatial scale defined by the stream network.
 
Figure 9. Drainage wings in NetMap allows for information in the terrestrial landscape to be reported to stream channels and then routed downstream, revealing patterns at any spatial scale defined by the channel network.
 
With drainage wings, entire watersheds are discretized into drainage wing facets (Figure 10). Using an average channel length of approximately 100 m, drainage wings average about 0.1 km2 in area. This allows for a finer scale analysis of watershed attributes such as road density (Figure 11).
Figure 10. An example of drainage wings are shown for a basin in California. Using channel segments that average 100 m in length, drainage wing areas are approximately 0.1 km2.
Figure 11. Calculating road density at the stream segment scale provides a much higher resolution mapping of variations in road density. For example, road density in the Clearwater Basin that is calculated at the scale of HUC 6th field subbasins ranged between 0.5 km/km2 and 4 km/km2 (top). In contrast road density measured at the stream segment scale ranged between zero and 100 km/km2 based on the use of drainage wings.
 
 

Made with help of Dr.Explain

Copyright TerrainWorks 2014