Difference between revisions of "Tutorials:12 Groundwater"

Latest revision as of 20:04, 6 August 2009

This tutorial is compatible with:

• WMS Version 8.2 and later
• GSSHA Version 3.0b and later
  
  

Disclaimer: GSSHA tutorial exercises do not represent real world conditions

1. In the 2D Grid Module select GSSHA™ | Open Project File.
2. Browse to the Groundwater folder.
3. Select the EauGalleGSSHA.prj file and select Open.

Your watershed should look similar to the one shown in the image below. This is a subbasin of the Eau Galle watershed in Wisconsin.

Figure 35. Eau Galle Watershed.

GSSHA™ Model Setup

Hydrologic Modeling Wizard

1. At the bottom of the WMS window, click on the Hydrologic Modeling Wizard icon to open the Modeling Wizard.
2. Click on Select Model in the Wizard Task Bar on the left.
3. Select GSSHA for the desired model, but DO NOT click on the Initialize Model Data button.
4. Click on Define Land Use and Soil Data in the Wizard Task Bar on the left.
5. In the Define Land Use and Soil Data window, click the file browse button next to Add shapefile and open the file labeled EauGalleLU.shp. Make sure the Type is set to Land Use
6. Do the same for the soil type data and open the file labeled EauGalleST.shp. Make sure the Type is set to Soil Type
7. Click the Create Coverages button to create coverages from the land use and soil type data.
8. Click Next, Next, and Finish to create the Land Use coverage
9. Click Next, Next, and Finish to create the Soil Type Coverage
10. Click Next in the Hydrologic Modeling Wizard
11. In the Hydrologic Computations frame, click the Compute Index Mapping Tables button.



To model groundwater in GSSHA™ we will use 3 index maps. A soil type index map, a land use index map, and a combined soil type and land use index map. We will create these index maps from inside the wizard.



12. In the GSSHA™ Maps dialog, make sure the input coverage says Land Use and the coverage attribute says Id.
13. Change the Index map name to LU
14. Click the Coverages->Index Map button. You should see a land use index map appear in the WMS canvas window.
15. Now change the input coverage to Soil Type and the coverage attribute to Texture
16. Change the Index map name to ST
17. Click the Coverages->Index Map button. You should see a soil type index map appear in the WMS canvas window
18. Change Input coverage (1) back to Land Use and set the Coverage attribute to Id
19. Toggle on the option for Input coverage (2) and select the Soil Type coverage. Change the Coverage attribute to Texture.
20. Change the Index map name to Combined and click the Coverages->Index Map button. You should see a new index map appear in the WMS canvas window.
21. Click Done on the GSSHA™ Maps dialog

GSSHA™ Mapping Tables

1. In the GSSHA™ Map Table Editor, click on the Roughness tab
2. Set the index map to be LU and click Generate IDs
3. Set the Row Crops roughness to be 0.24, the Grass roughness to be 0.165, and the Forest roughness to be 0.235.
4. Click on the Interception tab.
5. Click Yes to turn the interception option on in the job control
6. Set the index map to be LU and click Generate IDs
7. Fill in the Interception spreadsheet with values from the following table:



Table 11. Values for Eau Galle Interception Row Crops Grass Forest Storage capacity (mm) 1.143 1.0 0.0423 Interception coeff (0.0-1.0) 0.045 0.027 0.18



8. Click on the Retention tab.
9. Click Yes to turn the retention option on in the job control
10. Set the index map to be LU and click Generate IDs
11. Leave the retention values at 0 for all three land uses
12. Click on the Evapotranspiration tab.
13. Click Yes to turn the evapotranspiration option on in the job control
14. In the GSSHA™ Job Control window, change the Evapotranspiration type to Penman method
15. Toggle on the option for Seasonal resist. and click OK
16. Set the index map to be LU and click Generate IDs
17. Fill in the Evapotranspiration spreadsheet with values from the following table:



Table 12. Values for Eau Galle Evapotranspiration Row Crops Grass Forest Land-surface albedo 0.2 0.2 0.2 Vegetation Height (m) 0.5 0.25 10.0 Vegetation Radiation Coeff 0.5 0.30 0.15 Canopy stomatal resistance (s/m) 45.0 100.0 120.0



18. Click on the Infiltration tab.
19. Click Yes to turn the infiltration option on in the job control
20. In the GSSHA™ Job Control window, change the Infiltration type to Green & Ampt with soil moisture redistribution and click OK.
21. Set the index map to be Combined and click Generate IDs
22. Fill in the Infiltration spreadsheet with values from the following table:



Table 13. Values for Eau Galle Infiltration Silt loam Silt loam Loam Loam Coarse sand Silt loam Loam Coarse sand Land ID #21 Row Crops Land ID #33 Grass Land ID #21 Row Crops Land ID #33 Grass Land ID #21 Row Crops Land ID #43 Forest Land ID #43 Forest Land ID #33 Grass Hydraulic Conductivity (cm/hr) 0.048342 0.156576 0.606576 1.09 0.2001 0.156576 0.606576 12.0 Capillary head (cm) 8.34 16.68 8.89 11.01 4.95 16.68 8.89 4.95 Porosity (m^3/m^3) 0.501 0.501 0.463 0.453 0.437 0.501 0.463 0.437 Pore Distribution index (cm/cm) 0.234 0.234 0.252 0.378 0.234 0.234 0.252 0.694 Residual saturation (m^3/m^3) 0.015 0.015 0.027 0.041 0.02 0.015 0.027 0.02 Field capacity (m^3/m^3) 0.33 0.33 0.27 0.207 0.12 0.33 0.27 0.091 Wilting Point (m^3/m^3) 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12



23. Click on the Initial Moisture tab
24. Set the index map to be ST and click Generate IDs
25. Set the Coarse sand initial moisture to be 0.371, the Loam initial moisture to be 0.393, and the Silt loam initial moisture to be 0.42585.
26. Click Done to close the GSSHA™ Map Table Editor

Creating Groundwater Datasets

To add groundwater to the GSSHA™ model, we need to create five more datasets: a porosity map, a hydraulic conductivity map, an aquifer bottom map, a water table elevation map, and a boundary condition map. We'll set up the continuous maps outside of the Hydrologic Modeling Wizard.

Hydraulic Conductivity and Porosity

1. Click Close to close the Hydrologic Modeling Wizard.
2. In the 2D Grid Module select GSSHA™ | Maps...
3. Click on the tab labeled Continuous - Grid
4. Click on the Data Calculator... button
5. In the Expression field, enter 11.197
6. In the Result field, type HK, for hydraulic conductivity
7. Click the Compute button to create the continuous hydraulic conductivity map.
8. In the Expression field, enter 0.4
9. In the Result field, type Porosity
10. Click the Compute button to create the continuous porosity map.
11. Click Done to close the Data Calculator
12. Click Done to close the GSSHA™ Maps dialog

Aquifer Bottom

We will create the aquifer bottom map and the water table elevation map from x,y,z points that we have stored in a comma separated file.



1. In the WMS window, select File | Open
2. Open the file named aquifer_bot_elevs.csv
3. In the File Import Wizard, change the starting import row to 2
4. Make sure the Comma delimeter is toggled on and click Next>
5. Make sure the WMS data type is set to 2D Scatter Points
6. Make sure the first column is mapped to X, the second column to Y, and the third column to Data set
7. Click Finish to import the dataset as a 2D Scatter Point dataset



A new scatter point dataset named aquifer_bot_elevs should appear in the Project Explorer.




Figure 36. Aquifer Bottom Scatter Point Dataset


8. Right click on the aquifer_bot_elevs scatter set and choose Interpolate->...To Grid
9. Change the interpolated data set name to Aq_Bottom. DO NOT toggle on the option to Map elevations.
10. Click OK to create a grid from the aquifer bottom scatter points



You should now see a new continuous map named Aq_Bottom among the other continuous datasets.

Water Table Elevation

1. In the WMS window, select File | Open
2. Open the file named water_table_elevs.csv
3. In the File Import Wizard, change the starting import row to 2
4. Make sure the Comma delimeter is toggled on and click Next>
5. Make sure the WMS data type is set to 2D Scatter Points
6. Make sure the first column is mapped to X, the second column to Y, and the third column to Data set
7. Click Finish to import the dataset as a 2D Scatter Point dataset



A new scatter point dataset named water_table_elevs should appear in the Project Explorer.




Figure 37. Water Table Elevation Scatter Point Dataset


8. Right click on the water_table_elevs scatter set and choose Interpolate->...To Grid
9. Change the interpolated data set name to WTE. DO NOT toggle on the option to Map elevations.
10. Click OK to create a grid from the water table scatter points



You should now see a new continuous map named WTE among the other continuous datasets. Now we need to associate the new datasets with the Eau Galle GSSHA™ model.



11. Right click on the Continuous Maps folder underneath the EauGalleGSSHA folder in the 2D Grid Data section of the Project Explorer.
12. Select Assign->Aq_Bottom
13. Right click on the Continuous Maps folder again and select Assign->WTE



There should now be 5 continuous maps associated with the Eau Galle GSSHA™ model.

Groundwater Boundary Condition

1. Turn off the Scatter Point datasets in the Project Explorer
2. Click on the GSSHA™ coverage in the Project Explorer to make it the active coverage
3. Choose the Select Feature Line Branch tool and double click on the most downstream stream arc
4. In the All row (highlighted in yellow) of the Properties spreadsheet, change the Type to Trapezoidal Channel



For this example, we will assume the stream geometry is the same everywhere in the model.



5. Still in the All row, enter a Manning's n of 0.119, a Depth of 0.5, a Bottom Width of 1.0, and a Side slope of 4.2
6. For the Groundwater BC, select Flux River. This will assign a flux river boundary condition to all the stream cells that underlie the stream arcs.
7. Click OK to close the Properties dialog
8. Switch to the 2D Grid Module
9. With the 2D Grid Module selected, select GSSHA™ | Job Control...
10. Change the Channel routing computation scheme to Diffusive wave
11. Toggle on the Groundwater option and click Edit parameter...



When toggling on the Groundwater option, a groundwater boundary index map named Gw boundary is created automatically in the Project Explorer. Notice the stream cells representing the flux river boundary condition on the Gw boundary grid. Also, a no-flow boundary condition is assumed to exist around the perimeter of the watershed.

Groundwater Job Control

1. In the GSSHA Groundwater dialog, set the Aquifer... Data Set to Aq_Bottom, the Water... Data Set to WTE, the Hydra... Data Set to HK, and the Porosity Data Set to Porosity.
2. Set the Time Step to be 600
3. Make sure the LSOR direction is set to Vertical
4. Make sure the LSOR convergence is 0.00001
5. Make sure the Relaxation coefficient is 1.2
6. Make sure the Leakage rate is set to 0.0
7. Click OK to close the GSSHA Groundwater dialog
8. In the GSSHA™ Job Control window, make sure the Overland flow Computation method is set to ADE
9. Click OK to close the GSSHA™ Job Control Parameters dialog

Precipitation

Groundwater is most often used for long-term simulations. To finish setting up a long-term GSSHA™ simulation, we need to enter long-term precipitation data.

1. With the 2D Grid Module selected, select GSSHA™ | Precipitation
2. In the Rainfall event(s) combo box, change the type to Gage.
3. Click the Import Gage File... button
4. Open the file named rain_combine_events_may15.gag
5. Change the Multi-gage interpolation method to Thiessen polygons



You'll notice many storm events are now shown in the Project Explorer.



6. Click OK to close the GSSHA™ Precipitation dialog
7. Select GSSHA™ | Job Control...
8. Toggle on the option for Long term simulation and click Edit parameter...
9. Set the Latitude to be 44.81
10. Set the Longitude to be 267.83
11. Set the GMT at -6.0
12. Enter a Minimum event discharge of 0.1
13. Enter a Soil moisture depth of 0.25
14. Click the Browse button next to HMET Data File
15. Open the file named Wisconsin_HMet_formatted2002_may15.txt
16. Change the HMET Data Format to WES
17. Click OK to close the Continuous Simulation dialog
18. Click OK to close the GSSHA™ Job Control Parameters dialog

Saving and Running the GSSHA™ Model

We are now ready to save and run the GSSHA™ model.

1. With the 2D Grid Module selected, select GSSHA™ | Save Project File
2. Save the GSSHA™ Project as EauGalleGSSHARun1 and click Save
3. Select GSSHA™ | Run GSSHA™
4. Toggle off the option to Suppress screen printing and click OK

The model should take about an hour and a half to two hours to run, depending on the speed of your computer. If you would like to view the solution without waiting for the model to run to completion, click the Abort button, select GSSHA | Read Solution and open the solution for EauGalleGSSHASed.prj in the Sediment folder.

Viewing GSSHA™ Model Results

1. Once the model has finished running, make sure the Read solution on exit option is toggled on and click Close. It may take some time to read in the solution.
2. Once the solution file has been read in, choose the Select hydrograph tool and double click on the hydrograph icon near the watershed outlet
3. This will open up an outflow hydrograph plot which you can copy to the clipboard or export to Excel (Figure 38)



Figure 38. Eau Galle Groundwater Model Outflow Hydrograph.


4. You may also wish to explore the summary file or view depth contours
5. Though it will take some time, you may also wish to create a movie file of your rendered depth contours
6. You can tell GSSHA™ to output more than just the depth dataset. If you wish to view groundwater specific datasets such as infiltration rate or groundwater elevations, change the Output Control settings in the GSSHA™ Job Control dialog and rerun the model.

This concludes the GSSHA™ Groundwater Tutorial.

GSSHA Tutorials

GSSHA Tutorial Download Website