Water Vapour
Water Vapour
H2O is the only matter that can be present at ambient conditions under its three phases:
water vapour, liquid water and ice.
97.4% of the water on Earth is assembled in the oceans, 2.6 % are on the continent,
either on the surface or in the underground, and only 10-3 % are in the atmosphere.
Nevertheless, these 0.001% are the essential link in the global water cycle. Fig. 3.2 gives
some estimates on the reservoirs and fluxes involved.
The reservoirs (Fig. 3.2) are given in 103 km3, and estimate the water masses available
in the ocean (1.4 109 km3), in the glaciers (28 106 km3), in the surface water (lakes,
rivers, ..) (3 105 km3), in the ground water (8 106 km3), in the atmosphere (13 103km3)
and in the biosphere (103 km3). The arrows indicate the net fluxes, with their orders of
magnitude in 103 km3 yr-1.
The largest flux corresponds to the evaporation of water from the oceans where the
major part of the evaporation takes place in the tropical oceans. However, a significant
part of the evaporated water condenses, forms clouds and rains close to the evaporation
site. Thus, only 40 103 km3 yr-1 are transported away from the vapour sources,
condenses over the continent and waters, thus, the land. Even though the precipitation
flux over lands is only about 1/3 of the precipitation flux over the ocean, the watering is
sufficient to allow a diversified vegetation to develop. A significant part of the
precipitated water evaporates again, another part feeds the water reservoirs, and the
excess runs off and re-joins the oceans.
Viewing the evaporation processes on the molecular level, it occurs when the flow of
individual water molecules away from the water body exceeds the return flow to the
water from the air. The actual amount of evaporation that takes place is a rather
complex function of five different parameters:
    •    Energy availability: The flow of molecules away from the water body requires a
         certain amount of energy. This energy can be thermal, i.e. evaporation increases
         with increasing temperature and the photons of the solar radiation contribute an
         important amount to this energy.
A.I. Flossmann                              11
    Fig.3.2: The global water cycle; reservoirs are given in 103 km3 and fluxes in 103 km3 yr-1.
    •    Reduction in atmospheric pressure: The molecules can leave the water surface
         more easily if the counteracting atmospheric pressure force is smaller.
    •    Humidity: When the air layer above the water surface is relatively dry and far
         from the saturation state, the vapour gradient between water surface and
         atmosphere is significant. This steep gradient reduces the return flow of
         molecules to the water body. In calm conditions, however, vapour will rapidly
         build up in the air layer, saturate it and the gradient will rapidly disappear.
    •    Presence of wind: Strong winds will constantly remove the moist layer from the
         water surface and, thus, maintain a steep gradient of water vapour and high
         evaporation rates.
    •    Water availability: Through open surfaces, water is easily accessible and
         maximum evaporation rates can be achieved. If the water is present in the soil,
         the vegetation or animals and humans, evaporation is still taking place. However,
         it is regulated by pores or stomata, and is, thus, called transpiration.
Commonly, PET (potential evapotranspiration) is the evaporation rate, which will occur
from an open water surface, while AET (actual evapotranspiration) is the amount of
water vapour formed from “resisting” surfaces. PET=AET over the ocean, e.g., while
AET<<PET in the subtropical deserts, e.g.
As already suggested above, there is a limit to the quantity of water in the gaseous
phase. This limit becomes evident through a little experiment: a vacuum is generated
inside a box. Then, liquid water is introduced into the box. Initially, the entire liquid is
instantaneously vaporized. After a certain quantity, part of the water remains liquid,
and any new quantity of water added will not increase the amount of vapour generated.
If more molecules are vaporized, other molecules from the vapour will condense. This
causes a stable number of vapour molecules in the gas phase, resulting in a vapour
A.I. Flossmann                                 12
pressure called “saturation vapour pressure” esat. Only by increasing temperature, the
amount of vapour molecules, and the value of the saturation vapour pressure will
increase.
 Fig. 3.3: A vaporization in a vacuum will result in a saturation vapour pressure, as a function of
 temperature
In order to describe the equilibrium condition between two phases of the same matter,
the two variables of state: p (pressure) and T (temperature) are generally used. They
are convenient, as easily controlled. Thermodynamics proposes a number of state
functions to describe a process. The free enthalpy G(p,T) is most adapted to study the
problem of phase transitions, as a function of p and T.
The total free enthalpy G of a matter present in any two phases with masses M1 and M2,
resp. , can be described by the sum of the two specific free enthalpies g1(p,T) and g2(p,T)
of each phase:
𝐺 𝑝, 𝑇 = 𝑀! 𝑔! 𝑝, 𝑇 + 𝑀! 𝑔! 𝑝, 𝑇 = 𝑀! 𝑔! 𝑝, 𝑇 + 𝑀 − 𝑀! 𝑔! (𝑝, 𝑇)
                                𝐺(𝑝, 𝑇)
                  𝑔 𝑝, 𝑇, 𝑥 =           = 𝑥 𝑔! 𝑝, 𝑇 − 𝑔! (𝑝, 𝑇) + 𝑔! (𝑝, 𝑇)
                                  𝑀
The state of the system, thus, depends on T, p, and x. For the two phases to exist in
equilibrium at a fixed pressure and temperature, the specific free enthalpy g(p,T,x) has
to be at a minimum and obey:
                 !"
                          = 0 => 𝑔! 𝑝, 𝑇 = 𝑔! (𝑝, 𝑇)                                          (3.1)
                 !" !,!
This equality of the two specific free enthalpies confirms a relationship p=p(T) between
pressure and temperature for the equilibrium condition. For the condition of the H2O
water and vapour phases being in equilibrium, this means:
For a state T+ΔT and p+Δp, where the phases are still in equilibrium, Eq.(3.1) can be
written as
                      𝑔! 𝑝 + Δ𝑝, 𝑇 + Δ𝑇 = 𝑔! (𝑝 + Δ𝑝, 𝑇 + Δ𝑇)
A.I. Flossmann                                  13
or
                                  𝑔! 𝑝, 𝑇 + 𝑑𝑔! = 𝑔! 𝑝, 𝑇 + 𝑑𝑔!
or with Eq.(3.1)
                                        𝑑𝑔! 𝑝, 𝑇 = 𝑑𝑔! (𝑝, 𝑇)
Thermodynamics allows to identify the specific volume v and specific entropy s through:
                                     𝜕𝑔                   𝜕𝑔
                                              =𝑣 ;                  =−𝑠
                                     𝜕𝑝   !
                                                          𝜕𝑇    !
and yields for water with p=esat:
                                           𝑑𝑒 !"# 𝑠! − 𝑠!
                                                 =
                                            𝑑𝑇     𝑣! − 𝑣!
Introducing the specific latent heat of phase transition at constant pressure and
temperature:
                             𝑙!" =      𝛿𝑄!,! = 𝑇 𝑠! − 𝑠! = ℎ! − ℎ!
yields the equation of Clausius – Clapeyron:
                                          !! !"#             !!"
                                                    =                                         (3.2)
                                              !"          !(!! !!! )
All first order equilibrium states between two phases of the same matter are described
by a Clausius-Clapeyron equation of the type (3.2). The specific latent heat needs to be
adapted to the process in question, and is a function of temperature:
Table 3.1 summarizes different values of latent heat of H2O at T=0°C. They are defined
positive, when the phase transition takes place from the state with less entropy (e.g.
liquid: v1) to the state with more entropy (e.g. vapour: v2).
The temperature dependency of latent heat can be estimated, using the expressions for
specific enthalpy:
                               ℎ! 𝑇 = ℎ! 𝑇! + 𝑐!" 𝑇 − 𝑇!
                               ℎ! 𝑇 = ℎ! 𝑇! + 𝑐! 𝑇 − 𝑇!
                                ℎ! 𝑇 = ℎ! 𝑇! + 𝑐! 𝑇 − 𝑇!
where cw and ci are the heat capacities of liquid water and ice (cw=4218 J kg-1 K-1;
ci=2106 J kg-1 K-1) and cpv is the specific heat of vapour at constant pressure.
A.I. Flossmann                                       14
                             𝑙! 𝑇 = 𝑙! 𝑇! + (𝑐!" − 𝑐! )(𝑇 − 𝑇! )
                              𝑙! 𝑇 = 𝑙! 𝑇! + (𝑐!" − 𝑐! )(𝑇 − 𝑇! )
and                             𝑙! 𝑇 = 𝑙! 𝑇 − 𝑙! (𝑇)
                                                         𝑅! 𝑇
                                       𝑣! = 𝑣! =
                                                        𝑒 !"#,!
                                 𝑑                 𝑙!                                       (3.3)
                                   ln(𝑒 !"#,! ) =
                                𝑑𝑇                𝑅! 𝑇 !
Fig. 3.4 displays the dependence of the saturation vapour pressure with respect to the
temperature, as detailed in Eq.(3.3).
If water vapour is vaporized in the air, and not in a vacuum as assumed above, then the
derivation still holds. However, the pressure of water vapour becomes the partial
pressure, following Dalton’s law (compare Eq. 2.1).
The same reasoning applies to the equilibrium between the vapour and the ice phase,
described by
                      !                   !!
                        ln 𝑒 !"#,! = !         !
                     !"                  !!
Fig. 3.4: Saturation vapour pressure over liquid (upper light blue) and over ice (dark blue) as a
function of temperature
A.I. Flossmann                                     15
represented also in Fig.3.4. The equilibrium curve between ice and vapour below 0°C is
presented here together with the liquid supercooled curve.The concept of a saturation
state of water vapour in the atmosphere is essential for cloud physics and is used to
derive a certain number of humidity measures.
                           !
                                   𝑀!           𝑀!
                          𝑞 = !              ≈
                             𝑀 + 𝑀! + 𝑀! + 𝑀! 𝑀! + 𝑀!
When neglecting the presence of liquid water and ice, and using Eq. (2.1), yields:
                                           𝑒                𝑒
                         𝑞! = 0.622                ≈ 0.622
                                      𝑝 − 0.378 𝑒           𝑝
or
                                               𝑒          𝑒
                                  𝑟 = 0.622       ≈ 0.622
                                              𝑝−𝑒         𝑝
3.3.1.5 approximations
Combining the equations yields:
                                                𝑟
                                        𝑞! =
                                               1+𝑟
A.I. Flossmann                                16
or
                                               𝑞!
                                          𝑟=
                                             1 − 𝑞!
and for small values of r and qv (in the atmosphere the values for both vary between 0
and 0.05):
                                           𝑟 ≈ 𝑞!
In the literature, approximated formulas for the saturation vapour pressures can be
found. An example is the formula of Magnus, which gives:
                                                      !! (!!!! )
                       𝑒 !"# 𝑇 = 𝑒 !"# 𝑇! 10         !! !(!!!! )                   (3.4)
A.I. Flossmann                                 17
The appearance of a “gas constant” Rh that is a function of the variable vapour content r
is not very useful. Consequently, using Rd instead of Rh is more convenient which defines
a new temperature Tv called virtual temperature:
𝑅! 𝑇 = 𝑅! 𝑇!
                              1 + 1.608𝑟
                 𝑇! = = 𝑇 !              ! ≈ 𝑇(1 + 0.608𝑟)                             (3.7)
                                 1+𝑟
                                              𝑒
                                     𝑅𝐻 =             ∙ 100
                                            𝑒 !"#,!
The ratio of the partial pressures of the actual state of the air and the saturation state is
a measure smaller or equal than unity. Generally, this value is multiplied by 100 to
give values of relative humidity in %. Using the relationship between e and qv, yields
also:
                                      𝑟              𝑞!
                             𝑅𝐻 ≈ !"#,! ∙ 100 ≈ !"#,! ∙ 100
                                   𝑟               𝑞
Due to the different values of the saturation vapour pressure over liquid water and over
ice, the resulting relative humidities are not identical.
A.I. Flossmann                               18
corresponding temperature is lower that 0°C, the point is called frost point τi, as frost is
formed. Dew and frost point can be calculated using e.g. Magnus formula (3.4):
                                                           𝑑!
                           𝜏 = 𝑇! − 𝑏!                   𝑒         +1
                                                    lg !"#    − 𝑑!
                                                      𝑒    𝑇!
For bv and dv the appropriate values for water or ice have to be used.
Using eq. (3.2) for water and ice allows calculating the ratio of dew point temperature
and frost point temperature close to T0:
                             !!       ! ! !"#,!          !! (!! )
                                  =                  =              = 0.89            (3.8)
                             !        !   ! !"#,!        !! (!! )
This approximation is only valid close to T0. For other temperatures, the temperature
dependence of the latent heat needs to be taken into account.
where U is internal energy and H=U+pV is enthalpy. Q and W are the heat added to the
system and the volume work done to the system (δW=-pdV).
Applying 1st principle to the psychrometer, it can safely be assumed that the evaporation
process is isobaric (dp=0) and adiabatic (δQ=0). Then, the equation becomes dH=0 or
dh=0, using specific enthalpy.
Quantifying the evaporation process occurring between an initial state and a final state
in equilibrium allows calculating the variation of enthalpy.
The amount of vapour necessary results from liquid water present initially. The process
can thus be idealized by two separate steps: the humid air around the bulb decreases in
temperature from T to Tw, while a part of the liquid water experiences a phase transition
to the vapour. The amount of liquid water necessary corresponds to reaching the
saturation vapour pressure esat,w(Tw).
A.I. Flossmann                                        19
Fig.3.5: Schematic display of the psychrometer processes
The variation of the enthalpy can, thus, be written as the variation due to the 2 steps:
where cp is the specific heat at constant pressure and dqv is the amount of vapour
necessary to increase the qv(T) to qsat(Tw):
𝑐! 𝑇! − 𝑇 + 𝑙! 𝑞 !"# 𝑇! − 𝑞! 𝑇 =0
                                                          !       !.!""
                                              𝑇 − 𝑇! = !!          !
                                                                          𝑒𝑠𝑎𝑡 (𝑇𝑤 ) − 𝑒    (3.9)
                                                              !
an equation for the wet bulb depression results. This equation can be transformed into
an equation for the vapour pressure e:
                                                          𝑝 𝑐𝑝
                               𝑒 = 𝑒𝑠𝑎𝑡 𝑇𝑤 −          (𝑇 − 𝑇𝑤 )
                                             0.622 𝑙𝑣
Assuming cp as the specific heat of dry air and neglecting the specific heat of water
vapour:
                        !!,! ! ! !!!,! ! !       !!,! !!!,! !
                 𝑐! =                        =                    ≈ 𝑐!,! + 𝑐!,! 𝑟 ≈ 𝑐!"    (3.10)
                            ! ! !! !                !!!
with cpd=1.0035 J g-1 K-1 and cpv=2.080 J g-1 K-1, and assuming that l12≈l12(T0), a constant
can be defined:
                                     𝑐!"
                             𝐴! =            = 0.65 10!! 𝐾 !!
                                  0.622 𝑙!!
                                                𝑐!"
                                  𝐴! =                 = 0.57 10!! 𝐾 !!
                                             0.622 𝑙!!
A.I. Flossmann                                         20
with a “Sprung”-constant of Ci=1757 K. The ratio between the two constants gives again
a factor of 0.89, as already documented for the different between dew and frost point
(compare Eq.(3.6)). Thus,
                                                𝑝
                              𝑒 = 𝑒 !"#,! 𝑇! − (𝑇 − 𝑇! )
                                               𝐶!
and
                                               𝑝
                               𝑒 = 𝑒 !"#,! 𝑇! − (𝑇 − 𝑇! )
                                               𝐶!
Tables for the web bulb depression in terms of relative humidity and dew point can be
found in the appendix.
           Note :
           cvd=5/2 Rd = 717 J kg-1 K-1     => cpd-cvd=Rd (eq. of Mayer)
           cpd=7/2 Rd
If the air temperature is 20°C and the wet bulb depression is 5°C, then the relative
humidity is 58%. If there is no-wet-bulb depression, the air is saturated and the relative
humidity is 100%.
An instrument that is less accurate than a psychrometer is the hair hygrometer. As the
name implies, this instrument uses hair as the humidity sensor. Hair lengthens slightly
as the relative humidity increases and shrinks with humidity decrease. The typical
length change of hair is on the order of 2.5% for relative humidity ranging from 0% to
100%. Usually, a sheaf of blond human hair is connected to a pointer, allowing the
reading of RH in percent. When the pointer draws a curve on a clock-driven drum, the
instrument is called a hygrograph. It allows recording continuously the fluctuations in
relative humidity with time.
Chilled mirror dew point hygrometers are some of the most precise instruments
commonly available. These use a chilled mirror and optoelectronic mechanism to detect
condensation on the mirror surface. The temperature of the mirror is controlled by
electronic feedback to maintain a dynamic equilibrium between evaporation and
condensation on the mirror, thus closely measuring the dew point temperature.
For applications where cost, space, or fragility is relevant; other types of electronic
sensors are used, at the price of a lower accuracy. In capacitive humidity sensors, the
effect of humidity on the dielectric constant of a polymer or metal oxide material is
measured.
A.I. Flossmann                             21
3.4.2 Remote sensing measurement techniques
All remote sensing methods rely on the fact that water vapour scatters and absorbs
radiation at different wavelengths in so-called absorption-bands (compare chapter 14).
The amount of radiation absorbed is proportional to the amount of vapour in the
pathway between the radiative emission source and the receptor. Consequently, all
these methods give integrated quantities along the path that is covered.
Numerous remote sensing instruments are currently space-borne on Earth-orbiting
satellites. They use different wavelengths, that present different advantages and
disadvantages.
The Air Mass Corrected Differential Optical Absorption Spectroscopy (AMC-DOAS), e.g.
on GOME or SCIAMACHY, is a method to retrieve total water vapour column amounts
from spectral measurements in the visible wavelength region around 700 nm. Because
data in the visible spectral range are analysed, the AMC-DOAS method is only applicable
to measurements on the dayside and to (almost) cloud-free ground scenes. A significant
advantage of the AMC-DOAS method is that the derived water vapour columns do not
depend on additional external information, like a calibration using radio sonde data
which is often used in the microwave spectral region. The AMC-DOAS water vapour
columns therefore provide a completely independent data set.
Wavelength in the near IR are used also, e.g. by Moderate Resolution Imaging
Spectroradiometer (MODIS) on board of the satellites Terra and Aqua, e.g. Since the
MODIS water vapour retrieval relies on observations of water vapour attenuation of
near Infrared (IR) solar radiation reflected by surfaces and clouds, it is sensitive to the
presence of clouds. The frequency and the percentage of cloud free conditions at mid-
latitudes is only 15-30% on average, limiting the usefulness of this spectrometer in
certain regions.
Yet another method for estimating the amount of water vapour in the atmosphere is
linked to the exploitation of the signals propagating from GPS satellites to ground-based
A.I. Flossmann                              22
GPS receivers, which are delayed by atmospheric water vapour. This delay is
parameterized in terms of a time-varying zenith wet delay (ZWD), which is retrieved by
stochastic filtering of the GPS data. Given surface temperature and pressure readings at
the GPS receiver, the retrieved ZWD can be transformed into an estimate of the
integrated water vapor (IWV) overlying that receiver. The emerging GPS networks offer,
thus, also the possibility of observing the distribution of water vapour.
Yet, another method to obtain integrated information of water vapour is the use of
Raman lidars, equipped with a water vapour channel. The Raman scattering is a weak
scattering phenomenon of light from a molecular medium. In the Raman scattering
mechanism, also known as inelastic scattering process, the scattered radiation is shifted
from the exciting radiation by an amount that is unique to the scattering molecule. Since
the vibrational energy shift is different for each type of molecule, Raman scattering can
be used to sense specific constituents in the atmosphere. Raman scattering is a weaker
interaction than Rayleigh scattering, with a cross-section that is typically three orders of
magnitude smaller. Because of the weak nature of Raman scattering, lidars using Raman
principle have ben optimized for nighttime investigation of the atmosphere.
All techniques deployed on satellites are also used on the ground or on other airborne
platforms (planes, balloons, ..).
Exercises:
    1.) For a station where p=1000 hPa, T= 10°C and RH=70% calculate dew point
        temperature, partial pressure, mixing ratio, virtual temperature and saturation
        mixing ratio.
    2.) For a station where T=15°C, r=6.22 g kg-1 and RH=50% calculate dew point
        temperature, pressure and saturation vapour pressure.
    3.) Which three parameters allow to calculate all other humidity measures:
        a) pressure, dew point temperature and mixing ratio
        b) pressure, dew point temperature and saturation mixing ratio
        c) pressure, temperature and mixing ratio
        d) vapour pressure, mixing ratio and dew point temperature
        Which will be the dew point temperature and the relative humidity if the air
        inside the cylinder is heated to 20°C?
A.I. Flossmann                              23
        Then, the temperature in the cylinder is decreased to 1°C. Which will be the dew
        point temperature and the relative humidity?
A.I. Flossmann 24