3.9 KiB
2.7.5. Lateral Sub-surface Runoff¶
Lateral sub-surface runoff occurs when saturated soil moisture conditions exist within the soil column. Sub-surface runoff is
(2.7.108)¶\[q_{drai} = \Theta_{ice} K_{baseflow} tan \left( \beta \right) \Delta z_{sat}^{N_{baseflow}} \ ,\]
where \(K_{baseflow}\) is a calibration parameter, \(\beta\) is the topographic slope, the exponent \(N_{baseflow}\) = 1, and \(\Delta z_{sat}\) is the thickness of the saturated portion of the soil column.
The saturated thickness is
(2.7.109)¶\[\Delta z_{sat} = z_{bedrock} - z_{\nabla},\]
where the water table \(z_{\nabla}\) is determined by finding the first soil layer above the bedrock depth (section 2.2.2.2) in which the volumetric water content drops below a specified threshold. The default threshold is set to 0.9.
The specific yield, \(S_{y}\), which depends on the soil properties and the water table location, is derived by taking the difference between two equilibrium soil moisture profiles whose water tables differ by an infinitesimal amount
(2.7.110)¶\[S_{y} =\theta_{sat} \left(1-\left(1+\frac{z_{\nabla } }{\Psi _{sat} } \right)^{\frac{-1}{B} } \right)\]
where B is the Clapp-Hornberger exponent. Because \(S_{y}\) is a function of the soil properties, it results in water table dynamics that are consistent with the soil water fluxes described in section 2.7.3.
After the above calculations, two numerical adjustments are implemented to keep the liquid water content of each soil layer (\(w_{liq,\, i}\) ) within physical constraints of \(w_{liq}^{\min } \le w_{liq,\, i} \le \left(\theta_{sat,\, i} -\theta_{ice,\, i} \right)\Delta z_{i}\) where \(w_{liq}^{\min } =0.01\) (mm). First, beginning with the bottom soil layer \(i=N_{levsoi}\), any excess liquid water in each soil layer (\(w_{liq,\, i}^{excess} =w_{liq,\, i} -\left(\theta_{sat,\, i} -\theta_{ice,\, i} \right)\Delta z_{i} \ge 0\)) is successively added to the layer above. Any excess liquid water that remains after saturating the entire soil column is added to drainage \(q_{drai}\). Second, to prevent negative \(w_{liq,\, i}\), each layer is successively brought up to \(w_{liq,\, i} =w_{liq}^{\min }\) by taking the required amount of water from the layer below. If this results in \(w_{liq,\, N_{levsoi} } <w_{liq}^{\min }\), then the layers above are searched in succession for the required amount of water (\(w_{liq}^{\min } -w_{liq,\, N_{levsoi} }\) ) and removed from those layers subject to the constraint \(w_{liq,\, i} \ge w_{liq}^{\min }\). If sufficient water is not found, then the water is removed from \(W_{t}\) and \(q_{drai}\).
The soil surface layer liquid water and ice contents are then updated for dew \(q_{sdew}\), frost \(q_{frost}\), or sublimation \(q_{subl}\) (section 2.5.4) as
(2.7.111)¶\[w_{liq,\, 1}^{n+1} =w_{liq,\, 1}^{n} +q_{sdew} \Delta t\]
(2.7.112)¶\[w_{ice,\, 1}^{n+1} =w_{ice,\, 1}^{n} +q_{frost} \Delta t\]
(2.7.113)¶\[w_{ice,\, 1}^{n+1} =w_{ice,\, 1}^{n} -q_{subl} \Delta t.\]
Sublimation of ice is limited to the amount of ice available.