唐古拉山冬克玛底地区冰川变化遥感监测谯程骏（中国科学院寒区旱区环境与工程研究所寒旱区流域水文及应用生态实验室，甘肃兰州730000）摘要利用1973年MSS 、1992年TM 、2001年ETM +、2007年TM 数字遥感影像资料以及数字高程模型（DEM ），结合第一次冰川编目资料，通过地理信息系统和遥感图像处理技术提取了唐古拉山冬克玛底地区不同年份的冰川分布范围。
结果表明：研究区内的冰川面积1976年比1969年减小了0.84%，1992年冰川面积比1976年的减小了1.77%，2001年冰川面积比1992年减小了4.24%，2007年冰川面积比2002年减小了3.31%。
结合气象资料和有限的野外资料进行分析，该地区冰川退宿趋势同气温升高趋势基本一致。
1969年到2007年近40年时间里，冰川面积总共减小17.61km 2，占1969年面积的9.81%，冰储量减少5.75km 3，相当于水资源损失5.17km 3。
关键词唐古拉山；冰川变化；遥感；地理信息系统中图分类号P343.6文献标识码A 文章编号0517－6611（2010）14－07703－03Remote Sensing Monitoring of Glacier Changes in Dongkemadi Region of Tanggula MountainQIAO Cheng-jun （Laboratory of Watershed Hydrology and Applied Ecology in Cold and Arid Regions ，Cold and Arid Regions Environmental and Engineering Research Institute ，Chinese Academy of Sciences ，Lanzhou ，Gansu 730000）Abstract Based on the numeric remote sensing image data （MSS data in 1973，TM data in 1992，ETM +data in 2001，TM data in 2007）and nu-meric terrain model ，combined with the first glacier inventories data ，the glacier distribution range in Dongkemadi Region of Tanggula Mountain in different years were extracted by using GIS technology and remote sensing image processing technology.The results showed that the glaciers area in the study area had decreased by 0.84%，1.77%，4.24%，3.31%from 1969to 1976，1976to 1992，1992to 2001and 2001to bined with the meteorological data and limited field data ，the analysis results showed that glacier shrinkage in the study area was mainly subjected to continuous temperature rise.From 1969to 2007，the glacier area had decreased 17.61km 2，being 9.81%of that area in 1969.Theglacier volume had decreased 5.75km 3，which meant that water resources lost 5.17km 3.Key words Tanggula Mountain ；Glacier change ；Remote sensing ；Geographic information system基金项目国家重点基础研究发展计划资助项目（2007CB411502）；水利部公益性行业科研专项经费项目（2007SHZ1-46）资助。
作者简介谯程骏（1983－），男，四川什邡人，在读硕士，从事水文水资源以及冰川变化的研究。
E-mail ：qiao _cheng _jun @.鸣谢感谢美国马里兰大学和美国地质调查局提供MSS /TM /ETM +遥感影像数据以及DEM 数据；感谢中国气象局提供了气象数据。
收稿日期2010-03-02冰川是自然界重要而且极具潜力的淡水资源，它覆盖了地球10%的陆地面积，而且世界上有80%的淡水储存于冰川。
冰川是气候的产物，又是全球气候变化灵敏的指示器，对气候变暖高度敏感。
在气候变暖的影响下，冰川处于持续的负物质平衡而不断退缩［1］，作为内陆盆地水源以及外流水系源头的冰川，其在不同时间空间尺度的变化势必引起以冰川融水补给为主的河流水量的丰枯变化以及长期水量的不断减少，作为一种水资源，冰川还对区域生态环境的平衡与稳定起着至关重要的作用［2］。
青藏高原是中低纬度现代冰川较为发育的地区，20世纪90年代以来青藏高原冰川呈现退缩状态，但各地冰川变化具有明显的空间特征［1］。
由于大部分山地冰川地处偏远，并且数量众多，因此，卫星遥感技术被广泛应用于大尺度的冰川变化研究中。
陆地资源卫星（包括MSS ，TM 及ETM +）已经成为冰川学研究的主要数据源之一［3］。
综合利用遥感丰富的数据资源和地理信息系统技术高效的空间数据处理功能，可以有效地获取和分析冰川变化信息。
采用地理信息系统（GIS ）与遥感（RS ）技术，利用航空和卫星遥感资料，结合冰川编目资料和数字高程模型（DEM ），通过综合分析获得唐古拉山中段冬克玛底地区最近40年冰川变化情况。
最后结合气象观测数据分析了其对气候变化的响应。
该研究区冰川变化在指示青藏高原中部冰川变化与气候变化方面具有重要意义。
1研究区概况研究区域位于青藏高原腹地唐古拉山中段山区冬克玛底地区（图1），在青海省和西藏自治区的交界处，海拔均在5000m 以上，冰川面积共计179.43km 2（来源于1969年航测资料），冰川平均粒雪线海拔5530m ，发育的冰川主要为小型冰斗冰川和山谷冰川［4］，山谷冰川谷底为古冰川作用形成的平坦开阔的稀疏草地。
植被为单一的短草植物，生长最盛时高度约为3 5cm 。
年平均气温为－6.0ħ，气温年较差为24.9ħ，全年只有6 9月平均气温在0ħ以上。
年平均相对湿度为65%，降水集中于6 9月［5］。
研究区内的冬克玛底冰川是国内少数几条长期观测的冰川之一，已有的观测资料时间序列较长，有助于结合遥感资料共同分析该区域内冰川变化情况。
2数据源及处理方法笔者使用的数据如下：①中国冰川目录Ⅷ－长江水系，研究区冰川分布图来自1969年拍摄的航空相片，可以反映1969年的冰川分布情况。
②遥感数字影像资料为1976年MSS （Landsat Multispectral Scanner ，陆地资源卫星多光谱传感器，行列号为148，37），1992年和2007年TM （Landsat Themat-ic Mapper ，陆地资源卫星专题制图仪，行列号为137，37），2001年ETM +（Landsat Enhanced Thematic Mapper Plus ，陆地资源卫星增强型专题制图仪，行列号为137，37）。
③研究区的数字高程模型（DEM ，Digital Elevation Mode ，分辨率为30m ）。
④距离研究区最近的，有长时间序列的国家气象台站安多站气象资料。
冰雪在可见光及近红外波段（TM1 TM4）反射率都高，在这些波段中冰川的光谱值范围较宽，不利于冰川与背景及其他地物区别，而冰川在TM5中光谱值范围很窄，使用该波段有利于提取冰川信息；另外，水体会吸收中远红外波段，使冰雪在TM5波段光谱值降低，且明显低于冰川中冰碛的光安徽农业科学，Journal of Anhui Agri.Sci.2010，38（14）：7703－7705责任编辑胡剑胜责任校对况玲玲谱亮度值［6］。
因此，采用TM和ETM+5、4、3波段合成假彩色影像，便于冰川边界的提取。
对于MSS影像，则直接采用4、2、1波段合成假彩色影像。
ETM+的8波段为分辨率为15 m的全色波段，将该波段与合成的5、4、3假彩色影像融合，可大幅度提高冰川的目视解译精度。
由于地形会影响当地太阳入射角和高度角的入射角度，在影像中形成阴影或者地形起伏不明显，导致冰川信息提取精度降低，利用该地区的数字高程模型对数字遥感影像进行了地形纠正处理，纠正后的遥感图像清晰度提高，并且可根据地形起伏情况进一步确定冰川范围。
遥感数字影像资料和数字高程模型均采用统一的UTM46带投影和WGS84椭球体建立坐标系统。
遥感数字影像资料经上述方法处理后，以目视判读方法来识别各期冰川的边界，由于影像分辨率的局限性和人工目视判读不可避免存在误差，尽量将这种方法判读的冰川边界控制在一个像元之内。
最后在地理信息系统（GIS）支持下完成冰川范围的矢量化，得到1976、1992、2001和2007年4期冰川范围数据，将4期数据同冰川编目资料1969年的数据复合，得到5期资料的叠加数据集（图2），分析和比较不同时期的冰川变化。
注：左边为2007年Landsat TM影像，RGB：5、4、3。
Note：The left figure is Landsat TM image in2007，RGB of5，4，3.图1研究区地理位置Fig.1The geographical location of the studyarea图2冬克玛底地区各时期冰川分布叠加Fig.2The overlay of glacier distribution in Dongkemadi area in each period3结果分析通过GIS软件的相关量算（表1），可看出从20世纪60年代末到2007年的近40年时间内，研究区冰川面积减少了17.61km2，占1969年总面积的9.81%。
1969 1976年，冰川面积变化缓慢，而从1992 2007年，冰川面积减小速度加快。
与研究区相距不远的各拉丹冬冰川［7］在1969 2000年面积仅变化1.7%，可能是由于该研究区1969年冰川总面积只有格拉丹东地区冰川面积的1/5，冰川面积都相对较小，导致对气候变化的响应比较敏感，变化也较大。