WMS:GSSHA Groundwater: Difference between revisions

From XMS Wiki
Jump to navigationJump to search
No edit summary
 
 
(22 intermediate revisions by 3 users not shown)
Line 1: Line 1:
{{TOC right}}
Groundwater interaction in GSSHA models can only be performed in conjunction with one of the following infiltration methods:
Groundwater interaction in GSSHA models can only be performed in conjunction with one of the following infiltration methods:
* Green & Ampt with soil moisture redistribution
* Green & Ampt with soil moisture redistribution
Line 4: Line 5:


To model groundwater interaction in a GSSHA model:
To model groundwater interaction in a GSSHA model:
#Assign groundwater boundary conditions conceptually to feature points/nodes and arcs on the GSSHA coverage.  
#Assign groundwater boundary conditions conceptually to feature points/nodes and arcs on the GSSHA coverage.  
#Toggle on the Groundwater option in the [[WMS:GSSHA Job Control|Job Control]] and set the global parameters.
#Turn on the ''Groundwater'' option in the [[WMS:GSSHA Job Control|''GSSHA Job Control'' dialog]], and set the global parameters by clicking '''Edit parameters...''' to bring up the [[#GSSHA Groundwater Dialog|''GSSHA Groundwater'' dialog]].
#Toggle on the Sub-surface losses/gains option for GSSHA stream arcs in order to model groundwater interaction with the GSSHA channel network. (optional)
#Turn on ''Sub-surface losses/gains'' in the [[WMS:GSSHA Feature Arcs#Arcs Spreadsheet|GSSHA feature arcs ''Properties'' dialog]] in order to model groundwater interaction with the GSSHA channel network. This option is only available for trapezoidal channels. (optional)


==Groundwater Boundary Conditions==
==Modeling==
===Groundwater Boundary Conditions===
WMS can assign these boundary conditions to [[WMS:GSSHA Feature Nodes|feature points/nodes]]:
WMS can assign these boundary conditions to [[WMS:GSSHA Feature Nodes|feature points/nodes]]:
*Generic
*Generic – The cell is a regular infiltration cell, indicating no special boundary is desired.
*Constant head
*Constant head – The groundwater level in that cell remains constant for the duration of the simulation.
*Static well
*Static well &ndash;  The cell will have a source/sink term of constant rate, m<sup>3</sup>/s.
*Dynamic well
*Dynamic well &ndash;  A dynamic, temporally varying, pumping well exists in the cell.


The following groundwater boundary conditions are allowed for [[WMS:GSSHA Feature Arcs|feature arcs]]:
The following groundwater boundary conditions are allowed for [[WMS:GSSHA Feature Arcs|feature arcs]]:
*Generic
*Generic &ndash; Default boundary condition.
*No flow
*No flow &ndash; Usually created along the watershed boundaries or portions of the watershed boundaries. This typically represents a region where no lateral groundwater flow is permitted.
*Constant head
*Constant head &ndash; The groundwater level in that cell remains constant for the duration of the simulation.
*Flux river
*Flux river &ndash; Indicates that the cell contains a stream node and the flux between the stream node and saturated groundwater below the stream node be calculated based on Darcy’s law, as described by McDonald and Harbaugh (1988), allowing an exchange of water between the stream and groundwater during every groundwater update.
*Head river
*Head river &ndash; For the purposes of the saturated groundwater simulations, the cell will behave as a constant head cell and will not be updated during the saturated groundwater calculations.
 
WMS automatically generates groundwater boundary and well [[WMS:GSSHA Maps|index maps]], according to the groundwater boundary conditions conceptualized using feature points/nodes and arcs on the GSSHA coverage, when the ''Groundwater'' option is turned on in the [[WMS:GSSHA Job Control|''GSSHA Job Control'' dialog]]. These index maps appear in the "Index Maps" folder of the Project Explorer and are useful for visualizing boundary conditions and well locations and pumping rates.


WMS automatically generates Gw Boundary and Wells [[WMS:GSSHA Index Maps|index maps]], according to the groundwater boundary conditions conceptualized using feature points/nodes and arcs on the GSSHA coverage, when the Groundwater option is toggled on in the Job Control. These index maps will appear in the Index Maps folder of the Project Explorer and are useful for visualizing boundary conditions and well locations/pumping rates. If wells are defined then the Wells map table, which is not visible in the [[WMS:GSSHA Mapping Tables|Map Tables dialog]], is also populated when the Wells index map is created. Any changes that are made to either the boundary conditions or wells require that you regenerate the index maps (and wells map table) by right-clicking on either of the index maps in the Project Explorer and selecting the Regenerate command.
If wells are defined, the Wells map table—which is not visible in the [[WMS:GSSHA Mapping Tables|''Map Tables'' dialog]]—is also populated when the Wells index map is created. Any changes that are made to either the boundary conditions or wells require regenerating the index maps (and wells map table) by right-clicking on either of the index maps in the Project Explorer and selecting '''Regenerate''' .  


The groundwater boundary condition map and, if necessary, the well index map and map table are also written when the GSSHA .prj file is saved.
The groundwater boundary condition map and, if necessary, the well index map and map table are also written when the GSSHA *.prj file is saved.


==Global Parameters==
==GSSHA Groundwater Dialog==
[[Image:GsshaGroundwater.jpg]]
{{anchor|Global Parameters}}[[Image:GsshaGroundwater.jpg|thumb|350 px|''GSSHA Groundwater'' dialog]]
The ''GSSHA Groundwater'' dialog has a ''Groundwater'' section with two parts. The part on the left is for parameters that control the groundwater computations and includes these options:
*''Time step'' &ndash; A decimal number in seconds that specifies the groundwater model time step(s). Typically, the groundwater time step may be much larger than the overall model time step. Typical values are between "300" and "1200". The groundwater time step must be equal to or larger than the overall model time step, and the overall model time step must be an integer divisible into the groundwater model time step. During rainfall events, the groundwater model time step is temporarily changed to the overall model time step, and then changed back at the end of channel routing and overland flow routing. This switching is handled internally.
*''LSOR direction'' &ndash; Drop-down allowing "Vertical" or "Horizontal to be selected. This indicates the principal direction of flow is either in the vertical direction (along the y axis) or in the horizontal direction (along the x axis). The solution technique currently used to solve the 2D groundwater free surface flow equations is line successive over relaxation (LSOR). When LSOR is applied, the solution is either by rows or by columns. The solution should be aligned with the principal direction of saturated groundwater flow.
*''LSOR convergence'' &ndash; LSOR is an iterative method. The user-supplied convergence criterion is in meters of groundwater head or depth. The default value, 10<sup>-5</sup>, is typically sufficient for a good solution. Stricter criteria may lead to slower solutions or nonconvergence.
*''Relaxation coefficient'' &ndash; Used to project the current solution out into the future in an attempt to speed up convergence. A value of "1.0" (the default) indicates no projection. Values greater than "1.0" indicate that projections into the future are desired. Values up to "1.5" may speed the groundwater model solution. Typically, a value of about "1.2" provides the fastest solution. Increasing the relaxation coefficient to greater than "1.0" can result in nonconvergence of the solution. If this happens, the relaxation coefficient should be reduced. It is sometimes necessary to reduce the relaxation coefficient to less than "1.0" to obtain convergence of the solution.
*''Leakage rate'' &ndash; The leakage rate through the bottom of the aquifer in centimeters per hour. This is a uniform value that is applied to every cell in the grid. The default value is "0.0".
*''Aquifer cell size'' &ndash; The vertical cell size used with Richard's infiltration.


Enter parameters that control the groundwater computations in the GSSHA Groundwater dialog. The Aquifer cell size parameter is the vertical cell size used with Richard's infiltration. Specify continuous data sets used to define the aquifer bottom and water table. Hydraulic conductivity and porosity can be defined using continuous data sets or by assigning parameters in the Groundwater map table based on an index map generated using soil type data.
The part on the right is a spreadsheet with the following columns:
*''Name''
**''Aquifer bottom'' &ndash; Continuous dataset.
**''Water table'' &ndash; Continuous dataset.
**''Hydraulic conductivity'' &ndash; Continuous dataset or by assigning parameters in the Groundwater map table based on an index map generated using soil type data.
**''Porosity'' &ndash; Continuous dataset or by assigning parameters in the Groundwater map table based on an index map generated using soil type data.
*''Dataset'' &ndash; On each row is a drop-down listing all of the datasets in the project.


WMS writes all global groundwater parameters to the GSSHA .prj file. The groundwater map table (hydraulic conductivity and porosity) will be written to the .cmt file if Ids exist. Otherwise, hydraulic conductivity and porosity must be specified as continuous maps.
WMS writes all global groundwater parameters to the GSSHA *.prj file. The groundwater map table (hydraulic conductivity and porosity) will be written to the *.cmt file if IDs exist. Otherwise, hydraulic conductivity and porosity must be specified as continuous maps.


==Channels==
==Channels==
Interaction between the groundwater and the channel network is controlled for each link by toggling on the Sub-surface losses/gains option and specifying a sediment thickness and hydraulic conductivity. These values are written to the .cif file.
Interaction between the groundwater and the channel network is controlled for each link by toggling on the ''Sub-surface losses/gains'' option and specifying a sediment thickness and hydraulic conductivity. These values are written to the *.cif file.
 
==Related Topics==
* [[WMS:GSSHA Continuous Maps|GSSHA Continuous Maps]]
* [[WMS:GSSHA Feature Arcs|GSSHA Feature Arcs]]
* [[WMS:GSSHA Feature Nodes|GSSHA Feature Nodes]]
* [[WMS:GSSHA Index Maps|GSSHA Index Maps]]
* [[WMS:GSSHA Job Control|GSSHA Job Control]]
* [[WMS:GSSHA Mapping Tables|GSSHA Map Tables]]
 
{{GSSHAMain


|PLinks=
* [http://gsshawiki.com/index.php?title=Saturated_Flow:Modeling_Two-dimensional,_Saturated,_Lateral_Groundwater_Flow Groundwater]
* [http://gsshawiki.com/index.php?title=Richards_Equation:Modeling_the_Unsaturated_Zone_with_Richards_Equation Richards' Equation]


|GUMLinks=
{{GSSHA Wiki}}
* [http://gsshawiki.com/index.php?title=Project_File:Saturated_Groundwater_Flow_–_Optional Groundwater]
{{WMSMain}}
* [http://gsshawiki.com/index.php?title=Infiltration:Infiltration Infiltration]


}}
[[Category:GSSHA|Groundwater]]
[[Category:GSSHA Dialogs|Groundwater]]

Latest revision as of 16:13, 25 July 2024


Groundwater interaction in GSSHA models can only be performed in conjunction with one of the following infiltration methods:

  • Green & Ampt with soil moisture redistribution
  • Richards' equation

To model groundwater interaction in a GSSHA model:

  1. Assign groundwater boundary conditions conceptually to feature points/nodes and arcs on the GSSHA coverage.
  2. Turn on the Groundwater option in the GSSHA Job Control dialog, and set the global parameters by clicking Edit parameters... to bring up the GSSHA Groundwater dialog.
  3. Turn on Sub-surface losses/gains in the GSSHA feature arcs Properties dialog in order to model groundwater interaction with the GSSHA channel network. This option is only available for trapezoidal channels. (optional)

Modeling

Groundwater Boundary Conditions

WMS can assign these boundary conditions to feature points/nodes:

  • Generic – The cell is a regular infiltration cell, indicating no special boundary is desired.
  • Constant head – The groundwater level in that cell remains constant for the duration of the simulation.
  • Static well – The cell will have a source/sink term of constant rate, m3/s.
  • Dynamic well – A dynamic, temporally varying, pumping well exists in the cell.

The following groundwater boundary conditions are allowed for feature arcs:

  • Generic – Default boundary condition.
  • No flow – Usually created along the watershed boundaries or portions of the watershed boundaries. This typically represents a region where no lateral groundwater flow is permitted.
  • Constant head – The groundwater level in that cell remains constant for the duration of the simulation.
  • Flux river – Indicates that the cell contains a stream node and the flux between the stream node and saturated groundwater below the stream node be calculated based on Darcy’s law, as described by McDonald and Harbaugh (1988), allowing an exchange of water between the stream and groundwater during every groundwater update.
  • Head river – For the purposes of the saturated groundwater simulations, the cell will behave as a constant head cell and will not be updated during the saturated groundwater calculations.

WMS automatically generates groundwater boundary and well index maps, according to the groundwater boundary conditions conceptualized using feature points/nodes and arcs on the GSSHA coverage, when the Groundwater option is turned on in the GSSHA Job Control dialog. These index maps appear in the "Index Maps" folder of the Project Explorer and are useful for visualizing boundary conditions and well locations and pumping rates.

If wells are defined, the Wells map table—which is not visible in the Map Tables dialog—is also populated when the Wells index map is created. Any changes that are made to either the boundary conditions or wells require regenerating the index maps (and wells map table) by right-clicking on either of the index maps in the Project Explorer and selecting Regenerate .

The groundwater boundary condition map and, if necessary, the well index map and map table are also written when the GSSHA *.prj file is saved.

GSSHA Groundwater Dialog

GSSHA Groundwater dialog

The GSSHA Groundwater dialog has a Groundwater section with two parts. The part on the left is for parameters that control the groundwater computations and includes these options:

  • Time step – A decimal number in seconds that specifies the groundwater model time step(s). Typically, the groundwater time step may be much larger than the overall model time step. Typical values are between "300" and "1200". The groundwater time step must be equal to or larger than the overall model time step, and the overall model time step must be an integer divisible into the groundwater model time step. During rainfall events, the groundwater model time step is temporarily changed to the overall model time step, and then changed back at the end of channel routing and overland flow routing. This switching is handled internally.
  • LSOR direction – Drop-down allowing "Vertical" or "Horizontal to be selected. This indicates the principal direction of flow is either in the vertical direction (along the y axis) or in the horizontal direction (along the x axis). The solution technique currently used to solve the 2D groundwater free surface flow equations is line successive over relaxation (LSOR). When LSOR is applied, the solution is either by rows or by columns. The solution should be aligned with the principal direction of saturated groundwater flow.
  • LSOR convergence – LSOR is an iterative method. The user-supplied convergence criterion is in meters of groundwater head or depth. The default value, 10-5, is typically sufficient for a good solution. Stricter criteria may lead to slower solutions or nonconvergence.
  • Relaxation coefficient – Used to project the current solution out into the future in an attempt to speed up convergence. A value of "1.0" (the default) indicates no projection. Values greater than "1.0" indicate that projections into the future are desired. Values up to "1.5" may speed the groundwater model solution. Typically, a value of about "1.2" provides the fastest solution. Increasing the relaxation coefficient to greater than "1.0" can result in nonconvergence of the solution. If this happens, the relaxation coefficient should be reduced. It is sometimes necessary to reduce the relaxation coefficient to less than "1.0" to obtain convergence of the solution.
  • Leakage rate – The leakage rate through the bottom of the aquifer in centimeters per hour. This is a uniform value that is applied to every cell in the grid. The default value is "0.0".
  • Aquifer cell size – The vertical cell size used with Richard's infiltration.

The part on the right is a spreadsheet with the following columns:

  • Name
    • Aquifer bottom – Continuous dataset.
    • Water table – Continuous dataset.
    • Hydraulic conductivity – Continuous dataset or by assigning parameters in the Groundwater map table based on an index map generated using soil type data.
    • Porosity – Continuous dataset or by assigning parameters in the Groundwater map table based on an index map generated using soil type data.
  • Dataset – On each row is a drop-down listing all of the datasets in the project.

WMS writes all global groundwater parameters to the GSSHA *.prj file. The groundwater map table (hydraulic conductivity and porosity) will be written to the *.cmt file if IDs exist. Otherwise, hydraulic conductivity and porosity must be specified as continuous maps.

Channels

Interaction between the groundwater and the channel network is controlled for each link by toggling on the Sub-surface losses/gains option and specifying a sediment thickness and hydraulic conductivity. These values are written to the *.cif file.