next up previous
Next: About this document Up: The vertical diffusion and Previous: The numerical scheme used

The surface fluxes

The sensible and latent heat fluxes at the surface can be written in a generic form :

  equation586

Depending on the type of closure chosen by the land surface scheme the flux will be computed with variables at time t or t+1. For surface sensible heat ( tex2html_wrap_inline1355 ) flux the following correspondences are used :

  eqnarray599

where tex2html_wrap_inline1357 is the temperature at the first atmospheric level, p is the pressure at this level, tex2html_wrap_inline1361 pressure at some reference level, tex2html_wrap_inline1363 is surface temperature, tex2html_wrap_inline1365 is surface transfer coefficient for heat and moisture (Assumed to be same here.), tex2html_wrap_inline1367 is the density and tex2html_wrap_inline1369 is the wind speed at the first level. Equation 20 yields evaporation ( tex2html_wrap_inline1371 ) if the following definitions are used :

eqnarray619

where tex2html_wrap_inline1373 is the water vapor mixing ratio at the first atmospheric level, tex2html_wrap_inline1375 is the saturated humidity for the surface temperature and tex2html_wrap_inline1159 is the aridity coefficient. tex2html_wrap_inline1159 is defined as the ratio between actual and potential evaporation.

Computing these fluxes is the responsibility of the land surface scheme and our aim is to derive an equation for tex2html_wrap_inline1341 using tex2html_wrap_inline1383 as sole information from the surface. This ensures to the land surface scheme a maximum of freedom in the choice of the numerical scheme used for solving the surface energy balance equation. Inserting equation 20 into equation 15 yields :

 

Thus the only knowledge required for computing tex2html_wrap_inline1341 and iterating equation 16 from l=2 to the top are the old values of potential enthalpy and humidity at the lowest atmospheric level, the surface fluxes and finally the values of tex2html_wrap_inline1315 and tex2html_wrap_inline1317 for tex2html_wrap_inline1343 obtained during the descent.



POLCHER Jan
Fri Mar 6 16:09:11 MET 1998