Difference between revisions of "Process representation across tools"

From Hydromodel SA Wiki
Jump to navigation Jump to search
m
 
(58 intermediate revisions by the same user not shown)
Line 1: Line 1:
 +
<span id="Process rep - pg top"></span>
 +
The tables below summarise information about different process representation [[Terminology#algorithm anchor|algorithms]] across the set of modelling tools. They cover the inputs used to model the occurrence and rate (amount per model calculation timestep) of the given hydrological process (e.g. infiltration of water into the soil) as well as some general characteristics of the equations, particularly noting if there are thresholds involved and what these are. The tables can be used to understand what is being considered in the calculation of each process and what inputs each modelling tool needs. Not all the processes listed are represented explicitly by all the tools. This is highlighted in the tables below and an overview is also presented in the [[Modelling tool capability overview#Modelling tool capabilities summary|capabilities overview table]]. The nature of the algorithms used and their inputs are linked to the '''spatial and temporal scale''' at which the process is being represented in the given tool, described further for each one [[Model units & connections|here]].
 +
</br></br>
 +
In these tables '''thresholds''' refer to limits that determine when a process would start to occur or would stop occurring.  For example, the field capacity soil moisture of a soil is often used as the threshold soil moisture level for percolation of water downwards (to a lower soil layer or to recharge groundwater). If soil moisture is lower than this, no percolation will be calculated in the model. There may be multiple thresholds considered when modelling a process. For example, evapotranspiration of soil moisture may be assumed to stop once the atmospheric demand in the timestep has been met or once the soil moisture has reached wilting point level, even if demand has not yet been met.   
 +
</br></br>
 +
The inputs to a process algorithm can include:
 +
* input data (e.g., rainfall, using the rainfall input for a particular modelled unit, for a given timestep),
 +
* property parameter values that the user inputs in the model set-up (e.g., soil moisture content at saturation for a soil layer in a modelled unit)
 +
* states or water storage levels that the model calculated internally for the timestep, so not directly input by the user (e.g., soil moisture content in a certain soil layer at a given timestep)
 +
Different software tools may refer to equivalent inputs using different words and sometimes they require the user to input different specific property values, but end up calculating the same derived property. For example, on software tool may require porosity of soil to be input, another may require bulk density, and another may as for volumetric soil moisture at saturation, any of these, coupled with soil layer depth, can be used to calculate the maximum (or saturation) soil water storage volume. Where possible, similar terms have been used across tools in the tables below to highlight their conceptual similarities where these exist.  Details about some of the contrasting terminology used in the interfaces of different tools can be found [[Terminology#Hydrological process and parameter terms across tools|here]]. 
 +
</br>
  
−
''' PAGE IN PROGRESS - MORE COMING SOON! '''
 
−
 
−
building process representation overview tables + links to more specific pages about individual processes
 
−
</br>
 
  
 
== Canopy interception and evaporation (vs throughfall) ==
 
== Canopy interception and evaporation (vs throughfall) ==
Line 14: Line 21:
 
! scope="col" style="background: #CEE6F2; width:18%" |MIKE-SHE,<br/>[[Coverage of structural options within modelling tools#MIKE-SHE|semi-distributed &<br/> fully distributed]]
 
! scope="col" style="background: #CEE6F2; width:18%" |MIKE-SHE,<br/>[[Coverage of structural options within modelling tools#MIKE-SHE|semi-distributed &<br/> fully distributed]]
 
|-
 
|-
−
| style="vertical-align: top" |<br/>'''Algorithm <br/>inputs''' <br/>''(input data,<br/>parameters,<br/>model-calculated<br/>states/storages)''
+
| style="vertical-align: top" |<br/>'''Algorithm <br/>inputs''' <br/>''<small>(input data,<br/>model calculated<br/>states / storages,<br/>property parameters)</small>''
 
| style="background: #FFF5FA; vertical-align: top" |
 
| style="background: #FFF5FA; vertical-align: top" |
 
<br/>
 
<br/>
Line 72: Line 79:
 
timestep interception capacity ''(calculated)''  
 
timestep interception capacity ''(calculated)''  
 
|}
 
|}
 +
[[#Process rep - pg top|Back to top of page]]
 +
</br>
 
</br>
 
</br>
  
Line 85: Line 94:
 
! scope="col" style="background: #CEDAF2; width:15%" |MIKE-SHE,</br>[[Coverage of structural options within modelling tools#MIKE-SHE|fully-distributed, </br>more physical]]
 
! scope="col" style="background: #CEDAF2; width:15%" |MIKE-SHE,</br>[[Coverage of structural options within modelling tools#MIKE-SHE|fully-distributed, </br>more physical]]
 
|-
 
|-
−
| style="vertical-align: top" |<br/>'''Algorithm <br/>inputs''' <br/>''(input data,<br/>parameters,<br/>model-calculated<br/>states/storages)''
+
| style="vertical-align: top" |<br/>'''Algorithm <br/>inputs''' <br/>''<small>(input data,<br/>model calculated<br/>states / storages,<br/>property parameters)</small>''
 
| style="background: #FFF5FA; vertical-align: top" |
 
| style="background: #FFF5FA; vertical-align: top" |
 
<br/>
 
<br/>
Line 178: Line 187:
 
<br/>Saturation soil moisture  
 
<br/>Saturation soil moisture  
 
|}
 
|}
 +
[[#Process rep - pg top|Back to top of page]]
 +
</br>
 
</br>
 
</br>
  
 
== Surface runoff to channel network ==
 
== Surface runoff to channel network ==
−
'''Note on surface runoff vs ponding and surface storage:''' Water reaching the land surface which does not infiltrate into soil in a timestep, will either become surface runoff or remain as surface ponding/surface storage. In some models (e.g. Pitman models) all water that does not infiltrate becomes surface runoff. This makes sense when the timestep is long (e.g. a month). If the model represents surface storage, then in subsequent model timesteps this water can evaporate, infiltrate, and/or become surface runoff in the next timestep. Some models have relatively short timesteps (subdaily, daily) compared to the rate at which water would move across a modelled land unit (which could be large, rough, and/or flat, slowing the flow rate), explaining why some water will be considered 'surface storage' in one timestep and then 'surface runoff' in the next timestep in the model. MIKE-SHE can also consider that land surfaces can be very rough, or have dips that trap water, and so some amount surface water will not be able to runoff at all, regardless of a time delay. This water will both evaporate and infiltrate over time.   
+
'''Note on surface runoff vs ponding and surface storage:''' In the Pitman models, all water not calculated to infiltrate in a timestep, which is a month, is assumed to be runoff. For models with shorter timesteps, water reaching the land surface which does not infiltrate into soil in a timestep can become surface runoff or remain as surface ponding/surface storage. In subsequent model timesteps, water still on the land surface can evaporate, infiltrate, and/or become surface runoff in the next timestep. Some models have relatively short timesteps (subdaily, daily) compared to the rate at which water would move all the way across a modelled land unit (which could be large, rough, and/or flat, slowing the flow rate). This is why some water will be 'surface storage' in one timestep and then 'surface runoff' in the next timestep. MIKE-SHE can also consider that land surfaces can be very rough, or have dips that trap water, and so some amount surface water will not be able to run off at all (detention storage). This water will both evaporate and infiltrate over time.
 +
</br> ACRU4 differs notably, assuming the water that does not infiltrate into soil accounts for both surface runoff and interflow runoff.   
 
{| class="wikitable"
 
{| class="wikitable"
 
! scope="col" style="width:10%" | Algorithm <br/>description
 
! scope="col" style="width:10%" | Algorithm <br/>description
Line 191: Line 203:
 
! scope="col" style="background: #CEDAF2; width:15%" |MIKE-SHE,</br>[[Coverage of structural options within modelling tools#MIKE-SHE|fully-distributed, </br>more physical]]
 
! scope="col" style="background: #CEDAF2; width:15%" |MIKE-SHE,</br>[[Coverage of structural options within modelling tools#MIKE-SHE|fully-distributed, </br>more physical]]
 
|-
 
|-
−
| style="vertical-align: top" |<br/>'''Algorithm <br/>inputs''' <br/>''(input data,<br/>parameters,<br/>model-calculated<br/>states/storages)''
+
| style="vertical-align: top" |<br/>'''Algorithm <br/>inputs''' <br/>''<small>(input data,<br/>model calculated<br/>states / storages,<br/>property parameters)</small>''
 
| style="background: #FFF5FA; vertical-align: top" |
 
| style="background: #FFF5FA; vertical-align: top" |
 
<br/>
 
<br/>
 
* Surface water present* ''(state)''
 
* Surface water present* ''(state)''
 +
<br/>
 
<small> *Month timestep: all water not infiltrating becomes runoff </small>  
 
<small> *Month timestep: all water not infiltrating becomes runoff </small>  
 
<br/>
 
<br/>
Line 200: Line 213:
 
<br/>
 
<br/>
 
* Surface water present* ''(state)''
 
* Surface water present* ''(state)''
 +
<br/>
 
<small> *Month timestep: all water not infiltrating becomes runoff </small>  
 
<small> *Month timestep: all water not infiltrating becomes runoff </small>  
 
<br/>
 
<br/>
Line 206: Line 220:
 
* Surface water present* ''(state)''
 
* Surface water present* ''(state)''
 
* Lag coefficient*
 
* Lag coefficient*
 +
<br/>
 
<small> *Lag coefficient separates "quickflow"~surface runoff from "delayed-flow"~interflow. Quickflow portion reaches the channel on the same day generated (rain day). The rest is lagged over subsequent days </small>  
 
<small> *Lag coefficient separates "quickflow"~surface runoff from "delayed-flow"~interflow. Quickflow portion reaches the channel on the same day generated (rain day). The rest is lagged over subsequent days </small>  
 
| style="background: #FFFFF5; vertical-align: top"|
 
| style="background: #FFFFF5; vertical-align: top"|
Line 212: Line 227:
 
* '''''Surface path properties:'''''
 
* '''''Surface path properties:'''''
 
** Path length
 
** Path length
−
** Path slope
+
** Slope
−
** Path roughness (Manning's n)
+
** Roughness (Manning's n)
 
<br/>
 
<br/>
 
| style="background: #F5FCFF; vertical-align: top" |
 
| style="background: #F5FCFF; vertical-align: top" |
Line 221: Line 236:
 
** Detention storage
 
** Detention storage
 
** Path length
 
** Path length
−
** Path slope
+
** Slope
−
** Path roughness (Manning's n)
+
** Roughness (Manning's n)
 
<br/>
 
<br/>
 
| style="background: #F5F8FF; vertical-align: top" |
 
| style="background: #F5F8FF; vertical-align: top" |
Line 230: Line 245:
 
** Detention storage
 
** Detention storage
 
** Path length (from gridded topography)
 
** Path length (from gridded topography)
−
** Path slope(from gridded topography)
+
** Slope(from gridded topography)
−
** Path roughness (Manning's n)
+
** Roughness (Manning's n)
 
<br/>
 
<br/>
 
|-
 
|-
Line 262: Line 277:
 
<br/>Detention storage
 
<br/>Detention storage
 
|}
 
|}
 +
[[#Process rep - pg top|Back to top of page]]
 +
</br>
 
</br>
 
</br>
  
Line 275: Line 292:
 
! scope="col" style="background: #CEDAF2; width:15%" |MIKE-SHE,</br>[[Coverage of structural options within modelling tools#MIKE-SHE|fully-distributed, </br>more physical]]
 
! scope="col" style="background: #CEDAF2; width:15%" |MIKE-SHE,</br>[[Coverage of structural options within modelling tools#MIKE-SHE|fully-distributed, </br>more physical]]
 
|-
 
|-
−
| style="vertical-align: top" |<br/>'''Algorithm <br/>inputs''' <br/>''(input data,<br/>parameters,<br/>model-calculated<br/>states/storages)''
+
| style="vertical-align: top" |<br/>'''Algorithm <br/>inputs''' <br/>''<small>(input data,<br/>model calculated<br/>states / storages,<br/>property parameters)</small>''
 
| style="background: #FFF5FA; vertical-align: top" |
 
| style="background: #FFF5FA; vertical-align: top" |
 
<br/>
 
<br/>
Line 373: Line 390:
 
<br/>Wilting point SM
 
<br/>Wilting point SM
 
| style="background: #FFFFF5; vertical-align: top"|
 
| style="background: #FFFFF5; vertical-align: top"|
−
'''yes:'''
 
 
'''yes:'''
 
'''yes:'''
 
<br/>ET demand,
 
<br/>ET demand,
Line 386: Line 402:
 
<br/>Wilting point SM  
 
<br/>Wilting point SM  
 
|}
 
|}
 +
[[#Process rep - pg top|Back to top of page]]
 +
</br>
 
</br>
 
</br>
  
 
== Evapotranspiration (ET) from groundwater (GW) ==
 
== Evapotranspiration (ET) from groundwater (GW) ==
−
'''Note: More coverage of evapotranspiration related terminology and inputs across different tools can be found [[Terminology#Evapotranspiration terms across tools|here]]'''
+
'''Note:''' More coverage of evapotranspiration related terminology and inputs across different tools can be found '''[[Terminology#Evapotranspiration terms across tools|here]]'''
 
{| class="wikitable"
 
{| class="wikitable"
 
! scope="col" style="width:10%" | Algorithm <br/>description
 
! scope="col" style="width:10%" | Algorithm <br/>description
Line 399: Line 417:
 
! scope="col" style="background: #CEDAF2; width:15%" |MIKE-SHE,</br>[[Coverage of structural options within modelling tools#MIKE-SHE|fully-distributed, </br>more physical]]
 
! scope="col" style="background: #CEDAF2; width:15%" |MIKE-SHE,</br>[[Coverage of structural options within modelling tools#MIKE-SHE|fully-distributed, </br>more physical]]
 
|-
 
|-
−
| style="vertical-align: top" |<br/>'''Algorithm <br/>inputs''' <br/>''(input data,<br/>parameters,<br/>model-calculated<br/>states/storages)''
+
| style="vertical-align: top" |<br/>'''Algorithm <br/>inputs''' <br/>''<small>(input data,<br/>model calculated<br/>states / storages,<br/>property parameters)</small>''
 
| style="background: #FFF5FA; vertical-align: top" |
 
| style="background: #FFF5FA; vertical-align: top" |
 
<br/>
 
<br/>
Line 409: Line 427:
 
** "Crop" coefficient (vs pan)
 
** "Crop" coefficient (vs pan)
 
* '''''Aquifer properties:'''''
 
* '''''Aquifer properties:'''''
−
** Aquifer storage limit for ET
+
** Storage limit for ET
 
<br/>
 
<br/>
 
| style="background: #FFF7F5; vertical-align: top" |
 
| style="background: #FFF7F5; vertical-align: top" |
Line 421: Line 439:
 
** Pitman ET coefficient*  
 
** Pitman ET coefficient*  
 
* '''''Aquifer properties:'''''
 
* '''''Aquifer properties:'''''
−
** Aquifer storage limit for ET
+
** Storage limit for ET
 
<br/>
 
<br/>
 
| style="background: #F5FFF5; vertical-align: top" ; rowspan= "3" |
 
| style="background: #F5FFF5; vertical-align: top" ; rowspan= "3" |
Line 438: Line 456:
 
** Capillary rise rate coefficient
 
** Capillary rise rate coefficient
 
* '''''Aquifer properties:'''''
 
* '''''Aquifer properties:'''''
−
** Aquifer storage limit for capillary rise
+
** Storage limit for capillary rise
 
''+ all inputs for ET from soil''
 
''+ all inputs for ET from soil''
 
<br/>
 
<br/>
Line 448: Line 466:
 
* Aquifer storage ''(state)''  
 
* Aquifer storage ''(state)''  
 
* '''''Aquifer properties:'''''
 
* '''''Aquifer properties:'''''
−
** Aquifer storage limit for capillary rise
+
** Storage limit for capillary rise
 
''+ all inputs for [[#Evapotranspiration (ET) from soil moisture (SM)|ET from soil]]''
 
''+ all inputs for [[#Evapotranspiration (ET) from soil moisture (SM)|ET from soil]]''
 
<br/>
 
<br/>
Line 490: Line 508:
 
<br/>Water table depth vs root depth
 
<br/>Water table depth vs root depth
 
|}
 
|}
 +
[[#Process rep - pg top|Back to top of page]]
 +
</br>
 
</br>
 
</br>
  
 
== Interflow generation & routing to channel network ==
 
== Interflow generation & routing to channel network ==
 +
{| class="wikitable"
 +
! scope="col" style="width:10%" | Algorithm <br/>description
 +
! scope="col" style="background: #F2CEE0; width:15%" |WRSM-Pitman
 +
! scope="col" style="background: #F2D4CE; width:15%" |SPATSIM-Pitman
 +
! scope="col" style="background: #CEF2CE; width:15%" |ACRU4
 +
! scope="col" style="background: #F2F2CE; width:15%" |SWAT2012
 +
! scope="col" style="background: #CEE6F2; width:15%" |MIKE-SHE,</br>[[Coverage of structural options within modelling tools#MIKE-SHE|semi-distributed, </br>more conceptual]]
 +
! scope="col" style="background: #CEDAF2; width:15%" |MIKE-SHE,</br>[[Coverage of structural options within modelling tools#MIKE-SHE|fully-distributed, </br>more physical]]
 +
|-
 +
| style="vertical-align: top" |<br/>'''Algorithm <br/>inputs''' <br/>''<small>(input data,<br/>model calculated<br/>states / storages,<br/>property parameters)</small>''
 +
| style="background: #FFF5FA; vertical-align: top" |
 +
<br/>
 +
''Interflow from soil layer & from "percolation zone"''
 +
* Soil moisture ''(state)''
 +
* Percolation zone storage ''(state)''
 +
* Aquifer storage ''(state)''
 +
* '''''Soil properties:'''''
 +
** Saturation soil moisture
 +
** Field capacity SM
 +
** Interflow rate power coefficient
 +
** Max interflow rate
 +
** Interflow lag coefficient
 +
* '''''Percolation zone properties'''''
 +
** Max storage capacity
 +
* '''''Aquifer properties'''''
 +
** Max storage capacity
 +
<br/>
 +
| style="background: #FFF7F5; vertical-align: top" |
 +
<br/>
 +
''Interflow from soil layer''
 +
* Soil moisture ''(state)''
 +
* '''''Soil properties:'''''
 +
** Saturation soil moisture
 +
** Field capacity SM
 +
** Interflow rate power coefficient
 +
** Max interflow rate
 +
** Interflow lag coefficient
 +
<br/>
 +
| style="background: #F5FFF5; vertical-align: top"; rowspan= "3" |
 +
<br/>
 +
'''Interflow amount calculated in the same step as surface runoff, not as a separate process:'''
 +
</br>''total surface runoff + interflow runoff is calculated in the [[#Infiltration into soil moisture (vs surface runoff or surface ponding)|infiltration]] vs [[#Surface runoff to channel network|runoff]] step, then these are separated through lagging''
 +
* Surface water present* ''(state)''
 +
* Lag coefficient*
 +
<small> *Lag coefficient separates "quickflow"~surface runoff from "delayed-flow"~interflow. The lagged flow that doesn't reach the channel on the rain day, but instead reaches the channel over subsequent days can be considered interflow </small>
 +
<br/>
 +
| style="background: #FFFFF5; vertical-align: top"|
 +
<br/>
 +
''Interflow from soil layer''
 +
* Soil moisture ''(state)''
 +
* '''''Soil properties:'''''
 +
** Saturation soil moisture
 +
** Field capacity SM
 +
** Drainage slope
 +
** Saturated hydraulic conductivity (K<small>sat</small>)
 +
<br/>
 +
| style="background: #F5FCFF; vertical-align: top" |
 +
<br/>
 +
''Interflow from "interflow zone" (IZ)''
 +
</br>''to calc input into IZ:''
 +
* Macropore flow (Throughfall * by-pass coefficient)
 +
* Soil moisture ''(state)''
 +
* '''''Soil properties:'''''
 +
** Saturation soil moisture
 +
** Field capacity SM
 +
** Drainage slope
 +
** Saturated hydraulic conductivity (K<small>sat</small>)
 +
''to calc outflow from IZ:''
 +
* IZ storage ''(state)''
 +
* '''''IZ properties:'''''
 +
** Storage threshold for outflow
 +
** Lateral outflow rate constant
 +
<br/>
 +
| style="background: #F5F8FF; vertical-align: top"; rowspan= "3"  |
 +
<br/>
 +
'''Interflow not modelled as a separate process to "saturated zone" flow.'''
 +
</br> ''Lateral flow through a "saturated zone" material layer in the model which is not frequently saturated and/or is perched above more permanently saturated layers can be considered as interflow. Flows through different layers can be exported as separate model outputs.''
 +
</br></br>
 +
''See  [[#Aquifer exchanges with channel: aquifer outflow to channel and channel transmission loss|groundwater outflow to channel]]''
 +
<br/>
 +
|-
 +
| style="vertical-align: top" |'''Function <br/>type'''
 +
| style="background: #FFF5FA; vertical-align: top" |non-linear & threshold
 +
| style="background: #FFF7F5; vertical-align: top" |non-linear & threshold
 +
| style="background: #FFFFF5; vertical-align: top" |non-linear & threshold
 +
| style="background: #F5FCFF; vertical-align: top" |two step: non-linear & threshold, linear reservoir & threshold
 +
|-
 +
| style="vertical-align: top" |'''Thresholds'''
 +
| style="background: #FFF5FA; vertical-align: top" |
 +
'''yes:'''
 +
<br/>Field capacity SM,
 +
<br/>Percolation zone max storage + Aquifer max storage
 +
| style="background: #FFF7F5; vertical-align: top" |
 +
'''yes:'''
 +
<br/>Field capacity SM
 +
| style="background: #FFFFF5; vertical-align: top" |
 +
'''yes:'''
 +
<br/>Field capacity SM
 +
| style="background: #F5FCFF; vertical-align: top" |
 +
'''yes:'''
 +
<br/>Field capacity SM
 +
<br/>IZ storage threshold for outflow
 +
|}
 +
[[#Process rep - pg top|Back to top of page]]
 +
</br>
 +
</br>
  
 
== Aquifer recharge ==
 
== Aquifer recharge ==
 +
{| class="wikitable"
 +
! scope="col" style="width:10%" | Algorithm <br/>description
 +
! scope="col" style="background: #F2CEE0; width:15%" |WRSM-Pitman</br>(Sami GW)
 +
! scope="col" style="background: #F2D4CE; width:15%" |SPATSIM-Pitman</br>(Hughes GW)
 +
! scope="col" style="background: #CEF2CE; width:15%" |ACRU4
 +
! scope="col" style="background: #F2F2CE; width:15%" |SWAT2012
 +
! scope="col" style="background: #CEE6F2; width:15%" |MIKE-SHE,</br>[[Coverage of structural options within modelling tools#MIKE-SHE|semi-distributed, </br>more conceptual]]
 +
! scope="col" style="background: #CEDAF2; width:15%" |MIKE-SHE,</br>[[Coverage of structural options within modelling tools#MIKE-SHE|fully-distributed, </br>more physical]]
 +
|-
 +
| style="vertical-align: top" |<br/>'''Algorithm <br/>inputs''' <br/>''<small>(input data,<br/>model calculated<br/>states / storages,<br/>property parameters)</small>''
 +
| style="background: #FFF5FA; vertical-align: top" |
 +
<br/>
 +
''Recharge via "percolation zone" (vadose zone)''
 +
* Soil moisture ''(state)''
 +
* Vadose zone storage ''(state)''
 +
* Aquifer storage ''(state)''
 +
* '''''Soil properties:'''''
 +
** Saturation soil moisture
 +
** Field capacity SM
 +
** Percolation power coefficient
 +
** Max percolation rate
 +
*'''''Percolation zone properties:'''''
 +
** Recharge lag coefficient
 +
** Max storage capacity
 +
*'''''Aquifer properties:'''''
 +
** Max storage capacity
 +
<br/>
 +
| style="background: #FFF7F5; vertical-align: top" |
 +
<br/>
 +
''Recharge from soil column''
 +
* Soil moisture ''(state)''
 +
* Aquifer storage ''(state)''
 +
* '''''Soil properties:'''''
 +
** Saturation soil moisture
 +
** Field capacity SM
 +
** Percolation power coefficient
 +
** Max percolation rate
 +
*'''''Aquifer properties:'''''
 +
** Max storage capacity
 +
<br/>
 +
| style="background: #F5FFF5; vertical-align: top" |
 +
<br/>
 +
''Recharge from soil column*''
 +
* Soil moisture ''(state)''
 +
* '''''Soil properties:'''''
 +
** Saturation soil moisture
 +
** Field capacity SM
 +
** Percolation rate (~K<small>sat</small>)
 +
<br/>
 +
<small>*Aquifer storage is unlimited</small>
 +
<br/>
 +
| style="background: #FFFFF5; vertical-align: top"|
 +
<br/>
 +
''Recharge via vadose zone''
 +
* Soil moisture ''(state)''
 +
* Vadose zone storage ''(state)''
 +
* Aquifer storage ''(state)''
 +
* '''''Soil properties:'''''
 +
** Saturation soil moisture
 +
** Field capacity SM
 +
** Saturated hydraulic conductivity (K<small>sat</small>)
 +
*'''''Vadose zone (VZ) properties:'''''
 +
** Field capacity
 +
** Saturated hydraulic conductivity (K<small>sat</small>)
 +
*'''''Aquifer properties:'''''
 +
** Max storage capacity
 +
<br/>
 +
| style="background: #F5FCFF; vertical-align: top" |
 +
<br/>
 +
''Recharge via "interflow zone"''
 +
* Soil moisture ''(state)''
 +
* Interflow zone storage ''(state)''
 +
* Aquifer storage ''(state)''
 +
* '''''Soil properties:'''''
 +
** Saturation soil moisture
 +
** Field capacity SM
 +
** Saturated hydraulic conductivity (K<small>sat</small>)
 +
*'''''Interflow zone (IZ) properties:'''''
 +
**Storage threshold for outflow
 +
**Vertical outflow rate constant
 +
*'''''Aquifer properties:'''''
 +
** Max storage capacity
 +
<br/>
 +
| style="background: #F5F8FF; vertical-align: top" |
 +
<br/>
 +
''Recharge from soil column''
 +
* Soil moisture ''(state)''
 +
* Aquifer storage ''(state)''
 +
* '''''Soil properties:'''''
 +
** Saturation soil moisture
 +
** Moisture retention curve
 +
** Saturated hydraulic conductivity (K<small>sat</small>)
 +
*'''''Aquifer properties:'''''
 +
** Max storage capacity
 +
** Conductivity (K<small>sat</small>)
 +
<br/>
 +
|-
 +
| style="vertical-align: top" |'''Function <br/>type'''
 +
| style="background: #FFF5FA; vertical-align: top" |two step: non-linear & threshold, non-linear
 +
| style="background: #FFF7F5; vertical-align: top" |non-linear & threshold
 +
| style="background: #F5FFF5; vertical-align: top" |non-linear & threshold
 +
| style="background: #FFFFF5; vertical-align: top" |two step: non-linear & threshold, non-linear & threshold
 +
| style="background: #F5FCFF; vertical-align: top" |two step: non-linear & threshold, linear reservoir & threshold
 +
| style="background: #F5F8FF; vertical-align: top" |non-linear & threshold
 +
|-
 +
| style="vertical-align: top" |'''Thresholds'''
 +
| style="background: #FFF5FA; vertical-align: top" |
 +
'''yes:'''
 +
<br/>Field capacity SM,
 +
<br/>Vadose zone max storage
 +
<br/>Aquifer max storage
 +
| style="background: #FFF7F5; vertical-align: top" |
 +
'''yes:'''
 +
<br/>Field capacity SM,
 +
<br/>Aquifer max storage
 +
| style="background: #F5FFF5; vertical-align: top" |
 +
'''yes:'''
 +
<br/>Field capacity SM
 +
<br/><br/>
 +
<small>*Aquifer storage is unlimited</small>
 +
| style="background: #FFFFF5; vertical-align: top"|
 +
'''yes:'''
 +
<br/>Field capacity SM,
 +
<br/>Vadose zone field capacity
 +
<br/>Aquifer max storage
 +
| style="background: #F5FCFF; vertical-align: top" |
 +
'''yes:'''
 +
<br/>Field capacity SM,
 +
<br/>Aquifer max storage 
 +
| style="background: #F5F8FF; vertical-align: top" |
 +
'''yes:'''
 +
<br/>Field capacity SM,
 +
<br/>Aquifer max storage
 +
|}
 +
[[#Process rep - pg top|Back to top of page]]
 +
</br>
 +
</br>
  
−
== Groundwater (GW) flow ==
+
== Groundwater (GW) flow: lateral flow '''within''' the saturated zone ==
 +
'''Note:''' This refers to groundwater flow between modelled land units, such as grid cells, HRUs, or subcatchments, depending on the scale that groundwater stores and flows are modelled at (see more on scales of process representation for GW [[Model units & connections#Aquifers: spatial distribution, vertical layers, connections|here]])
 +
</br> For coverage of groundwater outflow into a river channel, see the table [[#Aquifer exchanges with channel: aquifer outflow to channel and channel transmission loss|below]]
 +
{| class="wikitable"
 +
! scope="col" style="width:10%" | Algorithm <br/>description
 +
! scope="col" style="background: #F2CEE0; width:15%" |WRSM-Pitman</br>(Sami GW)
 +
! scope="col" style="background: #F2D4CE; width:15%" |SPATSIM-Pitman</br>(Hughes GW)
 +
! scope="col" style="background: #CEF2CE; width:15%" |ACRU4
 +
! scope="col" style="background: #F2F2CE; width:15%" |SWAT2012
 +
! scope="col" style="background: #CEE6F2; width:15%" |MIKE-SHE,</br>[[Coverage of structural options within modelling tools#MIKE-SHE|semi-distributed, </br>more conceptual]]
 +
! scope="col" style="background: #CEDAF2; width:15%" |MIKE-SHE,</br>[[Coverage of structural options within modelling tools#MIKE-SHE|fully-distributed, </br>more physical]]
 +
|-
 +
| style="vertical-align: top" |<br/>'''Algorithm <br/>inputs''' <br/>''<small>(input data,<br/>model calculated<br/>states / storages,<br/>property parameters)</small>''
 +
| style="background: #FFF5FA; vertical-align: top" |
 +
<br/>
 +
''GW flow between subcatchments''
 +
* Aquifer storage ''(state)''
 +
*'''''Aquifer properties:'''''
 +
** Transmissivity
 +
** Storage threshold for flow
 +
** Regional GW slope
 +
<br/>
 +
| style="background: #FFF7F5; vertical-align: top" |
 +
<br/>
 +
''GW flow between subcatchments''
 +
* Aquifer storage ''(state)''
 +
*'''''Aquifer properties:'''''
 +
** Transmissivity
 +
** Storage threshold for flow
 +
** Regional GW slope'''''
 +
** Max storage capacity
 +
<br/>
 +
| style="background: #F5FFF5; vertical-align: top"; rowspan= "3" |
 +
<br/>
 +
'''GW flow between HRU aquifers is not modelled'''
 +
<br/><br/>
 +
'''''Special case, partial exception:''' upland HRU aquifer outflow routed to riparian HRU soil''
 +
<br/>
 +
| style="background: #FFFFF5; vertical-align: top"; rowspan= "3"|
 +
<br/>
 +
'''GW flow between subcatchment aquifers is not modelled'''
 +
<br/>
 +
| style="background: #F5FCFF; vertical-align: top"; rowspan= "3" |
 +
<br/>
 +
'''GW flow between aquifer storage units is not modelled'''
 +
<br/>
 +
| style="background: #F5F8FF; vertical-align: top" |
 +
<br/>
 +
''GW flow between grid cells (Darcy's Law)''
 +
* GW level gradient ''(state - compare levels in neighboring cells)''
 +
*'''''Aquifer properties:'''''
 +
** Specific yield
 +
** Saturated conductivity (K<small>sat</small>)
 +
<br/>
 +
|-
 +
| style="vertical-align: top" |'''Function <br/>type'''
 +
| style="background: #FFF5FA; vertical-align: top" |non-linear & threshold
 +
| style="background: #FFF7F5; vertical-align: top" |non-linear & threshold
 +
| style="background: #F5F8FF; vertical-align: top" |non-linear
 +
|-
 +
| style="vertical-align: top" |'''Thresholds'''
 +
| style="background: #FFF5FA; vertical-align: top" |
 +
'''yes:'''
 +
<br/>Aquifer storage threshold for flow
 +
| style="background: #FFF7F5; vertical-align: top" |
 +
'''yes:'''
 +
<br/>Aquifer storage threshold for flow
 +
| style="background: #F5F8FF; vertical-align: top" |
 +
n/a
 +
|}
 +
[[#Process rep - pg top|Back to top of page]]
 +
</br>
 +
</br>
  
−
== Aquifer exchange with channel: aquifer outflow or channel transmission losses ==
+
== Aquifer exchanges with channel: aquifer outflow to channel and channel transmission loss ==
 +
{| class="wikitable"
 +
! scope="col" style="width:10%" | Algorithm <br/>description
 +
! scope="col" style="background: #F2CEE0; width:15%" |WRSM-Pitman</br>(Sami GW)
 +
! scope="col" style="background: #F2D4CE; width:15%" |SPATSIM-Pitman</br>(Hughes GW)
 +
! scope="col" style="background: #CEF2CE; width:15%" |ACRU4
 +
! scope="col" style="background: #F2F2CE; width:15%" |SWAT2012
 +
! scope="col" style="background: #CEE6F2; width:15%" |MIKE-SHE,</br>[[Coverage of structural options within modelling tools#MIKE-SHE|semi-distributed, </br>more conceptual]]
 +
! scope="col" style="background: #CEDAF2; width:15%" |MIKE-SHE,</br>[[Coverage of structural options within modelling tools#MIKE-SHE|fully-distributed, </br>more physical]]
 +
|-
 +
| style="vertical-align: top" |<br/>'''Algorithm <br/>inputs''' <br/>''<small>(input data,<br/>model calculated<br/>states / storages,<br/>property parameters)</small>''
 +
| style="background: #FFF5FA; vertical-align: top" |
 +
<br/>
 +
''Dynamic two-way exchange within subcatchment <small>("runoff module")</small>:
 +
<small> aquifer-to-channel or channel-to-aquifer, depending on conditions (calculated state)''</small>
 +
''Additional one-way channel-to-aquifer representation within "channel module"''
 +
* River channel level ''(input)''
 +
* Subcatchment runoff / area (estimation of local channel store) ''(state)''
 +
* Aquifer storage/level -> GW slope ''(state)''
 +
* '''''Aquifer properties:'''''
 +
** Max flow rate for exchange
 +
** Exchange flow power coefficient
 +
* '''''Channel (module) properties:'''''
 +
** Transmission loss rate
 +
<br/>
 +
| style="background: #FFF7F5; vertical-align: top" |
 +
<br/>
 +
''Dynamic two-way exchange: <small>aquifer-to-channel or channel-to-aquifer, depending on conditions (calculated state)</small>''
 +
* River channel level ''(input)''
 +
* Aquifer storage/level -> GW slope ''(state)''
 +
* '''''Aquifer properties:'''''
 +
** Specific yield
 +
** Conductivity
 +
* '''''Channel properties:'''''
 +
** Length of channel
 +
** Subcatchment drainage density
 +
<br/>
 +
| style="background: #F5FFF5; vertical-align: top" |
 +
<br/>
 +
''One-way exchange: <small>aquifer-to-channel</small>''
 +
</br>''(Channel transmission loss not modelled)''
 +
* Aquifer storage ''(state)''
 +
* '''''Aquifer properties:'''''
 +
* Outflow proportion
 +
<br/>
 +
| style="background: #FFFFF5; vertical-align: top"|
 +
<br/>
 +
''Limited two-way exchange: <small>aquifer-to-channel or channel-to-'''"bank storage"'''(separate storage to the subcatchment's aquifer). Bank storage can release water to channel later. </small>''
 +
* Aquifer storage ''(state)''
 +
* Channel flow ''(state)''
 +
*'''''Aquifer properties:'''''
 +
** Storage limit for outflow
 +
** Specific yield
 +
** Outflow recession constant
 +
* '''''Channel properties:'''''
 +
* Shape - wetted perimeter
 +
* Bed conductivity 
 +
<br/>
 +
| style="background: #F5FCFF; vertical-align: top" |
 +
<br/>
 +
''One-way exchange per channel reach unit (can have multiple per subcatchment with different exchange directions): <small> each reach can have either aquifer-to-channel OR channel-to-aquifer, not switch between over time </small>''
 +
* Aquifer storage ''(state)''
 +
*'''''Aquifer properties:'''''
 +
** Storage threshold for outflow
 +
** Specific yield
 +
** Outflow recession constant
 +
* '''''Channel properties:'''''
 +
* Losing vs gaining specification
 +
* Shape - wetted perimeter
 +
* Bed conductivity 
 +
<br/>
 +
| style="background: #F5F8FF; vertical-align: top" |
 +
<br/>
 +
''Dynamic two-way exchange: <small>aquifer-to-channel or channel-to-aquifer, depending on conditions (calculated state) for each riparian grid cell and channel node pair</small>''
 +
* Water table height ''(state)''
 +
* Channel water level ''(state)''
 +
*'''''Aquifer properties:'''''
 +
** Specific yield
 +
** Conductivity (K<small>sat</small>)
 +
* '''''Channel properties:'''''
 +
* Shape - wetted perimeter
 +
* Bed conductivity
 +
<br/>
 +
|-
 +
| style="vertical-align: top" |'''Function <br/>type'''
 +
| style="background: #FFF5FA; vertical-align: top" |non-linear & threshold
 +
| style="background: #FFF7F5; vertical-align: top" |non-linear & threshold
 +
| style="background: #F5FFF5; vertical-align: top" |rate
 +
| style="background: #FFFFF5; vertical-align: top" |non-linear & threshold
 +
| style="background: #F5FCFF; vertical-align: top" |linear reservoir & threshold
 +
| style="background: #F5F8FF; vertical-align: top" |non-linear & threshold
 +
|-
 +
| style="vertical-align: top" |'''Thresholds'''
 +
| style="background: #FFF5FA; vertical-align: top" |
 +
'''yes - to switch flow direction:'''
 +
<br/>Aquifer water level vs River water level
 +
| style="background: #FFF7F5; vertical-align: top" |
 +
'''yes - to switch flow direction:'''
 +
<br/>Aquifer water level vs River water level
 +
| style="background: #F5FFF5; vertical-align: top" |
 +
no
 +
| style="background: #FFFFF5; vertical-align: top"|
 +
'''yes - to switch flow direction:'''
 +
<br/>Aquifer storage threshold for outflow
 +
| style="background: #F5FCFF; vertical-align: top" |
 +
'''yes - for aquifer-to-channel:'''
 +
<br/>Aquifer storage threshold for outflow
 +
| style="background: #F5F8FF; vertical-align: top" |
 +
'''yes - to switch flow direction:'''
 +
<br/>Aquifer water level vs River water level
 +
|}
 +
[[#Process rep - pg top|Back to top of page]]
 +
</br>

Latest revision as of 14:37, 1 December 2023

The tables below summarise information about different process representation algorithms across the set of modelling tools. They cover the inputs used to model the occurrence and rate (amount per model calculation timestep) of the given hydrological process (e.g. infiltration of water into the soil) as well as some general characteristics of the equations, particularly noting if there are thresholds involved and what these are. The tables can be used to understand what is being considered in the calculation of each process and what inputs each modelling tool needs. Not all the processes listed are represented explicitly by all the tools. This is highlighted in the tables below and an overview is also presented in the capabilities overview table. The nature of the algorithms used and their inputs are linked to the spatial and temporal scale at which the process is being represented in the given tool, described further for each one here.

In these tables thresholds refer to limits that determine when a process would start to occur or would stop occurring. For example, the field capacity soil moisture of a soil is often used as the threshold soil moisture level for percolation of water downwards (to a lower soil layer or to recharge groundwater). If soil moisture is lower than this, no percolation will be calculated in the model. There may be multiple thresholds considered when modelling a process. For example, evapotranspiration of soil moisture may be assumed to stop once the atmospheric demand in the timestep has been met or once the soil moisture has reached wilting point level, even if demand has not yet been met.

The inputs to a process algorithm can include:

  • input data (e.g., rainfall, using the rainfall input for a particular modelled unit, for a given timestep),
  • property parameter values that the user inputs in the model set-up (e.g., soil moisture content at saturation for a soil layer in a modelled unit)
  • states or water storage levels that the model calculated internally for the timestep, so not directly input by the user (e.g., soil moisture content in a certain soil layer at a given timestep)

Different software tools may refer to equivalent inputs using different words and sometimes they require the user to input different specific property values, but end up calculating the same derived property. For example, on software tool may require porosity of soil to be input, another may require bulk density, and another may as for volumetric soil moisture at saturation, any of these, coupled with soil layer depth, can be used to calculate the maximum (or saturation) soil water storage volume. Where possible, similar terms have been used across tools in the tables below to highlight their conceptual similarities where these exist. Details about some of the contrasting terminology used in the interfaces of different tools can be found here.


Canopy interception and evaporation (vs throughfall)

Algorithm
description 		WRSM-Pitman SPATSIM-Pitman ACRU4 SWAT2012 MIKE-SHE,
semi-distributed &
fully distributed

Algorithm
inputs
(input data,
model calculated
states / storages,
property parameters)


  • Rain
  • Cover properties:
    • max interception


  • Rain
  • Cover properties:
    • max interception


  • Rain
  • PET
  • Canopy storage (state)
  • Cover properties:
    • max interception (OR LAI)


Canopy interception is not explicitly modelled when using daily timestep modelling.
It is implicitly considered in the ‘initial abstraction’, see infiltration & surface runoff below


  • Rain
  • PET
  • Canopy storage (state)
  • Cover properties:
    • LAI
    • canopy interception coefficient


Function
type 		exponential & threshold exponential & threshold threshold threshold
Thresholds

yes:
timestep interception capacity (calculated)

yes:
timestep interception capacity (calculated)

yes:
timestep interception capacity (calculated)

yes:
timestep interception capacity (calculated)

Back to top of page

Infiltration into soil moisture (vs surface runoff or surface ponding)

Note: This excludes the case of a modelling unit (HRU, grid cell, area within a subcatchment, etc.) is designated as impervious. In this case, rain inputs would become surface runoff, potentially with a portion staying behind as surface storage/ponding if the area has attenuation specified (to represent roughness and flatness).

Algorithm
description 		WRSM-Pitman SPATSIM-Pitman ACRU4 SWAT2012 MIKE-SHE,
semi-distributed,
more conceptual 		MIKE-SHE,
fully-distributed,
more physical

Algorithm
inputs
(input data,
model calculated
states / storages,
property parameters)


  • Rain
  • Rainfall distribution factor
  • Soil moisture (state)
  • Soil properties:
    • Saturation soil moisture
    • Infiltration rate distribution (min, max)


  • Rain
  • Rainfall distribution factor
  • Soil moisture (state)
  • Soil properties:
    • Saturation soil moisture
    • Infiltration rate distribution (min, max)


  • Throughfall (rain-interception)
  • Soil moisture within "infiltration depth"(state)
  • Soil properties:
    • Saturation soil moisture
    • "Infiltration depth" parameter
    • Macropore by-pass (clay cracking)


  • Throughfall (rain-interception)
  • Soil moisture (state)
  • Soil properties:
    • Saturation soil moisture
    • Curve number (SCS-CN)
    • Macropore by-pass


  • Throughfall (rain-interception)
  • Soil moisture (state)
  • Soil properties:
    • Saturation soil moisture
    • Saturated hydraulic conductivity (Ksat),
    • Macropore by-pass


  • Throughfall (rain-interception)
  • Soil moisture (state)
  • Soil properties:
    • Saturation soil moisture
    • Soil moisture retention curve
    • Saturated hydraulic conductivity (Ksat),
    • Macropore by-pass


Function
type 		non-linear & threshold non-linear & threshold power & threshold power & threshold linear (rate) & threshold non-linear & threshold
Thresholds

yes:
Max infiltration rate,
Saturation soil moisture*
*Saturation excess becomes interflow, not surface runoff

yes:
Max infiltration rate,
Saturation soil moisture*
*Saturation excess becomes interflow, not surface runoff

yes:
Saturation soil moisture*
*Fraction of saturation excess becomes "delayed flow"~interflow, not "quickflow"~surface runoff

yes:
Max infiltration rate (calculated),
Saturation soil moisture

yes:
Max infiltration rate (uses Ksat),
Saturation soil moisture

yes:
Max infiltration rate (calculated),
Saturation soil moisture

Back to top of page

Surface runoff to channel network

Note on surface runoff vs ponding and surface storage: In the Pitman models, all water not calculated to infiltrate in a timestep, which is a month, is assumed to be runoff. For models with shorter timesteps, water reaching the land surface which does not infiltrate into soil in a timestep can become surface runoff or remain as surface ponding/surface storage. In subsequent model timesteps, water still on the land surface can evaporate, infiltrate, and/or become surface runoff in the next timestep. Some models have relatively short timesteps (subdaily, daily) compared to the rate at which water would move all the way across a modelled land unit (which could be large, rough, and/or flat, slowing the flow rate). This is why some water will be 'surface storage' in one timestep and then 'surface runoff' in the next timestep. MIKE-SHE can also consider that land surfaces can be very rough, or have dips that trap water, and so some amount surface water will not be able to run off at all (detention storage). This water will both evaporate and infiltrate over time.
ACRU4 differs notably, assuming the water that does not infiltrate into soil accounts for both surface runoff and interflow runoff.

Algorithm
description 		WRSM-Pitman SPATSIM-Pitman ACRU4 SWAT2012 MIKE-SHE,
semi-distributed,
more conceptual 		MIKE-SHE,
fully-distributed,
more physical

Algorithm
inputs
(input data,
model calculated
states / storages,
property parameters)


  • Surface water present* (state)


*Month timestep: all water not infiltrating becomes runoff


  • Surface water present* (state)


*Month timestep: all water not infiltrating becomes runoff


  • Surface water present* (state)
  • Lag coefficient*


*Lag coefficient separates "quickflow"~surface runoff from "delayed-flow"~interflow. Quickflow portion reaches the channel on the same day generated (rain day). The rest is lagged over subsequent days


  • Surface water present* (state)
  • Surface path properties:
    • Path length
    • Slope
    • Roughness (Manning's n)


  • Surface water present* (state)
  • Surface path properties:
    • Detention storage
    • Path length
    • Slope
    • Roughness (Manning's n)


  • Surface water present* (state)
  • Surface path properties:
    • Detention storage
    • Path length (from gridded topography)
    • Slope(from gridded topography)
    • Roughness (Manning's n)


Function
type 		(no transformation) (no transformation) (no transformation) non-linear non-linear & threshold non-linear & threshold
Thresholds

no*
*none after infiltration has been calculated

no*
*none after infiltration has been calculated

no*
*none after infiltration has been calculated

no*
*none after infiltration has been calculated

yes:
Detention storage

yes:
Detention storage

Back to top of page

Evapotranspiration (ET) from soil moisture (SM)

Note: More coverage of evapotranspiration related terminology and inputs across different tools can be found here

Algorithm
description 		WRSM-Pitman SPATSIM-Pitman ACRU4 SWAT2012 MIKE-SHE,
semi-distributed,
more conceptual 		MIKE-SHE,
fully-distributed,
more physical

Algorithm
inputs
(input data,
model calculated
states / storages,
property parameters)


  • S-Pan evaporation
  • Soil moisture (state)
  • Vegetation properties:
    • "Crop" coefficient (vs pan)
    • Pitman ET coefficient*
  • Soil properties:
    • Saturation soil moisture

*determines linear decline of ET with soil moisture decline


  • S-Pan evaporation
  • Soil moisture (state)
  • Vegetation properties:
    • "Crop" coefficient (vs pan)
    • Pitman ET coefficient*
  • Soil properties:
    • Saturation soil moisture

*determines linear decline of ET with soil moisture decline


  • A-Pan evaporation
  • Soil moisture (state)
  • Vegetation properties:
    • "Crop" coefficient (vs pan)
    • Root depth distribution
  • Soil properties:
    • Saturation soil moisture*
    • Field capacity SM
    • Wilting point SM

*ET assumed to also decline if soil gets close to saturation, waterlogging, unless wetland plants


  • Remaining ET demand (PET - Canopy evap.)
  • Soil moisture (state)
  • Vegetation properties:
    • LAI
    • Root depth (max)
    • Root distribution
    • Demand redistribution coefficient
  • Soil properties:
    • Field capacity SM
    • Wilting point SM


  • Remaining ET demand (PET - Canopy evap. - Ponded evap.)
  • Soil moisture (state)
  • Vegetation properties:
    • "Crop" coefficient
    • ET curtailment SM
    • Root depth (max)
    • Root distribution
  • Soil properties:
    • Field capacity SM
    • Wilting point SM


  • Remaining ET demand (PET - Canopy evap. - Ponded evap.)
  • Soil moisture (state)
  • Vegetation properties:
    • "Crop" coefficient
    • ET vs SM curve parameter
    • Root depth (max)
    • Root distribution
    • ET demand depth distribution
  • Soil properties:
    • Soil moisture retention curve


Function
type 		linear & threshold linear & threshold multi-part linear & threshold non-linear & threshold multi-part linear & threshold non-linear & threshold
Thresholds

yes:
ET demand,
Wilting point SM

yes:
ET demand,
Wilting point SM

yes:
ET demand,
Wilting point SM

yes:
ET demand,
Wilting point SM

yes:
ET demand,
Wilting point SM

yes:
ET demand,
Wilting point SM

Back to top of page

Evapotranspiration (ET) from groundwater (GW)

Note: More coverage of evapotranspiration related terminology and inputs across different tools can be found here

Algorithm
description 		WRSM-Pitman
(Sami GW) 		SPATSIM-Pitman
(Hughes GW) 		ACRU4 SWAT2012 MIKE-SHE,
semi-distributed,
more conceptual 		MIKE-SHE,
fully-distributed,
more physical

Algorithm
inputs
(input data,
model calculated
states / storages,
property parameters)


riparian area only

  • A-Pan evaporation
  • Riparian area
  • Aquifer storage (state)
  • Vegetation properties:
    • "Crop" coefficient (vs pan)
  • Aquifer properties:
    • Storage limit for ET


riparian area only

  • Remaining ET demand (PET - ET from soil)
  • Riparian area
  • Aquifer storage (state)
  • Vegetation properties:
    • "Crop" coefficient (vs pan)
    • Pitman ET coefficient*
  • Aquifer properties:
    • Storage limit for ET


riparian area only

Indirect representation:
Specified riparian zone HRU can receive a user-specified fraction of the aquifer outflow from connected upslope HRUs. This water is added to riparian HRU soil in the B-horizon, where it can be accessed for ET and can recharge GW.

No ET (or capillary rise) is simulated from the GW store of the riparian HRU.


represented through capillary rise into the soil driven by ET demand

  • PET
  • Aquifer storage (state)
  • Soil properties:
    • Capillary rise rate coefficient
  • Aquifer properties:
    • Storage limit for capillary rise

+ all inputs for ET from soil


riparian area only
represented through capillary rise into the soil driven by ET demand

  • Remaining ET demand (PET - Canopy evap. - Ponded evap. - ET from soil)
  • Aquifer storage (state)
  • Aquifer properties:
    • Storage limit for capillary rise

+ all inputs for ET from soil


capillary rise is calculated, but roots can also extend below the water table and withdraw directly

  • Remaining ET demand (PET - Canopy evap. - Ponded evap. - ET from soil)
  • Water table depth (state)
  • Aquifer properties:
    • Saturated hydraulic conductivity (Ksat)

+ vegetation property inputs for ET from soil

Function
type 		non-linear & threshold non-linear & threshold non-linear & threshold non-linear & threshold non-linear & threshold
Thresholds

yes:
ET demand,
Aquifer storage limit for ET

yes:
ET demand,
Aquifer storage limit for ET

yes:
ET demand,
Aquifer storage limit for capillary rise

yes:
ET demand,
Aquifer storage limit for capillary rise

yes:
ET demand,
Water table depth vs root depth

Back to top of page

Interflow generation & routing to channel network

Algorithm
description 		WRSM-Pitman SPATSIM-Pitman ACRU4 SWAT2012 MIKE-SHE,
semi-distributed,
more conceptual 		MIKE-SHE,
fully-distributed,
more physical

Algorithm
inputs
(input data,
model calculated
states / storages,
property parameters)


Interflow from soil layer & from "percolation zone"

  • Soil moisture (state)
  • Percolation zone storage (state)
  • Aquifer storage (state)
  • Soil properties:
    • Saturation soil moisture
    • Field capacity SM
    • Interflow rate power coefficient
    • Max interflow rate
    • Interflow lag coefficient
  • Percolation zone properties
    • Max storage capacity
  • Aquifer properties
    • Max storage capacity


Interflow from soil layer

  • Soil moisture (state)
  • Soil properties:
    • Saturation soil moisture
    • Field capacity SM
    • Interflow rate power coefficient
    • Max interflow rate
    • Interflow lag coefficient


Interflow amount calculated in the same step as surface runoff, not as a separate process:
total surface runoff + interflow runoff is calculated in the infiltration vs runoff step, then these are separated through lagging

  • Surface water present* (state)
  • Lag coefficient*

*Lag coefficient separates "quickflow"~surface runoff from "delayed-flow"~interflow. The lagged flow that doesn't reach the channel on the rain day, but instead reaches the channel over subsequent days can be considered interflow


Interflow from soil layer

  • Soil moisture (state)
  • Soil properties:
    • Saturation soil moisture
    • Field capacity SM
    • Drainage slope
    • Saturated hydraulic conductivity (Ksat)


Interflow from "interflow zone" (IZ)
to calc input into IZ:

  • Macropore flow (Throughfall * by-pass coefficient)
  • Soil moisture (state)
  • Soil properties:
    • Saturation soil moisture
    • Field capacity SM
    • Drainage slope
    • Saturated hydraulic conductivity (Ksat)

to calc outflow from IZ:

  • IZ storage (state)
  • IZ properties:
    • Storage threshold for outflow
    • Lateral outflow rate constant


Interflow not modelled as a separate process to "saturated zone" flow.
Lateral flow through a "saturated zone" material layer in the model which is not frequently saturated and/or is perched above more permanently saturated layers can be considered as interflow. Flows through different layers can be exported as separate model outputs.

See groundwater outflow to channel

Function
type 		non-linear & threshold non-linear & threshold non-linear & threshold two step: non-linear & threshold, linear reservoir & threshold
Thresholds

yes:
Field capacity SM,
Percolation zone max storage + Aquifer max storage

yes:
Field capacity SM

yes:
Field capacity SM

yes:
Field capacity SM
IZ storage threshold for outflow

Back to top of page

Aquifer recharge

Algorithm
description 		WRSM-Pitman
(Sami GW) 		SPATSIM-Pitman
(Hughes GW) 		ACRU4 SWAT2012 MIKE-SHE,
semi-distributed,
more conceptual 		MIKE-SHE,
fully-distributed,
more physical

Algorithm
inputs
(input data,
model calculated
states / storages,
property parameters)


Recharge via "percolation zone" (vadose zone)

  • Soil moisture (state)
  • Vadose zone storage (state)
  • Aquifer storage (state)
  • Soil properties:
    • Saturation soil moisture
    • Field capacity SM
    • Percolation power coefficient
    • Max percolation rate
  • Percolation zone properties:
    • Recharge lag coefficient
    • Max storage capacity
  • Aquifer properties:
    • Max storage capacity


Recharge from soil column

  • Soil moisture (state)
  • Aquifer storage (state)
  • Soil properties:
    • Saturation soil moisture
    • Field capacity SM
    • Percolation power coefficient
    • Max percolation rate
  • Aquifer properties:
    • Max storage capacity


Recharge from soil column*

  • Soil moisture (state)
  • Soil properties:
    • Saturation soil moisture
    • Field capacity SM
    • Percolation rate (~Ksat)


*Aquifer storage is unlimited


Recharge via vadose zone

  • Soil moisture (state)
  • Vadose zone storage (state)
  • Aquifer storage (state)
  • Soil properties:
    • Saturation soil moisture
    • Field capacity SM
    • Saturated hydraulic conductivity (Ksat)
  • Vadose zone (VZ) properties:
    • Field capacity
    • Saturated hydraulic conductivity (Ksat)
  • Aquifer properties:
    • Max storage capacity


Recharge via "interflow zone"

  • Soil moisture (state)
  • Interflow zone storage (state)
  • Aquifer storage (state)
  • Soil properties:
    • Saturation soil moisture
    • Field capacity SM
    • Saturated hydraulic conductivity (Ksat)
  • Interflow zone (IZ) properties:
    • Storage threshold for outflow
    • Vertical outflow rate constant
  • Aquifer properties:
    • Max storage capacity


Recharge from soil column

  • Soil moisture (state)
  • Aquifer storage (state)
  • Soil properties:
    • Saturation soil moisture
    • Moisture retention curve
    • Saturated hydraulic conductivity (Ksat)
  • Aquifer properties:
    • Max storage capacity
    • Conductivity (Ksat)


Function
type 		two step: non-linear & threshold, non-linear non-linear & threshold non-linear & threshold two step: non-linear & threshold, non-linear & threshold two step: non-linear & threshold, linear reservoir & threshold non-linear & threshold
Thresholds

yes:
Field capacity SM,
Vadose zone max storage
Aquifer max storage

yes:
Field capacity SM,
Aquifer max storage

yes:
Field capacity SM

*Aquifer storage is unlimited

yes:
Field capacity SM,
Vadose zone field capacity
Aquifer max storage

yes:
Field capacity SM,
Aquifer max storage

yes:
Field capacity SM,
Aquifer max storage

Back to top of page

Groundwater (GW) flow: lateral flow within the saturated zone

Note: This refers to groundwater flow between modelled land units, such as grid cells, HRUs, or subcatchments, depending on the scale that groundwater stores and flows are modelled at (see more on scales of process representation for GW here)
For coverage of groundwater outflow into a river channel, see the table below

Algorithm
description 		WRSM-Pitman
(Sami GW) 		SPATSIM-Pitman
(Hughes GW) 		ACRU4 SWAT2012 MIKE-SHE,
semi-distributed,
more conceptual 		MIKE-SHE,
fully-distributed,
more physical

Algorithm
inputs
(input data,
model calculated
states / storages,
property parameters)


GW flow between subcatchments

  • Aquifer storage (state)
  • Aquifer properties:
    • Transmissivity
    • Storage threshold for flow
    • Regional GW slope


GW flow between subcatchments

  • Aquifer storage (state)
  • Aquifer properties:
    • Transmissivity
    • Storage threshold for flow
    • Regional GW slope
    • Max storage capacity


GW flow between HRU aquifers is not modelled

Special case, partial exception: upland HRU aquifer outflow routed to riparian HRU soil


GW flow between subcatchment aquifers is not modelled


GW flow between aquifer storage units is not modelled


GW flow between grid cells (Darcy's Law)

  • GW level gradient (state - compare levels in neighboring cells)
  • Aquifer properties:
    • Specific yield
    • Saturated conductivity (Ksat)


Function
type 		non-linear & threshold non-linear & threshold non-linear
Thresholds

yes:
Aquifer storage threshold for flow

yes:
Aquifer storage threshold for flow

n/a

Back to top of page

Aquifer exchanges with channel: aquifer outflow to channel and channel transmission loss

Algorithm
description 		WRSM-Pitman
(Sami GW) 		SPATSIM-Pitman
(Hughes GW) 		ACRU4 SWAT2012 MIKE-SHE,
semi-distributed,
more conceptual 		MIKE-SHE,
fully-distributed,
more physical

Algorithm
inputs
(input data,
model calculated
states / storages,
property parameters)


Dynamic two-way exchange within subcatchment ("runoff module"): aquifer-to-channel or channel-to-aquifer, depending on conditions (calculated state) Additional one-way channel-to-aquifer representation within "channel module"

  • River channel level (input)
  • Subcatchment runoff / area (estimation of local channel store) (state)
  • Aquifer storage/level -> GW slope (state)
  • Aquifer properties:
    • Max flow rate for exchange
    • Exchange flow power coefficient
  • Channel (module) properties:
    • Transmission loss rate


Dynamic two-way exchange: aquifer-to-channel or channel-to-aquifer, depending on conditions (calculated state)

  • River channel level (input)
  • Aquifer storage/level -> GW slope (state)
  • Aquifer properties:
    • Specific yield
    • Conductivity
  • Channel properties:
    • Length of channel
    • Subcatchment drainage density


One-way exchange: aquifer-to-channel
(Channel transmission loss not modelled)

  • Aquifer storage (state)
  • Aquifer properties:
  • Outflow proportion


Limited two-way exchange: aquifer-to-channel or channel-to-"bank storage"(separate storage to the subcatchment's aquifer). Bank storage can release water to channel later.

  • Aquifer storage (state)
  • Channel flow (state)
  • Aquifer properties:
    • Storage limit for outflow
    • Specific yield
    • Outflow recession constant
  • Channel properties:
  • Shape - wetted perimeter
  • Bed conductivity


One-way exchange per channel reach unit (can have multiple per subcatchment with different exchange directions): each reach can have either aquifer-to-channel OR channel-to-aquifer, not switch between over time

  • Aquifer storage (state)
  • Aquifer properties:
    • Storage threshold for outflow
    • Specific yield
    • Outflow recession constant
  • Channel properties:
  • Losing vs gaining specification
  • Shape - wetted perimeter
  • Bed conductivity


Dynamic two-way exchange: aquifer-to-channel or channel-to-aquifer, depending on conditions (calculated state) for each riparian grid cell and channel node pair

  • Water table height (state)
  • Channel water level (state)
  • Aquifer properties:
    • Specific yield
    • Conductivity (Ksat)
  • Channel properties:
  • Shape - wetted perimeter
  • Bed conductivity


Function
type 		non-linear & threshold non-linear & threshold rate non-linear & threshold linear reservoir & threshold non-linear & threshold
Thresholds

yes - to switch flow direction:
Aquifer water level vs River water level

yes - to switch flow direction:
Aquifer water level vs River water level

no

yes - to switch flow direction:
Aquifer storage threshold for outflow

yes - for aquifer-to-channel:
Aquifer storage threshold for outflow

yes - to switch flow direction:
Aquifer water level vs River water level

Back to top of page

Retrieved from "https://hydromodel-sa-wiki.saeon.ac.za/index.php?title=Process_representation_across_tools&oldid=776"