Difference between revisions of "Surface Water Routing:Overland Flow Routing"

From Gsshawiki
Jump to: navigation, search
(5.2.2 Overland Flow Hydraulic Roughness)
 
(10 intermediate revisions by 2 users not shown)
Line 12: Line 12:
 
[[Image:Equation008.gif]]  (8)   
 
[[Image:Equation008.gif]]  (8)   
  
In addition to this original formulation in the CASC2D model, two additional methods of solving the equations have been added, an alternating direction explicit scheme (ADE) and an ADE scheme with an additional predictor-corrector step (ADE-PC) (Downer 2002a; Downer et al., 2000).  Both the ADE and ADE-PC methods employ the up-gradient difference technique, Equation 3, for flows in the upstream direction (Downer, 2002a).  Fluxes other than inter-cell fluxes, direct evaporation (DET), infiltration, exfiltration, are accounted for before overland routing is computed.
+
In addition to this original formulation in the CASC2D model, two additional methods of solving the equations have been added, an alternating direction explicit scheme (ADE) and an ADE scheme with an additional predictor-corrector step (ADE-PC) (Downer 2002a; Downer et al., 2000).  Both the ADE and ADE-PC methods employ the up-gradient difference technique, Equation 3, for flows in the upstream direction (Downer, 2002a). The explicit scheme is being deprecated in v8.0 and will not be available in future verions.  Fluxes other than inter-cell fluxes, direct evaporation (DET), infiltration, exfiltration, are accounted for before overland routing is computed.
  
In the ADE method, inter-cell flows are first calculated in the ''x'' direction according to Equation 6.  Depths in each row are updated based on the flows in the ''x'' direction:
+
In the ''ADE'' method, inter-cell flows are first calculated in the ''x'' direction according to Equation 6.  Depths in each row are updated based on the flows in the ''x'' direction:
  
 
[[Image:Equation009.gif]]    (9)
 
[[Image:Equation009.gif]]    (9)
Line 26: Line 26:
 
[[Image:Equation011.gif]]      (11)
 
[[Image:Equation011.gif]]      (11)
  
With the ADE-PC method additional steps are added to improve accuracy and stability.  As before, during each sweep, by rows or by columns, an estimate of heads is made based on the calculated flows, Equations 9 and 11.  Next, using the updated depths, updated estimates of flow are computed at the ''n''+''1'' time level
+
With the ''ADE-PC'' method additional steps are added to improve accuracy and stability.  As before, during each sweep, by rows or by columns, an estimate of heads is made based on the calculated flows, Equations 9 and 11.  Next, using the updated depths, updated estimates of flow are computed at the ''n''+''1'' time level
  
 
[[Image:Equation012.gif]]      (12)
 
[[Image:Equation012.gif]]      (12)
Line 36: Line 36:
 
These flows are then used to update the original depths, Equations 9 and 11.  This procedure is essentially the MacCormack method (MacCormack, 1969) except up-gradient differences are used in both the predictor and corrector steps.  A similar method was successfully implemented by Wang and Hjelmfelt (1998).   
 
These flows are then used to update the original depths, Equations 9 and 11.  This procedure is essentially the MacCormack method (MacCormack, 1969) except up-gradient differences are used in both the predictor and corrector steps.  A similar method was successfully implemented by Wang and Hjelmfelt (1998).   
  
In benchmark tests using the three methods: original explicit (EXPLICIT), ADE, and ADE-PC, in an contrived watershed consisting of two converging planes (open book), the Goodwin Creek Experimental Watershed (GCEW) (Senarath et al., 2000), and Poplar Creek (Downer et al., 2002a), the ADE and ADE-PC methods ran with significantly larger time steps (Downer et al., 2000).  Depending on the test case, time steps could be increased from 20% to 240% with commensurate decreases in simulation times.   
+
In benchmark tests using the three methods: original explicit ''EXPLICIT'', ''ADE'', and ''ADE-PC'', in an contrived watershed consisting of two converging planes (open book), the Goodwin Creek Experimental Watershed (GCEW) (Senarath et al., 2000), and Poplar Creek (Downer et al., 2002a), the ''ADE'' and ''ADE-PC'' methods ran with significantly larger time steps (Downer et al., 2000).  Depending on the test case, time steps could be increased from 20% to 240% with commensurate decreases in simulation times.   
  
The routing scheme is selected using the '''OVERTYPE''' card.  The default value is EXPLICIT.  The most efficient scheme to use depends on the particular watershed.  The ADE-PC scheme can generally handle rougher terrain and typically requires less smoothing of the DEM, but the additional computational steps result in greater computation time, unless use of the ADE-PC scheme permits substantially greater time steps than with one of the other two methods.  For smoothed DEMs or in watersheds with smoother terrain, the ADE and EXPLICIT methods usually can be employed, with a resulting savings in execution time.   
+
The routing scheme is selected using the '''OVERTYPE''' card.  The default value is ''ADE''.  The most efficient scheme to use depends on the particular watershed.  The ''ADE-PC'' scheme can generally handle rougher terrain and typically requires less smoothing of the DEM, but the additional computational steps result in greater computation time, unless use of the ''ADE-PC'' scheme permits substantially greater time steps than with one of the other two methods.  For smoothed DEMs or in watersheds with smoother terrain, the ''ADE'' and ''EXPLICIT'' methods usually can be employed, with a resulting savings in execution time.   
  
To improve stability, the timestep in both the EXPLICIT and ADE schemes is variable.  The time step is not variable in the ADEPC.  The addition of a variable timestep has allowed for increases in the overall model timestep when using either the EXPLICIT and ADE methods.  These new methods have not been benchmarked against the ADE-PC method.
+
To improve stability, the timestep in both the ''EXPLICIT'' and ''ADE'' schemes is variable.  The time step is not variable in the ''ADE-PC'' method.  The addition of a variable timestep has allows for increased overall model timestep when using either the ''EXPLICIT'' and ''ADE'' methods.  These new methods have not been benchmarked against the ''ADE-PC'' method.
 +
 
 +
in GSSHA version 8.0 the ''EXPLICIT'' option is deprecated and will not be available in future versions.
 +
 
 +
In GSSHA version 7.1.1 and higher, an inertial formulation of the flow equations, as described in Bates et al. (2010), is included as an option.  To use this option include the '''OVERLAND_MOMENTUM''' card in your project file along with ''ADE'' specified with the '''OVERTYPE''' card .  This option may be helpful for simulations where overland hydrodynamics are important, such as flooding due to a storm surge.  For these cases stability may be increased due the inclusion of inertial terms and a better accounting of the effects of friction.  The method may allow longer time steps and potentially reduce simulation times.  In addition, in GSSHA v7.1.1 and higher, a strict overland time step for the diffusive wave formulation, as described in Bates et al. (2010) has been included.  This option is specified by including the '''OVERLAND_STRICT_DT''' card in the project file.  This option also only works with the ''ADE'' method.  While use of this option may increase stability for high overland flows, the timestep is very restrictive compared to the normal criteria used for the ''ADE'' method.  If instability arises with use of the ''ADE'' method, it may be better to reduce the overall timestep, use the inertial formulation for flow '''OVERLAND_MOMENTUM''', or use the ''ADE-PC'' method, in lieu of the '''OVERLAND_STRICT_DT''' option.  For versions 7.14 and beyond, the user can control the time step coefficient for the momentum formulation by specifying a real value between 0.0 and 1.0 after the '''OVERLAND_MOMENTUM''' card.  The default value is 0.2.  Larger values increase speed, lower values increase stability.  The user is referred to Bates et al. (2010) for details:  Bates, P. D., Horritt, M. S., and T. J. Fewtrell (2010) A simple inertial formulation of the shallow water equations for efficient two-dimensional flood inundation modelling. ''J. Hydro.'' 387 (2010) 33-45. 
  
 
Special overland flow cell types that may defined and result in a special type of overland routing are wetlands described in the GSSHA primer here: http://www.gsshawiki.com/Wetlands:Wetlands
 
Special overland flow cell types that may defined and result in a special type of overland routing are wetlands described in the GSSHA primer here: http://www.gsshawiki.com/Wetlands:Wetlands
Line 109: Line 113:
 
n<sub>d</sub>=n<sub>0</sub>e<sup>-Bd</sup>
 
n<sub>d</sub>=n<sub>0</sub>e<sup>-Bd</sup>
  
where n<sub>d</sub> is the Manning roughness at depth (d), n<sub>0</sub> is the specified roughness in the '''ROUGHNESS''' mapping table, and B is an exponent specified in the '''ROUGHNESS''' mapping table.  The default value for B is zero, which results in a static value roughness with depth.  To modify this relationship, a value of B is specified for each roughness category in the '''ROUGHNESS''' table. The value is listed in the table after the specified roughness values, n<sub>0</sub>.  So that the table will look like
+
where ''n<sub>d</sub>'' is the Manning roughness at depth (d), ''n<sub>0</sub>'' is the specified roughness in the '''ROUGHNESS''' mapping table, and ''B'' is an exponent specified in the '''ROUGH_EXP''' mapping table.  The default value for ''B'' is zero, which results in a static value roughness with depth.  To modify this relationship, a value of ''B'' is specified for each roughness category in the '''ROUGH_EXP''' table. So that the tables will look like
  
 +
<pre>
 
ROUGHNESS "unif"
 
ROUGHNESS "unif"
|-
 
 
NUM_IDS 1
 
NUM_IDS 1
ID    DESCRIPTION1                            DESCRIPTION2                                  ROUGH EXPONENT
+
ID    DESCRIPTION1                            DESCRIPTION2                                  ROUGH  
1    Roughness ID                                                                       0.300000 0.5
+
1    Roughness ID                                                                           0.300000  
 +
 
 +
ROUGH_EXP "unif"
 +
NUM_IDS 1
 +
ID    DESCRIPTION1                            DESCRIPTION2                                  ROUGH_EXP
 +
1    Roughness ID                                                                          0.5
 +
</pre>
  
Positive values of B result in the roughness decreasing with depth, common; negative values result in roughness increasing with depth, uncommon.  A typical value is 0.5.  The exponent may require calibration.
+
Positive values of ''B'' result in the roughness decreasing with depth, common; negative values result in roughness increasing with depth, uncommon.  A typical value is 0.5.  The exponent may require calibration.
  
  
Line 131: Line 141:
 
==5.2.4 Specifying Initial Depths on the Watershed==
 
==5.2.4 Specifying Initial Depths on the Watershed==
  
Initial depths on the overland flow plane at the beginning of the simulation may be specified by use of the '''INITIAL_DEPTH''' project file card.  This card is used to specify an ASCII GRASS map containing initial overland depth values, m.  This feature is rarely used.
+
Initial depths on the overland flow plane at the beginning of the simulation may be specified by use of the '''OV_INIT_DEPTH''' project file card.  This card is used to specify an ASCII GRASS map containing initial overland depth values, m.  The initial depth can also be specified with the '''READ_OV_HOTSTART''' card.  A map of the final overland flow values from a simulation can be written out using the '''WRITE_OV_HOTSTART''' card.  This map can then be used to hotstart subsequent simulations.  A single value of initial water surface elevation (m) can be specified using '''INIT_ELEV_HEAD''' card.  The water surface elevation in every cell with an elevation less the specified elevation will be set to this value and the depths computed from the specified elevation and the cell elevation.
  
 
==5.2.5 Simulations without Channel Routing==
 
==5.2.5 Simulations without Channel Routing==

Latest revision as of 14:57, 13 October 2023

5.2.1 Overland Flow Routing Formulation

Overland flow in GSSHA employs the same methods described for 1-D channel routing, except the calculations are made in two dimensions. Flow is routed in two orthogonal directions in each grid cell during each time step. The watershed boundary represents a no flow boundary for the overland flow routing and when a grid cell lies on the watershed boundary, flow is not routed across the boundary. There is an option in GSSHA to impose head boundaries along the watershed boundary and compute flows across the overland boundary. In GSSHA, Δx = Δy. Inter-cell fluxes in the x and y directions, p and q, respectively, are computed in cell ij from the depth, dij, at the nth time level using the Manning equation for the head discharge relationship in the x and y directions, respectively, as

Equation006.gif (6)

Equation007.gif (7)

Depths in each cell are calculated at the n+1 time level based on the flows for each cell (Julien and Saghafian, 1991):

Equation008.gif (8)

In addition to this original formulation in the CASC2D model, two additional methods of solving the equations have been added, an alternating direction explicit scheme (ADE) and an ADE scheme with an additional predictor-corrector step (ADE-PC) (Downer 2002a; Downer et al., 2000). Both the ADE and ADE-PC methods employ the up-gradient difference technique, Equation 3, for flows in the upstream direction (Downer, 2002a). The explicit scheme is being deprecated in v8.0 and will not be available in future verions. Fluxes other than inter-cell fluxes, direct evaporation (DET), infiltration, exfiltration, are accounted for before overland routing is computed.

In the ADE method, inter-cell flows are first calculated in the x direction according to Equation 6. Depths in each row are updated based on the flows in the x direction:

Equation009.gif (9)

Inter-cell flows in the y direction are computed using the updated depths:

Equation010.gif (10)

Depths in each column are updated based on the flows in the y direction:

Equation011.gif (11)

With the ADE-PC method additional steps are added to improve accuracy and stability. As before, during each sweep, by rows or by columns, an estimate of heads is made based on the calculated flows, Equations 9 and 11. Next, using the updated depths, updated estimates of flow are computed at the n+1 time level

Equation012.gif (12)

The original flows and the updated flows are then averaged to come up with an estimate of flows for the time step:

Equation013.gif (13)

These flows are then used to update the original depths, Equations 9 and 11. This procedure is essentially the MacCormack method (MacCormack, 1969) except up-gradient differences are used in both the predictor and corrector steps. A similar method was successfully implemented by Wang and Hjelmfelt (1998).

In benchmark tests using the three methods: original explicit EXPLICIT, ADE, and ADE-PC, in an contrived watershed consisting of two converging planes (open book), the Goodwin Creek Experimental Watershed (GCEW) (Senarath et al., 2000), and Poplar Creek (Downer et al., 2002a), the ADE and ADE-PC methods ran with significantly larger time steps (Downer et al., 2000). Depending on the test case, time steps could be increased from 20% to 240% with commensurate decreases in simulation times.

The routing scheme is selected using the OVERTYPE card. The default value is ADE. The most efficient scheme to use depends on the particular watershed. The ADE-PC scheme can generally handle rougher terrain and typically requires less smoothing of the DEM, but the additional computational steps result in greater computation time, unless use of the ADE-PC scheme permits substantially greater time steps than with one of the other two methods. For smoothed DEMs or in watersheds with smoother terrain, the ADE and EXPLICIT methods usually can be employed, with a resulting savings in execution time.

To improve stability, the timestep in both the EXPLICIT and ADE schemes is variable. The time step is not variable in the ADE-PC method. The addition of a variable timestep has allows for increased overall model timestep when using either the EXPLICIT and ADE methods. These new methods have not been benchmarked against the ADE-PC method.

in GSSHA version 8.0 the EXPLICIT option is deprecated and will not be available in future versions.

In GSSHA version 7.1.1 and higher, an inertial formulation of the flow equations, as described in Bates et al. (2010), is included as an option. To use this option include the OVERLAND_MOMENTUM card in your project file along with ADE specified with the OVERTYPE card . This option may be helpful for simulations where overland hydrodynamics are important, such as flooding due to a storm surge. For these cases stability may be increased due the inclusion of inertial terms and a better accounting of the effects of friction. The method may allow longer time steps and potentially reduce simulation times. In addition, in GSSHA v7.1.1 and higher, a strict overland time step for the diffusive wave formulation, as described in Bates et al. (2010) has been included. This option is specified by including the OVERLAND_STRICT_DT card in the project file. This option also only works with the ADE method. While use of this option may increase stability for high overland flows, the timestep is very restrictive compared to the normal criteria used for the ADE method. If instability arises with use of the ADE method, it may be better to reduce the overall timestep, use the inertial formulation for flow OVERLAND_MOMENTUM, or use the ADE-PC method, in lieu of the OVERLAND_STRICT_DT option. For versions 7.14 and beyond, the user can control the time step coefficient for the momentum formulation by specifying a real value between 0.0 and 1.0 after the OVERLAND_MOMENTUM card. The default value is 0.2. Larger values increase speed, lower values increase stability. The user is referred to Bates et al. (2010) for details: Bates, P. D., Horritt, M. S., and T. J. Fewtrell (2010) A simple inertial formulation of the shallow water equations for efficient two-dimensional flood inundation modelling. J. Hydro. 387 (2010) 33-45.

Special overland flow cell types that may defined and result in a special type of overland routing are wetlands described in the GSSHA primer here: http://www.gsshawiki.com/Wetlands:Wetlands

And overland lakes, which are lakes not connected to the stream network. Overland lakes are specified with the OVERLAND_LAKE_MAP card that specifies the overland lake mask. The overland lake mask is a map with zeros and integer values, like the MASK_MAP except the integer values occur where overland lakes are desired. This mask map can be created in WMS by making a coverage in WMS for overland lakes and then assigning integer values to polygons that define the lakes, and then creating an index map from that WMS coverage. The significance of an overland lake is that no overland routing occurs within the overland lake cells. Other overland processes, such as infiltration, still occur.

5.2.2 Overland Flow Hydraulic Roughness

The GSSHA model requires that Manning roughness coefficients be assigned to every cell in the watershed mask. There are three ways to specify the hydraulic roughness of the overland flow planes in GSSHA. The first method is to apply a constant value over the entire watershed through the MANNING_N project file card. The second method is to use the MAPPING_TABLE to assign roughness coefficients using tabled values referenced to an index map. The third method is to produce a GRASS ASCII map of roughness coefficient, and provide the name of this map to GSSHA using the ROUGHNESS project file card. Table 8 provides typical values of the Manning roughness coefficient for overland flow over various surfaces:

Land Use or Cover Recommended
n-value
Range
Concrete or asphalt 0.011a 0.01-0.013a
0.05-0.15d
Developed/industrial 0.0137b -
Bare sand 0.01a 0.010-0.016a
Graveled surface 0.02a 0.012-0.03a
Bare clay-loam (eroded) 0.02a 0.012-0.033a
Gullied land - 0.320-0.357c
Bare field – no residue 0.05a 0.006-0.16a
Range (natural) 0.13a 0.01-0.32a
Range (clipped) 0.10a 0.02-0.24a
Grass and pasture - 0.05 – 0.15a
Pasture -
-
0.235-271c
0.30-0.40d
Clover - 0.08 – 0.25a
Small grain - 0.1 – 0.4a
Row crops - 0.07 – 0.2a
Cotton/soy - 0.246-0.261c
Grass (bluegrass sod) 0.45a 0.39-0.63a
0.30-0.50d
Short grass prairie 0.15a 0.10-0.20a,d
Dense grass 0.24a 0.17-0.30a
Bermuda grass 0.41a 0.30-0.48a
Lawns   0.40-0.50d
Forest 0.192b 0.184-198c
Sparsely vegetated 0.150b 0.05-0.13d
Dense Growth   0.40-0.50d

Table 8 - Values of overland flow roughness coefficient Notes: aEngman (1986), bDowner (2002b), cSenarath et al (2000), dHEC (1985)

These values should be considered guidelines. Manning roughness coefficients are typically assigned from literature values and then adjusted through calibration. Additional sources of literature values are Liong et al. (1989), Engman (1986), and Ree et al. (1977). If calibrated values differ significantly from published values, there may be appropriate justification.

Typical Manning roughness coefficient values for open channel flow are considerably smaller than overland flow values because of deeper flow depths in the channel. In cases where the overland flow may become very deep, such as simulating a tidal surge, the user may want to use depth varying overland roughness. GSSHA is currently formulated to calculate the overland roughness at any depth (d) using the formula

nd=n0e-Bd

where nd is the Manning roughness at depth (d), n0 is the specified roughness in the ROUGHNESS mapping table, and B is an exponent specified in the ROUGH_EXP mapping table. The default value for B is zero, which results in a static value roughness with depth. To modify this relationship, a value of B is specified for each roughness category in the ROUGH_EXP table. So that the tables will look like

ROUGHNESS "unif"
NUM_IDS 1
ID    DESCRIPTION1                            DESCRIPTION2                                   ROUGH   
1     Roughness ID                                                                           0.300000 

ROUGH_EXP "unif"
NUM_IDS 1
ID    DESCRIPTION1                            DESCRIPTION2                                   ROUGH_EXP
1     Roughness ID                                                                           0.5

Positive values of B result in the roughness decreasing with depth, common; negative values result in roughness increasing with depth, uncommon. A typical value is 0.5. The exponent may require calibration.


5.2.3 Runoff Retention

Natural land surfaces contain micro-topography, small depressions, that retain water prior to runoff. The water held in the grid cell, or retention storage, never becomes direct runoff and can only be removed from the land surface as infiltration or direct evaporation. In certain regions, the retention storage can be significant. Retention storage is input as a depth (mm) in each grid cell and may be optionally input to GSSHA as:

  • a uniform value using the RETENTION card,
  • as a table of values related to index maps by using the MAPPING_TABLE project card, and the RETEN_DEPTH card without specifying a file name or
  • as an ASCII GRASS map through the use of the RETEN_DEPTH project file card with a specified file name.

5.2.4 Specifying Initial Depths on the Watershed

Initial depths on the overland flow plane at the beginning of the simulation may be specified by use of the OV_INIT_DEPTH project file card. This card is used to specify an ASCII GRASS map containing initial overland depth values, m. The initial depth can also be specified with the READ_OV_HOTSTART card. A map of the final overland flow values from a simulation can be written out using the WRITE_OV_HOTSTART card. This map can then be used to hotstart subsequent simulations. A single value of initial water surface elevation (m) can be specified using INIT_ELEV_HEAD card. The water surface elevation in every cell with an elevation less the specified elevation will be set to this value and the depths computed from the specified elevation and the cell elevation.

5.2.5 Simulations without Channel Routing

It might be desirable to perform simulations without channel routing in very small watersheds lacking a defined channel network, or in the beginning stages of developing a GSSHA model of a watershed. It is always prudent to build a GSSHA model one process at a time. Getting just the overland flow portion of the model to run is always the first step in building a model. For these reasons, GSSHA may be used without channel routing. Because the channel network normally provides the outlet point, the overland cell containing the outlet must be specified during simulations without channel routing. Normally this grid cell will have the lowest elevation in the watershed. The row and column containing the outlet grid cell are specified using the OUTROW and OUTCOL cards, while the slope of the outlet grid cell is specified using the OUTSLOPE card. If channel routing is enabled through the inclusion of CHAN_EXPLIC in the project file, the OUTROW, OUTCOL, and OUTSLOPE cards are not required, and ignored if present.

GSSHA User's Manual

5 Surface Water Routing
5.1     Channel Routing
5.2     Overland Flow Routing
5.3     Channel Boundary Conditions
5.4     Overland Boundary Conditions
5.5     Embankments
5.6     Overland/Channel Interaction
5.7     Introducing Discharge/Constituent Hydrographs
5.8     Overland Routing with Snow
5.9     Overland Routing with BMPs