基于aser数据的黑河中下游水面蒸发量估算.docxVIP

基于aser数据的黑河中下游水面蒸发量估算 水的蒸发过程是自然界和生态系统的重要组成部分。在陆地、地区和流域尺度上，水蒸发损失的估计对监测、调查和管理水资源是重要的。Lenters研究了美国北部湖面蒸发量的变化趋势, 研究得出水面蒸发速率在时间尺度上变化很大, 季节变化最高可达18%。Alvarez研究了农业灌溉水库蒸发水量的损失量以及对水分储存的影响。Gianniou利用能量平衡法和一维水面温度涡度扩散模型的结合估算了水面的日蒸发量。Tanny利用Penman-Monteith模型模拟了水面蒸发量, 并和涡度相关法实测的水面蒸发量进行了对比。曾燕研究了分析出1960～2000年中国蒸发皿蒸发量呈明显下降趋势, 杜军研究得出藏北高原各站点蒸发皿蒸发量呈现显著的减少趋势。 利用气象资料估算的水面蒸发只能代表以点为基础的小范围的水面蒸发, 无法反应出水面蒸发的空间变化规律以及区域特性, 进而影响了区域尺度上水面蒸发变化规律的探求。近年来随着遥感的出现, 用遥感资料估算水面蒸发成为一个有效途径, 并且能够监测大范围水面蒸发的变化趋势。Bastawesy利用ASTER 数据和SPOT数据结合估算了埃及水面蒸发量的变化趋势。Wang利用MODIS数据估算了墨西哥的水面蒸发。国内一些学者利用遥感资料估算了区域蒸散发的时空变化规律, 如吴炳方、田辉等, 这些研究都把水面视为植被和土壤对待, 利用估算植被蒸散发模型估算水面蒸发。实际上, 水面蒸发有时要远远大于植被蒸散发。除此之外, 受遥感资料空间分辨率太低的影响, 估算非均匀下垫面蒸散发空间变化规律时产生了很大的误差。ASTER数据的空间分辨率最高可达15 m, 能够提高非均匀陆面蒸散发的估算精度。 基于以上问题的考虑, 本文利用ASTER遥感数据, 以黑河流域中下游为研究区域, 研究基于ASTER 数据估算水面蒸发量的精确性以及探讨ASTER 数据反演的反射率、水面温度、太阳净辐射以及水面蒸发量的空间变化范围以及不同时期的时间变化规律, 为西北干旱区水面蒸发规律的探究提供参考。 1 空气动力学参数 水面蒸发模型采用Penman-Monteith-Unsworth模型进行估算 λEi=Δ(Rni?Gi)+ρCp(es?ea)/ra(Δ+γ)λEi=Δ(Rni-Gi)+ρCp(es-ea)/ra(Δ+γ)(1) 式中, Rni为瞬时净辐射通量[W/ (m2·h) ];Gi为瞬时土壤热通量[W/ (m2·h) ];ρ、Cp分别为干空气的密度 (kg/m3) 和比热[1.013kJ/ (kg·℃) ];γ为干湿球常数 (kPa/℃) ;Δ为饱和水汽压与温度曲线的斜率 (kPa/℃) ;es为参照高度处的饱和水汽压 (kPa) ;ea为参考高度处的实际水汽压 (kPa) ;ra为空气动力学阻力 (s/m) ;λ为蒸发潜热 (MJ/kg) ;3.6为kJ/s转化为MJ/h的单位转换系数。 ra=4.72ln(zz0)2/(1+0.54u2)ra=4.72ln(zz0)2/(1+0.54u2)(2) 式中, u2为2m处的风速;z为参考高度 (m) ;z0为地表粗糙长度 (m) 。 1.1 地表反射率参数sr 由于辐射资料的观测站点较少, 因此辐射资料是通过遥感资料反演获得。由于遥感数据是观测日的瞬时数据, 本文通过遥感数据首先获得观测日的瞬时辐射数据 Rni= (1-α)Rs(in)+εsRL(in)-RL(out)(3) 式中,Rs(in)为入射太阳短波辐射[W/ (m2·h) ];α为地表反射率 (无量纲) ;RL(in)为入射的大气长波辐射[W/ (m2·h) ];RL(out)为反射的地表长波辐射[W/ (m2·h) ];εs为地表比辐射率 (无量纲) 。 地表反射率估算公式为 α=0.484α1+0.335α3-0.324α5+0.551α6+ 0.305α8-0.367α9-0.001 5 (4) 式中,α1、α3、α5、α6、α8、α9分别表示1、3、5、6、8、9波段的反射率。 1.2 土壤热态的瞬时值为 对于水体, 其瞬时土壤热通量Gi传输机制1～6月份和7～12月份有所不同, 用方程表示为 1.3 日蒸散与日土壤热通量之差之比对 蒸发比率不变法其基本假设是, 在能量平衡公式中, 能量通量组分的相对比例在白天稳定不变。根据蒸发比率不变法, 利用瞬时蒸散和瞬时太阳净辐射与瞬时土壤热通量之差之比等于日蒸散和日太阳净辐射与日土壤热通量之差之比来估算日蒸散 λEd=λEiRnd?GdRni?GiλEd=λEiRnd-GdRni-Gi(6) 式中, Rnd为日太阳净辐射[W/ (m2·d) ], 通过气象资料估算获得;Gd为日土壤热通量;λEi为瞬时水面蒸发量 (mm/h) ;λE