diff --git a/grass/hydrology.html b/grass/hydrology.html index 248ed07..f18c14b 100644 --- a/grass/hydrology.html +++ b/grass/hydrology.html @@ -1,10 +1,10 @@

Flow routing and watershed analysis

Resources:

@@ -42,24 +42,19 @@

Compute flow direction, flow accumulation and subwatersheds

-Extract more detailed streams from flow accumulation raster: +Extract more detailed streams from flow accumulation raster, +and convert basins and streams to vector format:


 r.mapcalc "streams_der_30m = if(abs(accum_10K) > 100, 1, null())"
-
- -

-Convert to vector format and display results -along with the official Wake county streams (red): - -


 r.to.vect -s basin_10K output=basin_10K type=area
 r.thin streams_der_30m output=streams_der_30m_t
 r.to.vect -s streams_der_30m_t out=streams_der_30m type=line
 

-Generate shaded map and display: +Generate relief shaded map of basins + and display extracted streams along with the official Wake county streams (red): How do the derived streams compare with the official stream map?
How can you modify the mapcalc expression to make stream origins fit better with the official stream map? @@ -70,14 +65,19 @@

Compute flow direction, flow accumulation and subwatersheds

d.vect basin_10K type=boundary d.rast lakes d.vect streams_der_30m color=blue -d.vect streams color=red d.out.file mystreams +d.vect streams color=red +d.out.file streams_compare +

Note that you can also use a MFD routing for r.watershed combined with r.stream.extract +for vectorized streams with topology of stream network. Read the manual pages to learn more. +

Create map of DEM depressions

-Depression filling is often necessary for certain flow routing algorithms but it can alter the elevation data significantly. +Depression filling is sometimes necessary for certain flow routing algorithms +but it can significantly alter the elevation data leading to srtificial straight channels. Find out how extensive the depressions are in our DEM.
-Note that r.watershed doesn't need any depression filling thanks to +Note that r.watershed tool doesn't need any depression filling thanks to its underlying algorithm which uses least cost path to get over depressions.

@@ -102,6 +102,10 @@ 

Create map of DEM depressions

d.out.file depressions
+What is the area mapped as depression in [ha] and [%] ? +(hint: use r.report with depr_bin)
+Optional: How would you compute the total volume of mapped depressions? +

Derive contributing area for a given outlet

Set region to the high resolution study area and zoom to it: