Personal tools

Difference between revisions of "Module PorousMediaProperties"

From MohidWiki

Jump to: navigation, search
Line 7: Line 7:
 
<math>J=-K\left ( \theta  \right )\left ( \frac{\partial h}{\partial \theta } + 1 \right )</math>
 
<math>J=-K\left ( \theta  \right )\left ( \frac{\partial h}{\partial \theta } + 1 \right )</math>
 
===Evapotranspiration===
 
===Evapotranspiration===
Some water may disappear from the soil because of the evaporation and transpiration processes, which become a sink in soil water profile. These two processes, currently named Evapotranspiration may be modeled using the Penmann Monteith equation:
 
  
:<math> \overset{\text{Energy flux rate}}{\lambda_v E=\frac{\Delta R_n  +  \rho_a c_p  \left(  \delta q  \right) g_a }{\Delta  + \gamma \left (    1 + g_a / g_s    \right)}}~ \iff ~  \overset{\text{Volume flux rate}}{ET_o=\frac{\Delta R_n  +  \rho_a c_p  \left(  \delta q  \right) g_a } { \left(  \Delta  + \gamma \left (    1 + g_a / g_s    \right)    \right) \lambda_v }}
+
Some water may disappear from the soil because of the evaporation and transpiration processes, which become a sink in soil water profile. These two processes, currently named Evapotranspiration may be modeled using the Penmann Monteith equation.
 +
 
 +
<math> \overset{\text{Energy flux rate}}{\lambda_v E=\frac{\Delta R_n  +  \rho_a c_p  \left(  \delta q  \right) g_a }{\Delta  + \gamma \left (    1 + g_a / g_s    \right)}}~ \iff ~  \overset{\text{Volume flux rate}}{ET_o=\frac{\Delta R_n  +  \rho_a c_p  \left(  \delta q  \right) g_a } { \left(  \Delta  + \gamma \left (    1 + g_a / g_s    \right)    \right) \lambda_v }}
 
</math>
 
</math>
  
Line 27: Line 28:
 
:''γ'' = Psychrometric constant  (''&gamma;'' ≈ 66 Pa K<sup>-1</sup>)
 
:''γ'' = Psychrometric constant  (''&gamma;'' ≈ 66 Pa K<sup>-1</sup>)
  
::In Figure below it can be seen that the actual transpiration and evaporation are calculated in Porous Media module of MohidLand. The actual evaporation, which happens only at the soil surface, is calculated based on a pressure head limit imposed by the user. It allows the model to not evaporate any surface water, even if it is available, when the soil is very far way from the saturation point.
+
In Figure below it can be seen that the actual transpiration and evaporation are calculated in Porous Media module of MohidLand. The actual evaporation, which happens only at the soil surface, is calculated based on a pressure head limit imposed by the user. It allows the model to not evaporate any surface water, even if it is available, when the soil is very far way from the saturation point.
 
[[Image:evapotranspiration fluxogram.jpg|thumb|center|400px|Evapotranspiration fluxogram in Mohid Land model]]
 
[[Image:evapotranspiration fluxogram.jpg|thumb|center|400px|Evapotranspiration fluxogram in Mohid Land model]]
  

Revision as of 17:22, 7 February 2011

Overview

Main Processes

Water flow

The water movement in soil will be dependent on the pressure gradients existing in the soil profile and also according to gravity. The equation that describes this motion is the Buckingham Darcy equation (Jury et al,1991) where J is the water velocity, h is the pressure head, θ is the water content and K is the hydraulic conductivity.

J=-K\left ( \theta  \right )\left ( \frac{\partial h}{\partial \theta } + 1 \right )

Evapotranspiration

Some water may disappear from the soil because of the evaporation and transpiration processes, which become a sink in soil water profile. These two processes, currently named Evapotranspiration may be modeled using the Penmann Monteith equation.

 \overset{\text{Energy flux rate}}{\lambda_v E=\frac{\Delta R_n   +   \rho_a c_p  \left(  \delta q  \right) g_a }{\Delta  + \gamma \left (    1 + g_a / g_s    \right)}}~ \iff ~  \overset{\text{Volume flux rate}}{ET_o=\frac{\Delta R_n   +   \rho_a c_p  \left(  \delta q  \right) g_a } { \left(   \Delta  + \gamma \left (    1 + g_a / g_s    \right)    \right) \lambda_v }}

λv = Latent heat of vaporization. Energy required per unit mass of water vaporized. (J/g)
Lv = Volumetric latent heat of vaporization. Energy required per water volume vaporized. (Lv = 2453 MJ m-3)
E = Mass water evapotranspiration rate (g s-1 m-2)
ETo = Water volume evapotranspired (m3 s-1 m-2)
Δ = Rate of change of saturation specific humidity with air temperature. (Pa K-1)
Rn = Net irradiance (W m-2), the external source of energy flux
cp = Specific heat capacity of air (J kg-1 K-1)
ρa = dry air density (kg m-3)
δe = vapor pressure deficit, or specific humidity (Pa)
ga = Hydraulic conductivity of air, atmospheric conductance (m s-1)
gs = Conductivity of stoma, surface conductance (m s-1)
γ = Psychrometric constant (γ ≈ 66 Pa K-1)

In Figure below it can be seen that the actual transpiration and evaporation are calculated in Porous Media module of MohidLand. The actual evaporation, which happens only at the soil surface, is calculated based on a pressure head limit imposed by the user. It allows the model to not evaporate any surface water, even if it is available, when the soil is very far way from the saturation point.

Evapotranspiration fluxogram in Mohid Land model

Other Features

Outputs

References

Data File

Keywords

Sample