全 文 ：第 !" 卷第 ## 期
!$$" 年 ## 月
生 态 学 报
%&’% (&)*)+,&% -,.,&%
/012 !"，.02 ##
.032，!$$"
!""#：$ $ %%%& ’()*)+,(-& (.
基金项目：国家自然科学基金资助项目（4$5"6$##）
收稿日期：!$$"7$575#；修订日期：!$$"7$87#8
作者简介：李登科（#895 :），男，陕西人，高级工程师，主要从事卫星遥感应用研究2 (7;<=1：1>?@#"##A B=C<2 D0;
/)0.1-",). ,"’2：’EF GH0IFDJ K OM JEF .
6,)+7-#!8：*, QFCR7SF，-FC=0H FCR=CFFH，;<=C1M FCR
陕北长城沿线风沙区植被指数变化及其与气候的关系
李登科#，郭T 铌!，何慧娟#
（#2陕西省农业遥感信息中心，西安T "#$$#4；!2中国气象局兰州干旱气象研究所，兰州T "5$$!$）
摘要：陕北长城沿线风沙区位于毛乌素沙漠东南部边沿，属毛乌素沙地向东南移动的最活跃地段，生态环境十分脆弱。使用
#8@# : !$$5年 !5<长时间序列的 .)%% U %VWW .Q/,数据、气候资料，分析了陕北长城沿线风沙区植被覆盖的历史演变及其与
气候因子的关系。结果表明：（#）陕北长城沿线风沙区植被覆盖状况 !5< 来尽管有波动起伏，但是整体在持续转好，年平均
.Q/,增加了 #$2 9!X。低覆盖率植被面积在减少，高覆盖率植被面积在增加。夏季的 .Q/,值最高、波动起伏最大，其次是秋
季；春、夏、秋三季的 .Q/,具有明显的上升趋势，季平均 .Q/,年增长率夏季最大，秋季次之；夏、秋季 .Q/,与年 .Q/,具有很
高的相关性，这两个季节的植被状况基本决定了全年的植被分布状况。.Q/,年变化曲线为单峰型，春季 .Q/,缓慢增加，秋季
.Q/,降低速度比较快。（!）年平均 .Q/,与温度的年际变化相关不明显，各季节 .Q/, 与温度相关也不明显。近年来长城沿
线风沙区的年降水量没有明显增加，而年平均 .Q/,线性增加趋势显著，降水量是引起 .Q/,年际波动的主要因子，非气候因素
是年平均 .Q/,线性增加的主要原因。降水量与 .Q/,存在着明显的年相关和隔季相关。年降水量与年 .Q/,的相关，冬季降
水量与春季 .Q/,的相关，春季降水量与夏季 .Q/,的相关，夏季降水量与秋季 .Q/,的相关性都非常高。（5）非气候因素中生
态保护和环境建设等人为措施，如植树造林、草原围栏封育等是导致植被显著增加的重要原因。
关键词：植被；.Q/,；气候；陕北长城沿线风沙区
文章编号：#$$$7$855（!$$"）##749!$7#$T 中图分类号：Y84@T 文献标识码：%
9’+’"-",). (!-.+’ -.1 ,": 7’*-",).:!,# %,"! (*,2-"’ ,. "!’ 7’+,). -*).+ "!’ ;7’-"
<-** ,. .)7"!’7. =!--.>,
*, QFCR7SF#，+Z) .=!，V( VN=7[N
! .$,+&+4+( */ 20&8 9(+(*0*5*-:，192，;#$<"*4 "5$$!$，1"&$#
!"#$ %"&’&()"$ *)+)"$，?@@A，?A（BB）：CD?@ E CD?F,
5G:"7-("：’EF K=C>7O0HC B
B0NJEF
_<0KNBN2 ,CJFH
KFHF F]<;=CF> OM NB=CR JEF ;0CJE1M +,__- （+10O<1 ,C3FCJ0HM _0>F1=CR
Q=LLFHFCDF /FRFJ
!""#：$ $ %%%& ’()*)+,(-& (.
-../-* 0-*/’1& 2!’ 3456 ,.(7’-1’8 1*)%*9 ,. 1#7,.+ -.8 8’(7’-1’8 7-#,8*9 ,. -/"/:.& 2!’ -../-* :’-. 3456 1!)%’8
#)1,",0’ ()77’*-",). %,"! #7’(,#,"-",).，;/" .)" %,"! "’:#’7-"/7’& 2!’7’ %’7’ 1,+.,<,(-." #)1,",0’ ()77’*-",).1 ;’"%’’. %,."’7
#7’(,#,"-",). -.8 1#7,.+ 3456，1#7,.+ #7’(,#,"-",). -.8 1/::’7 3456，-.8 1/::’7 #7’(,#,"-",). -.8 -/"/:. 3456& 2!’
()77’*-",).1 ;’"%’’. 1/::’7 -.8 -/"/:. #7’(,#,"-",).1 -.8 -../-* 3456 %’7’ ;)"! 1,+.,<,(-."*9 #)1,",0’& 2!’ 7’1/*"1
,.8,(-"’ "!-" 7-,.<-** ,1 - =’9 <-(")7 ,. 8’"’7:,.,.+ "!’ 8,1"7,;/",). )< 0’+’"-",). )0’7 "!’ 7’+,).& >../-* #7’(,#,"-",). 8,8 .)"
,.(7’-1’ 1,+.,<,(-."*9 8/7,.+ "!’ 1"/89 #’7,)8，;/" -../-* 3456 ().",./)/1*9 ,.(7’-1’8& 2!,1 ,:#*,’1 "!-" #7’(,#,"-",). %-1
"!’ :-,. 87,0’7 <)7 3456 ,."’7-../-* 0-7,-",).1 -.8 -."!7)#)+’.,( <-(")71 1/(! -1 7’<)7’1"-",). -.8 +7-11*-.8 7’()0’79 %’7’
"!’ (-/1’ <)7 "!’ ,.(7’-1,.+ 3456 "7’.8&
!"# $%&’(：0’+’"-",).；3456；(*,:-"’；3)7"!’7. ?!--.@,
陕北长城沿线风沙区位于陕西省最北部，毛乌素沙漠的东南部边缘（图 A），属毛乌素沙地向东南移动的
最活跃地段［A］。地处毛乌素沙地与黄土高原的过渡地带，地理位置为 BCDEFG H BIDBCG3、AFJD H AAFDK。地势
西北高、东南低，由西北向东南倾斜。海拔为 A FFF H A EFF : ，除白于山地外，大部分地区相对高差仅 BF H
EF:。以长城沿线为界，地貌上可分为北部风沙滩地，属温带半干旱、栗钙土干草原地带，地势平坦，沙丘绵
延。南部为黄土丘陵沟壑区，是陕北黄土高原的北缘，属温带半干旱半湿润、黑垆土森林草原带。地形破碎，
沟壑纵横，黄土层深厚，但结构疏松，侵蚀严重［L］。研究区行政范围包括府谷、神木、榆阳、横山、靖边和定边 C
县（区），总面积 BBIJL=:L，总人口 AMI& E 万人［B］。气候为温带干旱、半干旱向半湿润过渡气候。年平均降水
量 BFF H EFF:: ，降水少且多集中在 J H I 月份，日照时数达 LJNF H LICL!，全年蒸发量 LFFF H LEFF:: ，是降
雨量的 N H E 倍。无霜期 ABN H ACI8［N］。地带性植被从西北向东南由荒漠、荒漠草原向森林草原过渡［E］。在
经济活动方面是典型的农牧业交错带，随着煤炭、天然气的大规模开发，又成为农牧与工矿区的过渡带。多种
界质的叠加，致使生态环境呈现出明显的波动性、多样性和脆弱性［C］。
图 AO 陕北长城沿线风沙区示意图
O P,+& AO 2!’ %,.8Q;)7. 1-.89 7’+,). -*).+ "!’ R7’-" S-** ,. .)7"!’7.
?!--.@,O
近年来，该区域实施了种草种树、恢复和保护植被、
综合治理风蚀沙化的多项生态建设工程［J］，主要是扩
大林草比重，把护田林网、片林、人工饲草基地建设结合
起来，扩大地表绿色覆盖面积，削弱风沙活动强度，减少
生态系统中非生产性输出，调节生态系统能量和物质循
环，促进生态系统逐步平衡。限定旱作农业范围，把旱
作耕种区局限在地势开阔、水土条件较好的滩川地段。
改良天然草场，选择自然条件较好的河流沿岸及湖盆滩
地建立人工饲草基地，坚持以草定畜，实施有计划的轮
牧，实行草场适度利用，防止过牧引起草场退化。建立
绿洲防护体系，以扩大绿洲规模来分割、包围流沙，减弱
和遏止风蚀沙化过程。这些工程的实施，生态效益如
何，尤其是植被覆盖发生了怎样的变化，这是迫切需要
回答的问题。
与陕北长城沿线风沙区相关的区域植被变化与气候变化研究已有很多。马明国等［M］利用美国 KT? 探路
者（U-"!<,.8’7）提供的 AIMA H LFFA 年 LA- 3T>> $ >5VWW 3456 数字影像研究发现，中国西北植被覆盖在近
LA-里存在普遍退化的趋势。朴世龙等［I］利用相同的数据集研究发现植被退化的地区主要分布在西北地区
和青藏高原。李震等［AF］通过分析该数据集，得出中国西北地区大部分地区植被状况恶化，局部地区有所好
转；3456与降水存在明显的正相关关系，而与温度变化的关系并不明显。方精云等［AA］利用 R6XX?的 AIML H
AIII 年 3456 $ >5VWW数字影像，对中国近 LF-的植被覆盖变化的研究认为，中国西部地区的植被大都呈增
ALCNO AA期 O O O 李登科O 等：陕北长城沿线风沙区植被指数变化及其与气候的关系 O
!""#：$ $ %%%& ’()*)+,(-& (.
加趋势。李珍存等［/0］利用 /120 3 0445 年的 00-长时间系列的遥感 6789 数据，采用景观生态学方法分析得
出中国西北植被覆盖总体是向改善的趋势发展的，局部有退化现象。李月臣等［/5］基于 6:;; $ ;8<== 6789
数据、气候和土地覆盖数据对北方 /5 省（市、区）/120 3 /111 年植被动态变化及其与气候因子的关系进行了
分析，结果表明 /2-间研究区植被总体呈现增加趋势，植被变化与气温相关性显著而与降水无显著相关。张
莉等［/>］对 /12? 年和 /111 年两个时期的 @-.AB-" CD影像解译分析认为，土地沙漠化程度恶化的县主要是横
山县，趋于稳定或好转的区域是长城沿线以东的榆阳区、神木县等。这些研究因资料序列长度不同而得出不
同结论，甚至结论截然相反；或因研究范围太大，而不能具体涉及陕北长城沿线风沙区的植被变化情况和与气
候变化的关系，对于陕北长城沿线风沙区的研究资料时间序列短，所得的结论具有一定局限性。本文试图通
过遥感图像处理建立陕北长城沿线风沙区 /12/ 3 0445 年 6789时间序列，采用相关系数、变化斜率分析植被
变化情况及其与降水、温度变化的关系，期望得到陕北长城沿线风沙区近 05- 来的植被演变规律和主要驱动
因子。
!" 资料和方法
植被指数是卫星遥感中最具明确意义的指标，是基于植被叶绿素在 4& E1!F 处的强吸收，通过红外与近
红外波段的组合实现对植被信息状态的描述。在众多的植被指数中最为常用的是归一化植被指数 6789。它
能够敏感地反映出植被生长状况、生物物理化学性质及生态系统参数的变化，在一定程度上反映了像元所对
应区域的植被和土地覆盖类型的综合情况，因此，6789 常被直接或间接用于研究植被或土地覆盖变
化［/G 3 04］。H
本研究使用的资料是 /12/ 年 ? 月到 0445 年 /0 月 /GA 的最大值合成 6:;; $ ;8<== 6789，取自美国马
里兰大学 I9DDJ （I*)K-* 9.L’.")MN D)A’*,.+ -.A D-##,.+ J"OA,’B）2PF Q 2PF 分辨率的全球数据。这个数据
集中，对云、太阳高度角、仪器视场角、气溶胶的影响进行了校正。6789 的理论值在 R / 和到 / 之间，其中
R 4& /为水体，4 代表裸露地表，6789值从大于 4 到 / 代表了植被从疏到密的分布［0/，00］。气温资料和降水量
取自陕西气象档案馆。
图 0H 长城沿线风沙区年平均 6789曲线和线性趋势
H S,+& 0 H ;..O-* F’-. 6789 -.A ,"B "M’.A ,. "!’ %,.ATK)M. B-.AN
M’+,). -*).+ "!’ IM’-" U-** ,. .)M"!’M. J!--.V,
在 I9J支持下，将每个月两个 /GA 全球 6789 数据采用最大值合成法制作成月数据，逐月提取陕北长城
沿线风沙区及其所属各县（区）6789数据，统计月、春季（5 3 G 月份）、夏季（E 3 2 月份）、秋季（1 3 // 月份）、
冬季（/0 3 0 月份）、年平均值及不同级别 6789 的像元数。对于植被变化趋势的分析采用回归分析方法，即
把 6789值看作是一个时间的函数，对 6789与年份进行回归分析，得到一元线性回归方程，用来表示植被生
长和覆盖状况的变化趋势。如果方程的斜率为负值则认为植被覆盖度呈减小趋势，为正值则认为植被覆盖度
呈增加趋势。以相关系数通过的显著性检验水平 !判断植被覆盖减小或增加趋势明显程度，即 !*4& / 趋势
明显，!*4& 4G 趋势显著［1，/5］。长城沿线风沙区的年、
季温度、降水量，为该区域 E 县（区）年、季温度、降水量
的平均值。植被变化与气候因子之间关系的研究采用
相关分析方法，通过计算 6789 与降水、气温之间的相
关系数来表征植被与气候因子的相关性，判断 6789 与
气候因子的相关性也采用显著性检验方法。
#" 结果与分析
#& !" 空间平均 6789的变化特征
#& !& !" 年平均 6789的变化趋势
图 0 给出了长城沿线风沙区 /120 年到 0445 年的
年平均 6789 变化曲线和由线性回归得到的趋势线。
从年际变化来看，长城沿线风沙区的年平均 6789 在波
动起伏中增大，/120 年最差，为 4& /G?G；0440 年最好，
00E> H 生H 态H 学H 报H H H 0? 卷H
!""#：$ $ %%%& ’()*)+,(-& (.
为 /& 010/；22-平均 3456为 /& 0788。从线性趋势（图 0、表 0、表 2）来看，其线性回归方程的相关系数通过了
/& /0 显著水平检验，方程斜率的表明平均每年 3456 增加 7& 17 9 0/ :;，22- 间年平均 3456 上升了 0/& <2=。
趋势线滤去一些短期气候波动对 3456的影响，代表了 3456 随生态环境、长期气候变化及其综合效应的影
响，可以看出 01>2 ? 2//@ 年 3456总体上处于上升趋势。
表 !" 年、季平均 #$%&与年份的相关系数
’()*+ !" ,-..+*(/0-1 2-+33020+1/4 )+/5++1 #$%& (16 7+(.4
时间 A,B’
定边
4,.+C,-.
靖边
D,.+C,-.
横山
E’.+F!-.
榆阳
GH*,.
神木
I!’.BH
府谷
JH+H
长城沿线风沙区
A!’ %,.KLC)M.
F-.KN M’+,). -*).+
"!’ OM’-" P-** ,.
3)M"!’M. I!--.Q,
春 I#M,.+ /& 2;12 /& ;@>8!! : /& //2> /& 2><< /& @2@8 /& 82夏 IHBB’M /& 2210 /& ;70>!! /& 2>@/ /& <>@0@!
秋 RH"HB. /& ;2<@0!!! /& @87@ /& 8>/冬 P,."’M /& 2;<0 /& /<;@ : /& 0880 : /& /;70 : /& @0@8
年 R..H-* /& ;2;1!! /& <;8/!!! /& @@72 /& <>2;!!!! /& @>0;! /& @<>/! /& ;1;7!!!
S S !!!!! T /" //0，!!!! T /" /0，!!! T /" /8，!! T /" 0
表 8" 年、季节平均 #$%&与年份的回归方程斜率
’()*+ 8" 9*-:+ -3 *01+(. .+;.+440-1 +<=(/0-14 -3 #$%& /.+16
时间 A,B’
定边
4,.+C,-.
靖边
D,.+C,-.
横山
E’.+F!-.
榆阳
GH*,.
神木
I!’.BH
府谷
JH+H
长城沿线风沙区
A!’ %,.KLC)M.
F-.KN M’+,). -*).+
"!’ OM’-" P-** ,.
3)M"!’M. I!--.Q,
春 I#M,.+ /& ///281 /& ///;7@ /& /////2; /& ///220 /& ///2>> /& ///8>@ /& ///@/;
夏 IHBB’M /& //0/17 /& //01/; /& //022/ /& //22>7 /& //08@; /& //0700 /& //0<2<
秋 RH"HB. /& //00>/ /& //07;/ /& //000< /& //0;/< /& //088/ /& //0;>; /& //0;0@
冬 P,."’M /& ///;80 /& ///02/ : /& ///21/ : /& ///0// : /& ///70/ : /& ///;// : /& ///08/
年 R..H-* /& ///7;7 /& //0/81 /& ///80/ /& ///18@ /& ///<77 /& ///>;8 /& ///717
22-上升幅度（=）
6.(M’B’." ,. 22-
0/& 7< 0@& 1> <& >2 02& 72 >& 81 00& /< 0/& <2
图 @S 长城沿线风沙区季平均 3456年际变化曲线
S J,+& @ S I’-F).-* B’-. 3456 (!-.+’ ,. "!’ %,.KLC)M. F-.KN M’+,).
-*).+ "!’ OM’-" P-** ,. .)M"!’M. I!--.Q,
从各县（区）来看，年平均 3456变化趋势回归方程
的相关系数都为正值（表 0），方程的斜率也都为正值
（表 2），说明各个县（区）年平均 3456 变化趋势都是增
加。榆阳、靖边线性增大趋势尤其显著，它们的相关系
数是 /& <>2; 和 /& <;8/，分别通过了 /& //0 和 /& /0 显著
性水平检验。从各县（区）年回归方程的斜率可以看
出，靖边的年变化增量最大，为 0& /< 9 0/ :@，22-间年平
均 3456上升了 0@& 1>=；其余依次为榆阳、府谷、定边
和神木；横山最小，为 8& 0/ 9 0/ :;，22- 间年平均 3456
上升了 <& >2=，增加不显著。这与张莉等［0;］研究结论
有所不同，横山县植被覆盖没有出现明显的退化，这主
要是由于研究时段不同而致。
8& !& 8" 季平均 3456的变化趋势
图 @ 给出了长城沿线风沙区 01>2 ? 2//@ 年的各季
平均 3456变化曲线。从各季节平均 3456的年际变化
@2<;S 00期 S S S 李登科S 等：陕北长城沿线风沙区植被指数变化及其与气候的关系 S
!""#：$ $ %%%& ’()*)+,(-& (.
来看，长城沿线风沙区夏季的 /012平均值最高、波动起伏最大，其次是秋季。春季的 /012年际波动最小，冬
季虽然比春季波动大，但这主要是由积雪引起的。
表 3 中给出了长城沿线风沙区不同县（区）季平均 /012与年份相关系数。从整个陕北长城沿线风沙区
来看，秋季的相关系数最高为 4& 5367，通过了 4& 43 显著水平检验；春季和夏季的相关系数相近，通过了 4& 3 显
著水平检验；冬季相关系数最低，为负相关。说明陕北长城沿线风沙区秋季植被覆盖显著增加，春季、夏季明
显增加；冬季由于降雪增加而 /012 呈现减小趋势。从 /012 各季节变化趋势地域分布来看，冬季 /012 定
边、靖边即长城沿线风沙区西南部有所增加，横山、榆阳、神木、府谷即中部和东北部有所降低，但增加和降低
均不明显；春季 /012只有靖边、府谷显著增加，定边、榆阳、神木、府谷有所增加，横山有所降低；夏季定边、榆
阳显著增加，其余各县有所增加；秋季除横山有所增加外，其它县（区）均显著增加。
图 89 长城沿线风沙区月平均 /012年变化曲线
9 :,+& 89 66;<’-= ->’=-+’? @)."!*< /012 ,. "!’ %,.?;A)=. B-.?< =’+,).
-*).+ "!’ C=’-" D-** ,. .)="!’=. E!--.F,
!& "& #$ /012的季节变化
图 8 给出了陕北长城沿线风沙区 66-/012 的月平
均曲线。可以看出，/012从 8 月份开始增大，G 月份达
到极大值为 4& 6G57，H 月份开始下降，3 月份达到最低
值 4& 36I4，6 J 8 月份略有升高。春季 /012 缓慢增加，
秋季 /012降低速度比较快。
计算季节 /012与年 /012相关结果表明，春季、冬
季的 /012与年的 /012 是弱的正相关，而夏季、秋季
/012与年的 /012 相关系数分别达到了 4& G388 和
4K L574，显著性水平为 4& 443，说明夏季、秋季植被好，
全年的植被也会好的可信度达到了 HH& HM，即夏秋两
季的植被状况基本决定了陕北长城沿线风沙区全年的
植被分布状况。
!& "& %$ 不同等级 /012的面积变化
9 图 59 长城沿线风沙区 G 月份不同等级 /012 像元占总像元的百
分比
:,+& 59 N’=(’." )O #,F’*B %,"! ?,OO’=’." /012 *’>’* ,. PQ+QB" ,. "!’
%,.?;A)=. B-.?< =’+,). -*).+ "!’ C=’-" D-** ,. .)="!’=. E!--.F,
图 5 给出了从 3HG3 年到 644I 年每年 G 月长城沿
线风沙区不同级别的 /012 像元占整个区域总像元的
百分比。可以看出，/012在 4& 3 J 4& I 级别上像元所占
总像元的百分比是波动减少的，从 3HG3 J 3HG5 年的 5-
平均 L6& GM下降到 3HHH J 644I 年 5- 平均 87& GM；相
反，/012在 4& I J 4& 7 级别上所占总像元的百分比是波
动增加的，从 3HG3 J 3HG5 年的 5- 平均 6L& 6M上升到
3HHH J 644I 年 5- 平均 5I& 6M。说明长城沿线风沙区
低覆盖率植被面积在减少，高覆盖率植被面积在增加，
植被状况明显转好。
!& !$ /012的变化与气候变化的关系
!& !& "$ /012的变化与温度变化的关系
计算长城沿线风沙区年、季平均 /012 与温度的相
关系数可以得出年平均 /012 的年际变化与年平均温
度的年际变化相关不明显，各季节的平均 /012与相应季节的温度相关也不明显。由图 7 可知，年平均 /012
与温度年际波动没有一致性，尽管它们的线性变化趋势都是上升的，但由于各自上升的速度大小差异使得两
者之间没有明显的相关性。各季节 /012与温度的关系分析也获得类似的结果。例外的是春季 /012与夏季
温度显著相关，这种关系没有明确的生物学意义，可能是有偶然因素造成的。
8678 9 生9 态9 学9 报9 9 9 6L 卷9
!""#：$ $ %%%& ’()*)+,(-& (.
图 /0 长城沿线风沙区年 1234与温度的关系
0 5,+& / 0 6’*-",).7!,# 8’"%’’. -..9-* 1234 -.: "’;#’<-"9<’ ,. "!’
%,.:=8)<. 7-.:> <’+,). -*).+ "!’ ?<’-" @-** ,. .)<"!’<. A!--.B,
由表 C 可以得出，D 个数县年平均 1234 与年平均
温度的相关性不显著；只有榆阳、靖边的相关系数均分
别通过了 E& EF 和 E& F 显著水平的检验，说明年平均温
度的高低是影响榆阳植被覆盖变化的重要因子，年平均
温度对靖边植被覆盖有明显影响；各县季平均 1234 与
相应的平均温度的 GD 个相关系数中，GG 个未通过 !*
E& F 的显著水平检验，只有靖边春季、榆阳夏季的相关
系数通过了 E& F 显著水平的检验，说明绝大多数县
（区）各季节温度对植被覆盖变化的作用不显著。
!& !& !" 1234的变化与降水量变化的关系
表 D 给出了各县（区）年、季平均 1234与降水量的
相关系数。从表 D 可以看出，西南部的定边、靖边年平
均 1234与年降水量显著相关。冬季平均 1234与降水
量相关系数均为负值，靖边、横山分别通过了 E& EF 和
E& F 显著水平的检验，说明冬季降水多，即积雪多 1234
就低。春季靖边、神木的相关系数达到 E& EF 显著水平的检验，夏季靖边的相关系数通过 E& F 显著水平检验。
表 #" 各县（区）$%&’与温度的相关系数
()*+, #" -.//,+)01.2 3.,44131,205 *,06,,2 $%&’ )27 0,89,/)0:/, 12 7144,/,20 3.:201,5
时间 H,;’
定边温度
H’;#’<-"9<’
)I 2,.+8,-.
靖边温度
H’;#’<-"9<’
)I J,.+8,-.
横山温度
H’;#’<-"9<’
)I K’.+7!-.
榆阳温度
H’;#’<-"9<’
)I L9*,.
神木温度
H’;#’<-"9<’
)I A!’.;9
府谷温度
H’;#’<-"9<’
)I 59+9
春 A#<,.+ E& GF/M E& CNMD! E& E/OO E& FDNF E& FNPP E& GCMM
夏 A9;;’< Q E& EEGG E& FMGF Q E& FG/D E& CNFD! E& EGFC Q E& ENM/
秋 R9"9;. Q E& EE// E& FGD/ Q E& ENEE E& EPFM E& G/FE E& EMFD
冬 @,."’< E& GONO E& EOMN Q E& FEE/ Q E& FGDC Q E& GONO Q E& G/CO
年 R..9-* E& GEOP E& CONO! E& EP/E E& MEEC!!! E& GCFC E& FEE/
0 0 !!! ! S E" EF，! ! S E" F
表 ;" 各县（区）$%&’与降水量的相关系数
()*+, ;" -.//,+)01.2 3.,44131,205 *,06,,2 $%&’ )27 9/,31910)01.2 12 7144,/,20 3.:201,5
时间 H,;’
定边
T<’(,#,"-",).
)I 2,.+8,-.
靖边
T<’(,#,"-",).
)I J,.+8,-.
横山
T<’(,#,"-",).
)I K’.+7!-.
榆阳
T<’(,#,"-",).
)I L9*,.
神木
T<’(,#,"-",).
)I A!’.;9
府谷
T<’(,#,"-",).
)I 59+9
春 A#<,.+ E& FDOG E& DCFM!! E& GOCN E& GF/E E& DMOP!! E& CEFP
夏 A9;;’< E& GDPD E& CPM/! E& GGFE E& EN/C E& FF/E E& EGON
秋 R9"9;. E& EDFG Q E& FGNE Q E& GFGC Q E& GFOO Q E& EFCO Q E& FMDN
冬 @,."’< Q E& CDCP Q E& MFDD!!! Q E& C/NP! Q E& GDOP Q E& GFCG Q E& ENDG
年 R..9-* E& DNEN!! E& DDPE!! E& GMOP E& GPCP E& CDGN E& GCED
0 0 !!! ! S E" EF，!! ! S E" EM，! ! S E" F
表 <" 陕北长城沿线风沙区 $%&’与降水量的相关系数
()*+, <" -.//,+)01.2 3.,44131,205 *,06,,2 $%&’ )27 9/,31910)01.2 12 0=, 6127 > *./2 5)27? /,@1.2 )+.2@ 0=, A/,)0 B)++ 12 2./0=,/2 C=))2D1
时间 H,;’
春季降水量
A#<,.+ #<’(,#,"-",).
夏季降水量
A9;;’< #<’(,#,"-",).
秋季降水量
R9"9;. #<’(,#,"-",).
冬季降水量
@,."’< #<’(,#,"-",).
年降水量
R..9-* #<’(,#,"-",).
春 A#<,.+ E& CNMG E& EGFC Q E& EOPF E& DEDF! E& GE/C
夏 A9;;’< E& DEFC! E& FNFC Q E& GOFC E& GN/O E& GND/
秋 R9"9;. Q E& FCEN E& MFPN!!! Q E& FCEN E& GMOP E& CP/G!
冬 @,."’< E& EGNE E& GNPF Q E& GCFO Q E& CEMF E& F/MG
年 R..9-* E& G/FN E& DFOM Q E& CFPC E& GMCP E& DGCG!!
0 0 !!! ! S E" EF，!! ! S E" EM，! ! S E" F
MG/D0 FF期 0 0 0 李登科0 等：陕北长城沿线风沙区植被指数变化及其与气候的关系 0
!""#：$ $ %%%& ’()*)+,(-& (.
图 /0 长城沿线风沙区年 1234与年降水量的关系
0 5,+& / 0 6’*-",).7!,# 8’"%’’. -..9-* 1234 -.: #;’(,#,"-",). ,. "!’
%,.:<8);. 7-.:= ;’+,). -*).+ "!’ >;’-" ?-** ,. .);"!’;. @!--.A,
长城沿线风沙区年平均 1234 与降水量的显著相
关，它们之间的相关系数为 B& CDED，通过 B& BF 显著水
平检验（表 F）。图 / 是长城沿线风沙区年 1234 与降
水量的关系图。从图 / 可以看出，1234 的年际波动与
降水量的年际波动一致。年降水量多，该年的 1234 就
大，反之，1234就小。但 1234的年际变化趋势与降水
量的年际变化趋势有明显的不同。统计得到年降水量
与年份之间的相关系数仅为 B& B/CG，未通过 B& G 显著
水平检验，说明近 DD-来长城沿线风沙区的年降水量没
有明显增加。而从年平均 1234 线性增加趋势显著，说
明降水量是引起 1234年际波动的主要因子，非气候因
素是年平均 1234线性增加的主要原因。
拟合长城沿线风沙区年平均 1234 与年份线性回
归方程，再将残差与年降水量进行线性回归，所得的回
归方程的相关系数为 B& CCCH，通过了 B& BF 显著水平检
验。进一步证明了降水量是引起 1234 年际波动的主要因子，非气候因素是年平均 1234 线性增加的主要
原因。
长城沿线风沙区季平均 1234与季降水量年际变化存在滞后性。从表 F 可以看出，夏季 1234 与春季降
水量显著相关，它们之间的相关系数为 B& CBGE，通过 B& G 显著水平检验；秋季 1234 与夏季降水量显著相关，
它们之间的相关系数为 B& FGI/，通过 B& BG 显著水平检验。冬季 1234与春季降水量明显相关，它们之间的相
关系数为 B& CBCG，通过 B& G 显著水平检验。从图 H 可以看出，夏季 1234 与春季降水量年际变化波动具有一
致的规律，不同的是夏季 1234年际变化趋势是向上的，而春季降水量年际变化趋势方向相反。在图 I 中，秋
季 1234与夏季降水量年际变化波动也具有一致性，秋季 1234和夏季降水量年际变化趋势都是向上的，但秋
季 1234年际变化显著，而夏季降水量年际变化趋势不明显（夏季降水量与年份的相关系数仅为 B& BCBE，未通
图 H0 长城沿线风沙区夏季 1234与春季降水量的关系
0 5,+& H0 6’*-",).7!,# 8’"%’’. 79JJ’; 1234 -.: 7#;,.+ #;’(,#,"-",). ,.
"!’ %,.:<8);. 7-.:= ;’+,). -*).+ "!’ >;’-" ?-** ,. .);"!’;. @!--.A,
图 I0 长城沿线风沙区秋季 1234与夏季降水量的关系
0 5,+& I0 6’*-",).7!,# 8’"%’’. -9"9J. 1234 -.: 79JJ’; #;’(,#,"-",).
,. "!’ %,.:<8);. 7-.:= ;’+,). -*).+ "!’ >;’-" ?-** ,. .);"!’;. @!--.A,
KDKC 0 生0 态0 学0 报0 0 0 D/ 卷0
!""#：$ $ %%%& ’()*)+,(-& (.
过 /& 0 显著水平检验）；春季 1234与冬季降水量之间的关系与此类似。说明降水量对 1234 影响具有滞后
性，某季降水量是对下一季 1234具有正效应；季降水量是引起季 1234年际波动的主要因子，而非气候因素
是 1234线性增加的主要原因。
杨述河等［5］分析基于 67 影像解译土地利用数据，认为 089: ; 好转和恶化的状况，但总体上生态环境质量呈现改善趋向。从土地利用变化对生态环境影响的空间规律看，
北部风沙滩地区要比南部丘陵沟壑区生态环境改善的程度高。导致区域生态环境质量改善的主要驱动因子
包括退耕、还林、治沙种草与防护林建设（造林）、未利用地开垦等，其中治沙种草与造林对榆林市生态环境改
善的贡献率居首位。另外一些也研究发现，程度却有所好转，主要原因是实施了一系列积极的人为措施和政策，说明人类活动对现代土地沙漠化过程起
着重要的作用［0=，<>，<=］。榆林市是我国荒漠化治理比较有成效的地区，多年来开展了坚持不懈的沙漠化防治
和整治工作，具体措施有：限制放牧、鼓励圈养、减少山羊数量、增加优良品种；退耕还林、还草；建设三北防护
林；飞播造林、植草，节水灌溉等，这些措施对 影响［<:，<5］。
!" 结论
通过对长城沿线风沙区 0890 年到 年卫星遥感的 1234和地面气象观测数据的分析，可以得出如下
主要结论：
（0）长城沿线风沙区的植被状况 <<-来尽管有波动起伏，但是整体在持续转好。089< ; 年 1234 总
体上处于上升趋势，平均每年 1234增加 ?& 8? @ 0/ A=，<<-间年平均 1234上升了 0/& 5
长城沿线风沙区夏季的 1234平均值最高、波动起伏最大，其次是秋季；春、夏、秋三季的 1234 具有明显
的上升趋势，季 1234年增长率夏季最大，秋季次之；夏、秋季 1234 与年 1234 具有很高的相关性，这两个季
节的植被状况基本决定了全年的植被分布状况。1234 年变化曲线为单峰型。从 = 月份开始增大，9 月份达
到极大值为 /& <9:5，8 月份开始下降，0 月达到最低值 /& 0<>/。春季 1234 缓慢增加，秋季 1234 降低速度比
较快。
（<）年平均 1234 的年际变化与年平均温度的年际变化相关不明显，各季节的平均 1234 与相应季节的
温度相关也不明显。近 <<-来长城沿线风沙区的年降水量没有明显增加，而年平均 1234 线性增加趋势显
著，降水量是引起 1234年际波动的主要因子，非气候因素是年平均 1234 线性增加的主要原因。降水量与
1234存在着明显的年相关和隔季相关。年降水量与年 1234的相关，冬季降水量与春季 1234的相关，春季
降水量与夏季 1234的相关，夏季降水量与秋季 1234 的相关性都非常高。季降水量对 1234 影响具有滞后
性，某季降水量是对下一季 1234具有正效应；季降水量是引起季 1234年际波动的主要因子，而非气候因素
是 1234线性增加的主要原因。
（>）非气候因素中生态保护和环境建设等人为措施，如植树造林、草原围栏封育等是导致植被显著增加
的重要原因。
#$%$&$’($)：
［0 ］C DE D F，G!-.+ H F，IE J H，!" #$% K-EL’ )M 2’L’N",M,(-",). -.O FN’P’.",P’ 7’-LEN’L -+-,.L" ,. QE*,.& H)EN.-* )M 1)N"!%’L" I)N’L"NR S.,P’NL,"R，<0（<）：05 08&
［< ］C TE-. G H，U-.+ K H，U-.+ V W，!" #$% X’Y)"’ D’.L,.+ D"EOR ). "!’ V-.O 2’L’N",M,(-",). 6N’.O ,. QE*,. FN’M’("EN’ ,. 1)N"! D!-.Z,& K!,.’L’
H)EN.-* )M [()*)+R，0885，0:（:）：<> <5&
［> ］C U-.+ \ I，X’. G Q，6-. ] V& 3-*E’L )M ’()LRL"’Y L’NP,(’L -.O "!’,N (!-.+’ ,. "!’ N’+,).L -*).+ "!’ JN’-" %-** ,. .)N"! D!--.Z, FN)P,.(’& ^N,O
V-.O J’)+N-#!R， <=5&
［= ］C J-) \ W，U-.+ Q 7，U-.+ H W，!" #$% ^.-*RL,L ). 2’L’N",M,(-",). 2R.-Y,(L _-L’O ). X’Y)"’ D’.L,.+ -.O J4D ,. G).’ -*).+ "!’ JN’-" U-** ,.
?<5=C 00期 C C C 李登科C 等：陕北长城沿线风沙区植被指数变化及其与气候的关系 C
!""#：$ $ %%%& ’()*)+,(-& (.
/)0"!’0. 1!--.2, 30)4,.(’& 5)60.-* )7 8’9’0" :’9’-0(!，;<<=，;=（>）：?@ ??&
［= ］A B,6 C 1，5-D E& F0’.G H.-*D9,9 )7 B-.G 8’+0-G-",). ,. "!’ I).’ -*).+ "!’ E0’-" J-** ,. /)0"!’0. 1!--.2,& H("- E’)+0-#!,(- 1,.,(-，;<<;，=K（L）：
LL@ L=<&
［? ］A C-.+ 1 M，C-. M B，E6) B C& F!’ B-.G N9’ O!-.+’ -.G P"9 Q()R’.4,0).S’."-* Q77’("9 ,. F0-.9,",).-* H+0)R#-9")0-* :’+,). H O-9’ 1"6GD )7
C6*,. O,"D ,. /)0"!’0. 1!--.2, 30)4,.(’& 30)+0’99 P. E’)+0-#!D，;<
1!’.S6 O)6."D，H0,G B-.G E’)+0-#!D，;<<@，;?（L）：L<; L
5)60.-* )7 E*-(,)*)+D -.G E’)(0D)*)+D，;<<@，;=（;）：;@; ;@?&
［T ］A 3,-) 1 B，V-.+ 5 C& 8D.-S,( Z’+’"-",). O)4’0 O!-.+’ )4’0 "!’ B-9" >X C’-09 ,. O!,.-& W6-"’0.-0D 1(,’.(’9 ，;<<>，;>（L）：;TL @<;&
［><］A B, I，C-. V B，V-. U F& F!’ Z-0,-[,*,"D )7 /8ZP )4’0 /)0"!%’9" O!,.- -.G P"9 :’*-",). ") F’S#’0-"60’ -.G 30’(,#,"-",).& 5)60.-* )7 :’S)"’
1’.9,.+，;<<=，T（@）：@
［>>］A V-.+ 5 C，3,-) 1 B，M’ 5 1& Q.!-.(’S’." )7 Z’+’"-",). H(",4,"D ,. O!,.- ,. :’(’." ;< C’-09& 1(,’.(’ ,. O!,.-，1’0& O，;<<@，@@（?）：==L =?=&
［>;］A B, I O，Y- Y E，I!-.+ V& F!’ 8D.-S,( H.-*D9,9 )7 Z’+’"-",). 3-""’0. ,. "!’ /)0"!%’9" )7 O!,.-& :’S)"’ 1’.9,.+ F’(!.)*)+D -.G H##*,(-",).，
;<（L）：@@; @@K&
［>@］A B, C O，E).+ 3，B,6 O U& Z’+’"-",). O)4’0 O!-.+’9 -.G O)00’*-",). %,"! O*,S-",( V-(")09 ,. /)0"!’0. O!,.- 860,.+ >TX; >TTT& :’9)60(’9
1(,’.(’，;<
［>L］A I!-.+ B，J-.+ V C，C6’ B 3，!" #$% H 8D.-S,( 1"6GD )7 B-.G 8’9’0",7,(-",). ,. 8’9’0"RB)’99 F0-.9,",).-* I).’9 \-9’G ). :1 -.G EP1：H O-9’
1"6GD )7 C6*,. H0’-& H("- E’)9,(,’.",- 1,.,(-，;<
［>=］A F6(]’0 O 5 ，Z-.#0-’" O P ，1!-0S-. Y 5 ，!" #$& 1-"’**,"’ 0’S)"’ 9’.9,.+ )7 ")"-* !’0[-(’)69 [,)S-99 #0)G6(",). ,. "!’ 1’.’+-*’9’ 1-!’*：>TX<
>TXL& :’S)"’ 1’.9,.+ )7 Q.4,0).S’."，>TX=，>K：;@@ ;LT&
［>?］A E0-S’0 J，!" #$& O)S#-0,.+ E*)[-* Y)G’*9 )7 F’00’9"0,-* /’" 30,S-0D 30)G6(",4,"D （/33）：^ 4’04,’% -.G _’D :’96*"9& E*)[-* O!-.+’ \,)*)+D，
>TTT，=：> >=&
［>K］A O,!*-0 5，BD M，U,-) W& B-.G O)4’0 O*-99,7,(-",). %,"! HZM:: Y6*",(!-..’* O)S#)9,"’9 ,. /)0"!’0. Q.4,0).S’."9& :’S)"’ 1’.9,.+ )7
Q.4,0).S’."，>TT?，=X（>）：@? =>&
［>X］A 1’**’09 3 5，Y’’9). \ J，M-** V E，!" #$& :’S)"’ 1’.9,.+ )7 B-.G 1607-(’ 7)0 1"6G,’9 )7 E*)[-* O!-.+’：Y)G’*9R-*+)0,"!S9R’2#’0,S’."9% :’S)"’
1’.9,.+ )7 Q.4,0).S’."，>TT=，=>（>）：@ ;?&
［>T］A F6(]’0 O 5，F)%.9!’.G 5 : E& H70,(-. B-.GR()4’0 O*-99,7,(-",). N9,.+ 1-"’**,"’ 8-"-& 1(,’.(’，>TX=，;;K（L?X=）：@?T @K=&
［;<］A F)%.9!’.G 5 : E& E*)[-* 8-"- 1’"9 7)0 B-.G H##*,(-",).9 70)S "!’ HG4-.(’G Z’0D M,+! :’9)*6",). :-G,)S’"’0：-. P."0)G6(",).& P."’0.-",).-*
5)60.-* )7 :’S)"’ 1’.9,.+，>TTL，>=（>K）：@@>T @@@;&
［;>］A O!’. U E，B, 5 3，B, I 5，!" #$& Z’+’"-",). (!-.+’ ,. C-.(!, )7 /,.+2,- -.G ,"9 0’*-",).9!,# %,"! (*,S-"’ (!-.+’ ,. 0’(’." D’-09& H("- Q()*)+,(-
1,.,(-，;<=>? >=;;&
［;;］A B6 B，B, U，O!’.+ E 8& H.-*D9,9 ). "!’ 1’-9).-* 3!’.)*)+,(-* O!-0-("’0,9",(9 )7 "!’ M’,!’ :,4’0 \-9,. %,"! HZM:: /8ZP 8-"- 1’"& 5)60.-* )7
8’9’0" :’9’-0(!，;<<;，;;（;）：>XK >T>&
［;@］A /,’ M E，J-.+ Y，16. M，!" #$& 8D.-S,( H.-*D9,9 ). O!-.+’9 )7 8’9’0",7,(-",). ,. 8’9’0"R*)’99 \)6.G-0D \’*" 860,.+ "!’ B-9" >< C’-09 H O-9’
1"6GD ). C6*,. H0’- ,. /)0"! 1!--.2, 30)4,.(’& /)0"!%’9"’0. E’)*)+D ，;<<=，@X（;）：X? T@&
［;L］A E-) M 5，5,-.+ W E，M6) U \& :’S)"’ 1’.9,.+ :’9’-0(! )7 1-.GD 8’9’0",7,(-",). -*).+ "!’ E0’-" J-** ,. /)0"!’0. 1!--.2,& 5)60.-* )7 /-"60-*
:’9)60(’9，;<<=，;<（@）：LK> LK=&
［;=］A J-.+ U& 8,9(699,). ). 30)[*’S9 )7 Q()*)+,(-* Q.4,0).S’." O).9"06(",). ,. C6*,.& :’9’-0(! )7 1),* -.G J-"’0 O).9’04-",).，;<<@，><（L）：;KX
;KT&
［;?］A B6 I V，U,’ C 1，J-.+ 5 5& H.-*D9’9 ). O!-0-("’0,9",(9 -.G FD#,(-* Y)G’9 )7 Q()*)+,(-* -.G Q.4,0).S’."-* O).9"06(",). ,. /)0"! 1!--.2,&
:’9’-0(! )7 1),* -.G J-"’0 O).9’04-",).，;<<@，><（L）：=> =@&
［;K］A I!-) U \，5,- C Y& 8’9’0",7,(-",). -.G 30’4’.",4’ Y’-960’9 ,. C6*,.& 30)"’(",). V)0’9" 1(,’.(’ -.G F’(!.)*)+D，;<<=，（?）：L? LK&
参考文献：
［> ］A 苏世平，张继平，付广军，等&榆林沙区荒漠化成因及防治对策&西北林学院学报，;<（;）：>? ‘ >T&
［; ］A 全志杰，王成吉，王立宏，等&陕北榆林地区土地沙漠化动态遥感研究&生态学杂志，>TT?，>=（=）：;@ ‘ ;?&
［@ ］A 王晓峰，任志远，谭克龙&陕北长城沿线地区生态系统服务价值变化研究&干旱区地理，;<［L ］A 高小红，王一谋，王建华，等&陕北长城沿线地区 >TX? ‘ ;<<< 年沙漠化动态变化分析&中国沙漠，;<<=，;=（>）：?@ ‘ ??&
X;?L A 生A 态A 学A 报A A A ;K 卷A
!""#：$ $ %%%& ’()*)+,(-& (.
［/ ］0 刘彦随，1-2 3-)&陕北长城沿线地区土地退化态势分析&地理学报，4554，/6（7）：778 9 7/5&
［: ］0 杨述河，闫海利，郭丽英&北方农牧交错带土地利用变化及其生态环境效应———以陕北榆林市为例&地理科学进展，4557，48（:）：7; 9 //&
［6 ］0 王晓峰，任志远，黄青&农牧交错区县域土地利用变化及驱动力分析———以陕北神木县为例&干旱区地理，4558，4:（7）：754 9 756&
［< ］0 马明国，董立新，王雪梅&过去 4=-中国西北植被覆盖动态监测与模拟&冰川冻土，4558，4/（4）：484 9 48:&
［; ］0 朴世龙，方精云&最近 =< 年来中国植被覆盖的动态变化&第四纪研究，455=，4=（7）：4;7 9 854&
［=5］0 李震，阎福礼，范湘涛&中国西北地区 >?@A变化及其与温度和降水的关系&遥感学报，455/，;（8）：85< 9 8=8&
［==］0 方精云，朴世龙，贺金生，等&近 45 年来中国植被活动在增强&中国科学（B辑），4558，88（:）：//7 9 /:/&
［=4］0 李珍存，马明国，张峰，等& =;<4 9 4558 年中国西北地区植被动态变化格局分析&遥感技术应用，455:，4=（7）：884 9 886&
［=8］0 李月臣，宫鹏，刘春霞，等&北方 =8 省 =;<4 年 9 =;;; 年植被变化及其与气候因子的关系&资源科学，455:，4<（4）：=5; 9 ==6&
［=7］0 张莉，王飞跃，岳乐平，等&基于 CD和 3AD的沙漠E黄土过渡带土地沙漠化分布与变化研究———以陕北榆林地区为例& 西北地质，4557，4/
（=）：:8 9 ::&
［4=］0 陈晓光，李剑萍，李志军，等&宁夏盐池近年来植被与气候变化分析&生态学报，455:，4:（/）：=/=: 9 =/44&
［44］0 卢玲，李新，程国栋&利用 >FGG G@HCC植被指数数据集分析黑河流域季候特征&中国沙漠，4554，44（4）：=<6 9 =;=&
［48］0 聂浩刚，王岷，孙虎，等&近 =5 年来沙漠E黄土边界带荒漠化土地动态分析———以陕北榆林地区为例&西北地质，455/，8<（4）：<: 9 ;8&
［47］0 高会军，姜琦刚，霍晓斌&陕北长城沿线沙质荒漠化遥感研究&自然资源学报，455/，45（8）：76= 9 76/&
［4/］0 王晓&陕西榆林生态环境建设有关问题探讨&水土保持研究，4558，=5（7）：46< 9 46;&
［4:］0 卢宗凡，谢永生，王继军&陕北生态环境建设特点与典型模式分析&水土保持研究，4558，=5（7）：/= 9 /8&
［46］0 赵晓彬，贾艳梅&榆林沙区沙漠化及其防治措施&防护林科技，455/，（:）：7: 9 76&
;4:70 ==期 0 0 0 李登科0 等：陕北长城沿线风沙区植被指数变化及其与气候的关系 0