{"id":2171,"date":"2026-09-07T19:40:42","date_gmt":"2026-09-07T19:40:42","guid":{"rendered":"https:\/\/commons.mtholyoke.edu\/emarcano\/?page_id=2171"},"modified":"2026-09-21T17:48:05","modified_gmt":"2026-09-21T17:48:05","slug":"ex2-terrain","status":"publish","type":"page","link":"https:\/\/commons.mtholyoke.edu\/emarcano\/mhcgis\/geog323\/ex2-terrain\/","title":{"rendered":"Exercise 2 &#8211; Terrain analysis"},"content":{"rendered":"\n<h2 id=\"top\" class=\"wp-block-heading\"><strong>Terrain analysis<\/strong><\/h2>\n\n\n\n<h3 class=\"wp-block-heading\">Contents<\/h3>\n\n\n\n<ul class=\"wp-block-list\">\n<li><strong><a href=\"#part1\">Part 1: <strong>Watershed analysis<\/strong><\/a><\/strong>\n<ul class=\"wp-block-list\">\n<li><strong>Objective: <\/strong>To delineate the watersheds of the Dean and Roaring Brooks and the west branch of the Swift River, and calculate their respective areas<\/li>\n<\/ul>\n<\/li>\n\n\n\n<li><strong><a href=\"#part2\">Part 2: Stream ordering<\/a><\/strong>\n<ul class=\"wp-block-list\">\n<li><strong>Objective<\/strong>: To outline the streams in a watershed and classify them according to two methods<\/li>\n<\/ul>\n<\/li>\n\n\n\n<li><strong><a href=\"#part3\">Part 3: Viewshed<\/a><\/strong>\n<ul class=\"wp-block-list\">\n<li><strong>Objective<\/strong>: To determine how much of a watershed is visible from an observation tower<\/li>\n<\/ul>\n<\/li>\n<\/ul>\n\n\n\n<hr class=\"wp-block-separator has-text-color has-black-color has-alpha-channel-opacity has-black-background-color has-background is-style-wide\" \/>\n\n\n\n<p class=\"wp-block-paragraph\"><a href=\"https:\/\/commons.mtholyoke.edu\/emarcano\/geog323\">Return to Geog 323 page<\/a><\/p>\n\n\n\n<hr class=\"wp-block-separator has-text-color has-black-color has-alpha-channel-opacity has-black-background-color has-background is-style-wide\" \/>\n\n\n\n<h2 id=\"part1\" class=\"wp-block-heading\">Part 1: Watershed analysis<\/h2>\n\n\n\n<p class=\"wp-block-paragraph\"><em>Adapted from Spatialabs 2013<\/em><\/p>\n\n\n<div class=\"wp-block-image\">\n<figure class=\"alignright size-full is-resized\"><a href=\"https:\/\/commons.mtholyoke.edu\/emarcano\/wp-content\/uploads\/sites\/965\/2026\/03\/watersheds.png\"><img loading=\"lazy\" decoding=\"async\" width=\"486\" height=\"210\" src=\"https:\/\/commons.mtholyoke.edu\/emarcano\/wp-content\/uploads\/sites\/965\/2026\/03\/watersheds.png\" alt=\"watershed components\" class=\"wp-image-1727\" style=\"width:421px;height:auto\" srcset=\"https:\/\/commons.mtholyoke.edu\/emarcano\/wp-content\/uploads\/sites\/965\/2026\/03\/watersheds.png 486w, https:\/\/commons.mtholyoke.edu\/emarcano\/wp-content\/uploads\/sites\/965\/2026\/03\/watersheds-300x130.png 300w\" sizes=\"auto, (max-width: 486px) 100vw, 486px\" \/><\/a><\/figure>\n<\/div>\n\n\n<p class=\"wp-block-paragraph\">A watershed or drainage basin describes the area of land from which all surface waters (rain, melting snow or ice) converge at a single point. This point is the exit of the watershed and usually forms a stream or lake. Watersheds are important hydrologic and environmental units because all water within one watershed is connected. Thus, water pollution at one point in the watershed will affect the entire area downstream of that point and thus becomes an important issue in environmental restoration and protection. Since water always flows from the highest point to the lowest, watershed boundaries (drainage divides) are defined as the ridges separating flow.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\"><strong>The problem<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\">The geology dept. at MHC has an ongoing study on the surface and subsurface hydrology in the town of Shutesbury, located just west of the Quabbin reservoir.&nbsp; In this town the US Geological Survey (USGS) has installed 3 river gauges that track the surface water flow at Dean and Roaring Brooks, and in the west branch of the Swift River, a major source of water to the Quabbin.&nbsp; By determining the extent of the land area that contributes water to these streams, and eventually to the gauges, we can study different aspects of the hydrology of the area.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\"><strong>Objectives<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\">To delineate the watersheds of the Dean and Roaring Brooks, and the west branch of the Swift River, and calculate their respective areas.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\"><strong>Data description<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\">The dataset in <em>Watershed.gdb<\/em> in your working folder contains a 30m DEM and the streams of the state of Massachusetts, these will need to be clipped to our Area of Interest (<em>AOI<\/em>):<\/p>\n\n\n\n<ol style=\"list-style-type:lower-alpha\" class=\"wp-block-list\">\n<li>30m DEM of Massachusetts (<em>DEM_Mass_30m)<\/em> from MassGIS<\/li>\n\n\n\n<li>Streams (<em>HYDRO25K_ARC<\/em>) in Massachusetts from MassGIS<\/li>\n\n\n\n<li>Pour points (USGS Gauges), both vector (<em>Gauges<\/em>) and raster (<em>pour_points_raster<\/em>). The raster version is hard to see so the vector has been provided for visualization, the analysis will be done with the raster<\/li>\n\n\n\n<li><em>AOI<\/em> is a polygon that covers our area of interest (AOI), used to clip the DEM and streams to make them more manageable<\/li>\n<\/ol>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\" \/>\n\n\n\n<h3 class=\"wp-block-heading\"><strong>Task 0 \u2013 Preliminary <\/strong>tasks<\/h3>\n\n\n\n<ol class=\"wp-block-list\">\n<li>Load the data from the <em>Watershed<\/em> GeoDB and study it<\/li>\n\n\n\n<li>Clip the <em>HYDRO25K_ARC<\/em> using the <em>AOI<\/em>, call the result <em>HYDRO25K_ARC_AOI <\/em>(Hint: vector CLIP function)<\/li>\n\n\n\n<li>To Clip the DEM use the function <strong>Toolbox\\Spatial Analyst Tools\\ Extraction\\ Extract by Mask<\/strong><\/li>\n<\/ol>\n\n\n\n<ol style=\"list-style-type:lower-alpha\" class=\"wp-block-list\">\n<li>Input raster \u2013 <em>DEM_Mass_30m<\/em><\/li>\n\n\n\n<li>Mask data \u2013 <em>AOI<\/em><\/li>\n\n\n\n<li>Output \u2013 <em>DEM_AOI<\/em><\/li>\n<\/ol>\n\n\n\n<ol start=\"4\" class=\"wp-block-list\">\n<li>Create a hillshaded surface <strong>Toolbox\\Spatial Analyst Tools\\Surface\\Hillshade<\/strong><\/li>\n<\/ol>\n\n\n\n<ol style=\"list-style-type:lower-alpha\" class=\"wp-block-list\">\n<li>As input use <em>DEM_AO<\/em>I<\/li>\n\n\n\n<li>Output <em>Hillshade_AOI<\/em><\/li>\n\n\n\n<li>Leave the rest as default (but click on the help\/Info buttons to know what they are)<\/li>\n<\/ol>\n\n\n\n<p class=\"wp-block-paragraph\">Use this hill shaded theme to visualize your results. It can be made partly transparent and overlap over a nicely colored DEM for a good visual effect<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">When done turn off, better yet, remove, the statewide data. You don\u2019t want to keep drawing those big datasets, they slow down the software.<\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\" \/>\n\n\n\n<p class=\"wp-block-paragraph\">To delineate watersheds in GIS we require several steps including (detail instructions are shown below)<\/p>\n\n\n\n<ol style=\"list-style-type:upper-roman\" class=\"wp-block-list\">\n<li><strong>Filling sinks<\/strong> in the elevation model<\/li>\n\n\n\n<li>Determining the <strong>direction of water flow<\/strong> through each grid cell<\/li>\n\n\n\n<li>Identifying stream channels, <strong>flow accumulation<\/strong><\/li>\n\n\n\n<li>Determining the set of adjoining cells that contribute water to a particular pour point <strong>watershed<\/strong>.<\/li>\n<\/ol>\n\n\n\n<p class=\"wp-block-paragraph\">This process requires several input datasets derived from the DEM that need to be created first. These are: <em>filled sinks, flow direction, and water accumulation<\/em>.&nbsp; And we also need the location of the gauges or <em>pour points<\/em>.<\/p>\n\n\n\n<h3 class=\"wp-block-heading\">Task 1- Create a hydrologically correct DEM by filling any artificial holes (sinks) in the data.<\/h3>\n\n\n<div class=\"wp-block-image\">\n<figure class=\"alignright size-full is-resized\"><a href=\"https:\/\/commons.mtholyoke.edu\/emarcano\/wp-content\/uploads\/sites\/965\/2026\/03\/fill-sinks.png\"><img loading=\"lazy\" decoding=\"async\" width=\"332\" height=\"103\" src=\"https:\/\/commons.mtholyoke.edu\/emarcano\/wp-content\/uploads\/sites\/965\/2026\/03\/fill-sinks.png\" alt=\"\" class=\"wp-image-1728\" style=\"width:361px;height:auto\" srcset=\"https:\/\/commons.mtholyoke.edu\/emarcano\/wp-content\/uploads\/sites\/965\/2026\/03\/fill-sinks.png 332w, https:\/\/commons.mtholyoke.edu\/emarcano\/wp-content\/uploads\/sites\/965\/2026\/03\/fill-sinks-300x93.png 300w\" sizes=\"auto, (max-width: 332px) 100vw, 332px\" \/><\/a><\/figure>\n<\/div>\n\n\n<p class=\"wp-block-paragraph\">Sinks represent depressions (real or erroneous) in your data that impede water flow and therefore must be removed. &nbsp;The software fills the sinks by raising their elevations to equal neighboring cells so that water can continue to flow across them.<\/p>\n\n\n\n<ol start=\"5\" class=\"wp-block-list\">\n<li>Fill the sinks in the DEM<\/li>\n<\/ol>\n\n\n\n<ol style=\"list-style-type:lower-alpha\" class=\"wp-block-list\">\n<li><strong>Toolbox \u00bb Spatial Analyst Tools \u00bb Hydrology \u00bb Fill<\/strong><\/li>\n\n\n\n<li>Input surface raster &#8211; <em>DEM_AOI<\/em><\/li>\n\n\n\n<li>Output surface raster-<em> filled_DEM<\/em><\/li>\n<\/ol>\n\n\n\n<ol start=\"6\" class=\"wp-block-list\">\n<li>Assign the same color scheme to both <em>filled_DEM<\/em> and <em>DEM_AOI <\/em>and compare them (there are several color schemes suitable for elevation).<\/li>\n<\/ol>\n\n\n\n<p class=\"wp-block-paragraph\">Can you see any difference between the two DEM\u2019s?&nbsp; Don\u2019t worry, puny humans can\u2019t tell the difference between these but the computer can \ud83d\ude42<\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\" \/>\n\n\n\n<h3 class=\"wp-block-heading\">Task 2 &#8211; Determine which way water flows from cell to cell<\/h3>\n\n\n\n<ol start=\"7\" class=\"wp-block-list\">\n<li>Determine flow direction for your filled DEM<\/li>\n<\/ol>\n\n\n\n<ol start=\"1\" style=\"list-style-type:lower-alpha\" class=\"wp-block-list\">\n<li><strong><em>Spatial Analyst Tools &gt; Hydrology &gt; Flow Direction<\/em><\/strong><\/li>\n\n\n\n<li>Input here is <em>filled_DEM<\/em><\/li>\n\n\n\n<li>Output &#8211; <em>flow_dir<\/em><\/li>\n<\/ol>\n\n\n\n<p class=\"wp-block-paragraph\">Water always flows from high values (elevation) to low values.&nbsp; The flow direction layer looks strange but if you look at it carefully it will start to make sense.&nbsp;<\/p>\n\n\n<div class=\"wp-block-image\">\n<figure class=\"alignleft size-full is-resized\"><a href=\"https:\/\/commons.mtholyoke.edu\/emarcano\/wp-content\/uploads\/sites\/965\/2026\/03\/directions.png\"><img loading=\"lazy\" decoding=\"async\" width=\"251\" height=\"266\" src=\"https:\/\/commons.mtholyoke.edu\/emarcano\/wp-content\/uploads\/sites\/965\/2026\/03\/directions.png\" alt=\"\" class=\"wp-image-1742\" style=\"width:192px;height:auto\" \/><\/a><\/figure>\n<\/div>\n\n\n<p class=\"wp-block-paragraph\">The legend will have numbers in from 1 to 128 according to the figure on the right.&nbsp; Cells with values of 1 are cells flowing to the East, south flowing cells have a value of 4, 16 to the west, and 64 to the north.&nbsp; You could change the legend values using the symbology pane to make better sense of the image. See figure below for a suggested legend.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Look for areas of contrasting color changes and study the legend, you should be able to see where the flow of water changes from \u2018East to West\u2019 (i.e. west direction) to \u2018West to East\u2019 (east direction).&nbsp; See figure on the right for example and the figure below for a concept of water flow, remember the water will flow from higher to lower elevations.<\/p>\n\n\n\n<div style=\"height:4px\" aria-hidden=\"true\" class=\"wp-block-spacer\"><\/div>\n\n\n<div class=\"wp-block-image\">\n<figure class=\"aligncenter size-full\"><a href=\"https:\/\/commons.mtholyoke.edu\/emarcano\/wp-content\/uploads\/sites\/965\/2026\/09\/flow-directon.jpg\"><img loading=\"lazy\" decoding=\"async\" width=\"650\" height=\"312\" src=\"https:\/\/commons.mtholyoke.edu\/emarcano\/wp-content\/uploads\/sites\/965\/2026\/09\/flow-directon.jpg\" alt=\"\" class=\"wp-image-2217\" srcset=\"https:\/\/commons.mtholyoke.edu\/emarcano\/wp-content\/uploads\/sites\/965\/2026\/09\/flow-directon.jpg 650w, https:\/\/commons.mtholyoke.edu\/emarcano\/wp-content\/uploads\/sites\/965\/2026\/09\/flow-directon-300x144.jpg 300w\" sizes=\"auto, (max-width: 650px) 100vw, 650px\" \/><\/a><\/figure>\n<\/div>\n\n\n<div style=\"height:8px\" aria-hidden=\"true\" class=\"wp-block-spacer\"><\/div>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\" \/>\n\n\n\n<h3 class=\"wp-block-heading\">Task 3 &#8211; Determine where water accumulates<\/h3>\n\n\n\n<ol start=\"8\" class=\"wp-block-list\">\n<li>To determine where water accumulates and creates channels, go to:<\/li>\n<\/ol>\n\n\n\n<ol style=\"list-style-type:lower-alpha\" class=\"wp-block-list\">\n<li><strong>Spatial Analyst Tools &gt; Hydrology &gt; Flow Accumulation<\/strong>&nbsp;<\/li>\n\n\n\n<li>Use <em>flow_dir<\/em> as input<\/li>\n\n\n\n<li><em>flow_acc <\/em>as output.<\/li>\n\n\n\n<li>Output data type &#8211; Integer<\/li>\n<\/ol>\n\n\n\n<ol start=\"9\" class=\"wp-block-list\">\n<li>Invert the color scheme to show a water accumulation (streams) darker on a white background (Double click the color rectangle (legend) of <em>flow_acc<\/em> to go to the Symbology pane and tick Invert color box).<\/li>\n<\/ol>\n\n\n\n<p class=\"wp-block-paragraph\">You can just barely make out some of the stream channels. Every cell has a numerical value equal to the number of cells that \u201cprovide\u201d water to it (see illustrations)<\/p>\n\n\n<div class=\"wp-block-image\">\n<figure class=\"aligncenter size-full is-resized\"><a href=\"https:\/\/commons.mtholyoke.edu\/emarcano\/wp-content\/uploads\/sites\/965\/2026\/09\/flow-accumulation.jpg\"><img loading=\"lazy\" decoding=\"async\" width=\"909\" height=\"273\" src=\"https:\/\/commons.mtholyoke.edu\/emarcano\/wp-content\/uploads\/sites\/965\/2026\/09\/flow-accumulation.jpg\" alt=\"\" class=\"wp-image-2208\" style=\"aspect-ratio:3.3295896529989975;width:674px;height:auto\" srcset=\"https:\/\/commons.mtholyoke.edu\/emarcano\/wp-content\/uploads\/sites\/965\/2026\/09\/flow-accumulation.jpg 909w, https:\/\/commons.mtholyoke.edu\/emarcano\/wp-content\/uploads\/sites\/965\/2026\/09\/flow-accumulation-300x90.jpg 300w, https:\/\/commons.mtholyoke.edu\/emarcano\/wp-content\/uploads\/sites\/965\/2026\/09\/flow-accumulation-768x231.jpg 768w\" sizes=\"auto, (max-width: 909px) 100vw, 909px\" \/><\/a><\/figure>\n<\/div>\n\n\n<p class=\"wp-block-paragraph\">This raster layer could serve as a proxy for a hydrology layer (streams) if we don\u2019t have this as it shows areas where water accumulates over the terrain.&nbsp; The value of a cell reflects the number of cells, above it, that contribute water to that particular cell.&nbsp; You might need to reclassify the raster to \u201cdetermine\u201d what can be considered a stream and what not (a 0\/1 raster).&nbsp; Also, you might have to vectorize the result.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Using this raster, you can also do some extra hydrological analysis like stream orders.<\/p>\n\n\n\n<ol start=\"10\" class=\"wp-block-list\">\n<li>Turn on (and off) the <em>Hydro25K_Arc_AOI<\/em> layer and notice the similarities and differences.<\/li>\n\n\n\n<li>Also, click on different water accumulation pixels in the raster layer (turn off all layers except Flow Accumulation). Notice how the pixel values change in different parts of the streams. Darker pixels should have larger values than lighter ones.&nbsp; These numbers represent the numbers of cells that contributed water to that particular pixel where you clicked.<\/li>\n\n\n\n<li>Also verify that the theme <em>Gages<\/em> and <em>pour_points_raster<\/em> fall on top of the water accumulation channels.<\/li>\n<\/ol>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\" \/>\n\n\n\n<h3 class=\"wp-block-heading\">Task 4 &#8211; Determine watersheds<\/h3>\n\n\n\n<p class=\"wp-block-paragraph\">You are finally ready to delineate watersheds.<\/p>\n\n\n\n<ol start=\"13\" class=\"wp-block-list\">\n<li><strong>Spatial Analyst Tools &gt; Hydrology &gt; Watershed<\/strong><\/li>\n<\/ol>\n\n\n\n<ol style=\"list-style-type:lower-alpha\" class=\"wp-block-list\">\n<li><em>flow_dir<\/em> as your Input flow direction raster<\/li>\n\n\n\n<li><em>pour_points_raster<\/em> as Input raster or feature pour point data (check note 2 below)<\/li>\n\n\n\n<li>Name the output raster <em>Watersheds<\/em><\/li>\n\n\n\n<li>Disregard Pour point field<\/li>\n<\/ol>\n\n\n\n<p class=\"wp-block-paragraph\">This should create 3 watersheds, one for each gage site. You might want to add a \u201cName\u201d column in the attribute table of <em>Watershed<\/em> and add the corresponding name according to the info in <em>Gages<\/em> to facilitate answering the questions bellow.<\/p>\n\n\n\n<p class=\"has-pale-cyan-blue-background-color has-background wp-block-paragraph\"><strong>NOTE 1<\/strong>: In the Hydrology toolbox there\u2019s a function called \u201cBasin\u201d. This is similar to Watershed but it doesn\u2019t require pour points. It evaluates the whole area.&nbsp; Sometimes this is convenient but sometimes it gives strange results. Try it and compare the results.<\/p>\n\n\n\n<p class=\"has-pale-cyan-blue-background-color has-background wp-block-paragraph\"><strong>NOTE 2<\/strong>: In this exercise you were given the pour points as vectors (<em>Gauges) <\/em>and raster (<em>pour_points_raster).<\/em> Either version can be used in the Watershed function. The tricky part of using vector(s) is that the point(s) MUST fall on top of one of the rasters where water accumulates (despite the fact that Watershed doesn\u2019t use Water Accumulation). This condition is not always met just like you can see the difference between Water Accumulation and the stream layer, so the positions of the Gauges had to be adjusted to match the underlying raster for the purpose of this exercise.<\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\" \/>\n\n\n\n<h3 class=\"wp-block-heading\">Task 5 &#8211; Calculate area of watersheds<\/h3>\n\n\n\n<ol start=\"14\" class=\"wp-block-list\">\n<li>Open the attribute table of the watershed raster created and notice the count column in each watershed, this is the number of raster cells in each one.&nbsp; You multiply this number by the area of one cell (get the properties of the layer and find cell size and map units ).<\/li>\n<\/ol>\n\n\n\n<p class=\"has-light-green-cyan-background-color has-background wp-block-paragraph\"><em>Question 1. What is the area of one cell in meters square (m<sup>2<\/sup>)? _____<\/em><\/p>\n\n\n\n<p class=\"has-light-green-cyan-background-color has-background wp-block-paragraph\"><em>Question 2. What is the area of the Dean Brook watershed in square miles? _______<\/em><\/p>\n\n\n\n<p class=\"has-light-green-cyan-background-color has-background wp-block-paragraph\"><em>Question 3. What is the area of the Roaring Brook watershed in square miles? ________<\/em><\/p>\n\n\n\n<p class=\"has-light-green-cyan-background-color has-background wp-block-paragraph\"><em>Question 4. What is the area of the West branch of the Swift River watershed in square miles? ______<\/em><\/p>\n\n\n\n<ol start=\"15\" class=\"wp-block-list\">\n<li>Save the project for next part.<\/li>\n<\/ol>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity is-style-default\" \/>\n\n\n\n<p class=\"wp-block-paragraph\"><a href=\"#top\">Return to Top of page<\/a><\/p>\n\n\n\n<p class=\"wp-block-paragraph\"><a href=\"https:\/\/commons.mtholyoke.edu\/emarcano\/geog323\">Return to Geog 323 page<\/a><\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\" \/>\n\n\n\n<h2 id=\"part2\" class=\"wp-block-heading\">Part 2 \u2013 Stream Order Analysis<\/h2>\n\n\n\n<blockquote class=\"wp-block-quote is-layout-flow wp-block-quote-is-layout-flow\">\n<p class=\"wp-block-paragraph\">Stream ordering is a method of assigning a numeric order to links in a stream network. This order is a method for identifying and classifying types of streams based on their numbers of tributaries. Some characteristics of streams can be inferred by simply knowing their order.<br><br>For example, first-order streams are dominated by overland flow of water; they have no upstream concentrated flow. Because of this, they are most susceptible to non-point source pollution problems and can derive more benefit from wide riparian buffers than other areas of the watershed.<br><br>The Stream Order tool in ArcGIS (under Spatial Analyst\\Hydrology) has two methods you can use to assign orders. These are the methods proposed by Strahler (1957) and Shreve (1966).<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">&#8211; <a href=\"https:\/\/doc.esri.com\/en\/arcgis-pro\/latest\/tool-reference\/spatial-analyst\/how-stream-order-works.html\" target=\"_blank\" rel=\"noopener\">ArcGIS Pro 3.7 help<\/a>)<\/p>\n<\/blockquote>\n\n\n\n<figure class=\"wp-block-image size-full\"><a href=\"https:\/\/commons.mtholyoke.edu\/emarcano\/wp-content\/uploads\/sites\/965\/2026\/09\/stream-orders.jpg\"><img loading=\"lazy\" decoding=\"async\" width=\"712\" height=\"262\" src=\"https:\/\/commons.mtholyoke.edu\/emarcano\/wp-content\/uploads\/sites\/965\/2026\/09\/stream-orders.jpg\" alt=\"\" class=\"wp-image-2206\" srcset=\"https:\/\/commons.mtholyoke.edu\/emarcano\/wp-content\/uploads\/sites\/965\/2026\/09\/stream-orders.jpg 712w, https:\/\/commons.mtholyoke.edu\/emarcano\/wp-content\/uploads\/sites\/965\/2026\/09\/stream-orders-300x110.jpg 300w\" sizes=\"auto, (max-width: 712px) 100vw, 712px\" \/><\/a><\/figure>\n\n\n\n<div style=\"height:8px\" aria-hidden=\"true\" class=\"wp-block-spacer\"><\/div>\n\n\n\n<blockquote class=\"wp-block-quote is-layout-flow wp-block-quote-is-layout-flow\">\n<p class=\"wp-block-paragraph\">In both methods, the upstream stream segments, or exterior links, are always assigned an order of 1.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\"><br><strong>Strahler method<\/strong><br>In the Strahler method, all streams without any tributaries are assigned an order of 1 and are referred to as first order. The stream order increases when streams of the same order intersect. Therefore, the intersection of two first-order streams will create a second-order stream, the intersection of two second-order streams will create a third-order stream, and so on. The intersection of two streams of different orders, however, will not result in an increase in order. For example, the intersection of a first-order and second-order stream will not create a third-order stream but will retain the order of the highest ordered stream.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\"><strong>Shreve method<\/strong><br>As with the Strahler method, all exterior streams are assigned an order of 1. For all interior streams in the Shreve method, however, the orders are additive. For example, the intersection of two first-order streams creates a second-order stream, the intersection of a first-order and second-order stream creates a third-order stream, and the intersection of a second-order and third-order stream creates a fifth-order stream.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">&#8211; <a href=\"https:\/\/doc.esri.com\/en\/arcgis-pro\/latest\/tool-reference\/spatial-analyst\/how-stream-order-works.html\" target=\"_blank\" rel=\"noopener\">ArcGIS Pro 3.7 help<\/a>)<\/p>\n<\/blockquote>\n\n\n\n<h3 class=\"wp-block-heading\">Task 1 &#8211; Stream ordering<\/h3>\n\n\n\n<p class=\"wp-block-paragraph\">You have already most necessary data to compute stream order in our AOI (Flow direction and flow accumulation, both from <a href=\"#part1\">Part 1- watershed delineation<\/a>)<\/p>\n\n\n\n<ol class=\"wp-block-list\">\n<li>Open the project for Part 1.<\/li>\n<\/ol>\n\n\n\n<p class=\"wp-block-paragraph\">Before, we told you that the flow accumulation raster can work as a proxy for a river\/stream layer.&nbsp; So, let\u2019s prepare it by separating the cells where water accumulates from the cells with no accumulation. We\u2019ll use a value of accumulation of 100 to start separating the streams (the powers that be decided that an accumulation of less than 100 is not significant).<\/p>\n\n\n\n<ol start=\"2\" class=\"wp-block-list\">\n<li>Launch to the Reclassify tool (<strong>Spatial Analyst Tools &gt; Reclass &gt;<\/strong> <strong>Reclassify<\/strong>)<\/li>\n<\/ol>\n\n\n\n<ol style=\"list-style-type:lower-alpha\" class=\"wp-block-list\">\n<li>Use <em>flow_acc<\/em> as input<\/li>\n\n\n\n<li>Use Value as Reclass fiedl<\/li>\n\n\n\n<li>Click the Classify button and choose 2 classes<\/li>\n\n\n\n<li>In the upper values for the first class type 100  and leave the other upper value as it is (it&#8217;s the highest accumulation). Click OK<\/li>\n\n\n\n<li>Back at the first Reclassify dialog, enter the 0 as New values for the 0-100 range, and 1 for the second range, leave NoData untouched (you have to hit Enter after you change the numbers)<\/li>\n\n\n\n<li>Type <em>Streams_raster<\/em> as output and hit Run<\/li>\n<\/ol>\n\n\n\n<p class=\"wp-block-paragraph\">Now you should have a raster with <strong>only two values<\/strong>, 0 &amp; 1, showing where the water accumulates.&nbsp; Let\u2019s classify the streams next.<\/p>\n\n\n\n<ol start=\"3\" class=\"wp-block-list\">\n<li>Launch the Stream Order tool (<strong>Spatial Analyst &gt; Hydrology &gt; Stream Order<\/strong>)<\/li>\n<\/ol>\n\n\n\n<ol style=\"list-style-type:lower-alpha\" class=\"wp-block-list\">\n<li>Use <em>Streams_raster<\/em> as input<\/li>\n\n\n\n<li><em>flow_dir<\/em> as, you guessed it, the flow direction raster<\/li>\n\n\n\n<li>Call the output <em>Strahler_streams<\/em><\/li>\n\n\n\n<li>Leave STRAHLER as method of stream ordering.&nbsp; Click OK.<\/li>\n<\/ol>\n\n\n\n<p class=\"wp-block-paragraph\">Enjoy your multi color result.&nbsp; Compare the new stream layer with the <em>HYDRO25K_ARC_AOI<\/em> layer you have in the CP. Put the vector layer on top of the raster, in the CP, and turn on and off.&nbsp; Also add the <em>DEM_AOI<\/em>\/<em>Hillshade_AOI<\/em> with a suitable elevation symbology under those layers. Zoom in several places and pan around the scene to notice any differences. You should be able to see many artifacts resulting from the processing, a result of the shape of the surface (DEM), software misinterpretations.&nbsp; Also, many smaller channels were not digitized in the vector stream layer. The larger channels (higher flow accumulation) should be better than the small ones.<\/p>\n\n\n\n<ol start=\"4\" class=\"wp-block-list\">\n<li>Now run the Stream Order procedure using the Shreve method and compare those two results. Can you explain the difference in results?<\/li>\n<\/ol>\n\n\n\n<p class=\"has-background wp-block-paragraph\" style=\"background-color:#7bdcb5\"><em>Question 5.<\/em> <em>How many stream order values are in Strahler method? &nbsp;______<\/em><\/p>\n\n\n\n<p class=\"has-background wp-block-paragraph\" style=\"background-color:#7bdcb5\"><em>Question 6.<\/em> <em>How many stream order values are in Shreve method? &nbsp;______<\/em><\/p>\n\n\n\n<ol start=\"5\" class=\"wp-block-list\">\n<li>Open Attribute Table of <em>Strahler_streams<\/em>.<\/li>\n<\/ol>\n\n\n\n<p class=\"wp-block-paragraph\">In the attribute table you\u2019ll only see the classes and how many pixels there are in each stream classes (remember that this is still a raster layer).&nbsp; It doesn\u2019t tell us the number of individual streams (nor the length) in each class.&nbsp; For this we need to convert these to vectors.<\/p>\n\n\n\n<ol start=\"6\" class=\"wp-block-list\">\n<li>In the top tool box open <strong>Conversion Tools &gt; From Raster &gt; Raster to Polyline.<\/strong><\/li>\n<\/ol>\n\n\n\n<ol style=\"list-style-type:lower-alpha\" class=\"wp-block-list\">\n<li>As input raster enter <em>Strahler_streams<\/em><\/li>\n\n\n\n<li>Field should say \u201cValue\u201d<\/li>\n\n\n\n<li>Name the output <em>Strahler_streams_vector<\/em><\/li>\n\n\n\n<li>Make sure Simplify polylines box is ticked<\/li>\n\n\n\n<li>Background value Zero<\/li>\n<\/ol>\n\n\n\n<ol start=\"7\" class=\"wp-block-list\">\n<li>Compare the result, like you did before, with the <em>HYDRO25K_ARC_AOI<\/em> layer.<\/li>\n<\/ol>\n\n\n\n<p class=\"wp-block-paragraph\">You\u2019ll notice that the resulting layer is not as neat as the Hydro_arc layer.&nbsp; That\u2019s expected, your result comes from a raster-vector conversion and needs lots of refinement and editing, but for now it\u2019ll do just fine.&nbsp; Good job!<\/p>\n\n\n\n<ol start=\"8\" class=\"wp-block-list\">\n<li>Right click on <em>Strahler_streams_vector<\/em> layer in the CP and open the attribute table.&nbsp;<\/li>\n<\/ol>\n\n\n\n<p class=\"wp-block-paragraph\">Notice that there is now a column called <strong>\u201cgrid_code\u201d <\/strong>for each entry in the table (each line segment in the map).&nbsp; This is, as you might have guessed, the stream order. Please make a map that shows the streams with different symbols (you should know how to do this without instructions)<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Notice also that there is a column called <strong>Shape_Length<\/strong>, again this gives us the length of each segment (in map units, or meters in this case).<\/p>\n\n\n\n<ol start=\"9\" class=\"wp-block-list\">\n<li>In the attribute table of the streams vector layer, right click on the heading of the column called \u201cgrid_code\u201d and choose Summarize. This creates a table of summary statistics for each stream order.&nbsp; Right now, we only want the total length per stream order.<\/li>\n<\/ol>\n\n\n\n<ol style=\"list-style-type:lower-alpha\" class=\"wp-block-list\">\n<li>In the next dialog, under Statistics Fields select <strong>Shape_Length <\/strong>and for Statistics Type select <strong>Sum<\/strong>.<\/li>\n\n\n\n<li>Enter <em>Stream_order_summary<\/em> as output table. Click OK.<\/li>\n<\/ol>\n\n\n\n<p class=\"wp-block-paragraph\">At the bottom of the content pane (CP) there\u2019s a table, open it and observe the total length of each stream class.<\/p>\n\n\n\n<p class=\"has-background wp-block-paragraph\" style=\"background-color:#7bdcb5\"><em>Question 7.<\/em> <em>What is the total length of streams order 1 in miles? &nbsp;______<\/em><\/p>\n\n\n\n<p class=\"has-background wp-block-paragraph\" style=\"background-color:#7bdcb5\"><em>Question 8.<\/em> <em>What is the total length of streams order 6 in miles? ______<\/em><\/p>\n\n\n\n<p class=\"has-background wp-block-paragraph\" style=\"background-color:#7bdcb5\"><em>Question 9 (Challenge) What is the total length in miles of streams order 2 in the Dean Brook watershed? ____<\/em><br><br>Hint: You might want to extract that watershed and maybe vectorize it \ud83d\ude09<\/p>\n\n\n\n<p class=\"has-background wp-block-paragraph\" style=\"background-color:#7bdcb5\"><em>Question 10. What is the highest order of stream in Dean Brook watershed? <\/em>_______<\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\" \/>\n\n\n\n<h2 id=\"part3\" class=\"wp-block-heading\"><strong>Part 3 \u2013 Viewshed<\/strong><\/h2>\n\n\n\n<p class=\"wp-block-paragraph\">Now that we have delineated the extent of the Dean Brook watershed, the town of Shutesbury, MA, where most of this watershed is located, wants to establish an observation tower on top of nearby hill to monitor for forest fires.  They have selected the top of &#8220;January Hill&#8221;, northwest of Dean Brook watershed, where they&#8217;ll build a 10 meter high observation tower.<\/p>\n\n\n\n<ol class=\"wp-block-list\">\n<li>Add the point layer called <em>January_Hill<\/em><\/li>\n\n\n\n<li>Calculate the viewshed (i.e. visible areas) from this point (<strong>Spatial Analyst &gt; Surface &gt; Geodetic Viewshed)<\/strong><\/li>\n<\/ol>\n\n\n\n<ol style=\"list-style-type:lower-alpha\" class=\"wp-block-list\">\n<li>Input Raster <em>DEM_AOI<\/em><\/li>\n\n\n\n<li>Input point observer features <em>January_Hill<\/em><\/li>\n\n\n\n<li>Open Observer parameters and enter 10 <strong>meters <\/strong>for Observer offset (NOT &#8220;observer elevation&#8221;; don&#8217;t forget the unit!) and click Run<\/li>\n<\/ol>\n\n\n\n<p class=\"has-background wp-block-paragraph\" style=\"background-color:#7bdcb5\"><em>Question 11. (Challenge) What proportion, in percent, of the Dean Brook watershed is visible from top of the observation tower? <\/em>_____%<br><br>Hint: you know the area of the watershed (see question 2). You need to clip the viewshed raster to this watershed and calculate the area of visible areas.  Keep consistent units.<\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\" \/>\n\n\n\n<p class=\"wp-block-paragraph\"><a href=\"#top\">Return to Top of page<\/a><\/p>\n\n\n\n<p class=\"wp-block-paragraph\"><a href=\"https:\/\/commons.mtholyoke.edu\/emarcano\/geog323\">Return to Geog 323 page<\/a><\/p>\n\n\n\n<p class=\"wp-block-paragraph\"><\/p>\n","protected":false},"excerpt":{"rendered":"<p>Terrain analysis Contents Return to Geog 323 page Part 1: Watershed analysis Adapted from Spatialabs 2013 A watershed or drainage basin describes the area of land from which all surface waters (rain, melting snow or ice) converge at a single point. This point is the exit of the watershed and [&hellip;]<\/p>\n","protected":false},"author":793,"featured_media":0,"parent":2130,"menu_order":0,"comment_status":"closed","ping_status":"closed","template":"","meta":{"_monsterinsights_skip_tracking":false,"footnotes":""},"class_list":["post-2171","page","type-page","status-publish","hentry"],"_links":{"self":[{"href":"https:\/\/commons.mtholyoke.edu\/emarcano\/wp-json\/wp\/v2\/pages\/2171","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/commons.mtholyoke.edu\/emarcano\/wp-json\/wp\/v2\/pages"}],"about":[{"href":"https:\/\/commons.mtholyoke.edu\/emarcano\/wp-json\/wp\/v2\/types\/page"}],"author":[{"embeddable":true,"href":"https:\/\/commons.mtholyoke.edu\/emarcano\/wp-json\/wp\/v2\/users\/793"}],"replies":[{"embeddable":true,"href":"https:\/\/commons.mtholyoke.edu\/emarcano\/wp-json\/wp\/v2\/comments?post=2171"}],"version-history":[{"count":62,"href":"https:\/\/commons.mtholyoke.edu\/emarcano\/wp-json\/wp\/v2\/pages\/2171\/revisions"}],"predecessor-version":[{"id":2311,"href":"https:\/\/commons.mtholyoke.edu\/emarcano\/wp-json\/wp\/v2\/pages\/2171\/revisions\/2311"}],"up":[{"embeddable":true,"href":"https:\/\/commons.mtholyoke.edu\/emarcano\/wp-json\/wp\/v2\/pages\/2130"}],"wp:attachment":[{"href":"https:\/\/commons.mtholyoke.edu\/emarcano\/wp-json\/wp\/v2\/media?parent=2171"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}