实际蒸散发和潜在蒸散发

看了好多天蒸散发的资料,有所悟。
ET分为实际蒸散发ETa和潜在蒸散发ETp:
潜在蒸散发ETp包括好多种计算方法,ET0是ETp的一种,其基于P-M公式,是目前公认的比较准确的计算方法。但还有很多计算ETp的其他公式,如基于辐射计算的ETp、基于温度计算的ETp(见表1)。
ETc是作物实际蒸散发,由ET0乘以作物系数Kc得到,ETc主要的应用范围是农作物,如小麦、玉米、水稻、牧草等作物。除了作物类型的实际蒸散发,还有其他植被类型,如森林、灌从、草地等,因此ETc从属于ETp。
图

图2 张晓琳等: 五种潜在蒸散发公式在汉江流域的应用 表1来自张晓琳等: 五种潜在蒸散发公式在汉江流域的应用)

  • 2
    点赞
  • 16
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
根据提供的引用内容,MATLAB提供了多种方法来计算潜在蒸散。以下是两种常用的方法: 1. Penman-Monteith方法: ```matlab % 输入参数 T = 25; % 温度(摄氏度) RH = 60; % 相对湿度(百分比) Rs = 20; % 太阳辐射(MJ/m^2/day) u2 = 2; % 风速(m/s) % 计算潜在蒸散 ETo = evapotranspiration_penman_monteith(T, RH, Rs, u2); % 输出结果 disp(['潜在蒸散为:', num2str(ETo), ' mm/day']); function ETo = evapotranspiration_penman_monteith(T, RH, Rs, u2) % 常数 lambda = 2.45; % 水的潜热(MJ/kg) gamma = 0.067; % 斜率饱和汽压曲线(kPa/°C) % 饱和汽压 es = 0.6108 * exp(17.27 * T / (T + 237.3)); % 实际汽压 ea = es * RH / 100; % 风速修正 u2 = u2 * 3.6; % 将风速转换为km/h u2 = max(u2, 0.5); % 风速最小值为0.5 km/h % 潜在蒸散计算 ETo = (0.408 * delta * (Rs - Rn) + gamma * (900 / (T + 273)) * u2 * (es - ea)) / (delta + gamma * (1 + 0.34 * u2)); function delta = delta_calculation(T) % 斜率饱和汽压曲线计算 es = 0.6108 * exp(17.27 * T / (T + 237.3)); delta = (4098 * es) / ((T + 237.3) ^ 2); end function Rn = net_radiation(Rs, albedo) % 净辐射计算 Rn = (1 - albedo) * Rs; end end ``` 2. Hargreaves方法: ```matlab % 输入参数 Tmax = 30; % 最高温度(摄氏度) Tmin = 20; % 最低温度(摄氏度) Tmean = (Tmax + Tmin) / 2; % 平均温度(摄氏度) Ra = 25; % 太阳辐射(MJ/m^2/day) % 计算潜在蒸散 ETo = evapotranspiration_hargreaves(Tmean, Ra); % 输出结果 disp(['潜在蒸散为:', num2str(ETo), ' mm/day']); function ETo = evapotranspiration_hargreaves(Tmean, Ra) % 常数 a = 0.0023; b = 0.5; % 潜在蒸散计算 ETo = (0.0023 * Ra * (Tmean + 17.8) * (Tmean - 7.5) ^ 0.5); function Ra = extraterrestrial_radiation(latitude, day_of_year) % 太阳辐射计算 dr = 1 + 0.033 * cos(2 * pi / 365 * day_of_year); delta = 0.409 * sin(2 * pi / 365 * day_of_year - 1.39); phi = latitude * pi / 180; ws = acos(-tan(phi) * tan(delta)); Ra = (24 * 60 / pi) * 0.082 * dr * (ws * sin(phi) * sin(delta) + cos(phi) * cos(delta) * sin(ws)); end end ```
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值