Labs ERODE: Difference between revisions

Latest revision as of 16:33, 12 May 2016

Hydrology and Flow Routing

Classroom organization
If you have never used the Web Modeling Tool, learn how to use it here. The WMT allows you to set up simulations, but once you are ready to run them, you will need an account on the CSDMS supercomputer to submit your job. More information on getting an account can be found here HPCC Access. Note that this takes a few days!
This lab introduces you to raster-based flow routing; this basic concept of flow routing is a key process in many hydrological and sediment transport models. Whereas flow routing can be complicated, we here present a basic algorithm. This lab takes ~1 hr to complete.

Learning objectives

get familiar with the concept of flow routing in a digital elevation model/ grid
learn about the D8 method, contributing area and flow length in the drainage direction

Lab Notes
An important part of any landscape evolution model is the way it treats routing of water (and sediment) through the landscape. To give you a bit of a feel of the methodology, we will look at the ‘D8 method’ for determining the drainage directions in a raster grid. Assume you have a 3 by 3 grid. So if one has a source cell of interest in the middle, we assume that water and sediment drains in the direction of the steepest slope. This concept is also called ‘steepest descent’. The grid cell in the middle then has 8 neighboring cells, and depending on their elevations (z), and the distance between the cells, d, you can calculate the slope, S.

[math]\displaystyle{ S = \Delta z / d }[/math]

(1)

These surrounding cells are coded as shown in the figure here:

Question
In the hypothetical elevation grid below, the source cell and its 8 neighbors have defined elevation values.
Where does the water drain to? Select the correct cell and list its D8- grid value.

Question
Can you think of a theoretical situation when a D8 flow direction algorithm would not work? 
What would such a situation look like in a real DEM?

>> Open a new browser window and open the Web Modeling Tool WMT with hydrology components

>> We have two DEM datasets loaded into the WMT dataset for the D8 component. There is an example called 'HydrologyLab_D8_MountSopris' and one called 'HydrologyLabs_ D8Global_WindRiverRange'

>> Run the D8 component for either example and plot the contributing area grid (A) and the flow length (ds) resulting from the calculations.

@@ Line 1: / Line 1: @@
-==Landscape Evolution Modeling with ERODE==
+==Hydrology and Flow Routing==
 <br>
-This lab introduces you to the landscape evolution model ERODE. ERODE is a raster-based landscape evolution model; it evolves a landscape with a combination of hilslope erosion and transport and river channel erosion and transport.<br>
+<br>
+'''Classroom organization'''<br>
+If you have never used the Web Modeling Tool, learn how to use it [[WMT_tutorial|here]]. The WMT allows you to set up simulations, but once you are ready to run them, you will need an account on the CSDMS supercomputer to submit your job.
+More information on getting an account can be found here [[HPCC_Access| HPCC Access]]. Note that this takes a few days!<br>
+This lab introduces you to raster-based flow routing; this basic concept of flow routing is a key process in many hydrological and sediment transport models.
+Whereas flow routing can be complicated, we here present a basic algorithm. This lab takes ~1 hr to complete.
+<!-- It is also used in watershed or catchment delineation in GIS applications.
+it evolves a landscape with a combination of hillslope erosion and transport and river channel erosion and transport. Download an introductionary [[:File:lecture_ERODEbasics.pptx|presentation]] of the landscape evolution model Erode. <br> -->
+'''Learning objectives '''
-'''STEP1	Load the ERODE project'''
+* get familiar with the concept of flow routing in a digital elevation model/ grid
+* learn about the D8 method, contributing area and flow length in the drainage direction
+<br>
 <br><br>
-Once you are in the CSDMS Modeling Tool: <br>
->> Under the File Menu, choose ‘Open Project’.<br>
->> Go to the Terrestrial Group. Open Project: ERODE<br>
->> Drag the ERODEGlobal component in the driver palette<br>
->> Specify a unique working directory for each experiment, for example: /CMT_Output/ERODEexp1   (…etc for the next experiment).<br>
 <br>
+'''Lab Notes'''<br>
+An important part of any landscape evolution model is the way it treats routing of water (and sediment) through the landscape. To give you a bit of a feel of the methodology, we will look at the ‘D8 method’ for determining the drainage directions in a raster grid.
+Assume you have a 3 by 3 grid. So if one has a source cell of interest in the middle, we assume that water and sediment drains in the direction of the steepest slope. This concept is also called ‘steepest descent’. The grid cell in the middle then has 8 neighboring cells, and depending on their elevations (z), and the distance between the cells, d, you can calculate the slope, S.
+::::{|
+|width=500px|<math>S = \Delta z / d  </math>
+|width=50px align="right"|(1)
+|}
-'''Learn about Flow Routing with ERODE'''
+These surrounding cells are coded as shown in the figure here: <br><br>
-<br><br>
-An important part of ERODE is the way it treats routing of water (and sediment) through the landscape. To give you a bit of a feel of the methodology, we will look at the ‘D8 method’ for determining the drainage directions in a raster grid.
-Assume you have a 3 by 3 grid. So if one has a source cell of interest in the middle, we assume that water and sediment drains in the direction of the steepest slope. This concept is also called ‘steepest descent’. The grid cell in the middle then has 8 neighboring cells, and depending on their elevations (z), and the distance between the cells, d, you can calculate the slope, S. These surrounding cells are coded as shown in the figure here: <br><br>
 [[File:codingD8.jpg|200px]]
 <br>
@@ Line 28: / Line 40: @@
   What would such a situation look like in a real DEM?
 <br>
+>> Open a new browser window and open the Web Modeling Tool [https://csdms.colorado.edu/wmt-hydrology  WMT with hydrology components]<br>
+>> We have two DEM datasets loaded into the WMT dataset for the D8 component.
+There is an example called 'HydrologyLab_D8_MountSopris' and one called 'HydrologyLabs_ D8Global_WindRiverRange'<br>
-'''STEP2	Run and analyze an ERODE simulation'''
+>> Run the D8 component for either example and plot the contributing area grid (A) and the flow length (ds) resulting from the calculations.
-<br><br>
-ERODE uses an adaptive time-stepping scheme. Landscapes evolve over long periods of time, and to optimize how long a simulation takes one would like to have large timesteps. However, if the elevation differences created by erosion and sedimentation become too large, the numerical scheme can cause instabilities. One approach is to use adaptive time-stepping methods, where time-steps are large when relatively little happens, and time-steps are short when the landscape evolves more rapidly.
->> Start two short simulations of ERODE.  Use stop-method code 0, and use 1000 time-steps as your simulation duration.  Use for the first simulation a grid of 75 by 75 cells and for the second simulation a grid of 50 by 50 cells. Ask for output grids every 100 years (use the z-grid for elevations and the A-grid for contributing Area).
+[[File:ZoomedWindRivers_DEM.png| 400px]]
+[[File:ContributingArea_windRiverRange.png| 600px]]