祁连山是我国青藏高原东北部一个巨大的边缘山系,位于甘肃、青海两省交界地,由一系列近似西北—东南走向的平行山岭和山间盆地组成,海拔2500~5000m,降水300~500 mm,气候寒冷。许多海拔4000m以上的中、高山终年积雪,分布着现代冰川,是河西走廊的天然水库。3000m以下地区自然植被较好,有成片的天然林和灌丛,河谷、盆地水草丰盛是良好的天然牧场。
祁连山冰雪资源和森林资源是河西走廓经济建设的基础,是河西人民赖以生存的命脉。车克钧等[1]观测实验研究表明,祁连山水源涵养林不仅具有拦蓄降水、缩小温差、保持较高土壤湿度等作用,而且具有十分显著的消洪补枯、涵养水源、保持水土的水文生态作用。祁连山区位于我国气候敏感区和生态环境的脆弱带,由于人为和气候因素,自20世纪50年代以来祁连山森林覆盖率下降近50 %,草场退化,自然生态环境日趋恶化,为此国家于1986年成立祁连山国家自然保护区。
本文利用1989年和1998年NOAA卫星AVHRR资料提取出植被、积雪等信息,对祁连山自然保护区植被分布状况及其变化做了研究,同时研究了党河、疏勒河、北大河、黑河、西营河5个流域积雪面积年和年代变化特征,旨在为祁连山生态环境保护提供科学依据。
1 基本原理和资料 1.1 基本原理美国NOAA系列气象卫星携带改进的甚高分辨率辐射仪 (AVHRR) 有5个探测通道,其中通道1为可见光波段 (0.58~0.68 μm);通道2为近红外波段 (0.725~1.1 μm);通道3为中红外波段 (3.55~3.95 μm);通道4、5为热红外波段 (分别为10.3~11.3 μm,11.5~12.5 μm)[2],其星下点分辨率约为1.1 km2。
不同性质的地表有其自身特有的光谱分布特性。绿色植物的主要光谱反射特点是:绿叶中的叶绿素在0.5~0.7 μm的可见光谱段有2个强吸收谷,其反射率较小 (一段小于20 %);但在0.7~1.3 μm的近红外波段,由于叶肉海绵组织结构中有许多空腔,因此具有很大的反射表面,而且细胞内叶绿体呈水溶胶状态,其反射率较大,可达60 %。NOAA系列气象卫星AVHRR探测器通道1波长位于植物的可见光吸收区,通道2位于近红外强反射区,将这两个通道探测值组合起来便可以得到各种植被指数,很适合于监测植被[3]。
雪的光谱特点是在可见光和近红外波段 (即NOAA/AVHRR通道1、2) 具有高于陆地表面的反射率 (这点与云非常相近)[4],但在短波红外波段 (1.55~1.75 μm) 和 (2.10~2.35 μm) 雪吸收太阳辐射,反射率很低,而云在整个反射红外波段都具有较高的、均匀的光谱响应,因此在中短波红外波段 (1.55 μm或2.10 μm附近) 区分云和雪是可能的。雪的这一特点在NOAA/AVHRR通道3中也有反映[5~7],可作为区分雪和云的根据。
1.2 资料及处理本文资料由NOAA气象卫星资料、气象资料、祁连山国家自然保护区边界、祁连山流域边界资料组成。降水是影响西北地区植被生长的主要因素,为了减少气象条件对植被的影响,本文首先利用甘肃省55个气象站1986~1998年历年降水资料和农业气象资料进行年景分析,得出1989和1998年为相似年;再选祁连山区海拔2000m以上的10个站点,进行这两年降水偏差分析。从祁连山1989年和1998年1~8月降水量分析 (见表 1,注:卫星资料用7、8两月的资料,故考虑1~8月降水) 可以看到,除祁连山中部个别站点1998年降水较1989年略多外,祁连山区大部分地区1989年降水与1998年接近,降水因素对植被的影响可以忽略。
![]() |
表 1 祁连山区各站1989与1998年1~8月降水量偏差分析 |
对1989年NOAA-11和1998年NOAA-14资料做太阳高度角订正和临边变暗订正,对各时次各通道资料以青海湖为地物标志做定位处理。
将甘肃省林业勘察设计院规划的甘肃省祁连山国家级自然保护区总体设计图中祁连山保护区边界数字化后得到祁连山自然保护区边界;根据1:100万地形高度图读取祁连山党河、疏勒河、北大河、黑河、西营河流域边界。
2 计算方法(1) 植被本项目选归一化植被指数作为监测植被变化的指标,其计算公式为:
![]() |
(1) |
式中INDV为归一化植被指数,C1、C2分别为通道1、2的反射率值。INDV的特点是可以消除由于照明条件变化、倾斜、观测姿态所造成的各通道反射率的改变,基本上不依赖于卫星的观测角,适合于对植被进行动态监测[8]。由于本文所用资料涉及到NOAA-11和NOAA-14两颗卫星,考虑到卫星更换和AVHRR探测器通道1、2的衰减问题,我们在巴丹吉林沙漠选两个区域 (40.42°~40.47°N,101.60°~101.67°E和40.20°~40.27°N,101.62°~101.70°E),在黑河下游戈壁选一个区域 (41.15°~41.09°N,100.64°~100.70°E) 分析1989和1998年INDV变化情况 (表 2),可以看到1989和1998年7、8月这3个区域INDV没有明显差别,说明探测器的衰减对本文所用的归一化植被指数影响不大。
![]() |
表 2 沙漠戈壁地区1989年与1998年INDV对比 |
(2) 积雪根据通道1反射率高低和通道3和4亮温差的大小提取积雪信息。
![]() |
(2) |
![]() |
(3) |
式中C1为通道1反射率值,D34为通道3与通道4亮温差,T3、T4为通道3、4的亮温,a、b分别为通过大量实践得到区分云、雪的阈值。由于云移动性较大,采用几天资料相叠加合成即可排除大部分云干扰,详细计算方法参见文献[7]。
积雪面积的计算公式为
![]() |
(4) |
式中σ为某一区域中的积雪面积,σi为该区域中某一积雪像元面积 (考虑地球曲率影响),n为该区域中积雪像元总数,为了便于分析将各流域积雪面积除以流域总面积,得到积雪面积比率。
将卫星资料根据需要按上述公式计算出各时次植被指数、各流域积雪面积和积雪面积指数。
3 祁连山自然保护区植被变化特征 3.1 植被指数变化特征夏季是地表植被最丰富的阶段,也是监测植被的最佳时段。利用1989年和1998年7月NOAA/AVHRR资料分别计算出这两年的植被指数,根据我们实际工作经验并参考美国NESDIS方法,将植被指数按数值大小分成裸地、稀疏植被、较少植被、适中植被、茂密植被和很密植被6个等级,分别统计计算出整个祁连山自然保护区、保护区西段、保护区东段、大黄山等地区不同植被指数等级的面积及其变化。
图 1是祁连山自然保护区1989年与1998年夏季植被指数遥感监测图,从图中可以看到祁连山保护区主体部分为西北—东南向的狭长地带,山丹附近的不连续带将其分为东、西两段,大黄山、东大山、十八里铺为祁连山北部的几座独立的小山。山区植被指数基本自西向东降低,与降水量自东向西减少相对应,祁连山东段植被覆盖率和植被指数等级高于西段,保护区西端为植被指数最低区。比较这两年植被分布可以看到祁连山主体部分1998年植被指数大都较1989年低,特别是保护区西端祁丰管理站十年来植被指数下降最明显,裸地面积明显增加,而大黄山、东大山管理站1998年植被指数却较1989年有所增加。表 3是祁连山自然保护区十年来各种等级植被面积变化统计表,表中整个祁连山保护区十年来裸地面积增加20.7 %,很密植被减少了33.5 %,其中保护区西段裸地和稀疏植被面积增加30 %和12.9 %,较少植被、适中植被、茂密植被和很密植被分别下降11.5 %、12.6 %、13.1 %和23.6 %,东段裸地面积和很密植被分别下降6.6 %和38 %,稀疏植被、较少植被、适中植被、茂密植被分别增加52.5 %、6.6 %、20.7 %和14.8 %,表明祁连山区十年来植被指数在下降,植被退化严重。笔者曾到祁连山考察,山区深处仍有很多羊群,牧民骑摩托车放牧,可见祁连山区植被退化与过渡放牧有很大关系。大黄山、东大山十年来裸地面积分别减少14.1 %和82.8 %,很密植被增加85.2 %和100 %,植被生长良好。从甘肃省人大资源环境委员会了解到这些地方植被指数增加的主要原因是:近年来大黄山、东大山实行生态移民工程,将山中居民全部迁到山外,封山育林。
![]() |
|
图 1. 祁连山自然保护区1989年与1998年夏季植被指数遥感监测图 |
![]() |
表 3 1989、1998年祁连山自然保护区植被等级面积变化统计 |
3.2 针叶林、草地变化
通过实地考察了解到祁连山区水源涵养林由青海云杉、祁连山圆柏等常绿针叶林和灌木林组成,根据针叶林、灌木林和草的季节变化特点,冬季灌木林和草均已枯黄,植被指数较低,而针叶林仍然保持绿色,植被指数较高,因此用冬季的植被指数即可判识出针叶林。利用1989年和1998年1~2月资料提取出针叶林并分析其变化,取针叶林面积较大的肃南林场部分地区 (范围38.75°~38.81°N,99.6°~99.79°E),计算该区域1989年针叶林面积为103 km2,1998年为99.3 km2,减少3.6 %,与韩涛[9]用同期的TM资料计算结果对比,误差为13 %,表明10年来祁连山针叶林面积变化不大,植被指数减少的地区以灌木林和草为主。根据车克均等[1]10年间对不同森林类型林冠截留的观测表明,试验区雨季 (5~9月) 林冠层对降水的截留率:藓类云杉林为25.86 %;祁连圆柏林为28.95 %;灌木林为66.1 %。灌木林由于冠层枝叶表面积大且镶嵌分布紧密,截留降水量最大,灌木林的大幅度减少必将对祁连山生态环境产生负反馈作用。据我们对祁连山实地的调查,过渡放牧是祁连山植被退化的主要原因,建议各级政府加大对祁连山区的保护,退牧还林还草。
4 积雪变化特征 4.1 各流域积雪年变化特征将1997、1998年逐日卫星资料以10d为单位按文献[7]方法合成,得到逐旬积雪资料,分析积雪年变化特征 (由于本文计算山区积雪面积时没有考虑坡度坡向的影响,结果不免会受到一些影响)。图 2是党河、疏勒河、黑河、西营河流域1997、1998年旬积雪面积变化图,可以看到整个祁连山区积雪变化基本特点是年变化呈双峰形,主峰出现在春季 (3~6月),次峰出现在初冬 (11月),与陈乾等研究结果[10]接近。7~8月是积雪最少期,尽管祁连山降水集中在这个时段,但降水多以降雨为主,即使降雪也很快融化。分析各流域1997年和1998年旬积雪变化发现,夏季 (7~9月) 各流域积雪都很少,而多雪期从西到东出现的时段不尽相同,有各自的特点。西段的党河流域 (图 2a) 多雪期分别出现在11月中旬至次年1月和3~6月,2月为积雪较少期。中段的疏勒河、北大河、黑河流域 (图 2b、2c) 多雪期出现在3~6月和10~11月,与党河流域不同的是12月至次年2月为少雪期。西营河流域 (图 2d)1~6月和10~11月为多雪期。本文发现祁连山积雪要靠新雪不断补充才能维持,即使在冬季,如果没有降雪补充,山区积雪也会在一、两旬内迅速减少,因此有效开发祁连山区云水资源,实行祁连山区人工增雨雪十分必要。
![]() |
|
图 2. 祁连山1997、1998年旬积雪面积变化 (a) 党河流域,(b) 北大河流域,(c) 黑河流域,(d) 西营河流域 |
表 4给出了1998年各流域旬积雪相关系数,可以看到祁连山西段的党河流域积雪与相邻疏勒河流域相关较好,与北大河流域有一定的相关性但与其他流域相关不大;疏勒河与北大河有很好的相关,与黑河有一定的相关;黑河与北大河和西营河均有一定的相关性;表明祁连山东、中、西段积雪变化有各自的特点,存在一定差异性。
![]() |
表 4 1998年各流域旬积雪相关系数 |
4.2 1989年和1998年积雪变化特征
根据1989年和1998年积雪变化分析 (图略),1989与1998年春季积雪相差不大,但1998年秋、冬季积雪明显少于1989年;同时由于1998年气温明显偏高,1998年积雪面积变化幅度大于1989年。从1989年与1998年旬平均积雪比率及其变化统计表 (表 5) 可以看到,1998年各流域旬平均积雪面积均较1989年减小,自西向东减小幅度逐渐增大,党河流域减小2.17 %,西营河流域减小10.05 %。
![]() |
表 5 1989年与1998年旬平均积雪比率及其变化统计表 |
5 结语
(1) 祁连山自然保护区主体部分灌木林和草的退化是造成1998年植被指数较1989年有较大幅度下降的主要原因,通过生态移民可以使植被得以恢复,大黄山植被指数增加就是很好的例子,建议各级政府加大祁连山生态环境的保护力度。
(2) 祁连山区积雪变化基本特点是年变化呈双峰形,主峰出现在春季 (3~6月),次峰出现在初冬 (11月)。
(3) 各流域积雪年变化特征是:各流域夏季均为少雪期,但多雪期出现的时段不同。党河流域出现在11月中旬至次年1月和3~6月;疏勒河、北大河、黑河流域出现在3~6月和10~11月;西营河流域出现在1~6月和10~11月。
(4) 祁连山积雪要靠新雪不断补充才能维持,即使在冬季如果没有降雪补充,山区积雪也会在一、两旬内迅速减少,因此有效开发祁连山区云水资源,实行祁连山区人工增雨雪十分必要。
(5) 1989年和1998年积雪变化特征是:1998年秋冬季积雪明显少于1989年,同时由于1998年气温明显偏高,积雪面积变化振幅加大,1998年各流域旬平均积雪面积均较1989年减小,自西向东减小幅度逐渐增大,党河流域减小2.17 %,西营河流域减小10.05 %。
[1] | 车克钧, 傅辉恩, 王金叶. 祁连山水源林生态系统结构与功能的研究. 林业科学, 1998, 34, (5). |
[2] | Rao P K等编. 气象卫星--系统、资料及其在环境中的应用. 北京: 气象出版社, 1994. 140. |
[3] | 林培. 农业遥感. 北京: 农业出版社, 1990: 87. |
[4] | 斯韦恩PH, 戴韦SM. 遥感定量方法. 北京: 科学出版社, 1984: 160-169. |
[5] | 魏丽, 钟强. 中国西部云和积雪的多通道特征和识别. 高原气象, 1991, 10, (4): 371–380. |
[6] | 刘玉洁, 王丽波, 刘诚, 等. 卫星遥感雪灾监测与分析. 北京: 气象出版社, 1998. |
[7] | 郭铌, 倾继祖. 气象卫星资料对积雪的遥感监测与分析. 遥感技术与应用, 2000, 11, (4): 237–240. |
[8] | 吴艾笙, 钟强. 黑河实验区AVHRR反射率资料的各向异性订正. 高原气象, 1990, 9, (2): 136–144. |
[9] | 韩涛. 用TM卫星资料对祁连山部分地区进行针叶林、灌木林分类研究. 遥感技术与应用, 2002, 17, (6): 317–321. |
[10] | 陈乾, 陈添宇. 祁连山区季节性积雪资源的气候分析. 地理研究, 1991, 10, (4): 23–37. |