塔里木河下游应急输水前后生态变化遥感监测分析

1、塔里木河概况

　　塔里木河（以下简称塔河）是我国最长的内陆河，从它的源头叶尔羌河算起，全长2300km.主要支流阿克苏河、叶尔羌河、和田河汇合于肖夹克，肖夹克以下到台特玛湖为塔河干流，总长1321km[1].从阿拉尔至英巴扎为上游，长度447km；英巴扎至卡拉为中游，长度398km；卡拉至台特玛湖为下游，长度428km.塔河在卡拉以上的上中游河段，大体上沿着塔克拉玛干大沙漠的北部，由西向东；在卡拉以下的下游河段，则是沿着塔克拉玛干大沙漠的东边，由北向南，成为塔克拉玛干与其东的库木塔格大沙漠的分界线。

　　塔河干流两岸并不产水，其水源完全来自发源于四周山区的各支流。据1981-1993年的水文资料，塔河在阿拉尔站的年均水量是42.99亿m³，其中阿克苏河提供33.58亿m³，叶尔羌河0.23亿m³，和田河9.18亿m³.另外从孔雀河引入2.12亿m³，因此塔河干流的总水量为45.11亿m³[2].

　　2、塔河下游应急输水

　　塔河两岸过去没有农田，盛长胡杨、枝柳和草甸，形成乔灌草的绿色植被带，是塔里木盆地四周人工绿洲的生态屏障，下游的南北向河道两岸更是分隔塔克拉玛干和库木塔格两大沙漠的绿色走廊。沿岸分布着8个县市十多个团场，养育了新疆总人口的44％。随着人类活动的加剧及水资源的过度利用，自20世纪50年代以来，虽然各支流的天然来水没有减少，但由于源流区和干流上中游的用水不断增加，塔河在卡拉下泄的水量，从12亿m³（1957～1967）逐步减少到2.7亿m³（1986～1995）。塔河最终从它的归宿台特玛湖萎缩到了现在的大西海子水库，塔河下游大西海子水库以下河道于1972年断流，台特玛湖于1974年干涸，下游320km河道由此消失。

　　伴随着塔河的萎缩，河水的减少，断流正改变着塔河流域特别是下游的生态环境。目前上游依靠拦截塔河的大闸，将塔河来水几乎全部引入水库，又从孔雀河修渠引水接济，才得勉强维持，即使这样，水源仍不能得到保证，而且受下游地区的沙化威胁，生态环境日益恶化，植被衰退、土地盐碱化、沙尘肆虐、少雨干旱、沙进人退。生态恶化的持续已经引起了国家和新疆维吾尔自治区政府的高度重视，从1999年开始实施治理塔河的一系列重大项目，解决水资源的问题是抢救塔河下游不断恶化生态的根本所在。

　　1998年以来，博斯腾湖源流开都河来水一直偏丰，使得博斯腾湖水位长期居高不下，在湖水位不断升高的同时，湖面面积也在不断扩大，最大时增至1700多km²，周围及孔雀河沿岸部分地区不同程度的存在着防洪安全隐患，一部分沿湖地带常被湖水淹没。而为抢救塔河下游绿色走廊和生态系统正愁找不到水源，20世纪80年代以来为解决兵团几个团场的农业灌溉问题，建成了连通孔雀河和塔河干流的库塔干渠，这使得向塔河下游应急输水在水源和工程设施方面都成为可能。从2000年5月中旬开始到2003年6月，新疆自治区连续5次将博斯腾湖水引入塔河下游河道内（5次输水的详细资料见表1）。

　　表1　塔河5次应急输水一览

　　3、生态变化遥感监测分析

　　利用遥感等现代信息技术进行流域下游生态动态变化监测非常必要，通过卫星遥感手段提供准确的流域生态本底数据不仅事半功倍、快速有效，而且在一些无资料地区，对于大范围的变化监测又是不可替代的。为了解输水效果及其对塔河下游的生态影响，中国水利水电科学研究院遥感技术应用中心针对连续3年的5次应急输水，通过地面实地调查并结合遥感信息资料的分析应用，研究了塔河应急输水，恢复下游河道以后生态改善状况的变化。

　　研究资料主要包括Landsat7-ETM全波段图像和中巴资源卫星图像。选取1999年、2001年和2003年7月3个时相的图像数据进行分析。其中1999年代表应急输水前的情况，2001年代表输水过程中应急输水对生态环境的影响，2003年代表输水一段时间后（3年后）生态环境的情况。由于2003年5月，Landsat卫星7号星发生故障，不能接收图像，所以2003年的图像选用中巴资源卫星代替。

　　合成后的三幅假彩色影像一方面反映出在应急输水工程实施后，塔河干流下游河道两侧植被有所改善。另一方面，从2003年的影像上可以明显看出在部分原先干涸的河道中已经有水，沿河两侧也出现了很多大小不一的湖泊，反映了河流生态的变化。除了通过宏观上定性的观察以外，通过遥感影像的解译结果和植被指数等来定量探讨应急输水对生态环境的影响。

　　3.1　影像解译和分类　根据卫星遥感TM影像，进行土地利用分类解译，分类包括：耕地、林地、草地、水域、城乡、工矿、居民用地和未利用土地。其中林地包括有林地、灌木林和疏林地，草地包括高覆盖度草地、中覆盖度草地和低覆盖度草地，水域包括河渠、湖泊和水库，未利用土地包括沙地、戈壁等。

　　考虑到输水影响范围有限，经过充分分析，提取从大西海子到台特玛湖之间的塔河干流，然后作缓冲区分析，将其两侧各10km的范围作为影像解译的区间，解译结果形成专题图。从分类图可以看出，2001年与1999年的遥感影像相比，林地和草地的面积整体上变化不大，但沿着河道两侧局部范围内的林地面积有所增大（即局部绿色部分增多），2003年与2001年的遥感影像相比，林地和草地的面积整体上变化也不大，沿着河道两侧局部范围内的林地面积增加更为明显。

　　根据图1进行统计得到3个时相土地分类面积（表2），从中可以看出面积的变化，在塔河干流实施应急工程输水后生态环境有所改善，如输水处下游河道两侧10km范围内灌木林面积1999年为25.0km²，2001年增加到34.9km²，2003年增加到44.5km²，疏林地面积1999年为128.2km²，2001年增加到130.7km²，2003年增加到181.7km².

表2　3个时相土地分类面积变化（单位：km²）

　　3.2　植被指数及植被覆盖度分析　植被指数（Vegetation Index）指从多光谱遥感数据中提取的有关地球表面植被状况的定量数值。通常是用红波段（R）和近红外（IR）波段通过数学运算进行线性或非线性组合得到的数值，用以表征地表植被的数量分配和质量情况。常用的植被指数有很多种，如 RVI、GVI、DDVI、PVI、EVI及NDVI等。据前人的研究，经综合比较，选用NDVI（归一化植被指数）。它能更好地适应本地形区植被盖度稀疏、盖度差异悬殊的区域景观特点，兼之其应用最广，所以运算结果有较好的外延接轨前景，也容易从其他应用效果中取得补充和印证，为长期的脆弱生态环境监测创造了方便条件。

　　求取NDVI前，进行图像数据预处理，包括图像增强、几何配准、特征提取、统计分析等基础性工作。需特别指出的是，由于两种图像的地面、光谱分辨率不同以及动态分析的需要，两个时相的图像各数据层像素间需严格配准，误差必须小于一个像元，经过计算后NDVI的值域在-1.0～1.0，将该值转换成0～255[3].

　　以1999和2001年的NDVI图像为例，图像上NDVI值沿塔里河下游河道两侧距离近处值较高，反映了植被覆盖较好，而离河道距离较远处，植被覆盖较差，为沙地。

　　将植被指数应用于资源环境的监测和评价，最好是赋予NDVI值以相应的植被覆盖度含义，由于植被指数既反映了特定景观中群落面积同景观总面积的比例关系，也反映了植物群落层片结构的特点，即反映了植被的盖度分布，同时也间接反映了植物的生物量高低。所以把植物指数转化为植被盖度等级，实际上是对植被指数进行综合和简化，对于植被生态景观面积变化的定量评价更为直观。

　　为达到既不丢失信息少的灰度级，又便于等级划分的目的，先将植被指数直方图按高斯分布统一做直方图规定化图像增强处理，使二者辐射亮度值分布尽可能一致，利用Erdas遥感图像处理软件所提供的假彩色编辑功能在屏幕上进行彩色编码，就完成了植被盖度数字图像的编制，供植被生态景观格局质量变化的评价[4].

　　表3的4个等级植被盖度是根据国家《土地利用现状调查技术规程》，全国《草场资源调查技术规程》，《全国沙漠类型划分原则》的有关条款为指导，并结合干旱区绿洲植被特有的生态特征来划分的。

表3　植被指数灰度级与盖度级的对应关系

表4　塔河下游河道两侧10km范围内植被覆盖分类面积（单位:km²）