1、 南阳师范学院 20XX 届毕业生 毕业论文（设计） 题 目： 绿洲与周边荒漠地表温度差异的时序变化 完 成 人： 班 级： 学 制： 专 业： 地理信息系统 指导教师： 完成日期： 目 录 摘要 (1) 引言 (1) 1 文献综述 (2) 2 研究区域概况 (2) 3 研究区域概况研究方法 (3) 3.1 基于 MODIS 温度产品的 LST 获取 (3) 3.2 基于 Landsat TM 波段的地表温度（ LST）获取 (3) 4 结果与分析 (5) 4.1 LST 的时间变化 (6) 5 讨论 5.1 地表温度与植被状况的关系 (7) 5.2 绿 洲冷岛效应对农业气候的意义 (8) 6
2、 结 论 (8) 参考文献 (9) Abstract (10) 绿洲与周边荒漠的地表温度差异的时序变化 摘要： 本文以外流河绿洲卫宁平原为例，利用 MODIS 卫星遥感数据探索绿洲与其外围区域 LST 差异程度的时间变化规律。结果显示： 1）绿洲与外围区域间地表温度差异存在显著的时段是 4 月中旬至 10 月上旬；而此时段之外绿洲与外围区域间不存在显著的地表温度差异 ； 2）绿洲的“冷岛”结构造成了绿洲周围荒漠非常大势力差异，这会在绿洲和临近荒漠之间形成局地 势力环流。 关键词： 绿洲；地表温度； MODIS； landsat TM； 归一化植被指数 (NDVI) 引 言 地表温度 是地球表面
3、能量传输的重要因子，是地 -气间相互作用及能量交换的关键参数，在气候、水文、农业、生态以及环境等研究中重要意义 。 地表温度的时空变化在一定程度上影响了水分、热量等条件的变化，并进而对位于其中的人类社会以及生产生活产生了影响 。 MODIS 它包 含 36 个波段 ,分别针 对陆地、 海洋、水汽、气溶胶等来设计的 。 MODIS 的第 31 和 32 波 段 和 TM 的第 6 波段 常被用 来反演地表温度。 MODIS 地面分辨率分别为 250m、 500m、 1000m， 因此 MODIS适用于反映大区域 地表 温度的空间分布规律。 绿洲是自然界一种独特的地理区域类型 1,2，并在全球广泛
4、分布。它是人类文明发祥地的重要组成部分。绿洲是自然界唯一能抗御干旱气候环境的特殊生态环境系统之一。 由于绿洲对自然界水资源和光热能源的巧妙组合，使其拥有了得天独厚的生态技能和生产效率。 并且绿洲生态系统在对抗干旱气候环境的过程中，“绿洲效应”在某些方面也发挥了比较重要的作用。也正因为如此，我国西北地区面积不到5%的绿洲哺育着该地区的 95%以上的人口。 正是绿洲内有了人类 活动，逐步完善了绿洲的灌溉系统、改进了绿洲的植被类型、健全了绿洲的林网系统才使绿洲的结构更加合理生态能力更强。但人类活动对绿洲是双刃剑，它既创造了绿洲的辉煌也断送了绿洲的未来。譬如人类对绿洲相关区域环境 的干预（如山区涵养林
5、的破坏）以及全球范围内的人类活动引起的全球变暖所造成的山区冰雪融化和荒漠化扩大等。而通过对绿洲的小气候特征分析，发现存在一种“绿洲冷岛效应” 3， 4。而这种绿洲冷岛效应是由于绿洲下垫面热力性质与沙漠不同造成的。 本文尝试从以往研究成果中概括和提炼出绿洲在地理和气候效应方面的一些本质特性 。 近年来 ,随着遥感技术的快速发展 ,国内不少学者 在地表温度研究领域有了一定的成果并广泛的应用在一些研究领域，如 预测森林火灾、作物的长势、农田旱情监测、生态评估和城市环境监测等方面。 Gallo 等 人 分析了 NDVI 与利用遥感数据反演的地表温度之间的关系；戚鹏程 、覃志豪、白洁、孟宪红 等人利用
6、Landsat TM 的热红外数据反演得到地表温度来研究土地利用变化 与地表温度的关系 5-8 。高志强 、苏伟忠 等 人 利用遥感分类和定量反演技术，提取两期土地利用 /土地覆盖（ LUCC）信息和对应的地表温度信息 9,10 。 1 文献综 述 近年来 , 随着遥感技术的快速发展 ,国内不少学者利用 MODIS 等卫星遥感资料进行的 NDVI 和 LST 变化及应用方面的 研究领域有了一定的成果并广泛的应用在一些研究领域，如 预测森林火灾、作物的长势、农田旱情监测、生态评估和城市环境监测等方面 ，从而为防灾减灾提供了重要的参考依据 。 Gallo 等 人 分析了 NDVI 与利用遥感数据反
7、演的地表温度之间的关系； 戚鹏程 、覃志豪、白洁、孟宪红 等人利用 Landsat TM 的热红外数据反演得到地表温度来研究 土地利用变化与地表温度的关系 。高志强 、苏伟忠 等 人 利用遥感分类和定量反演技术，提 取两期土地利用 /土地覆盖（ LUCC）信息和对应的地表温度信息。 刘慧敏等对 MODIS与 AVHRR 数据的植被变化检测进行了比较。何彬方 11等利用 MODIS NDVI 对卫宁平原土地覆盖进行了分类研究。Price(1984)12最先把分裂窗算法 推广到陆地面积的反演，他通过引入比辐射率改正项来减少因陆地表面比辐射率变化而引起的误差。Beeker(1987)从理论上证明了用
8、分裂窗技术反演地表温度的可行性，并且第一次从理论上给出了使用分裂窗技术时的大气和比辐射率对地表温度反演的影响。由于自然界绿洲的千差万别，对绿 洲给出一个能高度概括其共性的科学定义在理论研究方面是很有重要意义的。有几个 学者从不同的角度给绿洲做了一些定义。比较有概括在干旱性的在我国有韩德麟给出的“绿洲是干旱地区人类对特殊自然环境的生态 、 景观。”在西方有 张强等人 给出的“绿洲是荒漠中的稠密植被的地方” 3 。 贾保全曾提出绿洲是在干旱气候条件下形成的在荒漠背景基质上以天然径流为依托 小尺度非地带性景观。 在 本文利用 MODIS遥感数据，对卫宁平原 地面温度的季节和空间分布特征进行了详尽分析
9、。 利用卫星资料进行的区域性 NVDI 和 LST 变化及其应用方面的研究已取得较 多进展，在农业生产和布局、生态系统、地表资源环境进行宏观、大范围、动态连续监测和研究中得到广泛应用，从而为减灾提供重要参考依据 13,14。 2 研究区域概况 卫 宁 平原 位于 北纬 35 14-39 23、 东经 104 17-1073 9之间，处在中国西部的黄河上游地区 ，展布于黑山峡至青铜峡间黄河干流的两侧 。 宁夏东邻陕西省，西部、北部接内蒙古自治区，南部与甘肃省相连。自古以来就是内接中原，西通西域，北连大漠， 与其东北的银川平原共同组成宁夏平原。长期以来，人类在此汲取黄河水灌溉农田，形成了人工绿洲（
10、下 简称卫 宁 绿洲），成为宁夏的重要 粮食产区 。卫 宁 绿洲西北侧为腾格里沙漠，分布着以格状沙丘为主的多种沙丘。 地貌复杂，山地迭起，盆地错落，大体可分为：黄土高原，鄂尔多斯台地，洪积冲积平原和六盘山、罗山、贺兰山南北中三段山地。平均海拔 1000m以上。按地表特征，还可分为南部暖温带平原地带，中部中温带半荒漠地带和北部中温带荒漠地带。全区从南向北表现出由流水地貌向风蚀地貌过渡的特征。 部分地表已发生沙化。气候为大陆兼沙漠性气候，据 1961-1990 年站点资料，年均日照时数 2878 h，年均气温 8.5 ， 10 积温 3319 ，年均降水量 仅 182 mm，年均蒸发量 达 300
11、1 mm。本文研究区包括卫 宁 绿洲及其外围区域（绿洲周边 15km的环形 区域）。 3 研究区域概况 研究方法 3.1 基于 MODIS 温度产品的 LST 获取 使用 MODIS 3级产品 （ MOD11A2） 获取 年内多时间点上的 LST 信息。从 2002 年的 MOD11A2 数据集中得到 “ 8 日 平均 合成 ” LST 图像47 幅。由于云对热红外成像过程的影响，各幅 MODIS LST 图像有不同程度的缺失。根据 MOD11A2 数据集中的质量控制图层（ Quality Control Layers）， 将缺失像元数占研究区总像元数 1/3以上的 MODIS LST图像排除
12、于研究之外。根据 “ 绿洲 ” 和 “外围区域” 矢量多边形的范围分别统计各时期绿洲和 外围 的 LST 平均值，同时从质量控制图层中分别统计相应误差的平均值 。 3.2 基于 Landsat TM 波段的地表温度（ LST）获取 获取 基于 Landsat TM6波段的地表温度（ LST） 的方法有多种，如普适性单通道算法、大气校正和单窗算法等。本文将采用覃志豪等以地表热辐射传导方程为基础，通过应用 Taloy公式对 Planck函数线性展开，对贡献相近的两个大气平均作用温度参数合 并，以及对辐射亮度 L6和亮度温度 T6进行线性回归等方式。经以往研究表明，单窗算法有较高的精度，在大气透射率
13、、大气平均作用温度和地表比辐射率估计有中等误差的情况下，基于单窗算法的地表温度反演误差约为1.2K6 。 信息源为 2009年 6月 25日的 TM影像的第 6波段以及 3、 4波段。首先，将原始数据中的第 6波段 DN值按下式进行转换，换算为大气层顶进入传感器的辐射亮度 L6 )/( 2 msrmW 6L = Gain DN + Bias （ 1） 式中， DN为 TM6波段的 DN值： Gain和 Bias为定标参数，从各期数据的头文件中获得。 用下式求得亮度温度 T6（ K）。 T6= 1260.56 / ln(607.76 / L6+1) （ 2） 利用下式求算地表温度 （） （ 3）
14、 式中 aT 为大气平均作用温度，在无实时探空资料时，可用地面附近气温的线性函 数近似 2 092621.00110.16 TT a （ 4） 式中 0T 为地面附近气温。 公式（ 3）中的 C和 D为中间变量 66C （ 5） )1(1)1( 666 D （ 6） 因此，只要知道参数 aT 、 6 和 6 ，即可用上述单窗算法推算任何像元的实际地表温度。（式中， 6 为大气透射率；为地表比辐射率。 6 、 6由下式计算 w1 1 5 3 6.00 3 1 4 1 2.16 （ 7） svvv PP )1(6 （ 8） 式中 w 为大气总水分密度， 1. 44 2/cmg ，可据近地面水汽压推
15、算，公式和参数见文献； v 为植被比辐射率； s 为 祼土比辐射率，分别取值为 0.986 和 0.970； vP 为植被覆盖度，可据 NDVI 推算。 NDVI=（ b4-b3） /（ b4+b3） （ 9） 式中， b4、 b3分别是 TM4、 TM3波段的 DN值。 NDVI值越大，地表越接近于完全的植被叶冠覆盖； NDVI值越小，越接近于完全裸地；而 DNVI介于植被于裸地之间时，则表明有一定比例的植被叶冠覆盖和一定比例的裸土。因此，可以 用以下公式确定各像元的植被覆盖度，即植被构成比例 vP , )/()( svsv N D V IN D V IN D V IN D V IP （ 1
16、0） 式中， vNDVI 和 sNDVI 分别是植被和裸土的 NDVI值。 4 结果与分析 4.1LST的时间变化 图 1 表 示 2002 年各时段（如 1月 1日 至 1 月 8 日，依次类推）绿洲和 外围区域 地表温度平均值，以及从 质量控制 图层 中 统计得到的误差平均值。从图中可见，绿洲与 外围区域 的 LST都表现为 “ 冬低夏高，春升秋降 ” 的规律。图 2中亦可见，从 4月中旬至 10 月上旬，绿洲与 外围区域 之间的 LST 差异 （下简称 LST 外围 -绿洲 ） 超过 MODIS LST产品误差的范围，说明这种差异是显著的，是自然规律的表现；而 1月至 4 月上旬以及 1
17、0 月中旬至 12 月，绿洲与 外围 之间的 LST 差异在MODIS LST 产品误差的范围之内，表明二者之 间的 LST 不存在显著差异 2002202402602803003203402002年1月1日2002年1月15日2002年1月29日2002年2月12日2002年2月26日2002年3月12日2002年3月26日2002年4月9日2002年4月23日2002年5月7日2002年5月21日2002年6月4日2002年6月18日2002年7月2日2002年7月16日2002年7月30日2002年8月13日2002年8月27日2002年9月10日2002年9月24日2002年10月8日
18、2002年10月22日2002年11月5日2002年11月19日2002年12月3日2002年12月17日日 期地表温度 K荒漠地表温度绿洲地表温度图 1 不同日期绿洲及其外围 区域 的地表温度 -2024681012地表温度差1 25 49 73 97 121 145 169 193 217 241 265 289 313 337 361日 期荒漠与绿洲的地表温度差图 2 绿洲 -外围区域地表温度差 5 讨论 5.1 地表温度 与植被状况的关系 以 卫 宁 绿洲及其外围区域 为研究区域， 基于 Landsat TM 影像 分别 获取 了 陆地表面温度（ LST）以及归一化植被指数 (NDVI
19、)，定量分析了 LST与 NDVI 的 关系。结果表明： LST 和 NDVI具有明显相反的变化趋势， LST 与 NDVI 平均值也具有显著的差异性 ； NDVI 越高，则LST越低。 通常认为， NDVI 的 季节变化 是十分 显著 的 ， 并 具有典型的地域特征，受当地气候影响，植被、农作物类型地域差异较大， 那么NDVI 的变化是否会对 LST 有明显的影响呢？干旱区的 植被覆盖程度的年内变化往往是很大的，其年际变化通常会很大，那么，这种明显的植被变化很可能会影响到地表温度的空间分异，从而影响绿洲的冷岛结构，进而影响绿洲冷岛效应。 285290295300305310315320-0.
20、20 0.00 0.20 0.40 0.60 0.80NVDI地表温度(K)图 3 地表温度 与 NDVI 之间的关系 5.2 绿洲冷岛效应 对 农业 气候 的 意义 西北干旱地区的沙漠在夏季长时间处在晴天或少云天气条件下，造成酷热和极干旱的环境，白天层结构不 稳定，湍流强烈发展，而在这时候绿洲冷岛效应内却处于逆温稳定层结，抑制湍流的发展，这是有利于植物生长的，因为绿洲 这种气温较低，湍流发展较弱的凉湿环境必然抑制植物的蒸腾过程和地面的蒸发过程，在干旱酷热的戈壁或沙漠之中蒸腾和蒸发的抑制，对于节约用水和植物或作物本身的生长都是极为有利的。就是说，只要在干旱的沙漠中绿洲的范围足够大，能形成绿洲冷
21、岛，那么就形成一种有利于植物生长的小气候，这一结论解释了沙漠中种树种草小片不易成活，大片易于成活的道理。自然这一结论有利于指导我们在沙漠或戈壁等干旱地区的种树种草活动和农业开发 。但这只是对冷岛 效应的一些初步研究结果，无论是实验还是理论都有待于进一步的深入研究。 6 结论 地表温度具有强烈空间异质性，单纯依靠地面实测难以对研究区的地表温度进行准确的分布式描述。而遥感资料在反映区域分异方面具有明显优势，可以 对 地表温度特征进行平面化描述，使研究者易于深入 观察绿洲与外围区域的地表温度分异情况 。 本文主要结论如下：1）绿洲与 外围区域 间 地表温度存在显著 差异 的时段是自 4 月中旬 起，
22、10 月上旬 止；而此时段之外 绿洲与 外围区域 间 不存在显著地表温度差异 。 2）土地开发利用或保护恢复措施会影响绿洲大小和形状， 从而影响冷岛结构；认为土地利用规划和管理中应将不同利用方式对冷岛结构的影响考虑在内。 参 考 文 献 1沈玉凌 . 绿洲 概念小议 J . 干旱区地理 , 1994, 17( 2) : 70-74. 2Gao Huajun. T he distribution and types of oasis in China J . AridLand Geography, 1987, 10( 4) : 23- 29. 高华君 . 我国绿洲分布与类型 J . 干旱区地理
23、, 1987, 10(4) : 23-29. 3张强，胡隐樵 ， Walter. 绿洲地理特征及其气候效应 J. 地球科学进展，2002， 17(4)： 477-486. 4苏从先，胡隐樵 . 绿洲和湖泊的冷岛效应 J. 科学通报， 1987， 10(3)：756-758. 5戚鹏程 ， 田国良 . 热红外遥感 M. 北京：电子工业出版社， 2006. 6覃志豪 , Zhang Minghua, Karnieli Arnon,等 . 用陆地卫星 TM6 数据演算地表温度的单窗算法 J. 地理学报， 2001， 56(4)： 456- 466. 7白洁，刘绍民，扈光 . 针对 TM/ETM+遥感数据的地表温度反演与验证 J. 农业工程学报， 2008， 24(9)： 148-154.
