Soil moisture calculation
All process-values which are calculated through soil moisture simulation are depicted in the following picture.
Hydrological Response Units (HRUs)
If Runoff is determined through soil moisture simulation, the HRU concept is applied. This results in a catchment being sub-sectioned into a various amount of hydrological homogenous areas.
The depicted assignment of land use and soil type is valid for each HRU. The resulting amount of water out of a HRU is applied at the outlet of the elemt. Therefore all HRUs discharge their Water with the same time delay, independent of their position within the catchment.
Soil moisture calculation is computationally intensive and therefore very time-consuming. This is especially true if there are many HRUs in each sub-catchment.
Providing information on land use is necessary to run the soil moisture calculation. E.g. the thickness of the rooted zone is determined by providing the root depth in the information regarding land use. Furthermore the parameters of land use are needed to calculate interception and transpiration. The parameters are:
- root depth
- soil coverage
- annual pattern of soil coverage
- leaf area index
- annual pattern of the leaf area index
These Haude-Coefficients scale the valid potential grass reference-evaporation for the time step by their ratio to the corresponding Haude-Coefficient for grass.
The soil moisture simulation is based on the non-linear calculation of the individual soil layers. The soil is divided into different soil layers. Each layer is calculated and (if existing) adjusted with the layers above or below. The parameters for the soil moisture calculation are the following physical soil properties:
- wilting point (
- field capacity (
- total pore volume (
- saturated hydraulic-conductivity (
- maximum infiltration capacity (
- maximum rate of capillary rise (
- Assignment to a soil type: sand, silt, clay
As few as one or as many as six soil layers are possible. A division into three layers has proven to yield the best results. That is why the entered layers are always distributed into three layers in BlueM internally.
- infiltration layer (standard thickness = 20 [cm])
- root layer (minimum thickness = 5 [cm])
- transport layer (minimum thickness = 5 [cm])
The new characteristic soil parameters for the internally used layers are calculated by weighing the original parameters of the original layers according to the original thickness of these layers.Saturated hydraulic conductivity is calculated according to conservation of continuity of the flow. For vertical Flow the velocity
v is to be constant within a layer for a given flow due to the continuity of the flow. Therefore the hydraulic gradient is no longer constant.
di= depth of each original layer [mm]
ki= saturated hydraulic conductivity of each original layer [mm/h]
kf,v= saturated hydraulic conductivity of the internally used layer [mm/h]
The Aggregation of the layers is depicted in Abbildung 38.
In BlueM the inflow to the interception reservoir (
QzuIC) is described as a linear function of the free interception reservoir.
kIC= 10.0 = Parameter to describe the fill rate of the interception reservoir [1/h] (Bug 409)
The maximum interception capacity for agricultural crops according to von Hoyningen-Huene (1983) is:
LAI= Leaf Area Index [-]
The rate of evaporation is assumed to be equal to the rate of potential evaporation (=
ETp). A simple differential equation can be derived due to the fact that inflow is equal to the amount of rain above the vegetation minus the rain which penetrates it :
Therefore the amount of penetrating rain can be determined through the subtraction of interception from rainfall onto the vegetation.
To simplify things it is assumed that this approach is valid for wooded areas, although this approach was derived out of studies for agricultural crops.
Furthermore the interception evaporation is calculated, which is subtracted from the potential evaporation which leads to a reduced potential evaporation remaining for the following processes of transpiration and evaporation.
Soil moisture calculation
The water balance equation for the soil layer is solved by using a stepwise linear approach for the processes effecting soil moisture, which are Infiltration, Evaporation, Transpiration, Percolation, Interflow and Capillary-Rise. The input value for evaporation and transpiration is determined through the reduced potential evaporation due to interception evaporation.
The equation to be solved is as follows:
θ(t)= current soil moisture
Inf(t)= infiltration into the soil
Evaakt(t)= actual evaporation
Transakt(t)= actual transpiration
Infiltration, Percolation, Evaporation, Transpiration, Interflow and Capillary-Rise are dependent on the current soil moisture. In the simulation this dependency is described by the following function curves.
kf= saturated hydraulic conductivity
nFK= available water capacity (
nFK = FK - WP)
WP= wilting point
FK= field capacity
GPV= total pore volume
I= gradient [-]
fPK= soil specific scaling factor for the percolation function
fEva= soil specific scaling factor for the evaporation function
fTP= soil specific scaling factor for the transpiration function
f1,Int= soil specific scaling factor for the interflow function
f2,Int= soil specific scaling factor for the interflow function
nInt= soil specific scaling factor for the interflow function
nTP= Krümmungsparameter der Transpirationsfunktion
Program parameters are calculated internally. The user only needs to supply the characteristic soil parameters
The simulation takes place with a newly developed Building-Block for the simulation of reservoirs, whichs process function are to be mapped through a stepwise linear approach.
The original approach according to Holtan (1961):
av= Vegetation parameter (between 0,1 and 1,0)
is used in a modified form in BlueM:
it needs to be noted that this function only describes the potential Infiltration. Actual Infiltration is limited to the available amount of rain. The modification is a results of the consideration that if it rains on a very dry soil the air within the soil cannot escape and thereby limits maximum infiltration. This leads to a stop of exponential increase in contrast to the original approach of Holtan. ("Flowerpot effect")
- → zu modifizieren in Ansatz nach van Genuchten
Interflow is realtivly independent of soil parameters and depends on soil moisture and slope of the individual HRU. (refer to Bug 28):
For the evaporation out of the soil the determined potential evaporation, which was adjusted to land use, is converted into evaporation for fallow land.
- Haude, W. (1954): Zur praktischen Bestimmung der aktuellen und potentiellen Evapotranspiration. – Mitteilungen des DWD, Bd. 8; Bad Kissingen
- Haude, W. (1955): Zur Bestimmung der Verdunstung auf möglichst einfache Weise. Mitteilungen des DWD, 2 (11), Bad Kissingen
- von Hoyningen-Huene, J. (1983): Die Interzeption des Niederschlages in landwirtschaftlichen Pflanzenbeständen. DVWK Schriften, Nr. 57, S. 1 - 53, PDF
- Holtan, H.N. (1961): A Concept for Infiltration Estimates in Watershed Engineering, U.S. Department of Agriculture, Agricultural Research Service, Bulletin 41-51, p. 25
- Ostrowski, M. (1992): Ein universeller Baustein zur Simulation hydrologischer Prozesse, Wasser und Boden, Heft 11 (PDF )
- Bear, J. (1988): Dynamics of fluids in porous media, American Elsevier Environmental Science Series (books.google.com)