< HydroGeoSphere

NOTE: Moisture contents referred to below and in section 2.5.3 are actually input in grok as saturations, which are less prone to error as they always vary between zero and 1, while moisture content varies between zero and porosity.

Default evapotranspiration values

Unless you modify the default values, all zones (and elements) in the ET domain will be assigned the default properties listed in Table 5.19:

Table 5.19: Default Properties for Evapotranspiration
ParameterValueUnit
NameDefault ET (grass)-
Canopy storage parameter0.04m
Initial interception storage0.04m
Transpiration fitting parameters:--
C10.5-
C20.0-
C31.0-
Transpiration limiting saturations:--
Wilting point0.05-
Field capacity0.2-
Oxic limit0.6-
Anoxic limit0.8-
Evaporation limiting saturations:--
Minimum0.5-
Maximum0.1-
Leaf area index1.0-
Root zone depth0.2m
Evaporation depth0.2m

For each instruction we will again indicate its scope (i.e. .grok, .etprops). Recall that if an instruction is used in the prefix.grok file, it will affect the current set of chosen zones, while in a properties (e.g. .etprops) file, it will only affect the named material of which it is a part.

Canopy storage parameter

Scope: .etprops

  1. cint_et Canopy storage parameter [L], in Equation 2.73.


Initial interception storage

Scope: .etprops

  1. init_sint_et initial canopy interception storage value [L], in Equation 3.59.


Transpiration fitting parameters

Scope: .etprops

  1. C_1 Coefficient [dimensionless] in Equation 2.75.
  2. C_2 Coefficient [dimensionless] in Equation 2.75.
  3. C_3 Coefficient [dimensionless] in Equation 2.77. By default, this coefficient is set to 1.0, which gives a linear ramping function whereas higher values would give higher order ramping functions.


Transpiration limiting saturations

Scope: .etprops

  1. thwp_et Saturation [dimensionless] at wilting point, in Equation 2.77.
  2. thfc_et Saturation [dimensionless] at field capacity, in Equation 2.77.
  3. tho_et Saturation [dimensionless] at oxic limit, in Equation 2.77.
  4. than_et Saturation [dimensionless] at anoxic limit, in Equation 2.77.


Evaporation limiting saturations

Scope: .etprops

  1. the2_et Saturation [dimensionless] below which evaporation is zero, in Equation 2.77.
  2. the1_et Saturation [dimensionless] above which full evaporation can occur, in Equation 2.77.


The following instruction can be used to modify the default value (1.0 for all time) of the leaf area index LAI:

Lai tables...End

Scope: .etprops

  1. time(1), lai(1) First entry.
  2. time(2), lai(2) Second entry.
  ...etc.
  n. time(n), lai(n) nth entry.

Causes grok to begin reading a group of time vs leaf area index instructions until it encounters an End instruction.

Paired values of time T and leaf are index LAI should be entered from earliest to latest time. The number of entries in the list are counted automatically to determine the table size.

Observed values of leaf are index [Scurlock et al., 2001] and maximum rooting depth [Canadell et al., 1996] for various terrestrial biomes are shown in Table 5.20:

Table 5.20: Observed Values of Leaf Area Index [Scurlock et al., 2001] and Maximum Rooting Depth [Canadell et al., 1996] for Various Terrestrial Biomes
Vegetative coverLeaf area indexMaximum rooting depth
Desert1.319.5
Tundra, circumpolar and alpine2.690.5
Wetlands, temperate and tropical6.34-
Grasslands, temperate2.502.6
Grasslands, tropical2.5015.0
Crops, temperate and tropical4.222.1
Shrubland, heath or Mediterranean-type vegetation-5.2
Forest:--
boreal deciduous broadleaf2.642.0
boreal evergreen needleleaf3.502.0
boreal/temperate deciduous needleleaf4.632.0
temperate deciduous broadleaf5.122.9
temperate evergreen needleleaf6.703.9
temperate evergreen broadleaf5.82-
tropical deciduous broadleaf3.923.7
tropical evergreen broadleaf4.907.3
Plantations (managed forests)8.72-
Overall mean5.23-


Root depth

Scope: .etprops

  1. root_depth_et Maximum root depth [L].

A normalized root depth function is mapped onto porous media elements above this maximum depth. Currently, four root depth functions are available; constant, linear, quadratic and cubic, as shown in Figure 5.16. These functions are defined such that the area under the normalized function is 1.0:

Figure 5.16: Normalized Root Depth Functions




By default, the linear form of the depth function is used. The following instructions are available for using the other forms:

Rdf constant function
Rdf quadratic decay function
Rdf cubic decay function

Evaporation depth

Scope: .etprops

  1. evap_depth_et Evaporation depth [L].

Evaporation as a function of depth is treated in a similar fashion as the root zone depth described above.



By default, the linear form of the evaporation function is used. The following instructions are available for using the other forms:

Edf constant function
Edf quadratic decay function
Edf cubic decay function

Echo et at point

Scope: .grok

  1. x1, y1 xy-coordinate.

This instruction finds the column of nodes that the given coordinate falls within and then writes the pertinent evapotranspiration information to the o.eco file.



For example:

ET properties for element column

ZONE: 1
ET MATERIAL: dense deciduous forest
    Canopy storage parameter = 0.0
    Transpiration fitting parameter C1              = 0.3
    Transpiration fitting parameter C2              = 0.2
    Transpiration fitting parameter C3              = 3.0E-6
    Moisture content at wilting point               = 0.32
    Moisture content at field capacity              = 0.2
    Moisture content at oxic limit                  = 0.76
    Moisture content at anoxic limit                = 0.9
    Moisture content at at energy limiting maximum  = 0.32
    Moisture content at at energy limiting minimum  = 0.2
    Initial interception storage = 0.0

Tabular data:
     Time           LAI
    0.00000       2.08000
   0.100000E+42   2.08000

Maximum evaporative zone depth  =   2.0
     Depth          EDF
    0.00000       2.00000
    1.00000       0.00000

Maximum root zone depth         =   4.0
    Depth          RDF
    0.00000       2.00000
    1.00000       0.00000

    ELEMENT    ET FLAG      DEPTH        EDF            RDF
    12889         T        0.00000     0.281045      0.146305
    11528         T       0.303244     0.234783      0.134740
    10167         T       0.606488     0.188522      0.123174
    8806          T       0.909732     0.142261      0.111609
    7445          T        1.21298     0.959994E-01  0.100044
    6084          T        1.51622     0.497381E-01  0.884783E-01
    4723          T        1.81946     0.765174E-02  0.769130E-01
    3362          T        2.1227      0.00000       0.653477E-01
    2001          T        2.4259      0.00000       0.537824E-01
    640           T        2.7292      0.00000       0.422171E-01
!                                     ------------  ------------
Totals                                 1.00000       0.942610


In this case, the element column falls in ET zone 1, and the default linear root depth and evaporation depth functions have been used and mapped onto it. All of the evaporative potential has been distributed to the element column, as indicated by a total EDF value of 1.0. However, only a portion (0.942610 or 94%) of the total RDF of 1.0 has been mapped onto it, because the maximum root zone depth (4 m) exceeds the total depth of the element column at this point (2.7292 m).

This article is issued from Wikibooks. The text is licensed under Creative Commons - Attribution - Sharealike. Additional terms may apply for the media files.