免费文献传递   相关文献

Analysis of vegetation dynamics using GIMMS NDVI3g in the Three-Rivers Headwater Region from1982 to 2012

三江源区近30年植被生长动态变化特征分析



全 文 :书犇犗犐:10.11686/犮狔狓犫2015111 犺狋狋狆://犮狔狓犫.犾狕狌.犲犱狌.犮狀
杜加强,贾尔恒·阿哈提,赵晨曦,房世峰,刘伟玲,阴俊齐,袁新杰,徐延达,舒俭民,何萍.三江源区近30年植被生长动态变化特征分析.草业
学报,2016,25(1):112.
DUJiaQiang,JIAERHENGAhati,ZHAOChenXi,FANGShiFeng,LIUWeiLing,YINJunQi,YUANXinJie,XUYanDa,SHUJianMin,
HEPing.AnalysisofvegetationdynamicsusingGIMMSNDVI3gintheThree-RiversHeadwaterRegionfrom1982to2012.ActaPrataculturae
Sinica,2016,25(1):112.
三江源区近30年植被生长动态变化特征分析
杜加强1,2,贾尔恒·阿哈提3,赵晨曦3,房世峰4,刘伟玲1,2,阴俊齐3,
袁新杰3,徐延达1,2,舒俭民1,2,何萍1,2
(1.中国环境科学研究院,北京100012;2.中国环境科学研究院环境基准与风险评估国家重点实验室,北京100012;3.新疆环境保护
科学研究院,新疆 乌鲁木齐830011;4.中国科学院地理科学与资源研究所资源与环境信息系统国家重点实验室,北京100101)
摘要:三江源区为长江、黄河、澜沧江提供了丰富的水资源,其地表植被状况与我国水安全和生态安全密切相关。
由于所使用的数据时间序列长度有限,以往研究缺乏近30年整体视角的连续时间序列研究,也很少分析多时段变
化趋势的持续性。为此,本文利用最新发布的1982-2012年的GIMMSNDVI3g数据集,采用多个嵌套时间序列
的分析方法,在多个空间尺度上,研究了近30年来三江源区各季节植被活动的变化趋势和过程。研究表明,1)三
江源区生长季和各季节NDVI在近30年均呈增强趋势,春季植被变化具有明显的阶段性,1998年前后NDVI分别
呈显著增加、减少趋势。2)生长季、春季和夏季NDVI的增加随时段延长而趋缓,尤其是春季减缓迅速甚至停滞,
而秋季植被活动持续增强。3)尽管区域尺度生长季、春季和夏季 NDVI变化趋缓,但在像元尺度显著变化的区域
范围却多呈增加趋势。4)不同植被类型显著变化的季节、时段存在差异。针叶林和灌丛NDVI主要在春季前几个
时段增加显著,荒漠、草甸的增加主要集中在春季、夏季前几个时段。5)各季节 NDVI在14个时段增加趋势持续
性的差异,引起各季节对生长季植被活动变化贡献的变化。总的趋势是,前几个时段,春季 NDVI变化较为剧烈,
夏季NDVI的贡献最大,而随时段延长,逐渐过渡到秋季NDVI变化幅度、变化量最大,夏季次之。
关键词:GIMMSNDVI3g;区域尺度;像元尺度;植被类型;季节;植被生长动态;三江源区  
犃狀犪犾狔狊犻狊狅犳狏犲犵犲狋犪狋犻狅狀犱狔狀犪犿犻犮狊狌狊犻狀犵犌犐犕犕犛犖犇犞犐3犵犻狀狋犺犲犜犺狉犲犲-犚犻狏犲狉狊犎犲犪犱
狑犪狋犲狉犚犲犵犻狅狀犳狉狅犿1982狋狅2012
DUJiaQiang1,2,JIAERHENGAhati3,ZHAOChenXi3,FANGShiFeng4,LIU WeiLing1,2,YINJunQi3,
YUANXinJie3,XUYanDa1,2,SHUJianMin1,2,HEPing1,2
1.犆犺犻狀犲狊犲犚犲狊犲犪狉犮犺犃犮犪犱犲犿狔狅犳犈狀狏犻狉狅狀犿犲狀狋犪犾犛犮犻犲狀犮犲狊,犅犲犻犼犻狀犵100012,犆犺犻狀犪;2.犛狋犪狋犲犓犲狔犔犪犫狅狉犪狋狅狉狔狅犳犈狀狏犻狉狅狀犿犲狀狋犪犾犆狉犻
狋犲狉犻犪犪狀犱犚犻狊犽犃狊狊犲狊狊犿犲狀狋,犆犺犻狀犲狊犲犚犲狊犲犪狉犮犺犃犮犪犱犲犿狔狅犳犈狀狏犻狉狅狀犿犲狀狋犪犾犛犮犻犲狀犮犲狊,犅犲犻犼犻狀犵100012,犆犺犻狀犪;3.犡犻狀犼犻犪狀犵犃犮犪犱犲犿狔
狅犳犈狀狏犻狉狅狀犿犲狀狋犪犾犘狉狅狋犲犮狋犻狅狀犛犮犻犲狀犮犲,犝狉狌犿狇犻830011,犆犺犻狀犪;4.犜犺犲犛狋犪狋犲犓犲狔犔犪犫狅狉犪狋狅狉狔狅犳犚犲狊狅狌狉犮犲狊犪狀犱犈狀狏犻狉狅狀犿犲狀狋犪犾犐狀
犳狅狉犿犪狋犻狅狀犛狔狊狋犲犿,犐狀狊狋犻狋狌狋犲狅犳犌犲狅犵狉犪狆犺犻犮犛犮犻犲狀犮犲狊犪狀犱 犖犪狋狌狉犪犾犚犲狊狅狌狉犮犲狊犚犲狊犲犪狉犮犺,犆犺犻狀犲狊犲犃犮犪犱犲犿狔狅犳犛犮犻犲狀犮犲狊,犅犲犻犼犻狀犵
100101,犆犺犻狀犪
犃犫狊狋狉犪犮狋:TheThree-RiversHeadwaterRegion(TRHR)referstothesourceareaoftheYangtzeRiver,
YelowRiverandLancangRiver,animportantwaterresource.Vegetationactivitiesinthisregionareclosely
第25卷 第1期
Vol.25,No.1
草 业 学 报
ACTAPRATACULTURAESINICA   
1-12
2016年1月
收稿日期:20150303;改回日期:20150430
基金项目:国家自然科学基金资助项目(41001055),国家环保公益性行业科研专项经费项目(2012090275),中央级公益性科研院所基本科研业
务专项(2012YSKY13)资助。
作者简介:杜加强(1980),男,山东临沂人,副研究员,博士。Email:dujq@craes.org.cn
通信作者Correspondingauthor.Email:heping@craes.org.cn
relatedtowaterandecologicalsecurityinChina.Asacriticalcomponentofterrestrialecosystems,vegetation
hasbecomeoneoftheimportantissuesinglobalclimatechange.Numerouscasestudieshavebeenconducted
onvegetationdynamictrendsoverrelativelyshortstudyperiodshowever,fewlongtermstudiesfromthe1980s
torecentyearshavebeencompletedbecausethenormalizeddifferencevegetationindex(NDVI)datasetslimit
ed.Additionalytherearefewquantitativeanalysesoftheconsistencyofvegetationtrendsovertime.Inthis
study,thelatestGIMMSNDVI3gdatasetsfrom1982to2012wereusedtoidentifyanyspatiotemporalpatterns
invegetationchangesoverthegrowingseasonandseasonalNDVIatregional,biomeandpixelscalesinthe
TRHR.Toexplorethechangeprocesslinearregressionoverprogressivelylongerperiodsfromtheinitialyear
ofthestudywasappliedtodetectvegetationgrowthchangesinTRHR.Attheregionalscale,althoughaposi
tivetrendofgrowingseasonandseasonalNDVIwerefoundduringpastthreedecades,thereweretwodistinct
periodswithoppositetrendsinspringNDVI,asignificantincreaseandasignificantdecreasebeforeandafter
1998,respectively.Asthestudyperiodincreases,theratesofNDVIincreaseinthespringandsummergrow
ingseasonmarkedlyreduced,butinautumntheNDVIincreased.However,significantNDVIincreaseswere
foundinmostseasons.TheseasonandperiodwithsignificanttrendinNDVIdifferedamongvegetationtypes.
SignificantincreasesinNDVIwereobservedinspringduringthefirstfewperiodsforconiferousforestand
shrubandinspringandsummerduringthefirstfewperiodsfordesertandmeadowareas.Thedifferencein
persistenceofNDVItrendscausedvariationinthecontributionofseasonalvegetationactivitytoplantgrowth.
Duringthefirstfewperiods,thelargestNDVIincreaseandthelargestcontributionweregeneralyfoundin
springandinsummer,respectively.Duringthelastfewperiods,boththelargestNDVIincreaseandcontribu
tionwerefoundinautumn.ItcanbeinferredthatincreasesinNDVIduringthegrowingseasonwasmainlydue
toboththeearlierstartofvegetationgrowthandacceleratedvegetationactivityduringpreviousperiodsand
fromboththeacceleratedvegetationgrowthandadelayintheendofthegrowingseasonduringlatterperiods.
犓犲狔狑狅狉犱狊:GIMMSNDVI3g;regionalscale;pixelscale;vegetationtypes;seasonalNDVI;vegetationactivi
ty;theThreeRiverHeadwaterRegion
三江源地区地处青藏高原腹地,是长江、黄河、澜沧江的发源地,是我国淡水资源的重要补给地,也是亚洲、北
半球乃至全球气候变化的敏感区和重要启动区[12]。近年来,区域的气候变化以波动中显著增温、降水量增加为
主[1,34],这些变化的幅度比周边区域更为剧烈[56],对生物圈产生了深远影响。随着社会经济的发展,人类活动在
逐渐增强[7],冬虫夏草采挖、滥垦乱伐、矿产资源开采[8]等不合理的人类活动破坏当地的草地生态系统。为此,
2003年成立了国家级三江源自然保护区,总面积达15.23km2,占三江源区总面积的42%;2005年,实施了《三江
源生态保护与建设总体规划》,投资75亿元启动了生态保护与生态建设工程。2011年11月,国务院批准实施
《青海三江源国家生态保护综合试验区总体方案》,继续巩固和维护高原生态屏障功能[9]。研究气候变化和人类
活动的叠加作用下,三江源地区植被生长如何变化,对于合理开展区域草地生态保护与恢复建设,具有重要的现
实指导意义。
为此,广大学者开展了大量相关研究,但受限于数据时间序列长度以及研究方法上的差异,出现了研究结果
不一致甚至截然相反的问题[2]。现有利用归一化差分植被指数NDVI(NDVI是植被活动的有效表征指标之一)
的研究主要存在以下两点不足:一是缺乏30年视角的长时间序列研究,包含2006年之后的研究基本没有分析
2000年之前的植被状况[2,8,1012];利用1982-2006年GIMMSNDVI数据集的研究又大多缺少2006年之后植被
生长的评估[6,1314];二是研究多在单一时段,研究时段总体的线性趋势并不能完全准确地描述植被生长变化的实
际动态模式[1517]。此外,关注植被活动变化趋势持续性[1819]的研究相对较少。植被生长是一个对外界环境、干
扰不断适应的动态平衡过程,植被活动强度具有地域性和相对性,需要从时间序列角度进行研究。长时间序列的
NDVI数据是了解植被过去、监测现状、迎接未来挑战的强有力工具[2021]。2014年,最新的AVHRRNDVI3g数
2 ACTAPRATACULTURAESINICA(2016) Vol.25,No.1
据集发布[22],并开展了数据质量检查工作,结果显示数据质量适合开展区域植被动态监测[2324]。有关该数据集
在三江源区的应用尚未见到报道。因此,本文采用1982-2012年的GIMMSNDVI3g数据集,分别在区域尺度、
植被类型尺度、像元尺度,研究了三江源区生长季、春季、夏季、秋季植被活动在1982-1999年、1982-2000年、
……、1982-2012年14个时段的变化特征,以期深入了解三江源区草地生态系统变化过程与特征。
1 材料与方法
1.1 研究区概况
三江源区地处青藏高原腹地、青海省南部,地理位置为北纬31°39′-36°12′,东经89°45′-102°23′,总面积为
36.3万km2。三江源区多年平均气温变化范围在-5.6~3.8℃,其中最热月7月平均气温为6.4~13.2℃,最冷
月1月为-6.6~-13.8℃。年平均降水量262.2~772.8mm,其中6-9月降水量约占全年降水量的75%,而
夜雨量比例则达55%~66%。年蒸发量在730~1700mm之间。日照百分率为50%~65%,年日照时数2300~
2900h,三江源区以山地地貌为主,地势总体呈现西高东低的变化态势,平均海拔约4500m。高山草甸和高寒草
原是三江源地区主要植被类型。土壤类型以高山草甸土为主,沼泽化草甸土也较为普遍,冻土层极为发育。
区域分布有我国面积最大的江河源区和海拔最高的天然湿地,生态系统水源涵养、气候调节等服务功能十分
重要。由于生态效应的开放性、共享性和外溢性,加上特殊的地理位置和丰富的自然资源,使三江源区成为我国
青藏高原生态安全屏障的重要组成部分,关系到全国的生态安全和中华民族的长远发展。由于自然条件严酷、气
候恶劣,三江源地区生态系统结构相对简单,系统内物质、能量和信息流动缓慢,生态系统十分脆弱和敏感。
1.2 数据来源与处理
GIMMSNDVI3g数据集来源于 NASA 戈达德航天中心,合成时段均为15d,空间分辨率为0.083°×
0.083°,时间跨度是1981-2012年。采用地理信息系统软件,NDVI数据经过子集提取、图像镶嵌、裁剪数据格
式转换、投影转换及质量检验等预处理过程,形成三江源区GIMMSNDVI3g数据集。采用最大值合成方法得到
月尺度GIMMSNDVI数据,以进一步去除云的影响,并减少月内物候循环的影响[25]。参考相关研究[1920,2627],
采用1982-2012年多年生长季NDVI平均值0.1作为植被阈值,剔除NDVI值小于0.1的区域,排除非植被因
素的影响。
1∶1000000矢量化植被类型图来自中国植被类型图[28],将三江源区的植被主要分为针叶林、灌丛、荒漠、草
原、草甸5大类。
1.3 研究方法
为反映植被变化的年际和年内特征,采用生长季(4-10月)、春季(4-5月)、夏季(6-8月)和秋季(9-10
月)NDVI合成值来表征植被生长,各季节分别为时段内月NDVI的平均值。在区域尺度和像元尺度对NDVI与
年份进行最小二乘法回归分析,得到回归方程的斜率和Pearson相关系数,分别用来表示植被生长的变化率和变
化趋势。以NDVI变化率与多年平均的比值作为相对变化率[29],反映相对于多年平均值的变化速率和变化强烈
程度。以季节NDVI变化率与各季节月份个数的乘积计算NDVI变化总量,并以此反映各季节NDVI变化对生
长季NDVI变化的贡献。
为探讨三江源地区植被变化的动态过程,分别在1982-1999年、1982-2000年、……、1982-2012年14个
时间段上计算NDVI的变化趋势。时段划分主要是考虑到时段的总长度以及目前已有研究采用的时段状况。
显著增加、显著减少区域面积在14个时段的变化趋势、强度,采用其与1999-2012年年份的Pearson相关系数、
斜率来计算。根据显著性检验结果,将变化趋势分为如下3个等级:极显著(犘<0.01);显著(犘<0.05);不显著
(犘>0.05)。
2 结果与分析
2.1 区域尺度
生长季、各季节NDVI在14个时段均呈明显的波动增加趋势(图1),但显著增加的时段相对较少。生长季
3第25卷第1期 草业学报2016年
NDVI仅在前4个时段(1982-1999年、1982-2000年、1982-2001年、1982-2002年)、1982-2011年增加趋
势显著;春季植被活动仅在第一个时段呈显著增加,夏季植被活动仅在第二个时段显著增加。春季植被活动的变
化明显分为两个阶段:1982-1998年显著增加(犚2=0.36,狀=19),1998-2012年显著减小(犚2=0.34,狀=15)。
5年滑动平均显示,各季节NDVI均呈波动增加趋势,其中,夏季波动幅度较大,生长季和秋季NDVI增加显著
(犚2=0.18,犚2=0.39,狀=27)。
图1 各季节和生长季区域平均犖犇犞犐变化趋势
犉犻犵.1 犐狀狋犲狉犪狀狀狌犪犾狏犪狉犻犪狋犻狅狀狊狅犳狊犲犪狊狅狀犪犾犖犇犞犐犻狀狋犺犲犜犺狉犲犲犚犻狏犲狉犎犲犪犱狑犪狋犲狉犚犲犵犻狅狀
 
  生长季和夏季NDVI在14个时段的变化率呈先
图2 生长季和季节犖犇犞犐在14个时段的变化率
犉犻犵.2 犞犪狉犻犪狋犻狅狀狊狅犳狊犲犪狊狅狀犪犾犖犇犞犐犻狀狋犲狉犪狀狀狌犪犾
狊犾狅狆犲犱狌狉犻狀犵14狆犲狉犻狅犱狊
实心代表趋势显著Solidsymbolfor犘<0.05.
快速下降,而后在1982-2009年时段之后有所上升
(图2),总体均呈随时段长度延长极显著减小趋势
(犚2=0.70,犚2=0.52,狀=14)。春季植被活动则呈
随时段延长其 NDVI变化率持续快速下降的趋势
(犚2=0.94,狀=14,犘<0.01),而秋季NDVI变化率
在14个时段中则呈波动中显著增加趋势(犚2=0.38,
狀=14)。
春季NDVI变化率(斜率)在1、5~8时段均为3
个季节中最大的,2~4时段夏季NDVI变化率最大;1
和5时段夏季NDVI变化率排名第二,2~4时段排名
第二的为春季,6~8时段为秋季。9~14时段,秋季
NDVI变化率最大,9~11时段其次为春季,11~14时
段为夏季。相应地,在前8个时段,春季NDVI相对增加率均大于夏季和秋季,在随后的3个时段中,秋季NDVI
相对增加率超过夏季,成为3个季节中NDVI相对增加率最大的季节;在最后3个时段,春季NDVI相对增加率
在3个季节中最小,而秋季依然最大。
生长季NDVI在14个时段相对增加率的最大、最小值分别发生在1982-1999年、1982-2009年(表1)。春
季NDVI的相对增加率呈不断降低趋势,从1982-1999年的最大值持续减少到1982-2012年的最小值。夏季
4 ACTAPRATACULTURAESINICA(2016) Vol.25,No.1
NDVI相对增加率最大、最小值分别发生在1982-2000年和1982-2008年。秋季NDVI的相对增加率基本呈
逐渐增加趋势,最大、最小值分别发生在1982-2010年和1982-2000年。前5个时段,在3个季节中,春季ND
VI相对增加率最大,其次为夏季;6~11时段,春季NDVI相对增加率仍然最大,秋季则次之;12~14时段,秋季、
夏季NDVI相对增加率排名前两位。从3个季节的增加总量来看,前5个时段夏季最大,其次为春季;6~8时
段,春季最大,其次为秋季;9~11时段秋季最大,春季次之;12~14时段夏季和秋季较大。
表1 生长季和各季节犖犇犞犐在14个时段的相对增加率和增加总量
犜犪犫犾犲1 犜犺犲狉犪狋犲狊犪狀犱狋犺犲狋狅狋犪犾犪犿狅狌狀狋狊狅犳犖犇犞犐犮犺犪狀犵犲犻狀犪犾狊犲犪狊狅狀狊犱狌狉犻狀犵犳狅狌狉狋犲犲狀狆犲狉犻狅犱狊
时段
Periods
相对增加率 ThepercentNDVIchanges(%/a)
生长季Growingseason 春季Spring 夏季Summer 秋季Autumn
增加总量 ThetotalamountofNDVIchange(×10-2/a)
春季Spring 夏季Summer 秋季Autumn
1982-1999 0.24 0.60 0.22 0.05 0.18 0.23 0.02
1982-2000 0.23 0.43 0.29 0.03 0.13 0.29 0.01
1982-2001 0.22 0.38 0.19 0.15 0.12 0.19 0.08
1982-2002 0.18 0.35 0.18 0.09 0.11 0.18 0.04
1982-2003 0.14 0.30 0.10 0.07 0.09 0.10 0.04
1982-2004 0.11 0.29 0.04 0.10 0.09 0.04 0.05
1982-2005 0.11 0.25 0.02 0.13 0.08 0.02 0.07
1982-2006 0.09 0.21 0.03 0.10 0.06 0.03 0.05
1982-2007 0.09 0.17 0.02 0.12 0.05 0.02 0.06
1982-2008 0.06 0.14 0.01 0.09 0.04 0.01 0.04
1982-2009 0.06 0.11 0.02 0.07 0.03 0.02 0.03
1982-2010 0.11 0.10 0.08 0.16 0.03 0.08 0.08
1982-2011 0.11 0.08 0.09 0.16 0.02 0.09 0.08
1982-2012 0.09 0.04 0.06 0.16 0.01 0.06 0.08
2.2 植被类型尺度
除针叶林秋季所有时段、草原和荒漠春季后几个时段、草甸秋季前几个时段NDVI呈减少趋势外,其他植被
类型在各季节和所有时段NDVI均呈增加趋势。其中,针叶林在生长季、春季的前几个时段NDVI增加趋势显
著;灌丛仅在春季部分时段增加趋势显著;荒漠则仅在夏季前几个时段显著增加(表2)。与区域尺度的结果类
似,随时段延长,生长季和夏季NDVI的增加率仅在后几个时段有所增大,在靠前的大部分时段呈随时段延长而
减小趋势;春季多呈持续减小趋势,而秋季则多呈显著增大趋势。
2.3 像元尺度
NDVI变化趋势的统计结果显示(表3),生长季、春季和夏季全部时段,秋季绝大多数时段,NDVI呈升高趋
势的区域均大于50%。除秋季前几个时段之外,生长季和各季节NDVI显著提高的区域面积总体上也远大于显
著降低的区域。
除秋季NDVI与年份相关系数大于0的像元随时段长度延长呈极显著增加外(犚2=0.56,狀=14),其他3个
季节呈增加趋势的像元均减少,春季减少趋势达到了0.01的显著性水平(犚2=0.83,狀=14),夏季减少趋势达到
0.05的显著性水平(犚2=0.35,狀=14),生长季减少趋势不显著。随时段延长,各季节NDVI呈显著降低区域的
范围变化分为两类:秋季呈不显著的减少趋势,而生长季、春季和夏季均呈极显著增加趋势(犚2=0.65,犚2=
0.83,犚2=0.58,狀=14),尤其是春季,1982-1999年仅有0.13%的区域NDVI呈显著减少,而1982-2012年
则快速增加到10.93%。生长季NDVI呈显著提高区域随时段延长而不显著增加,夏季此类区域则呈略有增加
后快速减少,后两个时段又快速扩大的趋势,总体呈不显著的减少;春季和秋季NDVI呈显著提高区域随时段延
长而极显著增加(犚2=0.58,犚2=0.89,狀=14),尤其是秋季增加速率较快。
5第25卷第1期 草业学报2016年
表2 不同植被类型14个时段的犖犇犞犐年际变化率
犜犪犫犾犲2 犛狋犪狋犻狊狋犻犮狊狅犳狊犾狅狆犲犻狀犵狉狅狑犻狀犵狊犲犪狊狅狀犪狀犱狊犲犪狊狅狀犪犾犖犇犞犐犳狅狉犱犻犳犳犲狉犲狀狋犫犻狅犿犲狊 ×10-4/a
类型
Types
1982-
1999
1982-
2000
1982-
2001
1982-
2002
1982-
2003
1982-
2004
1982-
2005
1982-
2006
1982-
2007
1982-
2008
1982-
2009
1982-
2010
1982-
2011
1982-
2012
生长季 Growingseason
NF 9.49 8.74 8.73 7.93 6.23 6.38 6.33 5.11 4.39 3.49 3.49 5.19 5.05 3.18
SH 3.48 4.95 5.42 4.76 3.45 2.93 3.16 2.72 2.84 1.85 1.49 3.14 2.98 2.17
DE 6.60 7.25 7.34 6.13 4.79 3.93 3.40 2.64 2.43 1.50 1.12 2.05 2.26 1.59
ST 5.24 5.64 5.96 4.86 3.55 2.91 3.30 2.80 2.42 1.60 1.37 2.81 2.73 1.98
ME 5.03 5.57 5.93 5.18 3.81 3.34 3.83 3.27 3.06 2.08 2.04 4.12 3.83 2.86
春季Spring
NF 28.3020.76 17.66 16.26 15.45 14.76 13.81 12.20 10.53 8.43 9.30 9.27 7.71 6.42
SH 9.17 6.58 6.38 5.76 5.76 6.52 5.61 5.32 5.67 5.35 4.45 4.56 3.69 2.20
DE 6.73 5.47 5.78 4.90 4.24 4.39 3.40 2.89 2.59 2.11 0.99 0.44 0.11 -0.77
ST 7.73 5.05 4.46 3.69 3.46 3.32 3.55 3.13 2.34 1.98 1.00 0.92 0.36 -0.31
ME 12.20 8.09 6.72 6.22 6.20 6.75 5.97 5.43 5.22 4.83 4.52 4.97 3.78 2.17
夏季Summer
NF 11.70 17.09 11.45 11.95 8.14 5.56 4.87 4.81 2.97 0.93 1.95 5.38 6.47 3.25
SH 8.00 13.90 10.00 9.58 4.99 2.02 1.09 1.74 0.30 -0.85 -0.02 2.75 3.23 2.51
DE 11.52 14.4011.5710.43 7.16 5.31 3.87 3.61 2.65 2.19 2.39 3.84 4.50 3.45
ST 8.39 12.09 9.25 9.05 5.72 3.02 2.19 2.38 1.40 0.53 1.46 3.59 4.20 2.85
ME 9.56 15.33 10.63 10.09 5.16 1.98 1.23 1.80 0.50 -0.82 0.33 3.28 3.76 3.05
秋季 Autumn
NF -11.52 -11.64 -2.92 -4.43 -4.89 -1.18 0.33 -1.68 -0.34 1.10 -0.78 0.93 0.97 -0.13
SH -6.73 -5.62 -0.14 -1.07 -0.40 0.27 2.79 1.09 2.55 1.05 0.04 2.10 2.03 1.82
DE 1.54 1.88 4.67 3.07 2.97 2.08 2.93 1.41 2.06 0.19 0.00 1.86 2.15 2.09
ST -0.40 -0.23 4.18 1.85 1.48 2.40 4.15 2.89 3.50 2.29 1.66 3.91 3.63 3.39
ME -6.68 -6.71 0.44 -0.76 0.07 1.29 4.28 2.58 3.45 2.23 1.26 4.11 3.95 3.36
 NF:针叶林 Needleforest,SH:灌木林Shrub,DE:荒漠Desert,ST:草原Steppe,ME:草甸 Meadow.犘<0.05,犘<0.01(Pearson相关
线性回归的统计显著性StatisticalsignificanceoflinearregressionofPearsoncorrelation).
生长季植被活动显著增强的像元1982-1999年主要集中分布在研究区域东部和西北部,随着时段延长,尤
其是2006年之后,植被活动显著增强的区域较为分散;显著减少的区域零星分布于研究区域(图3)。春季植被
活动显著增加的区域前几个时段主要集中分布在研究区域东南部以及囊谦县,随着时段的延长,显著增加区域先
是有所减少,后在研究区域中部迅速增加,至1982-2012年,显著增加区域主要分布在中部的曲麻莱县、治多县
东南部、囊谦县等区域;显著减少区域主要集中在三江源区的西北部以及兴海县、同德县、泽库县部分地区(图
4)。夏季植被NDVI呈显著增加的趋势先是主要集中在研究区域东部玛沁县、河南县、甘德县部分地区以及三
江源区西北部(1982-1999年、1982-2002年),随后显著增加区域快速萎缩,并分散,而1982-2012年显著增加
区域主要分布在研究区域西部和南部部分地区,显著减少的区域则主要位于三江源区东部的黄河源区。1982-
1999年,秋季植被活动显著增加的区域较小,主要集中分布在三江源区的东部,随着时段延长,西部地区也有部
分区域显著增加;1982-1999年显著减少区域主要分布在中部地区,随后散布在三江源区。
1982-1999年、1982-2002年、1982-2006年、1982-2006年和1982-2012年共5个时段生长季、各季节
NDVI变化量见图5、图6(春季和秋季图略),(-0.001,0]和(0,0.001]区间像元之和占总像元数的比例见表4。
6 ACTAPRATACULTURAESINICA(2016) Vol.25,No.1
图3 生长季犖犇犞犐变化趋势及其显著性
犉犻犵.3 犛狆犪狋犻犪犾犱犻狊狋狉犻犫狌狋犻狅狀狅犳犖犇犞犐狋狉犲狀犱狊犪狀犱犻狋狊狊犻犵狀犻犳犻犮犪狀狋犾犲狏犲犾犻狀犵狉狅狑犻狀犵狊犲犪狊狅狀犱狌狉犻狀犵狋犺犲犳犻狏犲狆犲狉犻狅犱狊
 
生长季和各季节NDVI变化量的数值分布随时段延长的变化趋势总体一致,即变化量前期基本呈正态分
布,(0,0.001]区间或(-0.001,0]区间(秋季)的像元比例最大,并以其为中心,左右两侧共有2~3个区间的像
元比例超过10%;随着年份增加,(0,0.001]区间或(-0.001,0]区间的像元比例逐渐增加,而其他变化量绝对
值相对较大的区域则逐渐减少,至1982-2012年,(0,0.001]区间和(-0.001,0]区间的像元比例之和在80%
左右,且其他NDVI变化量区间像元比例均小于15%。
3 讨论
3.1 植被活动变化趋势
尽管不同研究的起始时段不同、采用方法存在差异,但三江源区植被生长总体趋于改善的结果是基本一致
的[2,6,1011,13,3031],本文的结果也与此一致。尽管2005年三江源区开始的生态保护与建设工程对促进三江源区植
被恢复具有重要作用[2,6,10],但若要定量生态工程对植被改善的贡献率,需要详细的工程地点、规模、实施过程等
方面的资料来验证。春季NDVI增加趋势变缓,可能与春季物候提前趋势变弱[26,3234]、草地退化[35]、春季大气污
染[36]等因素有关。
7第25卷第1期 草业学报2016年
图4 春季犖犇犞犐变化趋势及其显著性
犉犻犵.4 犛狆犪狋犻犪犾犱犻狊狋狉犻犫狌狋犻狅狀狅犳犖犇犞犐狋狉犲狀犱狊犪狀犱犻狋狊狊犻犵狀犻犳犻犮犪狀狋犾犲狏犲犾犻狀狊狆狉犻狀犵犱狌狉犻狀犵狋犺犲犳犻狏犲狆犲狉犻狅犱狊
 
1982-1999年春季是3个季节中增加量和增加率最大的季节,其次是夏季,与文献[29]青藏高原的结果一
致。但随着春季NDVI增速的放缓、夏季NDVI增速在前几个时段的大幅增加和秋季NDVI增速的持续增加,
到了1982-2012年,相对大小的排序变为秋季、夏季和春季,反映了生长季NDVI增加在各季节之间分配的变
化。
季节NDVI对生长季NDVI增加的贡献随时段长度增加而发生变化。从最初的夏季、春季NDVI变化对生
长季NDVI变化贡献率最大,逐渐过渡到秋季、夏季贡献率最大。
季节之间NDVI变化趋势的差异表明,三江源区前期(1982-2000s)植被NDVI的提高可能主要依赖于生长
季开始提前和生长加速,而最近十年生长季结束延迟和生长加速对NDVI增加的贡献更大。春季NDVI的变化
趋势与相关研究得出的1990s中期之前返青期呈提前趋势,之后呈推迟趋势[26,3334]可能导致的结果是一致的。
秋季NDVI持续增加的趋势与青藏高原生长季结束一再延迟[3738]的结果一致。NDVI变化趋势持续性季节之间
的差异,可能与不同季节气候变化趋势强弱不同、草地退化导致生长季起始时间发生等有关。
8 ACTAPRATACULTURAESINICA(2016) Vol.25,No.1
表3 生长季和各季节犖犇犞犐变化趋势统计
犜犪犫犾犲3 犜犺犲狆狉狅狆狅狉狋犻狅狀狅犳狆犻狓犲犾狊犮犺犪狉犪犮狋犲狉犻狕犲犱犫狔犱犻犳犳犲狉犲狀狋狊狋犪狋犻狊狋犻犮犪犾狊犻犵狀犻犳犻犮犪狀犮犲犾犲狏犲犾狊犱狌狉犻狀犵犳狅狌狉狋犲犲狀狆犲狉犻狅犱狊 %
时段Periods
生长季Growingseason
P SN SP
春季Spring
P SN SP
夏季Summer
P SN SP
秋季 Autumn
P SN SP
1982-1999 69.06 1.34 12.07 84.35 0.13 13.78 74.40 0.88 13.85 40.74 4.23 1.32
1982-2000 76.83 0.46 14.83 77.59 0.31 7.23 88.79 0.38 24.19 41.68 3.94 1.32
1982-2001 80.46 0.34 17.11 80.59 0.27 8.80 84.40 0.40 15.14 56.96 1.15 2.64
1982-2002 78.14 0.50 15.46 77.65 0.42 8.98 84.02 0.36 15.94 51.64 1.82 2.35
1982-2003 73.61 0.84 11.43 77.17 0.80 10.97 73.86 0.59 7.73 51.77 1.38 2.43
1982-2004 70.49 1.45 11.35 79.87 0.86 14.70 63.90 1.40 4.69 52.10 2.37 3.50
1982-2005 73.53 1.30 13.61 77.13 1.21 14.97 60.06 1.93 4.63 61.86 1.42 5.68
1982-2006 69.61 1.72 12.92 76.38 1.76 16.36 61.84 1.65 5.47 54.72 2.05 4.75
1982-2007 68.65 2.03 12.98 74.41 2.79 19.25 55.56 2.12 4.06 58.60 1.95 6.53
1982-2008 62.68 2.85 10.20 72.23 3.73 21.21 49.30 3.02 3.46 51.69 2.89 5.72
1982-2009 61.95 2.95 9.70 65.03 5.95 17.99 55.58 2.47 4.77 49.57 2.53 4.57
1982-2010 75.60 1.55 17.34 64.36 6.95 20.86 69.57 1.26 9.70 63.81 0.88 8.04
1982-2011 75.12 2.03 19.23 61.24 8.98 18.81 72.19 1.57 14.09 65.93 1.09 8.38
1982-2012 69.03 2.93 17.38 55.25 10.93 14.72 67.06 2.76 14.18 65.19 1.38 9.30
 P为NDVI变化趋势大于0的像元比例;SN为变化趋势小于0且显著性水平达到0.05的像元比例;SP为变化趋势大于0且显著性水平达到0.05
的像元比例。Pstandfor犚>0,SNstandfor犚<0andstatisticalsignificanceoflinearregressionofPearsoncorrelationless0.05,SPstandfor犚>
0andstatisticalsignificanceoflinearregressionofPearsoncorrelationless0.05.
图5 部分时段生长季犖犇犞犐变化量
数值分布
犉犻犵.5 犜犺犲狀狌犿犲狉犻犮犪犾犱犻狊狋狉犻犫狌狋犻狅狀狅犳犖犇犞犐
狊犾狅狆犲犻狀犵狉狅狑犻狀犵狊犲犪狊狅狀犱狌狉犻狀犵狊狅犿犲狆犲狉犻狅犱狊
   A:≤-0.002;B:(-0.001,0.002];C:
(-0.001,0];D:(0,0.001];E:(0.001,
0.002];F:(0.002,0.004];G:>0.004.
下同Thesamebelow.
3.2 嵌套时段对结果的影响
随时段延长,显著减少、显著增加区域的增加,部分原因是否来自于统计学上关于稳态系统序列长度的推论,
即随着年份逐渐加入、时间序列逐渐延长,理论上变化趋势几乎确定更小,除非在加入的年份中发生较大的一些
变化,即统计学上的“可塑性时间单元问题”(themodifiabletemporalunitproblem,MTUP)[3940]。
9第25卷第1期 草业学报2016年
图6 部分时段夏季犖犇犞犐变化量数值分布
犉犻犵.6 犜犺犲狀狌犿犲狉犻犮犪犾犱犻狊狋狉犻犫狌狋犻狅狀狅犳犖犇犞犐
狊犾狅狆犲犻狀狊狌犿犿犲狉犱狌狉犻狀犵狊狅犿犲狆犲狉犻狅犱狊
 
  植被活动趋势研究需要 NDVI长期时间序列数
据,且基于不同的开始年,结果可能不同。本文发展的
固定起始年份、不断增加结束年份计算趋势的方法,
对于评价植被的长期变化趋势、识别生态系统退化是
一种适宜的方法:固定的年份可以作为一个恒定的参
考、基准,以趋势变化的持续性评判植被活动状况。尽
管采用时段不断延长计算NDVI趋势的方法,可能会
引起统计上的质疑,对区域尺度NDVI增加量减少有
一定影响,但为了反映植被NDVI变化趋势的持续性
和鲁棒性(Robustness),采用嵌套的时间序列来评价
表4 犖犇犞犐变化量位于(-0.001,0]和
(0,0.001]区间像元比例
犜犪犫犾犲4 犘犻狓犲犾犳狉犪犮狋犻狅狀狑犻狋犺犖犇犞犐狊犾狅狆犲犻狀狋犺犲狉犪狀犵犲狅犳
(-0.001,0]犪狀犱(0,0.001] %
季节
Seasons
1982-
1999
1982-
2002
1982-
2006
1982-
2009
1982-
2012
生长季Growingseason70.47 77.91 86.58 92.65 92.29
春季Spring 61.86 76.86 78.83 80.92 85.84
夏季Summer 47.18 51.16 77.99 85.21 79.73
秋季 Autumn 52.25 64.46 71.48 78.43 79.98
植被的变化趋势仍不失为一种可供选择的有效方法。
4 结论
本文利用最新发布的1982-2012年的GIMMSNDVI3g数据集,采用固定起始年份、结束年份逐渐增加的
多个嵌套时间序列的分析方法,在区域、植被类型和像元3个空间尺度上,对近30年来三江源区生长季、春季、夏
季和秋季植被活动变化趋势和过程进行了分析。主要结论如下:
1)三江源区生长季和春夏秋各季节植被活动在近30年均呈增强趋势,1982-2012年NDVI增加趋势分别
为每年0.0002,0.0001,0.0002和0.0004。春季植被变化具有明显的阶段性,1998年前后分别呈显著增加、减少
趋势,变化率分别为每年0.0011和-0.0008。
2)随时段延长,生长季和夏季NDVI的增加随时段延长而趋缓,春季增加趋势迅速减弱接近停滞,而秋季植
被活动则持续增强。生长季和春季植被显著改善的区域由1982-1999年主要集中在研究区域东南部、玉树州东
南部,逐渐向三江源区全区域扩展。
3)尽管区域尺度生长季、春季和夏季植被活动趋缓,但在像元尺度显著变化的区域范围却多呈增加趋势,尤
其是呈显著减少的区域,呈极显著扩大趋势。区域尺度春季 NDVI增长的持续趋缓主要是由于像元尺度呈减
少、显著减少的范围快速扩大。秋季NDVI增长率随时段延长增加,则主要是由显著增加的范围快速扩大和显
著减少范围减小所致。
01 ACTAPRATACULTURAESINICA(2016) Vol.25,No.1
4)不同植被类型显著变化的季节、时段存在差异。各季节草原在所有时段变化均不显著,针叶林和灌丛
NDVI主要在春季前几个时段增加显著,荒漠、草甸主要在春季、夏季前几个时段增加幅度较大。
5)各季节NDVI在14个时段增加趋势持续性的差异,引起各季节对生长季植被活动增加贡献的变化。前几
个时段,NDVI变化幅度较大的季节是春季,对生长季NDVI变化贡献最大的是夏季;随时段延长,逐渐过渡到秋
季NDVI变化幅度、变化量均最大,夏季次之。
犚犲犳犲狉犲狀犮犲狊:
[1] ShaoQQ,FanJW.ComprehensiveMonitoringandAssessmentofEcosystemintheThreeRiverHeadwaterRegion[M].
Beijing:SciencePress,2012.
[2] LiuHX,LiuGH,FuBJ.ResponseofvegetationtoclimatechangeandhumanactivitybasedonNDVIintheThreeRiver
Headwatersregion.ActaEcologicaSinica,2011,31(19):54955504.
[3] HuangL,LiuJY,ShaoQQ,犲狋犪犾.ChanginginlandlakesrespondingtoclimatewarminginNortheasternTibetanPlateau.
ClimaticChange,2011,109(34):479502.
[4] YangK,YeBS,ZhouDG,犲狋犪犾.ResponseofhydrologicalcycletorecentclimatechangesintheTibetanPlateau.Climate
Change,2011,109(34):517534.
[5] QinJ,YangK,LiangSL,犲狋犪犾.ThealtitudinaldependenceofrecentrapidwarmingovertheTibetanPlateau.Climatic
Change,2009,97(12):321327.
[6] XuWX,GuS,ZhaoXQ,犲狋犪犾.HighpositivecorrelationbetweensoiltemperatureandNDVIfrom1982to2006inalpine
meadowoftheThreeRiverSourceRegionontheQinghaiTibetanPlateau.InternationalJournalofAppliedEarthObservation
andGeoinformation,2011,13:528535.
[7] LiuJY,XuXL,ShaoQQ.ThespatialandtemporalcharacteristicsofgrasslanddegradationintheThreeRiverHeadwaters
RegioninQinghaiProvince.ActaGeographicaSinica,2008,63(4):364376.
[8] ZhangC,LiQ,LiZF.InfluenceofhumanactivitiesonvariationofvegetationcoverintheThreeRiverSourceRegion.China
Population,ResourcesandEnvironment,2014,24(5):139144.
[9] SunHL,ZhengD,YaoTD,犲狋犪犾.ProtectionandconstructionofthenationalecologicalsecurityshelterzoneonTibetan
Plateau.ActaGeographicaSinica,2012,67(1):312.
[10] LiuXF,RenZY,LinZH,犲狋犪犾.ThespatialtemporalchangesofvegetationcoverageintheThreeRiverHeadwatersRe
gioninrecent12years.ActaGeographicaSinica,2013,68(7):897908.
[11] XiaoT,WangCZ,FengM,犲狋犪犾.DynamiccharacteristicofvegetationcoverageintheThreeRiverSourceRegionfrom
2000to2011.ActaAgrestiaSinica,2014,22(1):3945.
[12] WangSH,JiaSF,LvAF.TherelationshipbetweenNDVIandresidentialsitesacross“ThreeRiverSource”Area.Re
sourcesScience,2012,34(11):20452050.
[13] HuMQ,MaoF,SunH,犲狋犪犾.Studyofnormalizeddifferencevegetationindexvariationanditscorrelationwithclimatefac
torsinthethreeriversourceregion.InternationalJournalofAppliedEarthObservationandGeoinformation,2011,13(1):
2433.
[14] QianS,FuY,PanFF.ClimatechangetendencyandgrasslandvegetationresponseduringthegrowthseasoninThreeRiver
SourceRegion.ScienceChinaEarthSciences,2010,40(10):14391445.
[15] BeckPS,GoetzSJ.Sateliteobservationsofhighnorthernlatitudevegetationproductivitychangesbetween1982and2008:
ecologicalvariabilityandregionaldifferences.EnvironmentalResearchLetters,2011,6(4):55015511.
[16] WesselsKJ,PrinceSD,MalherbeJ,犲狋犪犾.Canhumaninducedlanddegradationbedistinguishedfromtheeffectsofrainfal
variability:AcasestudyinSouthAfrica.JournalofAridEnvironments,2007,68(2):271297.
[17] DeJR,VerbesseltJ,SchaepmanME,犲狋犪犾.Trendchangesinglobalgreeningandbrowning:contributionofshortterm
trendstolongertermchange.GlobalChangBiology,2012,18(2):642655.
[18] PengSS,ChenAP,XuL,犲狋犪犾.RecentchangeofvegetationgrowthtrendinChina.EnvironmentalResearchLetters,
2011,6(4):40274039.
[19] MohammatA,WangXH,XuXT,犲狋犪犾.DroughtandspringcoolinginducedrecentdecreaseinvegetationgrowthinInner
Asia.AgriculturalandForestMeteorology,2013,178179(9):2130.
[20] PiaoSL,WangXH,CiaisP,犲狋犪犾.ChangesinsatelitederivedvegetationgrowthtrendintemperateandborealEurasia
from1982to2006.GlobalChangeBiology,2011,17(10):32283239.
[21] vanLeeuwenW,OrrBJ,MarshSE,犲狋犪犾.MultisensorNDVIdatacontinuity:Uncertaintiesandimplicationsforvegeta
tionmonitoringapplications.RemoteSensingofEnvironment,2006,100:6781.
11第25卷第1期 草业学报2016年
[22] PinzonJ,TuckerC.Anonstationary1981-2012AVHRRNDVI3gtimeseries.RemoteSensing,2014,6(8):69296960.
[23] WangJ,DongJ,LiuJ,犲狋犪犾.ComparisonofgrossprimaryproductivityderivedfromGIMMSNDVI3g,GIMMS,andMO
DISinSoutheastAsia.RemoteSensing,2014,6(3):21082133.
[24] ZengFW,ColatzG,PinzonJ,犲狋犪犾.Evaluatingandquantifyingtheclimatedriveninterannualvariabilityinglobalinventory
modelingandmappingstudies(GIMMS)normalizeddifferencevegetationindex(NDVI3g)atglobalscales.RemoteSensing,
2013,5(8):39183950.
[25] FensholtR,ProudSR.Evaluationofearthobservationbasedgloballongtermvegetationtrends———ComparingGIMMSand
MODISglobalNDVItimeseries.RemoteSensingofEnvironment,2012,119(16):131147.
[26] PiaoSL,CuiMD,ChenAP,犲狋犪犾.Altitudeandtemperaturedependenceofchangeinthespringvegetationgreenupdate
from1982to2006intheQinghaiXizangPlateau.AgriculturalandForestMeteorology,2011,151(12):15991608.
[27] ZhangGL,ZhangYJ,DongJW,犲狋犪犾.GreenupdatesintheTibetanPlateauhavecontinuouslyadvancedfrom1982to
2011.ProceedingoftheNationalAcademyofScience,2013,110(11):43094314.
[28] EditorialCommitteeforVegetationMapofChina.VegetationAtlasofChina[M].Beijing:SciencePress,2001.
[29] YangYH,PiaoSL.VariationsingrasslandvegetationcoverinrelationtoclimaticfactorsontheTibetanPlateau.Journalof
PlantEcology,2006,30(1):18.
[30] XiaoT,LiuJY,ShaoQQ.AsimulationonchangesinvegetationproductivityinThreeRiverSourcesNatureReserve,
QinghaiProvinceoverpast20years.JournalofGeoinformationScience,2009,11(5):557565.
[31] FanJW,ShaoQQ,LiuJY,犲狋犪犾.DynamicchangesofgrasslandyieldinThreeRiverHeadwaterRegionfrom1988to
2005.ActaAgrestiaSinica,2010,18(1):511.
[32] YuHY,LuedelingE,XuJC.WinterandspringwarmingresultindelayedspringphenologyontheTibetanPlateau.Pro
ceedingoftheNationalAcademyofScience,2010,107(51):2215122156.
[33] YuHY,XuJC,ErickO,犲狋犪犾.SeasonalresponseofgrasslandstoclimatechangeontheTibetanPlateau.PLoSONE,
2012,11(7):e49230.
[34] ShenMG,TangYH,ChenJ,犲狋犪犾.Influenceoftemperatureandprecipitationbeforethegrowingseasononspringphenol
ogyingrasslandsofthecentralandeasternQinghaiTibetanPlateau.AgriculturalandForestMeteorology,2011,151(12):
17111722.
[35] ChenH,ZhuQA,WuN,犲狋犪犾.DelayedspringphenologyontheTibetanPlateaumayalsobeattributabletootherfactors
thanwinterandspringwarming.ProceedingoftheNationalAcademyofScience,2011,108(19):E93.
[36] YiSH,ZhouZY.IncreasingcontaminationmighthavedelayedspringphenologyontheTibetanPlateau.Proceedingofthe
NationalAcademyofScience,2011,108(19):E94.
[37] JeongSJ,HoCH,Gim HJ,犲狋犪犾.Phenologyshiftsatstartvs.endofgrowingseasonintemperatevegetationoverthe
NorthernHemispherefortheperiod1982-2008.GlobalChangeBiology,2011,17(7):23852399.
[38] CheML,ChenBZ,InesJohnL,犲狋犪犾.Spatialandtemporalvariationsintheenddateofthevegetationgrowingseason
throughouttheQinghaiTibetanPlateaufrom1982to2011.AgriculturalandForestMeteorology,2014,189190:8190.
[39] deJongR,deBruinS.Lineartrendsinseasonalvegetationtimeseriesandthemodifiabletemporalunitproblem.Biogeo
sciences,2012,9:7177.
[40] ChengT,AdepejuM.Modifiabletemporalunitproblem(MTUP)anditseffectonspacetimeclusterdetection.PLoSONE,
2014,9(6):e100465.
参考文献:
[1] 邵全琴,樊江文.三江源区生态系统综合监测与评估[M].北京:科学出版社,2012.
[2] 李辉霞,刘国华,傅伯杰.基于NDVI的三江源地区植被生长对气候变化和人类活动的响应研究.生态学报,2011,31(19):
54955504.
[7] 刘纪远,徐新良,邵全琴.近30年来青海三江源地区草地退化的时空特征.地理学报,2008,63(4):364376.
[8] 张罛,李强,李忠峰.三江源地区人类活动对植被覆盖度的影响.中国人口·资源与环境,2014,24(5):139144.
[9] 孙鸿烈,郑度,姚檀栋,等.青藏高原国家生态安全屏障保护与建设.地理学报,2012,67(1):312.
[10] 刘宪锋,任志远,林志慧,等.2000-2011年三江源区植被覆盖时空变化特征.地理学报,2013,68(7):897908.
[11] 肖桐,王昌佐,冯敏,等.2000-2011年青海三江源地区草地覆盖度的动态变化特征.草地学报,2014,22(1):3945.
[12] 王素慧,贾绍凤,吕爱锋.三江源地区植被盖度与居民点的关系研究.资源科学,2012,34(11):20452050.
[28] 中国植被图编辑委员会.1∶1000000中国植被图集[M].北京:科学出版社,2001.
[29] 杨元合,朴世龙.青藏高原草地植被覆盖变化及其与气候因子的关系.植物生态学报,2006,30(1):18.
[30] 肖桐,刘纪远,邵全琴.近20年青海三江源自然保护区植被生产力变化模拟.地球信息科学学报,2009,11(5):557565.
[31] 樊江文,邵全琴,刘纪远,等.1988-2005年三江源草地产草量变化动态分析.草地学报,2010,18(1):511.
21 ACTAPRATACULTURAESINICA(2016) Vol.25,No.1