【精品硕博论文-水利水电工程】基于卫星遥感地表温度与植被指数特征空间的陆面蒸散发估算(002)

基于卫星遥感地表温度与植被指数特征空间的陆面蒸散发估算1遥感蒸散发2采用的模型3研究区域及采用的数据4具体的数据处理5结果验证6存在的问题综述遥感蒸散发蒸散发既是地表水分循环的重要组成部分，同时也是地表能量平衡方程的主要项。遥感方法以其能够获知大范围地表特征信息的优势，为较准确估算地表蒸散发提供了可能，从而受到人们的日益重视采用的模型本文采用的是SEBAL模型。SEBAL是基于陆面能量平衡方程的一种模型，是仅仅利用遥感影像的可见光、近红外以及热红外波段以及少量的气象数据，比如风速、气温，来计算地表净辐射量、土壤热通量、感热通量以及用于蒸散发计算的瞬时潜热通量。RN地表净辐射，为地表净得短波和长波辐射之和H下垫面到大气的感热通量G下垫面地表的热通量，即下垫面土壤中的热交换，土壤热通量LE潜热通量，即下垫面与大气间水汽的热交换S、PH分别为作物光合作用和生物量增加的能量，比重小，可以忽略研究区域及采用的数据研究区域安徽省蚌埠市五道沟试验站，地处东经11721，北纬3309，是河海大学水资源开发教育部重点实验室，其气象数据精度能达到10MIN的精度（如蒸发数据、气温数据等）采用的数据本文采用的数据是2013年2月11号发射的LANDSAT8数据，LANDSAT8携带了2个传感器，包括OLI（陆地成像仪）和TIRS（热红外传感器）OLI与ETM对比OLIETM序号波段M空间分辨率M序号波段M空间分辨率M10433045330204500515301045005153030525060030205250605304063006803030630069030508450885304077509003061560166030515501750307210023003072090235030805000680158052009001591360139030蓝色波段BAND104330453M主要应用海岸带观测，短波红外波段BAND913601390M包括水汽强吸收特征可用于云检测。具体的数据处理现有的地表温度反演算法大致有以下四种大气校正法、单通道算法、劈窗算法和多波段算法。大气校正法，需要大气实时剖面数据，数据获取比较困难，而多波段算法适合于多个热红外波段的数据，所需参数多，运算复杂，反演难度较大单通道算法适用于一个热红外波段，而劈窗算法适用于两个热红外波段，因而本文采用单通道以及劈窗算法进行温度的演算。名称优缺点优点不足JIMNEZMUNOZ和SOBRINO单通道算法1既考虑了地表比辐射率的影响,也考虑了大气辐射的影响2反演过程所需要的大气参数仅为大气水分含量1模型的普适性没有得到验证。2模型的过程参数依赖过多经验性公式。劈窗算法1既考虑了地表比辐射率的影响,也考虑了大气辐射的影响2反演过程所需要的大气参数比传统的辐射传导方程法要少得多,仅需要近地表气温和大气水分含量这两个参数1模型的普适性没有得到验证。2模型的过程参数依赖过多经验性公式。3基于的假定条件较多。数据处理预处理流程为纬度，可由MODIS一级产品数据MOD03获得；S为坡度，R为坡向，均可由DEM得到；和分别代表太阳赤纬和太阳时角相关能量值的计算为斯一玻常数其中G8，OBS为008M处2块土壤热流板观测值的平均值，T，T分别为TCAV平均热电偶观测得到的平均土壤温度及其时间分辨率30MIN；Z为ZR之上土壤厚度，为008M（可以根据实际情况进行调整）TCAV平均土壤热电阻法感热通量H其中Z1和Z2分别代表近地表的两个高度01米和2米，U为摩阻速度MS1。L0和L0分别代表的是大气稳定及不稳定两种状况，G为重力加速度，98MS2。通过一个迭代计算过程对空气动力学阻力RA进行稳定度修正，需要通过多次5次以上重复计算,得到稳定的显热通量结果验证地表温度日期日序地表温度（K）劈窗算法（K）单通道算法（K）201324630454300853044420132783033430287304592014121300822970229740相比较而言，单通道算法的计算值与实际值比较接近，而单通道算法应用的是波段10，原因是因为波段11的噪声比较大；因此应用了2个波段的劈窗算法结果有些偏差。由于受到云层，以及LANDSAT8时间分辨率的影响。自2013年7月至2014年5月仅仅选择出3期可用图像。DOY246劈窗算法与单通道算法温度对比空间分布图劈窗算法与单通道算法温度对比空间分布图DOY278劈窗算法与单通道算法温度对比空间分布图DOY121辐射通量RN日期过境时间日序RN计算RN实测绝对误差均方根误差2013090310450024657683254927927553759168201310051045002785338935550262113344660420140501104500121651696601005506912569577平均3312635475土壤热通量日期过境时间日序G计算G测量值绝对误差均方根误差201309031045002467880857430213784570272013100510450027871700716590041000220140501104500121633006332800280001平均值714912343潜热通量日期日序LE实测LE计算绝对误差均方根误差201309032463555884118005621331598472013100527813208821811002407557959620140501121443685394027496582465880平均4331545480潜热通量与日蒸散发量对比潜热通量与日蒸散发空间分布图存在的问题误差来源显热通量计算中冷热点选取的好坏直接影响着最终精度，冷热点的选择均需要一定的技能和经验性，主观性强，存在较大的人为误差利用能量平衡方程,潜热通量为净辐射与土壤热通量、感热通量的差值,而通量又可以间接利用遥感反演的参数值和模型参数求值。由于参数值大都需要通过遥感数据获得,同时地表温度、比辐射率、植被指数等一些重要参数至今仍是遥感反演的瓶颈,所以方程应用存在一定的适用范围和误差。能量不闭合的问题是观测手段的缺陷日期日序H实测LE实测RN实测G测量值能量闭合度LEH/RNG20130903246129809433874985492791574300095262013100527818003511320882557726372900006438201405011216678993341795625980593620007532仅靠3期图像的计算结果，并不能可靠的说明SEBAL模型在五道沟水文水资源试验站的适用性，因此还需要引入其它时间分辨率高的遥感影像来分析，例如HJ1（国产环境小卫星），MODIS等。遥感估算蒸散发，是水文模型中模拟地表水循环的一部分，水文模型需要长时间序列的