免费文献传递   相关文献

Spatiotemporal variation of grassland vegetation NDVI in the middle and upper reaches of the Hei River and its response to climatic factors

黑河中上游草地NDVI时空变化规律及其对气候因子的响应分析



全 文 :书黑河中上游草地犖犇犞犐时空变化规律及其
对气候因子的响应分析
周伟1,王倩2,章超斌1,李建龙1
(1.南京大学生命科学学院,江苏 南京210093;2.甘肃民族师范学院历史文化系,甘肃 合作747000)
摘要:植被动态监测以及与气候变化的耦合关系是陆地生态系统研究的热点,而草地 NDVI时空动态以及影响其
变化的关键气候因子分析,对于草地覆盖变化的动态机理研究具有重要意义。基于遥感和GIS技术,利用黑河流
域1999-2007年SPOT-NDVI数据和气象站点数据,分析了黑河中上游草地植被NDVI的时空变化特征以及与
气候因子的相关性。结果表明,黑河中上游草地植被NDVI在1999-2007年间呈明显增加趋势;草地植被 NDVI
增加和减少的面积分别占草地总面积的71.53%和19.43%;增加的区域主要分布在托来山、走廊南山、冷龙岭和
河西走廊中段的典型草原和高寒草甸草地,减少的区域主要分布在托勒南山、走廊南山和冷龙岭的荒漠草地向非
草地过渡区域;黑河中上游不同草地植被NDVI的年内变化均呈单峰型,7月份达最大值;草地NDVI与气温呈极
显著正相关(狉=0.942,犘<0.01),与降水量呈显著正相关(狉=0.922,犘<0.05),NDVI与相对湿度呈不显著正相
关,与日照时数呈负相关性,并且不同草地类型牧草生长的限制因子不同。近年来伴随西北地区气候向暖湿转型,
黑河中上游草地植被覆盖明显增加,气温是影响研究区牧草生长的关键气候因子。
关键词:天然草地;SPOT-NDVI;时空动态;气候因子;黑河中上游
中图分类号:S812.1  文献标识码:A  文章编号:10045759(2013)01013810
  全球气候变化与陆地生态系统是当前全球变化研究的重要内容,气候变化对陆地生态系统的影响及其反馈
一直是其中的研究焦点之一[1]。植被是陆地生态系统的主体[2],也是连接大气、土壤和水分等自然要素的“纽
带”[3]。植被作为环境变化重要的指示标志,植被指数的变化能够揭示环境的演化、变迁,植被指数遥感数据具有
覆盖范围大,周期短的优点,可为大面积植被动态监测提供技术支撑。归一化植被指数(NDVI,normalized
differencevegetationindex)是表征地表植被覆盖和对光合有效辐射吸收能力的一个较好指标[4],被广泛用于生
物量、叶面积指数和植被生产力格局的估算和草地动态监测[5]。草地作为我国生态系统最大的天然屏障[6,7],具
有防风固沙、调节气候、保持水土和促进生态平衡的重要作用[8]。
NDVI在高植被覆盖地区存在过饱和现象,而对植被稀疏地区的植被变化尤其敏感。NDVI在草地植被研
究方面被广泛应用,Piao等[9]分析了草地植被NDVI动态变化以及对气候变化响应。王新欣等[10]利用 MODIS
NDVI对天山北坡中段草地进行生物量估测,建立了生物量动态估测模型。卫亚星等[11]利用 MODISNDVI数
据对青海省的草场进行分级,分析了其草地状况;毛飞等[12]利用美国国家海洋和大气局的高级甚高分辨率辐射
仪获得的NDVI(NOAA/AVHRRNDVI,nationaloceanicatmosphericadminister/advancedveryhighresolu
tionradiometerNDVI)数据对藏北那曲地区草地植被进行了分类。戴声佩等[13]利用SPOTNDVI分析了祁连
山地区草地植被的时空变化趋势。齐述华等[14]利用 NOAA/AVHRR数据分析了青藏高原中东部植被长势对
气候因子的响应,表明气候因子是影响植被覆盖的一个重要因子。辜智慧等[15]利用NOAA/AVHRR数据分析
了锡林郭勒草原植被覆盖变化以及对气候因子的响应。王军邦等[16]利用 MODISNDVI对内蒙古中部植被进行
动态监测,并且发现降水量引起的NDVI波动变化大于气温。
138-147
2013年2月
   草 业 学 报   
   ACTAPRATACULTURAESINICA   
第22卷 第1期
Vol.22,No.1
收稿日期:20120109;改回日期:20120329
基金项目:国家重点基础研究发展计划(973计划)项目(2010CB950702),国家自然科学基金项目(41271361)和国家863计划专题项目
(2007AA10Z231)资助。
作者简介:周伟(1985),女,山东泰安人,在读博士。Email:zhouw866@163.com
通讯作者。Email:jli2008@nju.edu.cn
黑河中上游位于祁连山北坡和河西走廊中段,处于气候变化的敏感区和生态环境脆弱区[17],是研究陆地生
态系统对气候变化响应机制的理想场所。草地是该流域的重要植被类型,近年来随着气候变化以及人类活动的
影响,黑河流域草地覆盖发生明显变化,因此黑河流域草地植被变化与气候因子的关系研究,将有助于理解流域
植被变化对气候变化的响应特征,对分析草地生态系统在气候变暖背景下的变化趋势具有重要的理论意义。目
前关于黑河流域草地植被覆盖变化的研究较少;并且在分析植被NDVI与气候因子间的相关性时,多采用气象
站点上的植被NDVI和气象因子的点状数据进行分析,分析结果不能综合反映面状植被覆盖变化对气候变化的
响应特征。鉴于此,本研究利用由SPOT4搭载的Vegetation传感器获取的SPOTVegetationNDVI数据分析
了黑河中上游草地植被NDVI时空变化特征,并采用每种草地类型的NDVI和气候因子的平均值进行相关性分
析,进一步探讨影响牧草生长的关键气候因子,以期为该区域应对气候变化以及草地资源的保护和管理提供理论
依据。
1 材料与方法
1.1 研究区概况
黑河流域是我国西北干旱区第二大内陆河流域(38°~42°N,98°~101°E),面积约13×104km2。研究区域
位于黑河上中游的祁连山地和河西走廊中段(图1)。黑河上中游地区自然条件复杂,水热条件差异大,具有典型
大陆性气候特征,由于受东南季风影响的强弱、水热条件的再分配以及复杂多样的地形变化等多因素的综合作
用,发育了多种多样的地貌、土壤类型和植被类型,草地类型呈现明显的垂直地带性分异[18,19]。黑河中上游的天
然草地分布于海拔1300~4100m的祁连山区和祁连山与龙首山、合黎山之间的山前倾斜平原,草地占研究区
总面积的73%,从低海拔到高海拔依次分布着沼泽草地、低湿地草甸、干荒漠草地、山地荒漠草地,草原化荒漠草
地、荒漠化草原、山地草原、山地草甸草原、山地草甸、高寒草原、高寒草甸[20]。
图1 研究区位置图
犉犻犵.1 犔狅犮犪狋犻狅狀狅犳狊狋狌犱狔犪狉犲犪
1.2 数据来源与处理
本研究所采用的数据为1999-2007年的SPOTVGTNDVI数据集,来源于国家自然科学基金委员会“中
国西部环境与生态科学数据中心”(http://westdc.westgis.ac.cn),该数据是通过对SPOT 源数据进行大气校
正、辐射校正、几何校正和拉伸后生成10d最大化合成 NDVI数据,其空间分辨率为1km。土地覆盖数据来源
于 GLC2000项目开发的基于SPOT4遥感数据的全球土地覆盖数据中国子集,在 ArcGIS9.3软件中利用公
931第22卷第1期 草业学报2013年
式:NDVI=DN×0.004-0.1将DN值转换为标准 NDVI值,然后利用黑河中上游矢量边界进行剪裁,得到黑河
中上游NDVI数据。
气象数据由国家气象信息中心资料室提供,包括7个气象站点1999-2007年逐月平均气温、降水量、相对湿
度和日照时数。在ArcGIS9.3中使用反距离权重(inversedistanceweighted,IDW)方法进行插值,得到研究区
气象栅格数据。NDVI和气象数据均投影为 Albers投影。
1.3 研究方法
1.3.1 NDVI数据处理 采用 ArcGIS空间分析技术和数理统计方法,对研究区草地植被NDVI数据进行统计
和分析。采用国际上惯用的最大值合成 MVC(maximumvaluecomposites)法对每月10d合成 NDVI进行处
理,即图像中每一个像元值用该月最大 NDVI值代替,该处理可以减少大气的云、颗粒、阴影、视角以及太阳高度
角的影响[21]。使用常用的累积平均法和均值法合成草地年累积 NDVI值,NDVI年均值,生长期的NDVI年均
值和月均值。
1.3.2 NDVI变化趋势线分析 趋势线分析法能模拟每个栅格的变化趋势[22],反映不同时期植被覆盖变化趋
势的空间特征,本研究采用趋势线分析模拟1999-2007年黑河中上游草地NDVI变化趋势,计算公式如下:
θslope=
狀×∑

犼=1
犼×犖犇犞犐犼-∑

犼=1
犼∑

犼=1
犖犇犞犐犼
狀×∑

犼=1
犼2- ∑

犼=1
( )犼 2
式中,狀为监测年数;犖犇犞犐犼为第犼年犖犇犞犐平均值;θslope为趋势线的斜率,其中θslope>0,说明犖犇犞犐在狀年间
的变化趋势是增加的,反之则是减少。根据θslope的变化范围[13,23],定义显著减少(θslope<-0.006),轻度减少
(-0.006<θslope<-0.001),基本不变(-0.001<θslope<0.001),轻度增加(0.001<θslope<0.006)和显著增加
(θslope>0.006)5个变化区间,并统计每个变化区间的面积及其面积百分比。
1.3.3 草地NDVI的空间变化量分析 本研究使用影像差异法分析草地植被NDVI空间变化量[24]。影像差异
是将2个不同时相的影像在完全配准情况下,从一个影像中将另一个影像的信息提取出来。它通常用于提取多
时相图像中随时间而变化的信息,常用于动态监测。研究中将2007年与1999年的2期NDVI影像相减,将结果
像元值大于0的记作增加,表示1999-2007年黑河流域草地植被NDVI增加,反之为减少。
1.3.4 草地NDVI与气候因子的相关性分析 根据GLC2000全球土地覆盖数据中的中国子集,并结合甘肃省
草地资源[18,25]调查资料,将研究区草地分为4类,荒漠草地(GLC2000中的类型11-荒漠草地)、平原草地
(GLC2000中的类型10-平原草地)、高寒草甸草地(GLC2000中的类型8-高山亚高山草甸和类型22-高山亚
高山草地)、典型草原(GLC2000中的类型12-草甸)。研究区气象站点分布不均匀并且每个站点附近的草地类
型不唯一,为了避免以气象站点的气候因子和NDVI的相关性系数代替面状草地类型的NDVI与气候因子相关
性,本研究根据研究区气象因子的插值数据和NDVI数据,提取不同草地类型多年月平均气温、降水量、相对湿
度、日照时数和月平均NDVI数据,分别计算NDVI与各气候因子间的相关性。5-9月黑河流域气候雨热同期,
处于植被生长季,因此本研究选择了草地生长期5-9月的NDVI数据分析植被NDVI与各气候因子间的相关
性。相关性系数公式如下:
狉=


犻=1
(狓犻-珚狓)(狔犻-珔狔)


犻=1
(狓犻-珚狓)2(狔犻-珔狔)槡 2
式中,狓犻是5-9月第犻月的植被NDVI,珚狓是5-9月NDVI的平均值;狔犻 是5-9月第犻月的气候因子值,珔狔是
5-9月气候因子的平均值。并对相关系数进行显著性检验。
1.3.5 数据统计分析 草地NDVI的时间变化特征、空间分布格局、草地NDVI变化趋势系数等数据统计分析
均在ArcGIS9.3中完成,通过SPSS16.0进行草地NDVI与气候因子的相关性分析和显著性检验,数据分析和
图表制作在Excel中完成。
041 ACTAPRATACULTURAESINICA(2013) Vol.22,No.1
2 结果与分析
2.1 黑河中上游草地植被NDVI空间分布特征分析
黑河中上游草地植被NDVI空间分布具有明显的地区差异性(图2和3),整体呈现东南部多,西北部少的格
局。位于黑河上游及其支流的托来山、走廊南山、冷龙岭和张掖市是典型草原覆盖区域,草地覆盖高,其中祁连
县、肃南县东南部和张掖市南部地区植被的年累计NDVI值较高,介于4~5。位于托勒南山北坡和走廊南山两
侧的高寒草甸草原覆盖区域以及酒泉和张掖市南部的典型草原覆盖区,NDVI值次之,介于3~4。从高寒草甸
草原向荒漠草原过渡的区域草地植被 NDVI介于2~3。祁连县和肃南县西北部以及张掖市北部的荒漠草地
NDVI较低,为1~2。位于高台县南部荒漠的草地NDVI值最低,小于1。
图2 黑河中上游草地类型分布
犉犻犵.2 犌狉犪狊狊犾犪狀犱狋狔狆犲狊犱犻狊狋狉犻犫狌狋犻狅狀犻狀狋犺犲犿犻犱犱犾犲犪狀犱狌狆狆犲狉狉犲犪犮犺犲狊狅犳犎犲犻犚犻狏犲狉
图3 1999-2007年黑河中上游草地植被年累计犖犇犞犐空间分布
犉犻犵.3 犜犺犲犱犻狊狋狉犻犫狌狋犻狅狀狅犳犮狌犿狌犾犪狋犻狏犲犖犇犞犐犳狉狅犿1999狋狅2007犻狀狋犺犲犿犻犱犱犾犲犪狀犱狌狆狆犲狉狉犲犪犮犺犲狊狅犳犎犲犻犚犻狏犲狉
141第22卷第1期 草业学报2013年
2.2 黑河中上游草地植被NDVI变化时间特征分析
2.2.1 NDVI年际变化特征 通过计算得到1999-2007年4种不同草地类型和全部草地植被年均NDVI和生
长季(5-9月)年均NDVI值(图4)。年均值NDVI和生长季的平均NDVI变化趋势基本一致,都呈现波动上升
趋势,这与程瑛等[26]的研究结论一致,由于祁连山中西段和黑河流域处于气候从暖干向暖湿转型的显著转型
区[27],降水量和融雪量增加,草地覆盖度增加。典型草原植被NDVI最高,平原草地次之,高寒草甸草地的ND
VI较低,荒漠草地最低,全部草地植被NDVI值介于平原草地和高寒草甸草地之间。年均NDVI最高值出现在
2007年,达到0.241,最低值出现在2001年,为0.202,生长季 NDVI最高值在2007年,达到0.381,最低值在
2001年,为0.301。
图4 黑河中上游不同草地植被犖犇犞犐年际变化
犉犻犵.4 犐狀狋犲狉犪狀狀狌犪犾犮犺犪狀犵犲狊狅犳犖犇犞犐犳狅狉犱犻犳犳犲狉犲狀狋犵狉犪狊狊犾犪狀犱狋狔狆犲狊犻狀狋犺犲犿犻犱犱犾犲犪狀犱狌狆狆犲狉狉犲犪犮犺犲狊狅犳犎犲犻犚犻狏犲狉
a:年均变化 Theannualchangeincluding12months;b:生长季变化(5-9月份)Representstheannualchangefrom MaytoSeptember.
2.2.2 NDVI年内变化特征 黑河中上游全部草地、
图5 黑河中上游不同草地植被犖犇犞犐年内变化特征
犉犻犵.5 犜犺犲犿狅狀狋犺犾狔犮犺犪狀犵犲狊狅犳犖犇犞犐犻狀犱犻犳犳犲狉犲狀狋犵狉犪狊狊犾犪狀犱狋狔狆犲狊
典型草原、平原草地、高寒草甸草地和荒漠草地的
NDVI年内月均变化曲线均呈单峰型(图5),年最大
NDVI值均出现在7月份,分别为0.430,0.569,
0.503,0.352,0.182。5-9月是牧草生长季,NDVI
变化明显,草地NDVI从5月份开始增加,6月份增加
迅速,7月份达到最大值,8月份 NDVI开始下降,10
月份牧草基本停止生长,每年的1-4月和11-12月
份草地NDVI变化不大。这主要受气候特征的影响,
黑河中上游地区5-9月雨热同期,并且在7-8月份
达到最佳,牧草生长茂盛,从11月到翌年2月牧草处
于枯黄期,植被覆盖较低,3-4月草地处于返青期,植
被生长缓慢,NDVI值也较低。其中荒漠草地 NDVI
的变化幅度明显低于其他草地类型。
2.3 黑河中上游草地植被NDVI变化的空间特征分析
2.3.1 黑河中上游草地植被NDVI变化空间分布特征 利用趋势分析方法模拟了研究区草地植被NDVI随时
间的变化趋势,并且趋势系数大于0.001时,植被 NDVI在1999-2007间处于增加趋势,当趋势系数小于
-0.001时植被覆盖度降低(表1,图6)。黑河中上游草地植被NDVI总体呈增加趋势,草地覆盖趋于改善,并且
241 ACTAPRATACULTURAESINICA(2013) Vol.22,No.1
存在明显的地区差异性。71.53%的天然草地植被NDVI增加,NDVI增加的面积为28418km2,NDVI增加的
区域主要分布在托来山、走廊南山、冷龙岭和张掖市、山丹县南部海拔在1800~4100m的山前倾斜平原和河西
走廊中段的天然草地地区。这些地区处于高山亚高山地区,受人类活动影响较小,且随着全球气候变暖,降水量
和冰雪融水增多[28,29],使得黑河流域草地植被覆盖度增大。草地 NDVI减少的区域主要分布在研究区托勒南
山、走廊南山和冷龙岭的荒漠草地向非草地过度区域,NDVI减少的面积为7717km2,占黑河中上游草地总面积
的19.43%。由于该地区受人类活动影响较大,过渡放牧等导致天然草地退化明显[30],草地覆盖降低。1999-
2007年只有9.04%的草地NDVI基本不变。
黑河中上游不同草地类型的NDVI变化存在明显差异,典型草原、平原草地、高寒草甸草地和荒漠草地植被
NDVI增加的面积分别为14778,515,8095和5300km2,分别占黑河中上游总草地面积的36.52%,1.30%,
20.37%,13.34%。典型草原、平原草地、高寒草甸草地和荒漠草地植被 NDVI分别减少3143,22,3664,887
km2,分别占研究区总草地面积的7.92%,0.06%,9.22%,2.23%。其中典型草原植被NDVI增加最明显,高寒
草甸草地增加面积次之,平原草地最低。
表1 1999-2007年黑河中上游草地植被犖犇犞犐变化统计
犜犪犫犾犲1 犜犺犲狊狋犪狋犻狊狋犻犮狊狅犳犵狉犪狊狊犾犪狀犱犖犇犞犐犮犺犪狀犵犲犻狀狋犺犲犿犻犱犱犾犲犪狀犱狌狆狆犲狉狉犲犪犮犺犲狊狅犳犎犲犻犚犻狏犲狉
趋势系数
Trend
index
变化等级
Change
level
全部草地
Totalgrassland
面积
Area
(km2)
比例
Ratio
(%)
典型草原
Typicalgrassland
面积
Area
(km2)
比例
Ratio
(%)
平原草地
Plaingrassland
面积
Area
(km2)
比例
Ratio
(%)
高寒草甸草地
Alpinemeadowsteppe
面积
Area
(km2)
比例
Ratio
(%)
荒漠草地
Desertgrassland
面积
Area
(km2)
比例
Ratio
(%)
<-0.006 ED 2085 5.25 643 1.62 3 0.008 1339 3.37 99 0.249
-0.006~-0.001MD 5632 14.18 2500 6.30 19 0.048 2325 5.85 788 1.984
-0.001~0.001 BI 3592 9.04 1631 4.18 27 0.068 1220 3.07 686 1.727
0.001~0.006 MI 16899 42.54 7927 19.95 167 0.420 4920 12.38 3885 9.779
>0.006 EI 11519 28.99 6851 16.57 348 0.876 3175 7.99 1415 3.562
 注:ED表示显著减少,MD轻度减少,BI基本不变,MI轻度增加,EI显著增加。
 Note:EDrepresentextremelydecrease,MDrepresentmilddecrease,BIrepresentbasicalyinvariability,MIrepresentmildincrease,EIrepresentex
tremelyincrease.
2.3.2 黑河中上游草地植被NDVI空间变化量 通过2007年和1999年2期年均NDVI值进行影像差异分析,
从中提取得到1999-2007年研究区域草地NDVI变化量分布(图7)。NDVI变化的分布特征与NDVI多年变
化趋势分布格局基本相同。NDVI增加量最大的区域分布于河西走廊的山丹军马场、张掖市东部和高台县西南
部的典型草原;典型草原、平原草地、高寒草甸草地、荒漠草地的最大增加量分别为0.187,0.273,0.186,0.211。
NDVI减少量最多的地区分布在托勒南山、走廊南山和托莱山和冷龙岭的高寒草甸草地和荒漠草原向非草地过
度区域,典型草原、平原草地、高寒草甸草地、荒漠草地的最大减少量分别为0.105,0.127,0.094,0.067。
2.4 黑河中上游草地植被NDVI与气候因子关系分析
影响牧草生长的气候因子很多,本研究选择了月平均气温、月平均降水量、相对湿度和日照时数。分析4种
草地类型1999-2007年牧草主要生长期(5-9月)月平均NDVI与上述4个气候因子的相关性(表2)。研究区
草地植被NDVI与月平均气温、月平均降水量、相对湿度的相关性系数分别为0.942(犘<0.01),0.922(犘<
0.05),0.762,其中与平均气温呈极显著正相关关系,并且NDVI与气温的相关系数大于降水量,这与郭铌等[31]
研究结论一致。由于黑河流域处于西北内陆干旱区,降水量较少,对海拔较高的区域气温较低将不利于植被的生
长;而近年来随着西北地区气候向暖湿转型[29],气温升高、降水量增加使植被生长加速,NDVI增大。草地NDVI
与日照时数呈负相关关系(狉=-0.291),因为研究区日照充足,当日照条件得到基本满足情况下,多余的日照反
而对牧草生长不利,这是由于日照越多,蒸散越大,最终导致NDVI下降[32]。
341第22卷第1期 草业学报2013年
图6 1999-2007年黑河中上游草地犖犇犞犐变化趋势空间分布
犉犻犵.6 犛狆犪狋犻犪犾犱犻狊狋狉犻犫狌狋犻狅狀狅犳犖犇犞犐犮犺犪狀犵犲狊狋狉犲狀犱犳狉狅犿1999狋狅2007
图7 1999-2007年黑河中上游草地植被犖犇犞犐变化量
犉犻犵.7 犛狆犪狋犻犪犾犱犻狊狋狉犻犫狌狋犻狅狀狅犳犖犇犞犐狏犪狉犻犪狋犻狅狀犳狉狅犿1999狋狅2007犻狀狋犺犲犿犻犱犱犾犲犪狀犱狌狆狆犲狉狉犲犪犮犺犲狊狅犳犎犲犻犚犻狏犲狉
典型草原NDVI与月平均气温、月平均降水量、相对湿度的相关性系数分别为0.907(犘<0.05),0.936(犘<
0.05),0.761;与日照时数相关系数为-0.308;牧草生长受降水量和气温的影响较大。
平原草地NDVI与月平均气温、月平均降水量、相对湿度的相关性系数分别为0.980,0.902,0.764,与平均
气温的相关性最高,通过0.01显著性水平检验;与日照时数相关性为-0.256。
高寒草甸草地NDVI与月平均气温呈显著性正相关关系(狉=0.934,犘<0.05),与降水量和相对湿度相关系
数分别为0.833,0.504,与日照时数相关性为-0.062。说明在海拔高降水较丰富的高寒草甸草地,气温升高有
利于植被生长,草地NDVI增大。相对湿度对草地NDVI影响不大。
荒漠草地NDVI与月平均气温、月平均降水量、相对湿度的相关性系数分别为0.862,0.931(犘<0.05),
441 ACTAPRATACULTURAESINICA(2013) Vol.22,No.1
0.82,与日照时数为-0.389。由于荒漠草地降水量较低,牧草生长容易受到干旱胁迫,因此水分是限制荒漠草地
植被生长的主要因子。
表2 黑河中上游草地植被犖犇犞犐与气候因子的相关性
犜犪犫犾犲2 犆狅狉狉犲犾犪狋犻狅狀犫犲狋狑犲犲狀犵狉犪狊狊犾犪狀犱狏犲犵犲狋犪狋犻狅狀犖犇犞犐犪狀犱犮犾犻犿犪狋犲犳犪犮狋狅狉狊犻狀狋犺犲犿犻犱犱犾犲犪狀犱狌狆狆犲狉狉犲犪犮犺犲狊狅犳犎犲犻犚犻狏犲狉
NDVI平均值
AverageNDVI
平均气温
Averagetemperature
平均降水量
Averageprecipitation
相对湿度
Relativehumidity
日照时数
Sunshineduration
全部草地Totalgrassland 0.942 0.922 0.762 -0.291
典型草原Typicalgrassland 0.907 0.936 0.761 -0.308
平原草地Plaingrassland 0.980 0.902 0.764 -0.256
高寒草甸草地Alpinemeadow 0.934 0.833 0.504 -0.062
荒漠草地Desertgrassland 0.862 0.931 0.820 -0.389
 注:表示极显著相关(犘<0.01);表示显著相关(犘<0.05)。
 Note:representextremelysignificantcorrelationat犘<0.01;representsignificantcorrelationat犘<0.05.
3 讨论与结论
植被被认为是综合反映生态环境变化的敏感指示器。有研究表明近20年来中国植被活动在增强[33]。黑河
流域处于生态环境脆弱区和气候变化的敏感区,是研究植被响应气候变化的理想区域。草地作为黑河中上游的
主要植被类型,在1999-2007近10年间,其草地植被覆盖变化显著,不同草地植被的年均NDVI和生长期ND
VI均呈增加趋势,2007年NDVI值最高,2001年最低。分析其动态变化的原因,可能存在两方面因素:一是近年
来随着祁连山自然保护区的成立,一系列生态保护政策的实施,如封山育林、退耕还林还草工程等,促使了该区域
植被覆盖的增加;另一方面,受西风环流的影响,该研究区降水量呈现出明显增加趋势[34]。本研究表明,不同草
地类型的NDVI年内变化呈单峰曲线型,5月份开始增大,7月份达到最高值,10月份牧草开始枯黄,从11月到
翌年4月份草地NDVI值较低且基本不变,NDVI年内变化特征与马明国等[23]和程瑛等[26]结论一致,这种变化
特征主要受该地区气候条件以及牧草自身生长特性等因素的影响。
黑河中上游草地植被NDVI空间分布呈现明显的地区差异性,整体呈现东南部的典型草原和高寒草甸草地
较高,西北部荒漠草地较低的格局。1999-2007年草地植被NDVI整体呈增加趋势,71.53%的天然草地植被覆
盖在增加,增加的区域主要分布在托来山、走廊南山、冷龙岭和张掖市、山丹县南部海拔1800~4100m的山前
倾斜平原和河西走廊中段的典型草原和高寒草甸草地;这些地区处于高山亚高山地区,受人类活动影响较小,且
随着全球气候变化,该区域气候由暖干向暖湿转型[27],降水量和冰雪融水增多[13,26]。典型草原、平原草地、高寒
草甸草地、荒漠草地的NDVI最大增加量分别为0.187,0.273,0.186,0.211。减少的区域主要分布在研究区托
勒南山、走廊南山和冷龙岭的荒漠草地向非草地过渡区域,由于该地区受人类活动影响较大,天然草地退化明显,
草地覆盖降低[30]。减少面积占黑河中上游草地总面积的19.43%;该区域典型草原、平原草地、高寒草甸草地、荒
漠草地的NDVI最大减少量分别为0.105,0.127,0.094,0.067。
气温、降水、日照和空气湿度与植物的生长密切相关。本研究中,黑河中上游牧草生长受气温和降水量的影
响较为明显,并且表现为NDVI与气温的相关系数大于降水,该结论与郭铌等[31]和毛飞[32]的研究一致,说明热
量充足,水资源充沛的条件有利于牧草生长;相对湿度对草地覆盖度的影响较小,日照时数与 NDVI呈负相关
性,说明日照时数越长,光照越多,植被蒸腾和土壤水分蒸发量就越大,从而导致空气相对湿度下降,结果不利于
土壤保墒和牧草生长[32]。对于不同草地类型,其牧草生长的限制因子不同。其中高寒草甸草地NDVI与气温的
相关系数大于降水,由于高寒草甸地区降水量相对丰富,对植被生长而言,热量条件较降水更为重要[31]。而对于
荒漠草地,由于其热量相对充足,水分则成为牧草生长的最主要限制因子。
541第22卷第1期 草业学报2013年
致谢:感谢国家自然科学基金委员会中国西部环境与生态科学数据中心提供SPOTVGT-NDVI数据集;同时
感谢中国气象局气象信息中心资料室提供的气象数据。
参考文献:
[1] 中国国家自然科学基金委员会生命科学部,中国科学院上海文献情报中心.全球变化与陆地生态系统-IGBP的核心计划之
一[A].中国国家自然科学基金委员会生命科学部,中国科学院上海文献情报中心编.全球变化与生态系统[M].上海:上海
科学技术出版社,1994:6295.
[2] 索玉霞,王正兴,刘闯,等.中亚地区1982年至2002年植被指数与气温和降水的相关性分析[J].资源科学,2009,31(8):
14221429.
[3] 孙红雨,王常耀,牛铮.中国植被覆盖变化及其与气候因子关系———基于NOAA时间序列数据[J].遥感学报,1998,2(3):
204210.
[4] AsrarG,FuchsM,KanemaseET,犲狋犪犾.Estimatingabsorbedphotosyntheticradiationandleafareaindexfromspectralre
flectanceinwheat[J].AgronomyJournal,1984,76:300306.
[5] PettoreliN.UsingthesatelitederivedNDVItoassessecologicalresponsestoenvironmentalchange[J].TrendsinEcology
&Evolution,2005,20:504510.
[6] 丁一汇,石永怡.气候变化高端访谈[M].北京:气象出版社,2004:126129.
[7] 任继周,梁天刚,林慧龙,等.草地对全球气候变化的响应及其碳汇潜势研究[J].草业学报,2011,20(2):122.
[8] 查勇,Jaygao,倪绍祥.国际草地资源遥感研究新进展[J].地理科学进展,2003,22(6):607617.
[9] PiaoSL,MohammatA,FangJY,犲狋犪犾.NDVIbasedincreaseingrowthoftemperategrasslandsanditsresponsestoclimate
changesinChina[J].GlobalEnvironmentChangehumanandPolicyDimensions,2006,16:340348.
[10] 王新欣,朱进忠,范艳敏,等.基于 MOODISNDVI的天山北坡中段草地动态估产模型研究[J].草业科学,2009,26(7):
2427.
[11] 卫亚星,王莉雯,刘闯.基于 MODISNDVI时序数据的青海省草地分级[J].资源科学,2008,30(5):688693.
[12] 毛飞,侯英雨,唐世浩,等.基于近20年遥感数据的藏北草地分类及其动态变化[J].应用生态学报,2007,18(8):1745
1750.
[13] 戴声佩,张勃,王海军,等.基于SPOTNDVI的祁连山草地植被覆盖时空变化趋势分析[J].地理科学进展,2010,29(9):
10751080.
[14] 齐述华,王长耀,牛铮,等.利用NDVI时间序列数据分析植被长势对气候因子的响应[J].地理科学进展,2004,23(3):
9199.
[15] 辜智慧,陈晋,史培军,等.锡林郭勒草原1983~1999年 NDVI逐旬变化量与气象因子的相关分析[J].植物生态学报,
2005,29(5):753765.
[16] 王军邦,陶健,李贵才,等.内蒙古中部 MODIS植被动态监测分析[J].地球信息科学学报,2010,12(6):835842.
[17] 刘晶,刘学录,王哲锋.祁连山东段景观格局变化及其驱动因子研究[J].草业学报,2011,20(6):2633.
[18] 甘肃草原总站.甘肃草地资源[M].兰州:甘肃科学技术出版社,1999:56120.
[19] 黄德青,于兰,张耀生,等.祁连山北坡天然草地地上生物量及其与土壤水分关系的比较研究[J].草业学报,2011,20(3):
2027.
[20] 柴来智,张和平,李瑜年.甘肃省张掖地区草地资源[M].中国草业科学,1988,(专辑):4987.
[21] 顾娟,李新,黄春林.NDVI时间序列数据集重建方法述评[J].遥感技术与应用,2006,21(4):391395.
[22] StowD,DaeschnerS,HopeA,犲狋犪犾.Variabilityoftheseasonalyintegratednormalizeddifferencevegetationindexacross
thenorthslopeofAlaskainthe1990s[J].InternationalJournalofRemoteSensing,2003,24:11111117.
[23] 马明国,董立新,王雪梅.过去21a中国西北植被覆盖动态监测与模拟[J].冰川冻土,2003,25(2):232236.
[24] 章明,张培松,刘洪斌,等.基于SPOTVEGETATION数据的海南岛年际植被变化研究[J].西南大学学报(自然科学
版),2009,31(3):148153.
[25] 赵忠,何毅,李青,等.肃南裕固族自治县草原资源调查[J].草业学报,2010,19(6):231247.
[26] 程瑛,徐殿祥,郭铌.近20年来祁连山区植被变化特征分析[J].干旱区研究,2008,25(6):772777.
641 ACTAPRATACULTURAESINICA(2013) Vol.22,No.1
[27] 施雅风,沈永平,李栋梁,等.中国西北气候由暖干向暖湿转型的特征和趋势探讨[J].第四纪研究,2003,23(2):152
164.
[28] 贾文雄,何元庆,李宗省,等.祁连山区气候变化的区域差异特征及突变分析[J].地理学报,2008,63(3):257269.
[29] 张杰,李栋梁.祁连山及黑河流域降雨量的分布特征分析[J].高原气象,2004,23(1):8188.
[30] 郭铌,杨兰芳,李民轩.利用气象卫星资料研究祁连山区植被和积雪变化[J].应用气象学报,2003,14(6):700707.
[31] 郭铌,朱燕君,王介民,等.近22年来西北不同类型植被 NDVI变化与气候因子的关系[J].植物生态学报,2008,32(2):
319327.
[32] 毛飞,卢志光,张佳华,等.近20年藏北地区AVHRRNDVI与气候因子的关系[J].生态学报,2007,27(8):31983205.
[33] 方精云,朴世龙,贺金生,等.近20年来中国植被活动在增强[J].中国科学C辑,2003,33(6):554565.
[34] 王胜,王亚明,王润元,等.近10a来祁连山地区气候及水资源研究现状[J].干旱气象,2007,23(3):8287.
犛狆犪狋犻狅狋犲犿狆狅狉犪犾狏犪狉犻犪狋犻狅狀狅犳犵狉犪狊狊犾犪狀犱狏犲犵犲狋犪狋犻狅狀犖犇犞犐犻狀狋犺犲犿犻犱犱犾犲犪狀犱狌狆狆犲狉狉犲犪犮犺犲狊狅犳
狋犺犲犎犲犻犚犻狏犲狉犪狀犱犻狋狊狉犲狊狆狅狀狊犲狋狅犮犾犻犿犪狋犻犮犳犪犮狋狅狉狊
ZHOU Wei1,WANGQian2,ZHANGChaobin1,LIJianlong1
(1.SchoolofLifeScience,NanjingUniversity,Nanjing210093,China;2.DepartmentofHistory
andCulture,GansuNormalUniversityforNationalities,Hezuo747000,China)
犃犫狊狋狉犪犮狋:Vegetationdynamicsanditscoupledrelationwiththeclimateisahotspotinthestudyofhowterres
trialecosystemsrespondtoclimatesystems.AnalysisofgrasslandNDVI(normalizeddifferencevegetationin
dex)spatiotemporaldynamicsandexplorationofkeyclimatefactorshasimportantsignificanceforresearchof
grasslandcoverchange.BasedontheRS,GISandSPOTNDVIdataandmonthlymeteorologicaldatafrom
1999to2007,thespatiotemporalvariationcharacteristicsofgrasslandvegetationNDVIinthemiddleandupper
reachesoftheHeiRiverandthecorrelationbetweenNDVIandfourclimatefactorswereanalyzed.1)Thean
nualaverageNDVIandgrowingseasonNDVIincreasedinthemiddleandupperreachesoftheHeiRiver.The
areaofgrasslandNDVIthatincreasedaccountedfor71.53%ofthetotalgrasslandarea,whilethatwhichde
creasedaccountedfor19.43%.ThegrasslandvegetationNDVIintypicalsteppe,alpinesteppeandmeadowof
theTuolaiMountain,SouthernMountainofHexicorridor,LenglongRidgeandthemiddleoftheHexicorridor
increased,whilethedecreasedregionswerelocatedinthezonebetweendesertgrasslandandnongrasslandof
theTuolaiMountain,SouthernMountainoftheHexicorridorandLenglongRidge.2)Thegrasslandvegeta
tionNDVIoftheHeiRivershowsasinglepeakcurvewithintheyearlychangesofdifferentgrasslandtypes,
andthevalueofNDVIreachedapeakinJuly.3)Therewasanextremelysignificantpositivecorrelationbe
tweengrasslandNDVIandaveragetemperature(犚=0.942,犘<0.01),andasignificantcorrelationwithaver
ageprecipitation(犚=0.922,犘<0.05).ThecoefficientofcorrelationbetweengrasslandNDVIandrelativehu
miditywaspositivebutnotsignificant,whilethecorrelationwithsunshinehourswasnegative.Thelimiting
factorwasdifferentinthefourgrasslandtypes.Inrecentyears,withtheshiftfromwarmdrytowarmhumid
climateinNorthwestChina,grasslandvegetationcoverageinthemiddleandupperreachesoftheHeiRiver
significantlyincreasedandtemperaturewastheimportantfactorforgrassgrowth.
犓犲狔狑狅狉犱狊:naturalgrassland;SPOT-NDVI;spatiotemporaldynamics;climatefactors;themiddleandupper
reachesofHeiRiver
741第22卷第1期 草业学报2013年