Welcome to pyTSEB’s documentation!


The Two Source Energy Balance (TSEB) model computes the turbulent fluxes for two layers, soil and vegetation, with the interaction between the layers being allowed. Therefore, following an analogy between the flux transport and the Ohm’s Law for transport of electricity, the network of sensible H and latent \(\lambda E\) heat flux transfers can be thought of as being in series and represented as in:

TSEB resistance network in series

In the TSEB scheme \(\lambda E\) is usually estimated as a residual of the surface energy balance:

\[\begin{split}\lambda{}E_{S}&\approx R_{n,S}-G-H_{S}\\ \lambda{}E_{C}&\approx R_{n,C}-H_{C}\end{split}\]


\[\begin{split}H&=H_C+H_S\\ &=\rho_{air} C_{p}\frac{T_{AC}-T_{A}}{r_{a}}\\ H_{C}&=\rho_{air} C_{p}\frac{T_{C}-T_{AC}}{r_{x}}\\ H_{S}&=\rho_{air} C_{p}\frac{T_{S}-T_{AC}}{r_{s}}\end{split}\]

In the above equations \(\rho_{air}\) is the density of air, \(C_{p}\) is the heat capacity of air at constant pressure, \(T_{AC}\) is the air temperature at the canopy interface:


\(T_C\) and \(T_S\) are related to the directional radiometric temperature \(T_{rad}\left(\theta\right)\) by

\[\sigma T_{rad}^4\left(\theta\right)=f_c\left(\theta\right)\sigma\,T_{C}^4+\left[1-f_{c}\left(\theta\right)\right]\sigma\,T_{S}^4\]

with \(f_c\left(\theta\right)\) representing the fraction of vegetation observed by a sensor pointing at a zenith angle \(\theta\)


and \(\kappa_{be}\left(\theta\right)=\frac{\sqrt{\chi^2+\tan^2\theta}}{\chi+1.774\left(\chi+1.182\right)^{-0.733}}\) being the extinction coefficient of a canopy with a leaf angle distribution function defined by the Cambpell 1990 \(\chi\) parameter.

\(r_{a}\), \(r_{x}\) and \(r_{s}\) are the resistances to heat and momentum transfer in the surface layer, from the leaf canopy and from the soil surface respectively.

Indices and tables