Penman Monteith公式

Monteith在1981改善Penman公式为Penman Monteith公式为以下形式：^[6]

${\displaystyle ET_{0}={\Delta (R_{n}-G)+\rho _{a}c_{p}{e_{s}-e_{a} \over r_{a}} \over \Delta +\gamma (1+{r_{s} \over r_{a}})}}$ （式1）

其中：

R_n为全天空日射量(MJ m^-2 day^-1)

G为地表热通量(MJ m^-2 day^-1)

ρ_a为空气密度

c_p为空气比热

e_s 为饱和蒸汽压，可由Goff-Gratch方程式计算

e_a 为实测蒸气压，即为相对湿度乘上e_s

r_a为空气动力学阻力

r_s为植披所造成的传输阻力

Δ为水蒸气饱和曲线对温度作图之斜率（ kPa/ °C），可用式2计算^[7]：

${\displaystyle \Delta ={4099e_{a} \over (T+237.3)^{2}}}$ ，温度单位为°C（式2）

γ：干湿度常数（psychrometric constant）可用式3计算^[7]：

${\displaystyle \gamma =0.00163{P \over \lambda }}$ ，P为大气压力（kPa），${\displaystyle \lambda }$ 为水汽化潜热2.453 （MJ /kg）（式3）

r_a的计算较为繁复，需要风速、地表粗糙长度等资料；r_s则受叶面积指数(LAI)、叶片形状等因素影响，因此这两个参数随不同地表、不同作物而有差异，这使得Penman-Monteith方程式之泛用性受限。而后联合国粮食及农业组织为推广此公式适用于所有地区，故选定参考作物之实验数据，并以以下假设简化参数^[8]：

(1)  假设参考作物高度为0.12 （m），且其叶面反射率为0.23^[9]。

(2)      令空气动力学阻力r_a = 208/u₂ s m^-1

(3)      令植披所造成的传输阻力r_s = 70 s m^-1

故将Penman-Monteith方程式简化如下：

${\displaystyle ET_{0}={0.408\Delta (R_{n}-G)+{900 \over T-273}u_{2}(e_{s}-e_{a}) \over \Delta +\gamma (1+0.34u_{2})}}$ （式4）

其中：

ET_o reference evapotranspiration (mm day^-1)

R_n 全天空日射量(MJ m^-2 day^-1)

G 地表热通量(MJ m^-2 day^-1)

T 日平均气温(°C)

u₂ 离地表两米之风速(m s^-1)

e_s 饱和蒸汽压(kPa)，可由Goff-Gratch方程式计算

e_a 实测蒸气压(kPa)，即为相对湿度乘上e_s

e_s - e_a 饱和蒸汽压差 (kPa)

Δ 水蒸气饱和曲线对温度作图之斜率 (kPa °C^-1)

γ 干湿度常数 (kPa °C^-1)

式4为联合国粮食及农业组织所开发之作物模拟软件 AquaCrop^[10]所采用。在 AquaCrop中，各种形式的作物系数（K _c ）解释了建模的特定植被与参考蒸散量（RET或ET ₀ ）标准之间的差异。应力系数（K _s ）解释了由于环境应力而导致的ET降低（例如，土壤饱和度降低了根区O ₂ ，低土壤湿度会引起枯萎，空气污染和盐分化）。

Priestley-Taylor公式是Penman-Monteith公式的另一简化版本，其将空气非饱和水蒸气和对流项简化为单个经验常数α。 PT模型可以表示如下：

${\displaystyle ET_{0}pt=\alpha {\Delta (R_{n}-G) \over \lambda _{u}(\Delta +\gamma )}1000}$ （式5）

其中：

ETpt：Priestley-Taylor蒸散量，单位为mm day-1

α ：解释蒸汽压差和电阻值的经验常数[-]

通常，开放水域的α为1.26，但取值范围很广，从小于1（潮湿条件）到几乎2（干旱条件）^[11]。