Changes between Version 6 and Version 7 of doc/tec/lsm


Ignore:
Timestamp:
Apr 4, 2016 2:54:19 PM (9 years ago)
Author:
maronga
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • doc/tec/lsm

    v6 v7  
    1010{{{
    1111#!Latex
    12 \begin{equation}
     12\begin{equation*}
    1313   C_0 \dfrac{dT_0}{dt} = R_\mathrm{n} - H - LE - G
    14 \end{equation}
     14\end{equation*}
    1515}}}
    1616where ''C'',,0,, and ''T'',,0,, are the heat capacity and radiative temperature of the surface skin layer, respectively. ''R'',,n,,, ''H'', ''LE'', and ''G'' are the net radiation, sensible heat flux, latent heat flux, and ground (soil) heat flux at the surface, respectively. ''H'' is calculated as
    1717{{{
    1818#!Latex
    19 \begin{equation}
     19\begin{equation*}
    2020  H = - \rho\ c_\mathrm{p}\ \dfrac{1}{r_\mathrm{a}} ( \theta_1 - \theta_0 )
    21 \end{equation}
     21\end{equation*}
    2222}}}
    2323where ''ρ'' is the density of the air, ''c'',,p,, = 1005 J kg^-1^ K^-1^$ is the specific heat at constant pressure, ''r'',,a,, is the aerodynamic resistance, and ''θ'',,0,, and ''θ'',,1,, are the potential temperature at the surface and at the first grid level above the surface, respectively. ''r'',,a,, is calculated via Monin-Obukhov similarity theory, based on roughness lengths for heat and momentum and the assumption of a constant flux layer between the surface and the first grid level.
     24
     25''G'' is parametrized as (Duynkerke 1999)
     26{{{
     27#!Latex
     28\begin{equation*}
     29  G = \Lambda ( T_0 - T_{\mathrm{soil},1} )
     30\end{equation*}
     31}}}
     32with ''Λ'' being the heat conductivity between skin layer and the soil, and ''T'',,soil,1,, being the temperature of the uppermost soil layer. The latent heat flux is calculated as
     33{{{
     34#!Latex
     35\begin{equation*}
     36  LE = - \rho\ l_\mathrm{v}\ \dfrac{1}{r_\mathrm{a} + r_\mathrm{s}} ( q_{\mathrm{v},1} - q_{\mathrm{v,sat}}(T_0) )\;.
     37\end{equation*}
     38}}}
     39Here, ''l'',,v,, = 2.5 * 10^6^ J kg^-1^ is the latent heat of vaporisation, ''r'',,s,, is the surface resistance, ''q'',,v,1,, is the specific humidity at first grid level, and ''q'',,v,sat,, is the saturation specific humidity at temperature ''T'',,0,,.
     40
     41All equations above are solved locally for each surface element of the LES grid. Each element can consist of both patches of bare soil, vegetation, and a liquid water reservoir, which is the interception water stored on plants and soil from precipitation. Therefore, an additional equation is solved for the liquid water reservoir. ''LE'' is then calculated for each of the three components (bare soil, vegetation, liquid water). The resistances are calculated following Jarvis (1976).
     42
     43''C'',,0,, is set to zero, the energy balance is solved implicitly by linearising ''q'',,v,sat,,.
    2444
    2545== Soil model ==