基于能量平衡原理的黄河源区逐时陆面蒸散发量估算.docxVIP

基于能量平衡原理的黄河源区逐时陆面蒸散发量估算 遥感估算的蒸发分布是获得区域规模的主要手段。国内外很多学者在遥感反演蒸散发精度方面做了大量的工作,陈添宇等利用卫星资料和地表能量平衡算法估算了中国西北地区东部的日蒸散发;刘朝顺等利用卫星资料和改进的遥感模型反演了山东省月蒸散发的变化情况;Bastiaanssen等利用能量平衡方程改善了遥感模型的精度;Su等通过结合地面观测资料,发展了不同尺度的遥感估算蒸散发的模型。但以上方法基本是利用极轨气象卫星的资料,估算晴天无云条件下蒸散发,对于有云情况都进行回避。本文基于陆表能量平衡系统(Surface Energy Balance System,SEBS)模型,利用极轨气象卫星高空间分辨率和静止气象卫星高时间分辨率的优势,结合黄河源区气象观测资料,实现一天中包括有云覆盖的(无降水过程)不同时刻黄河源区陆面蒸散发量的实时监测。 1 陆面蒸发面积的估计理论 1.1 地表反照率法 Rn=(1-α)Q+L↓-L↑ (1) 式中α为地表反照率,由FY-2D卫星遥感观测资料确定;Q为太阳总辐射;L↓为向下的长波辐射;L↑为向上的长波辐射。 1.1.1 采用自适应滤波方案计算晴天条件下的净辐射 根据SEBS模型估算晴天条件下净辐射空间分布。 1.1.2 无云区地面的温度 黑体亮度温度TBB可表示有云区云顶或无云区地面的温度。在云区,一般TBB≤0℃,可以根据云顶亮温的演变推断天气系统的强度以及可能引起的相应天气现象。 (1) 通过率r，可以表示环境衰减对太阳辐射的衰减作用 τ=QEQ0(2)τ=QEQ0(2) 式中QE为太阳辐射经大气衰减到达地面的太阳辐射能;Q0为大气顶的太阳总辐射。 (2) 传感器的能量 卫星传感器接收到的能量分为两部分:一部分是太阳辐射未经过大气衰减的能量,一部分是经过大气衰减后到达地面又反射的能量(图1)。 αTOAQ0=Q0ταsτ+(1-τ)Q0(3) 式中αTOAQ0为传感器接收到的能量;Q0τ αsτ为经过大气衰减到地面,又经过地面反射和衰减之后的能量;(1-τ)Q0为在云顶就被反射的能量。由式(3)可以求解出大气透过率的公式为 τ=1?1?4αs(1?αTOA)√2αs(4)τ=1-1-4αs(1-αΤΟA)2αs(4) 式中αTOA和αs都可以由卫星资料观测获得。 通过式(4)可以求出逐时的透过率,阴天条件下到达地面的太阳总辐射为 QE=(1-α)Q0τ(5) 需要说明,通过以上公式得到的大气透过率,并不是完整的大气对太阳辐射的透过率,表示的是从大气层外直到地面的辐射传输路径上,地气系统总的等效消光作用。 1.2 不同覆盖物参数 G0=Rn[Γc+(1-fv)(Γs-Γc)] (6) 式中Γc为植被覆盖区参数;Γs为裸土区参数;植被覆盖度fv用经验公式进行计算。 1.3 计算热敏度的方案 本文采用迭代法计算感热通量,具体计算方法详见文献。 1.4 土壤水分充分供应的温湿度4. 蒸发比是蒸散发量和有效能量(净辐射与土壤热通量之差)的比值。可表示为 Λ=LERn?G(7)Λ=LERn-G(7) 式中Λ为蒸发比;LE为潜热通量;Rn为净辐射通量;G为感热通量。 在土壤水分亏缺的干燥地表环境下,由于无可用土壤水分供给蒸散发,即 λEdry=Rn-G0-Hdry≡0 或Hdry=Rn-G0(8) 式中Hdry为干燥地表环境下的感热通量。而在土壤水分充分供应的湿润地表环境下,为 λEwet=Rn-G0-Hwet或Hwet=Rn-G0-λEwet(9) 式中Hwet、Ewet分别为湿润地表环境下的感热通量和潜热通量。则相对蒸发比为 Λr=λEλEwet=1?λEwet?λEλEwet(10)Λr=λEλEwet=1-λEwet-λEλEwet(10) 由式(7)～式(10),蒸发比最终表示为 Λ=λERn?G0=ΛrλEwetRn?G0(11)Λ=λERn-G0=ΛrλEwetRn-G0(11) 土壤水分充分供应的湿润地表环境下的潜热通量和感热通量通过结合Monteith提出的 Penman-Monteith公式得出。 因为蒸发比全天基本保持稳定,若已知某一时刻的净辐射和土壤热通量,就可以求解出该时刻的潜热通量,通过汽化潜热系数,即可算出蒸散发量。 2 地面气象观测资料及预处理 黄河源区(95°30′E～103°30′E,31°30′N～36°30′N)位于青藏高原东北部,面积约12.2万km2,横跨青海、四川、甘肃三省,绝大部分海拔为3 000～5 000 m, 该地区陆面蒸散发量的大小与黄河水资源变化有密切的关系。黄河源区年平均气温为1.2 ℃,为典型的高原大陆性气候。图2为黄河源区土地利用分布情况,区域内植被类型以高原灌丛草甸为主,间以零星灌木丛。左上