Difference between revisions of "Tutorials:8 Infiltration"

From Gsshawiki
Jump to: navigation, search
(Setting up the index map)
 
(No difference)

Latest revision as of 21:42, 2 May 2009

This tutorial is compatible with:

  • WMS Version 8.1
  • GSSHA Version 3.0b and later

Disclaimer: GSSHA tutorial exercises do not represent real world conditions

Infiltration is a key loss mechanism in a watershed; no watershed model is complete without it. GSSHA™ has four different infiltration models. During this tutorial you will set up the inputs needed for the Green & Ampt with Soil Moisture Redistribution model.

If you are starting at this tutorial,

  1. In the 2D Grid Module Icon 2DGrid.png select GSSHA™ | Open Project File.
  2. Browse to the Finished_tutorial\streams_embank folder.
  3. Select the streams_embank.prj file and select Open.

Index Map Setup

Importing the Shapefile as a GIS Layer

The first step that needs to be done is to set up an index map that describes the spatial variation in parameters needed by the infiltration model. We shall use a soil type shapefile from the SSURGO database of the area to create a soil type index map. Before we do anything with the shapefile in WMS, though, let’s look at it so that we know what we are working with.

  1. Launch Microsoft® Excel®.
  2. Select File | Open.
  3. Browse to the Soils folder of the GSSHA™ tutorial data.
  4. Using the control key, select both the soil_clip.dbf and IDs_for_WMS.dbf files. Select open.
  5. If soil_clip.dbf is not the active worksheet change to it by selecting Window | soil_clip.dbf.

  6. This is the attribute data that comes in the shape file. Notice that there are four parameters, ‘areasymbol’, ‘spatialver’, ‘musym’, and ‘mukey’. What WMS is looking for is a parameter that will serve as the soil type ID. There is not a suitable one in this file.


  7. Switch to the IDs_for_WMS.dbf worksheet by selecting Window | IDs_for_WMS.dbf.

  8. This dbf (dbase IV) file has been created to provide the information that WMS needs for the soil_clip shapefile to be useful. Usually when setting up a GSSHA™ model the soils data base file needs to be manipulated so that it has a number that indicates a soil classification. Notice that the previous symbols are there, although the number of lines of data is greatly reduced. Additionally four parameters have been created, ‘newid’, ‘classifica’, ‘erosion’, and ‘descriptio’. Looking at the ‘classifica’ and ‘descriptio’ parameters we can see that the soils were grouped according to soil type classification.

    Table 8. Soils type classifications.
    Classification Description
    1 Urban Land complex
    2 Silt Loam
    3 Silty clay Loam
    4 Fine sand
    5 Silty Clay
    6 Water

    First we want to bring in the soil_clip shapefile as a GIS layer. You can close Excel® and switch back to WMS.


  9. Select Edit | Current Coordinates...
  10. Select Set Projection.
  11. Ensure the coordinate system is set to UTM, NAD83, Meters, Zone 16.
  12. Select OK.
  13. Select OK.
  14. Right-click on GIS Layers in the data tree.
  15. Select Add Shapefile Data.

  16. Browse to the Soils folder and open the soil_clip.shp shapefile. This shapefile is added as a GIS layer; however it is not “officially” part of your project. It can now be converted into a WMS coverage in a much simpler manner, however before we convert it to a WMS coverage, we want to join the IDs_for_WMS.dbf file to it.


  17. Right-click on soil_clip.shp in the data tree.
  18. Select Join Table To Layer.
  19. Open the IDs_for_WMS.dbf file.
  20. Under Shapefile Join Field select mukey.
  21. Under table data field select classifica.
  22. Select OK.
  23. Right-click on soil_clip.shp in the data tree.
  24. Select Open Attribute Table.
  25. Make sure that the classifica field was added to the table and hit Ok.

  26. The next step is to create the soil type coverage in WMS to receive the polygons.


  27. Right-click on Coverages in the data tree.
  28. Select New Coverage…
  29. Change the Coverage Type to Soil Type.
  30. Hit OK.

  31. Now we can convert the shapefile to a coverage.


  32. Make sure the Soil Type coverage is the current active coverage by clicking on it in the Coverage list.
  33. Select the soil_clip.shp GIS layer by clicking on it. This will change the active module to the GIS model Icon GIS Module.png.
  34. In the GIS Module select, Mapping | Shapes -> Feature Objects.
  35. Select Yes for use all shapes in visible shapefiles.
  36. Select Next.
  37. Scroll over to the Classifica column.
  38. In the drop-down box that says ‘not mapped’ change it to be ‘SCS Soil type’
  39. Hit Next and then Finish.
  40. Wait for WMS to convert the shapefile into the coverage.

Cleaning up the Soil Type Coverage

Let’s visualize the soil type coverage.

  1. In the 2D Grid Module Icon 2DGrid.png select Display | Display Options.
  2. Turn off all of the 2D Grid options except the boundary.
  3. Click OK.
  4. In the Data Tree, uncheck the land use and GSSHA™ coverages.
  5. In the Data Tree, select the Soil Type coverage.
  6. Select Display | Display Options.
  7. Switch to the Map Data tab.
  8. Under the Polygons field, turn on the Color Fill Polygons.
  9. Select Soil Type Display Options.
  10. Set up the colors and patterns to make them more visible. Pick colors similar to the following image. Select OK when you are done.
  11. In the Points/Nodes area, turn off the Points/Nodes and the Vertices.
  12. In the Legends field, turn on the Soil Type legend.
  13. Select OK.

  14. Figure 31. Initial soil types in the watershed.

    Based on the earlier description of the soil types, soil type 1 is classified as a reworked residential area. Soil type 6 is classified as water. In order to assign soil properties to soil types 1 and 6 we will edit the polygons with these soil types and change them to have the the soil type of the neighboring polygons.

  15. Using the Pan, Zoom, and Select Polygon tools change the soil polygon IDs for Soil ID 1 (cyan) and Soil ID 6 (blue) to be what their neighbors are. This is accomplished by double-clicking on the desired polygon. This will bring up the Soil type mapping dialog as shown in the figure below. The soil type ID is changed in the WMS soil ID field. Click ‘Apply’
  16. Repeat step 14 for all the soil polygons of type 1 and 6.

  17. Figure 32. Changing a soils type mapping.

Setting up the index map

Now we can make an index map out of the soil type coverage. The index map will be used to describe the spatial variability of the infiltration parameters for the simulation.

  1. Switch to the 2D Grid Module Icon 2DGrid.png.
  2. Select GSSHA™ | Maps…
  3. For the coverage use the Soil Type coverage. Set the Coverage attribute to be ID. Do not use a second GIS data source.
  4. Change the result name to Soil Type
  5. Click on Coverages -> Index Map
  6. Select Done.

  7. Now we have an index map of the soils shapefile. Now we need to turn on the infiltration parameters in the Job Control.


  8. In the 2D Grid module Icon 2DGrid.png select GSSHA™ | Job Control.
  9. In the Job Control dialog change the infiltration option to Green + Ampt with Soil Moisture Redistribution.
  10. Click OK.

  11. Now we can set up the mapping tables for infiltration.


  12. In the 2D Grid Module Icon 2DGrid.png select GSSHA™ | Map Tables….
  13. Select the infiltration tab.
  14. In the Using index map box choose Soil Type.
  15. Select Generate IDs.

  16. The Generate IDs from Map button created three IDs, 2, 3, and 5. You will recall from when we created the land use index map that these three IDs came from the polygon IDs 2, 3, and 5. So the soil type index map ID #2 represents silt loam soils; the index map #3 represents silty clay loam; and the ID #5 represents silt.


  17. Using the following table, enter the values for each parameter.
  18. Table 9. Soil properties.
    ID Description Hyd. Conductivity (cm/hr) Capillary Head (cm) Porosity (m3/m3) Pore Dist. Index (cm/cm) Residual Saturation (m3/m3) Field Capacity (m3/m3)
    2 Silt Loam 0.68 16.68 0.501 0.234 0.015 0.33
    3 Silty Clay Loam 0.20 27.30 0.471 0.177 0.040 0.366
    5 Silty Clay 0.10 29.22 0.479 0.150 0.056 0.387

  19. Select the Initial Moisture tab in the process window.
  20. In the Using index map combo box select the uniform map.
  21. Select Generate IDs.
  22. Enter a value of 0.3 for the initial moisture.
  23. Click Done.

  24. You can now save the model and run.

Visualization

Since no water ran off we can guess that it all infiltratated. But let’s check to make sure.

  1. In the data tree, under the solution that was just read in, double-click on the summary file.
  2. Verify that there are no mass balance errors, check the amount that it rained and the amount that infiltrated, and close the file.

  3. Since we need more water lets increase the precipitation.


  4. In the 2D Grid Module select GSSHA™ | Precipitation.
  5. Increase the precipitation rate to 25 (mm/hr).
  6. Select OK.

  7. Let’s also turn on the infiltration output options.


  8. Select GSSHA™ | Job Control
  9. Select Output Control...
  10. Turn on Cumulative Infiltration Depth and Infiltration Rate.
  11. Select OK.
  12. Select OK.

  13. Save and run the simulation. You can now use the same techniques from Chapter 2 to visualize the infiltration data sets.

    One of the more interesting movies you can make is of the Infiltration Depth data set (To turn this dataset on, open the GSSHA Job Control window, click on Output Control, and toggle on the Cumulative Infiltration depth dataset. The two different soil types saturate at different rates and make for a pretty impressive movie.


GSSHA Tutorials

GSSHA Tutorial Download Website