...

Finally, you will calculate flow accumulation in each cell from the cells upstream, as illustrated below. Flow accumulation is dependent on the flow direction you calculated previously.

1. Return to the Geoprocessing pane.
2. At the top left of the Geoprocessing pane, click the Back button.
3. In the search box, type "flow".
4. Click the Flow Accumulation tool.
5. For ‘Input flow direction raster’, select the FDR layer.
6. For ‘Output accumulation raster’, rename the raster from “FlowAcc_FDR1” to “FAC”.
7. Ensure your Geoprocessing pane appears as shown below and click Run.
   Image Removed
   
   Image Added
   
   
8. In the Contents pane, uncheck Turn off the Flowlines_StatePlane layer to better see the FAC layer.
9. Right-click the FAC layer and click the select Symbology.
10. For ‘Symbology’, use drop-down menu to select‘Primary symbology’, select Classify.
11. In the ‘Method’ drop-down menu to For ‘Method’, select Equal Interval.
12. Use the ‘Classes’ drop-down menu to select8 classes.
13. For ‘Classes’, select 8.
14. In the 'Classes' tableFor ‘Class Breaks’, double-click the numbers in the Upper Value column  and column and typethe following list: 100, 300, 1000, 3000, 10000, 30000, 100000. (Leave the final max value, 351408 as is.)
15. Use the ‘Color Ramp’ drop-down menu to select the multipart color schemeFor ‘Color scheme’, scroll down to the very bottom and select the Yellow-Green-Blue (Continuous) color scheme, as shown below.
    
    Image Modified

Image Modified

Visually trace along the streams until you find the darkest path exiting the watershed. Notice that, according to the model you have just generated, the Buffalo-San Jacinto basin drains north into Lake Houston, instead of south southeast into Galveston Bay. Look for the exact point on the map where this divergence is created. It is circled in the image below. A problem is that the topography in that area is extremely flat and, rather than a single well-defined channel, there are wide bodies of water connecting the bayous to Galveston Bay. Calculating the flow accumulation more accurately would require further editing the base raster so that the wide flat bodies of water are not treated as narrow flowlines. This example also illustrates how you must always double-check the results of your models to ensure they make sense. For the purposes of learning the remaining geoprocessing tools in this lab, you will continue with the existing flow accumulation model, as if it was correct.

Image Removed

Zoom into the northern outlet point, as indicated on the image above, so that you can see the individual pixels.

Image Removed

1. Zoom in to the southeastern outlet point, as indicated below, so that you can see the individual pixels.
   
   

Image Added

Creating and editing feature classes

You will now create a new point feature class to store the outlet point you just identified.

1. Click At the bottom of the Symbology pane, click the Catalog tab.
2. Right-click the HydrologyLab  the HydrologyLab.gdb geodatabase and select New > Feature Class….
3. For ‘Feature Class Name’‘Name’, type “Outlet”.
4. Use the ‘Geometry Type’ drop-down menu to For ‘Feature Class Type’, select Point.
5. For Coordinate System, select Current Map.
6. Click Run.
7. At the bottom, click Next.
8. Click Next again.
9. For 'Current XY', ensure NAD 1983 StatePlane Texas S Central FIPS 4204 (US Feet) is selected.
10. Click Finish.
11. In the Catalog pane, drag the newly created Outlet feature class into the Lab3 map view.
12. On the ribbonOn the Standard toolbar, click the Edit tab and click the Create button.
    Image Removed
    Image Added
    
    
13. On the Create Features pane, select the Outlet layer and click Point. Click -> button.
    
    
    
    
14. In the Map Displaymap view, click near the center of the outlet pixel to draw a point there.
    
    
    
    
15. On the Editor toolbarribbon, click the Save button.
    
    Image Added
    
    
16. When asked if you want to save your all edits, click Yes.
17. Close the Editor toolbar and the Create Features tab.

Delineating watersheds

Now you will delineate the watershed that flows to the outlet point.

Image Removed

the outlet point.

Image Added

1. At the bottom of the Create Features pane, click the Geoprocessing tab.
2. At the top left of the Geoprocessing pane, click the Back button.
3. In the search box, type "watershed".
4. Click the Watershed Watershed (Spatial Analyst Tools) tool.
5. For ‘Input D8 flow direction raster’, select the FDR layer raster.
6. For ‘Input raster or feature pour point data’ select the Outlet layer.
7. For ‘Output raster’, rename the raster from “Watersh_FDR1” to “Watershed”.
8. Ensure your Geoprocessing pane appears as shown below and click Run.
   Image Removed
   Image Added
   
   
9. In the Contents pane, right-click the Watershed layer and select Zoom To Layer.

You can toggle the Subbasin_StatePlane layer on and off. Notice that the watershed boundary you have created Image RemovedNotice that this is very close, but does not exactly match the boundaries of the subbasin you were originally provided from NHDPlus.

Image Added

1. In the Contents pane, turn off all layers beneath the Watershed layer, so that only the cells within the newly defined subbasin are visible.

...

1. At the top left of the Geoprocessing pane, click the Back button.
2. In the search box, type "calculator".
3. Click the Raster Calculator Calculator (Spatial Analyst Tools) tool.
4. Using the same clicking technique as before, enter the following equation: ("FAC" > 3000) & ("Watershed" > 0).
5. For ‘Output raster’, rename the raster from “fac_rasterca” to “Streams”.
6. Ensure your Geoprocessing pane appears as shown below and click Run.
   Image Removed
   Image Added
   

Again you should have a binary raster, where flowline cells are assigned a value of 1 and all other cells are assigned a value of 0.

...

1. At the top left of the Geoprocessing pane, click the Back button.
2. In the search box, type "stream".
3. Click the Stream Link tool.
4. For ‘Input stream raster’, select the Streams layer raster.
5. For ‘Input flow direction raster’, select the FDR layer raster.
6. For ‘Output raster’, rename “StreamL_Stre1” to “StreamLinks”.
7. Ensure your Geoprocessing pane appears as shown below and click Run.
   Image Removed
   Image Added
   
   
8. Symbolize the StreamLinks layer using unique values.
9. In the Contents pane, collapse the StreamLinks symbology.

...

Now you will delineate the catchment basin associated with each stream link. These catchment basins are similar to the subwatersheds you downloaded from NHDPlus.

1. Return to the Geoprocessing pane.
2. At the top left of the Geoprocessing pane, click the Back button.
3. In the search box, type "watershed".
4. Click the Watershed Watershed (Spatial Analyst Tools) tool.
5. For ‘Input D8 flow direction raster’, select the FDR layer raster.
6. For ‘Input raster or feature pour point data’ select the StreamLinks layer.
7. For ‘Output raster’, rename “Watersh_FDR2” to “CatchmentsCatchmentsRaster”.
8. Ensure your Geoprocessing pane appears as shown below and click Run.
   Image Removed
   Image Added
   
   
9. Symbolize the Catchments CatchmentsRaster layer using unique values.Image Removed
   
10. In the Contents pane, collapse the Catchments CatchmentsRaster symbology.

Converting rasters to features

Now you will convert the Catchments CatchmentsRaster raster into a polygon feature class.

1. Return to the Geoprocessing pane.
2. At the top left of the Geoprocessing pane, click the Back button.
3. In the search box, type "raster to polygon".
4. Click the Raster to Polygon tool.
5. For ‘Input raster’, select the Catchments layer CatchmentsRaster layer.
6. For ‘Output polyline features’, rename “RatsterT_Catchme1” to “CatchmentPolyCatchmentsPoly”.
7. Uncheck Simplify polylinespolygons.
8. Ensure your Geoprocessing pane appears as shown below and click Run.
   Image Removed
   Image Added
   
   
9. Symbolize the CatchmentPoly layer using unique values with ID as the value field. Make sure to click the Add all values button.
   Image Removed
   Image Added
   
   
10. In the Contents pane, collapse the CatchmentPoly symbology.

Now you will also convert the raster StreamLinks into a polyline feature class. While you could use the Raster to Polyline tool, you will instead use a tool in the Hydrology toolset, designed specifically for converting raster stream links into polyline features. The difference between the two methods is illustrated below.

...

Image Added

1. Return to the Geoprocessing pane.
2. At the top left of the Geoprocessing pane, click the Back button.
3. In the search box, type "stream".
4. Click the Stream to Feature tool.
5. For ‘Input stream raster’, select the StreamLinks layer raster.
6. For ‘Input flow direction raster’, select the FDR layer raster.
7. For ‘Output polyline features’, rename the feature class from “StreamT_StreamL1” to “DrainageLines”.
8. Uncheck Simplify polylines.
9. Ensure your Geoprocessing pane appears as shown below and click Run.
   Image Removed
   Image Added
   
   

Image RemovedImage Added

Calculating Strahler stream order

...

1. At the top left of the Geoprocessing pane, click the Back button.
2. In the search box, type "stream".
3. Click the Stream Order tool.
4. For ‘Input stream raster’, select the Streams layer raster.
5. For ‘Input flow direction raster’, select the FDR layer raster.
6. For ‘Output raster’, rename the raster from “StreamO_Stre1” to “StrahlerOrder”.
7. Ensure your Geoprocessing pane appears as shown below and click Run.
   Image Removed
   Image Added
   
   
8. In the Contents pane, turn off all layers except for the new StrahlerOrder layer, so that its symbology is more apparent.
   
   

...

1. At the top left of the Geoprocessing pane, click the Back button.
2. In the search box, type "zonal statistics".
3. Click the Zonal Statistics as Table toolTable (Spatial Analyst) tool.
4. For ‘Input raster or feature zone data’, select the StreamLinks layer from the drop-down menu.
5. Use the For ‘Zone field’ drop-down menu to , select the Value field.
6. For ‘Input value raster’, select the StrahlerOrder layer from the drop-down menu raster.
7. For ‘Output table’, rename the table from “ZonalSt_StreamL1” to “OrderTable”.
8. Ensure your Geoprocessing pane appears as shown below and click Run.
   Image Removed
   Image Added

Adding a new field to an attribute table

...

1. Open the DrainageLines attribute table.
2. At the top of the DrainageLines attribute table, click the Add Field button.
   Image Removed
   
   Image Added
   
   
3. In the bottom row of the table, for ‘ Field NameFor ‘Name:’, type “StrahlerOrder”.
4. Use the ‘Type‘Data Type:’ drop-down menu to select Short Integer for short integer.
5. Ensure your ‘Add Field’ window Fields table view appears as shown below and click Save on top of  the main tab.
   Image Removed
   Image Removed.
   
   Image Added
   
   
6. On the ribbon, under the Fields tab, click the Save button.
   
   Image Added
   
   
7. Close the Fields table view.

Notice Scroll all the way to the right and notice the new StrahlerOrder field has been added to the end of the DrainageLines attribute table. All of the values are currently null, but will later be populated with the Strahler Order integers from the OrderTable.

...

You are now ready to join the OrderTable to the DrainageLines layer based on their common identification fields.

1. Click At the top right of the DrainageLines attribute table, click the Table Options button and select Joins and Relates > Add Join….
   Image Removed
   Image Added
   
   
2. For ‘Input Join Field’, select the grid_code field.
3. For ‘Join Table’, select the OrderTable table.
4. For ‘Output Join Field’, select the VALUE Value field.
5. Ensure that your Geoprocessing pane Add Join window appears as shown below and clickRun OK.
   Image Removed
   Image Added
   
   

Scroll across the table and notice that most of the statistics (Min, Max, Mean, Majority, Minority, Median) are the same and all contain the Strahler order integer, so it doesn’t matter which field you copy, though you will choose the MIN field for this lab.

1. In the attribute table, right-click on the StrahlerOrder field header and select Calculate Field.
2. Scroll down the list of ‘Fields:’ ‘Fields’ and double-click the OrderTable.MINORITY field MIN field and click RunOK.
3. Scroll down the newly populated StrahlerOrder field and notice the values of 1, 2, and 3 corresponding to Strahler order value. You

Since you have copied the Strahler order into a field in the original DrainageLines attribute table, you no longer need the tabular join.

1. Again, click the Table Options button and select Joins and Relates > Remove Join(s) > > Remove All Joins.
2. When asked if you are sure you want to remove all joins, click Yes.
3. Symbolize the DrainageLines layer using graduated symbols based on the StrahlerOrder field.
4. In the Contents pane, turn on only the Outlet, DrainageLines, and CatchmentPoly, layers and turn offall other layers.
   Image Removed
   Image Added

FOR MAP LAYOUT TO BE TURNED IN

Create an 8.5 x 11 layout showing the vector catchment basins along with the stream links symbolized using graduated symbols according to their Strahler order and the outlet point that you used.

1. An exported map of Create an 8.5 x 11 layout showing the DEMRecon layer turned on.
2. Create an 8.5 x 11 layout showing the vector catchment basins along with the stream links symbolized using graduated symbols according to their Strahler order and the outlet point that you used.