HydroGeoSphere/Evapotranspiration

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:

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], $$c_{int}$$ in Equation 2.73.
 * &bull; &bull; &bull;

Initial interception storage
Scope: .etprops
 * 1) init_sint_et initial canopy interception storage value [L], $$S^0_{int}$$ in Equation 3.59.
 * &bull; &bull; &bull;

Transpiration fitting parameters
Scope: .etprops
 * 1) C_1 Coefficient $$C_1$$ [dimensionless] in Equation 2.75.
 * 2) C_2 Coefficient $$C_2$$ [dimensionless] in Equation 2.75.
 * 3) C_3 Coefficient $$C_3$$ [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.
 * &bull; &bull; &bull;

Transpiration limiting saturations
Scope: .etprops
 * 1) thwp_et Saturation [dimensionless] at wilting point, $${\theta}_{wp}$$ in Equation 2.77.
 * 2) thfc_et Saturation [dimensionless] at field capacity, $${\theta}_{fc}$$ in Equation 2.77.
 * 3) tho_et Saturation [dimensionless] at oxic limit, $${\theta}_o$$ in Equation 2.77.
 * 4) than_et Saturation [dimensionless] at anoxic limit, $${\theta}_{an}$$ in Equation 2.77.
 * &bull; &bull; &bull;

Evaporation limiting saturations
<tt>Scope: .etprops</tt>
 * 1) the2_et Saturation [dimensionless] below which evaporation is zero, $${\theta}_{e2}$$ in Equation 2.77.
 * 2) the1_et Saturation [dimensionless] above which full evaporation can occur, $${\theta}_{e1}$$ in Equation 2.77.
 * &bull; &bull; &bull;

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
<tt>Scope: .etprops</tt>
 * 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:


 * &bull; &bull; &bull;

Root depth
<tt>Scope: .etprops</tt>
 * 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:




 * &bull; &bull; &bull;

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
<tt>Scope: .etprops</tt>
 * 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.
 * &bull; &bull; &bull;

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
<tt>Scope: .grok</tt>
 * 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 <tt>o.eco</tt> file.
 * &bull; &bull; &bull;

For example:

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).