Exercise 2:

Basic Terrain Analysis and Flow Path Mapping

Due Date: Feb. 10

Submit your report file in your personal working directory in the GEOG591 course directory. If your ONYEN ID is smith, name the report file as ex1_smith_report.doc. In your report, put your histograms and maps which are produced by each process.

Datasets:

  • Raster Dataset

1)btl_20ft_dem - 20ft Interpolated LIDAR Elevation Data of Battle Creek, Chapel Hill

  • Vector Dataset

1)btl_middle_subwatershed – Middle Battle Creek subwatershed

2) btl_upper_subwatershed – Upper Battle Creek subwatershed

3) ch_imperviousness_surface– impervious surface map of Chapel Hill

4) ch_hydroarc – stream network in the area

All datasets of this exercise are in \\isis.unc.edu\html\courses\2010spring\geog\591\001\exercises_2012\ex2.

Preparation

  1. Create a directory named ex2 under your personal working directory, and copy all files and data in exercises/ex2 to that directory.
  1. Open ArcMap.
  1. Addbtl_20ft_demraster data on your ArcMap.
  1. Overlay all vector data.

Basic Terrain Analysis

  1. Compute slope and aspect from btl_20ft_demraster dataset.
  1. Make histograms for btl_20ft_dem and aspect map in Middle and Upper Battle Creek subwatersheds by selecting each subwatershed and using the zonal histogram function in Spatial Analyst.

-Paste all four histograms to your reports.

Q1. Differentiate between terrain characteristics in each of the subwatershed areas based on the DEM and aspect distributions. How are DEM and aspect distributed in each of the two subwatersheds?

Flow Direction Map

  1. Make a D8 flow direction map

-Go to ArcToolbox Spatial Analyst Tools  Hydrology  Flow Direction

-Name your results as “FlowDir_btl”

Flow Accumulation Map

  1. Make a D8 flow specific catchment area map

-Go to ArcToolbox Spatial Analyst Tools  Hydrology  Flow Accumulation

-Name your results as “FlowAcc_btl”

Stream Network Map with threshold 1000

  1. Derive a finite stream network with an accumulation threshold in raster calculator

-Go to Spatial Analyst Map Algebra Raster Calculator

-Double click “FlowAcc_btl”, “ >”, and type 1000.

-Hit “Evaluate”. You can get binary stream network map, which is temporary file with name “Calculation”

-Make this file permanent; right click “Calculation” file  Data  Make Permanent  Name it as “stream_1000”

Filling Pits

  1. Fill pits

-Go to ArcToolbox Spatial Analyst Tools  Hydrology  Fill

-Name your result as “fill_dem_20ft.”

  1. Redo all processes from 7 to 9 with the pit-filled DEM (fill_dem_20ft).

Note: Make sure that you are using different raster output names for this process (e.g. fill_FlowDir).

  1. Repeat the identification of stream channels with a threshold of 500

Q2. Compare the stream network and flow direction maps between the original and filled DEM. What is the major difference between them?

  1. Use the Map Algebra to detect locations of all pits by subtracting the original and pit removed DEM.

-Make the temporary result file permanent (like step 9) and name it “pit_loc.”

Q3. Comment on the spatial distribution of depression cells. Where can you find major pits?How does filling pits affect the flow paths of this area? Why?

Note: Overlay vector files given in this exercises to locate the depression cells.

Q4. How many hectares is the 1000 and the 500 grid cell threshold that defines the stream heads?

Q5. What range of drainage area (in hectares) appears to define a stream head in this area?

Feel free to ask any question to Yuri Kim (.)