利用分裂窗算法反演陆地-气象科学.pdfVIP

维普资讯 第 22卷 第 4期 气 象 科 学 Vo1．22，NO．4 2002年 12月 SCIENTIA M ETEOROLOGICA SINICA Dec．．2002 利用分裂窗算法反演陆地 表面温度的研究进展① 俞 宏 石汉青 (解放军理工大学气象学院 。南京 211101) 摘 要 本文对LST反演的分裂窗算法作了回颐和 比较 。LST反演的分裂窗 算法最一般的形式是Ts—T+A(T 一T)+B，其 中T 是陆地表面温度 ，T 和 T 分别是 AVHRR分裂窗通道 4和 5的亮温度 ，A和 B是由大气、观测角及地表辐射 率所决定的系数，不同的算法其系数的确定有很大的不同。地表 比辐射率的确定是 LST反演 中的主要难点，目前也有多种方法来确定地表 比辐射率值 。最后本文还指 出了目前在 LST反演中存在的局限性 。 关键词 分裂窗算法 反演 陆地表面温度 分类号 S165 引 言 陆地表面温度 (LST)是一个重要的地球物理参数，它在地一气间的物质与能量交换的过程 中起着重要的作用 。同时，它也是地球表面能量平衡和温室效应的一个很好的指标 ，是区域和 全球尺度地表物理过程 中的一个关键因子 ，陆地表面温度在气候、水文 、生态、环境、生地化学 等研究中都有重要应用n ]。分裂窗反演法在针对某些 已知地表 比辐射率的地表进行局地反 演时，可达到较高的反演精度 ，且有较高的反演效率，是 目前应用非常广泛的陆地表面温度反 演法 。关于陆地表面温度的分裂窗反演算法有很多，本文将对这些算法进行较为全面的评述 。 1 陆地表面温度反演的分裂窗算法 由于NOAA／AVHRR的分裂窗通道 4(10．3m一11．3m)和通道 5(11．5m一12．5m) 的资料容易获得，因此在陆地表面温度的计算中被广泛地应用[6-lo]，近年来 ，EOS／MODIS通 道 31(1o．78m一11．28#m)和通道 32(11．77m一12．27#m)的资料也被用来进行陆地表面温 度 的分裂窗反演 ]，此外，也有利用 GMS5／VISSR两个红外分裂窗通道反演陆地表面温度的 报道n 。近二十年来 已发展了利用AVHRR通道 4和通道 5反演 LST 的多种分裂窗算法 。尽 管这些算法在计算方式上有很大不同，但根据它们的表达形式，基本上可分为两类 。下面就多 种不同分裂窗算法进行 比较，以对它们各 自的特点及之间的差异有一个全面的了解 。 ① 收稿 日期：2002～O1—17；修改稿 日期 ：2002一O4～17 维普资讯 4期 俞 宏等 ：利用分裂窗算法反演陆地表面温度的研究进展 495 1．1 分裂窗算法的一般形式 研究反演算法的基础是辐射传输方程。在 1om～14p．m的热红外波段 ，假定：①晴空大气 条件 ；②地球表面为郎伯体 ；③大气满足局地热力平衡条件 ；④忽略大气分子和气溶胶的散射， 则卫星观测到 的辐射表示为 ： L (丁 ，)= ￡r() ·B (丁 )+ L ()+ (1一 ￡) () ·L (1) 式中左边第一项是卫星观测辐射，T 是波长为 入时的亮温值 ， 是波长入下的地表 比辐射率， (e)是在观测角 0方 向上从地面到卫星高度 Z之 间的大气透过率 ， (T)是波长 入、地表温度 T 所相应的黑体辐射率， (e)是到达卫星的在 0方 向的大气 向上发射辐射 ， 是大气 向下 发射辐射 。分裂窗反演算法的研制就是基于上述传输方程 。 假设 T 、T 分别为AVHRR通道 4和通道 5的亮温值 ，那么分裂窗算法的一般形式可用 下式来表示 ： Ts— T4+ A (丁4