Personal tools

Difference between revisions of "Module PorousMediaProperties"

From MohidWiki

Jump to: navigation, search
(Replacing page with '==Overview== ==Main Processes== ==Other Features== ==Outputs== ==References == ==Data File == ===Keywords=== ===Sample=== Category:MOHID Land')
Line 2: Line 2:
  
 
==Main Processes==
 
==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.
 
  
<math>J=-K\left ( \theta  \right )\left ( \frac{\partial h}{\partial \theta } + 1 \right )</math>
+
==Other Features==
===Water retention===
 
The shape of water retention curves can be characterized by several models, one of them known as the van Genuchten model:
 
:<math>\theta(h) = \theta_r + \frac{\theta_s - \theta_r}{\left[ 1+(\alpha |h|)^n \right]^{1-1/n}}</math>
 
where
 
:<math>\theta(h)</math> is the water retention curve [L<sup>3</sup>L<sup>−3</sup>];
 
:<math>|h|</math> is suction pressure ([L<sup>−1</sup>] or cm of water);
 
:<math>\theta_s</math> saturated water content [L<sup>3</sup>L<sup>−3</sup>];
 
:<math>\theta_r</math> residual water content [L<sup>3</sup>L<sup>−3</sup>];
 
:<math>\alpha</math> is related to the inverse of the air entry suction, <math>\alpha >0</math> ([L<sup>−1</sup>], or cm<sup>−1</sup>); and,
 
:<math>n</math> is a measure of the pore-size distribution, <math>n>1</math> (dimensionless).
 
===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 }}
 
</math>
 
 
 
:''λ''<sub>v</sub> =  Latent heat of vaporization. Energy required per unit mass of water vaporized.  (J/g)
 
:''L''<sub>v</sub> = Volumetric latent heat of vaporization. Energy required per water volume vaporized. (''L''<sub>v</sub> = 2453 MJ m<sup>-3</sup>)
 
 
 
:''E'' = Mass water evapotranspiration rate (g s<sup>-1</sup> m<sup>-2</sup>)
 
:''ET''<sub>o</sub> = Water volume evapotranspired (m<sup>3</sup> s<sup>-1</sup> m<sup>-2</sup>)
 
 
 
:&Delta; =    Rate of change of saturation specific humidity with air temperature.  (Pa K<sup>-1</sup>)
 
:''R''<sub>n</sub> = Net irradiance (W m<sup>-2</sup>), the external source of energy flux
 
:''c''<sub>p</sub> = Specific heat capacity of air (J kg<sup>-1</sup> K<sup>-1</sup>)
 
:''ρ''<sub>a</sub> = dry air density (kg m<sup>-3</sup>)
 
:δ''e'' = vapor pressure deficit, or specific humidity (Pa)
 
:''g''<sub>a</sub> = Hydraulic conductivity of air, atmospheric conductance (m s<sup>-1</sup>)
 
:''g''<sub>s</sub> = Conductivity of stoma, surface conductance  (m s<sup>-1</sup>)
 
:''γ'' = Psychrometric constant  (''&gamma;'' ≈ 66 Pa K<sup>-1</sup>)
 
 
 
All of these calculations about potential evapotranspiration are made in module Basin in MohiLand model. However, not all of the potential water that can be evaporated or transpirated will be in fact removed from the soil. The water that will really leave the soil through these processes is calculated in the Porous Media module.
 
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]]
 
  
==Other Features==
 
 
==Outputs==
 
==Outputs==
  
Line 52: Line 14:
  
 
===Sample===
 
===Sample===
 +
 +
[[Category:MOHID Land]]

Revision as of 18:00, 7 February 2011

Overview

Main Processes

Other Features

Outputs

References

Data File

Keywords

Sample