Demystifying Evaporation Rates: The Physics, Mechanics, and Mathematical Models (Part 1)

Introduction: The Ubiquity and Complexity of Evaporation

Evaporation is one of the most fundamental phase-change processes in nature and industry. Whether observing a puddle evaporating on a sunlit pavement, managing industrial chemical drying vats, or designing massive agricultural irrigation systems, understanding and accurately predicting the rate of evaporation is critical.

At its core, evaporation is the process by which molecules at the liquid-gas interface gain sufficient kinetic energy to overcome intermolecular cohesive forces and escape into the surrounding atmosphere as vapor. While the concept sounds deceptively simple, quantifying the exact rate at which this occurs is a complex multiphysics problem. It requires integrating thermodynamics, fluid mechanics, meteorology, and boundary-layer theory.

There is no single "universal" formula that applies universally to every scenario. Instead, scientists, engineers, and hydrologists rely on a hierarchy of specialized formulas tailored to specific conditions—ranging from simple empirical mass-transfer equations like Dalton’s Law to highly sophisticated energy-balance models like the Penman-Monteith equation.

In this first part of our comprehensive expert guide, we will explore the fundamental driving forces of evaporation, dissect the classic mass-transfer models, and analyze how variables such as vapor pressure, wind speed, and temperature dictate the speed of phase transition.

The Primary Driving Force: Vapor Pressure and Deficit

To understand any evaporation formula, one must first master the concept of vapor pressure. Evaporation is not merely driven by heat; it is fundamentally driven by a concentration gradient—specifically, the difference between the moisture concentration at the evaporating surface and the moisture concentration in the bulk atmosphere above it.

Saturation Vapor Pressure ($e_s$)

Every liquid surface maintains an invisible layer of saturated air directly in contact with it. The pressure exerted by water vapor molecules in this saturated boundary layer is known as the saturation vapor pressure ($e_s$). This value is heavily dependent on temperature. As temperature increases, the kinetic energy of the water molecules rises exponentially, allowing more molecules to escape into the vapor phase, which drastically increases $e_s$.

The relationship between temperature and saturation vapor pressure can be approximated using the Tetens equation:

Where:

• $e_s$ is the saturation vapor pressure in kilopascals ($\text{kPa}$).

• $T$ is the air temperature in degrees Celsius (°C).

Actual Vapor Pressure ($e_a$)

Conversely, the actual vapor pressure ($e_a$) represents the partial pressure exerted by the water vapor already present in the surrounding ambient air.

Vapor Pressure Deficit (VPD)

The true engine driving evaporation is the Vapor Pressure Deficit (VPD), defined mathematically as the difference between the saturation vapor pressure at the surface and the actual vapor pressure of the ambient air:

If $e_s>e_a$, evaporation occurs because molecules are leaving the liquid faster than they are condensing back into it. If $e_s=e_a$, the system is in dynamic equilibrium (100% relative humidity relative to the surface temperature), and the net evaporation rate drops to zero. If $e_a>e_s$, condensation occurs instead.

The Historical Foundation: Dalton’s Law of Evaporation

The earliest systematic approach to quantifying evaporation rates was formulated by the English chemist and physicist John Dalton in 1802. Dalton recognized that evaporation is directly proportional to the difference between the vapor pressure of the water surface and that of the overlying air, and that air movement accelerates this process.

The Mathematical Formulation

Dalton expressed the rate of evaporation ($E$) through the following classic mass-transfer relationship:

Where:

• $E$ is the rate of evaporation (typically expressed in depth per unit time, such as $\text{mm/day}$ or ${\text{kg/m}}^2\text{s}$).

• $C$ is an empirical mass-transfer coefficient that accounts for atmospheric pressure, geometry, and surface roughness.

• $e_s$ is the saturation vapor pressure at the water surface temperature.

• $e_a$ is the vapor pressure of the ambient air at a specified height.

Expanding Dalton’s Law with Wind Speed

Dalton observed that stagnant air quickly becomes saturated, slowing evaporation to a crawl. To account for this, modern adaptations of Dalton’s formula incorporate a wind function, $f(u)$, which scales the mass-transfer coefficient based on wind velocity:

Where:

• $u_z$ is the wind speed measured at height $z$ (usually in meters per second, $\text{m/s}$).

• $a$ and $b$ are empirical constants derived experimentally for specific water bodies (such as lakes, reservoirs, or drying pans).

While Dalton’s approach is foundational, it remains limited because it assumes constant surface temperature and does not explicitly account for net radiation or sensible heat transfer from the atmosphere.

Aerodynamic Mass Transfer and Boundary Layer Theory

To build upon Dalton's empirical foundation, fluid dynamics introduces the concept of the atmospheric boundary layer. When air moves across an evaporating surface (like a wet soil bed, a crop canopy, or a lake), a thin layer of air—known as the vapor boundary layer—forms immediately above the interface.

The Role of Turbulence

Within this boundary layer, moisture transport occurs through two distinct mechanisms:

1. Molecular Diffusion: Movement of water molecules along a concentration gradient driven by random thermal motion.

2. Eddy Diffusion (Turbulence): Bulk transport of moisture packets driven by swirling air currents and wind shear.

Because molecular diffusion is extremely slow, the rate-limiting step of evaporation in open environments is almost entirely governed by turbulent transport through the boundary layer.

The Aerodynamic Evaporation Equation

Using turbulent transfer coefficients and logarithmic wind profile laws, meteorologists express the aerodynamic evaporation rate ($E_a$) as:

Where:

• ${\rho }_a$ is the density of the air (${\text{kg/m}}^3$).

• $ϵ$ is the ratio of the molecular weight of water vapor to dry air ($\approx 0.622$).

• $P$ is the atmospheric pressure ($\text{kPa}$).

• $r_a$ is the aerodynamic resistance ($\text{s/m}$), which quantifies how efficiently turbulence transports vapor away from the surface.

The aerodynamic resistance $r_a$ is inversely proportional to wind speed and surface roughness. It is mathematically modeled as:

Where $k$ is von Kármán's constant ($\approx 0.41$), $z_m$ and $z_h$ are measurement heights for wind and humidity, $d$ is the zero-plane displacement height, and $z_{om}$ and $z_{oh}$ are roughness lengths for momentum and heat/vapor transfer.

Summary of Part 1

In this first part, we have established that calculating the rate of evaporation requires analyzing the vapor pressure deficit ($e_s-e_a$) as the fundamental thermodynamic driver. We explored:

• How saturation vapor pressure varies non-linearly with temperature via the Tetens equation.

• How John Dalton pioneered the connection between vapor pressure gradients and evaporation rates.

• How wind speed and boundary layer turbulence modify mass transfer through aerodynamic resistance.

In the upcoming second part of this guide, we will bridge mass transfer with energy conservation, examining the Energy Balance Method and the industry-standard Penman-Monteith Equation, which unites aerodynamic principles with net radiation to yield precise evaporation rates across hydrology and agronomy.