FLUXES AT INTERFACES
How is this done you ask??? Let's discuss fluxes at interfaces.
Fluxes at Interfaces
The vertical flux (Fz), of any variable
z, is assumed to be driven by the difference in z
across the interface. Thus: 
Fz = -UT (ztop - zbottom) (1)
where UT represents a transport velocity across the interface, while ztop and zbottom are the values just above and below the boundary.
Hence, we can write the following equations for the surface fluxes of momentum, heat, and moisture.
(u'w')s = -CD M (Uair - Uground) = -CD M (Uair) (1)
(v'w')s = -CD M (Vair - Vground) = -CD M (Vair) (2)
where M is the mean horizontal wind speed and CD is the bulk drag coefficient. CD, M , and Uair should all be calculated at the same height, usually 10 meters is used.
For heat and moisture:
(w'q')s = -CH M (qair - qground) (3)
(w'q')s = -CE M (qair - qground) (4)
where CH and CE are the bulk transfer coefficients for heat and moisture, respectively. Typical values for these coefficients are 1x10-3 to 5 x 10-3 (they are dimensionless).
Notice from (3) and (4) that to calculate surface heat and moisture fluxes, one needs the values of theta and mixing ratio at the ground, within the top 1 mm of soil or sea. These ground values of theta and mixing ratio are called the skin values.
How in the world does an NWP model figure out skin values of temperature and moisture????? Answer: many of the newer, more sophisticated models have soil-layer models that are coupled to the atmospheric model.
typical soil layer models contain a few (on order of 5) vertical layers of soil and describe the transfer and storage of heat and moisture within the ground.
They also contain detailed information about the land use and type over the model domain.
Now you can understand the need for a soil-layer model for NWP.
Let's discuss the basics of a soil layer model in more detail