摘 要 :地表水分是监测土地退化的一个重要指标,是气候、水文、生态、农业等领域的主要参数。选择陕北黄土高原地区作为研究区域,首先采用基于植被覆盖特征的植被状态指数(Vegetation Condition Index,简称VCI)和基于地表温度特征(Land Surface Temperature,简称Ts)的温度状态指数(Temperature Condition Index, 简称TCI)分别评价了区域地表水分状况的季节变化。在此基础上分析了植被指数与地表温度特征线性关系的季节变化规律,计算了基于两者经验关系的地表干燥度指数(Temperature-Vegetation Dryness Index,TVDI)。该指数对Ts/NDVI特征空间的生态特征的解释,对土壤和作物的水分含量具有综合的指示意义。文中利用该指数综合评价了研究区域地表水分状况的时空分布差异,进一步对VCI、TCI与TVDI相关关系的季节变化进行比较分析,并结合气候因子进行了相关验证,从而对不同指数的应用范围做出判定。研究结果表明,单独采用TCI或VCI表征地表水分会受到明显的季节影响,而TVDI能在不同季节综合体现植被覆盖和地表温度特征对地表水分的影响,从而能较好的反映区域地表水分状况的时空变化特征。各区域的TVDI值季节分布上皆为4~7月份高于10~翌年1月份,但各区TVDI值的季节变化则存在显著不同,而各流域内部TVDI值的空间变异性也存在季节差异,其中在10月份较为显著。
Abstract:Vegetation coverage and surface temperature are important parameters to describe characters of land covers, and can provide information on vegetation and moisture conditions at the surface. The present paper aims at demonstrating how Terra/MODIS data may be used to estimate spatial patterns of soil moisture, a key variable in distributed hydrological models. The basic approach is to interpret the so-called Ts/NDVI space in terms of surface soil moisture status. Three kinds of index, such as vegetation condition index (VCI), temperature condition index (TCI) and a simplified land surface dryness index (Temperature Vegetation Dryness Index, TVDI) based on an empirical parameterisation of the relationship between surface temperature and vegetation index are suggested. These indexes are related to soil moisture and, in comparison to existing interpretations of the Ts/NDVI space, are conceptually and computationally straightforward. These indexes are based on satellite-derived information only, and the potential for application is therefore large.
Using 36 Terra/MODIS images from 2003 of the semiarid North Shaanxi Loess Plateau, NDVI (MOD11) and Ts (MOD13) were investigated as indicators of vegetation abundance and land surface temperature. To examine soil moisture conditions at the surface, the spatial pattern and temporal evolution in VCI, TCI and TVDI was analysed with Ts and NDVI. The results were: (1) Land surface moisture represented by VCI or TCI separately could not eliminate the effects caused from seasonal change of Ts and NDVI. However, TVDI, derived from analysis of the temperature/vegetation index space and created by integrated analyses of these two kinds of data, helped to determine main principles of the temporal and spatial variation of surface soil moisture. (2) From temporal evolution of TVDI, it could be inferred that the trend in TVDI is high values in the dry season (spring or autumn) and low values in the rainy season (summer or winter). (3) Spatial evolution of TVDI indicates that surface soil moisture of each region has its own seasonal rhythm. Mean TVDI in the northern region increased rapidly from winter to spring and decreased from summer to autumn in south; in the middle and west, the increment of TVDI from winter to spring was almost equal to decrement from summer to autumn. (4) Similarly, spatial variation of surface soil moisture in watersheds was seasonal. Typically, in October the variation in TVDI was greatest in most watersheds, especially in north.
全 文 :第 !" 卷第 ## 期
!" 年 ## 月 生 态 学 报 %&’% (&)*)+,&% -,.,&% /012 !",.02 ## .032,!
"
!""#:
%%%& ’()*)+,(-& (.
基金项目:国家科技支撑计划资助项目(!
45
#);国家自然科学基金重点研究资助项目
(7
489!:)
收稿日期:!
";$!;$:;修订日期:!
";#
;!9 作者简介:李正国(#<: =),男,江西上饶人,博士生,主要从事景观生态与土地利用的研究2 (;>?@1:1AB2 #!8C !482 DEF
!通讯作者 &0GGEHI0DJ@DB ?KFL0G2 (;>?@1:M1N?DBC KGO?D2 IPK2 EJK2 QD
致谢:感谢国家农业部资源遥感与数字农业重点开放实验室提供 ’(RR% S T)6,-数据产品
/)0.1-",). ,"’2:’LE IG0UEQF N?H V@D?DQ@?11M HKII0GFEJ OM .?F@0D?1 -Q@EDQE ?DJ ’EQLD010BM -KII0GF WG0BG?> 0V &L@D?(.02 !
45
#);.?F@0D?1 .?FKG?1 -Q@EDQE YEM X0KDJ?F@0D 0V &L@D? (.02 7
489!:)
3’(’,4’1 1-"’:!
";$!;$:;5((’#"’1 1-"’:!
";#
;!9 6,)+7-#!8:*, ZLEDB;+K0,IL2 62 Q?DJ@J?FE2 >?@D1M EDB?BEJ @D 1?DJHQ?IE EQ010BM ?DJ 1?DJ KHE2 (;>?@1:1AB2 #!8C !482 DEF 基于植被温度特征的黄土高原地表水分季节变化
李正国#,!,王仰麟!,!,吴健生!,张小飞8
(#2农业部资源遥感与数字农业重点开放实验室,北京[ #
$:#;!2 北京大学 城市与环境学院,北京[ #
:"#;
82中国科学院生态环境研究中心,北京[ #
$:9) 摘要:地表水分是监测土地退化的一个重要指标,是气候、水文、生态、农业等领域的主要参数。选择陕北黄土高原地区作为研 究区域,首先采用基于植被覆盖特征的植被状态指数(/EBEF?F@0D &0DJ@F@0D ,DJE\,简称 /&,)和基于地表温度特征(*?DJ -KGV?QE ’E>IEG?FKGE,简称 ’H)的温度状态指数(’E>IEG?FKGE &0DJ@F@0D ,DJE\,简称 ’&,)分别评价了区域地表水分状况的季节变化。在此 基础上分析了植被指数与地表温度特征线性关系的季节变化规律,计算了基于两者经验关系的地表干燥度指数(’E>IEG?FKGE; /EBEF?F@0D 6GMDEHH ,DJE\,’/6,)。该指数对 ’H S .6/,特征空间的生态特征的解释,对土壤和作物的水分含量具有综合的指示意 义。文中利用该指数综合评价了研究区域地表水分状况的时空分布差异,进一步对 /&,、’&, 与 ’/6, 相关关系的季节变化进 行比较分析,并结合气候因子进行了相关验证,从而对不同指数的应用范围做出判定。研究结果表明,单独采用 ’&, 或 /&, 表 征地表水分会受到明显的季节影响,而 ’/6,能在不同季节综合体现植被覆盖和地表温度特征对地表水分的影响,从而能较好 的反映区域地表水分状况的时空变化特征。各区域的 ’/6,值季节分布上皆为 7 = " 月份高于 #$ =翌年 # 月份,但各区 ’/6, 值的季节变化则存在显著不同,而各流域内部 ’/6,值的空间变异性也存在季节差异,其中在 #$ 月份较为显著。 关键词:’/6,;地表水分;黄土高原 文章编号:#
;<88(!
")##;7948;#8[ 中图分类号:]#7<[ 文献标识码:% 9."7-:-..0-* ;07<-(’ ;),* 2),;"07’ (!-.+’ =-;’1 ). 4’+’"-",). > "’2#’7-"07’ (!-7-("’7,;",(; ,. ?)’;; @*-"’-0 -7’- *, ZLEDB;+K0#,!,^ %.+ _?DB;*@D!,!,^‘ a@?D;-LEDB!,&b%.+ bH@?0;XE@8 # !"# $%&’(%)’(# ’* +",’-(.", +"/’)" 0"1,213 %14 5232)%6 73(2.-6)-(" ’* 8212,)(# ’* 73(2.-6)-(" (897),:"2;213 #
:#,<=21% ! <’66"3" ’* >(&%1 ? @1A2(’1/"1)%6 0.2"1.",,B"C213 >12A"(,2)#,:"2;213 #:"#,<=21% 8D +","%(.= <"1)"( *’( @.’E@1A2(’1/"1)%6 0.2"1.",,<=21"," 7.%4"/# ’* 0.2"1.",,BD 9D :’F !:"#,:"2;213 #$$:9,<=21%
!"#$ %"&’&()"$ *)+)"$,ABBC,AC(DD):EFGH I EFCF,
5=;"7-(":/EBEF?F@0D Q03EG?BE ?DJ HKGV?QE FE>IEG?FKGE ?GE @>I0GF?DF I?G?>EFEGH F0 JEHQG@OE QL?G?QFEGH 0V 1?DJ Q03EGH,?DJ
Q?D IG03@JE @DV0G>?F@0D 0D 3EBEF?F@0D ?DJ >0@HFKGE Q0DJ@F@0DH ?F FLE HKGV?QE2 ’LE IGEHEDF I?IEG ?@>H ?F JE>0DHFG?F@DB L0N
’EGG? S T)6,- J?F? >?M OE KHEJ F0 EHF@>?FE HI?F@?1 I?FFEGDH 0V H0@1 >0@HFKGE,? PEM 3?G@?O1E @D J@HFG@OKFEJ LMJG010B@Q?1
>0JE1H2 ’LE O?H@Q ?IIG0?QL @H F0 @DFEGIGEF FLE H0;Q?11EJ ’H S .6/, HI?QE @D FEG>H 0V HKGV?QE H0@1 >0@HFKGE HF?FKH2 ’LGEE P@DJH
0V @DJE\,HKQL ?H 3EBEF?F@0D Q0DJ@F@0D @DJE\ (/&,),FE>IEG?FKGE Q0DJ@F@0D @DJE\ (’&,)?DJ ? H@>I1@V@EJ 1?DJ HKGV?QE
!""#:$ $ %%%& ’()*)+,(-& (.
/01.’22 ,./’3 (4’5#’0-"60’ 7’+’"-",). 801.’22 9./’3,4789):-2’/ ). -. ’5#,0,(-* #-0-5’"’0,2-",). ); "!’ 0’*-",).2!,#
:’"%’’. 260;-(’ "’5#’0-"60’ -./ <’+’"-",). ,./’3 -0’ 26++’2"’/& 4!’2’ ,./’3’2 -0’ 0’*-"’/ ") 2),* 5),2"60’ -./,,.
()5#-0,2). ") ’3,2",.+ ,."’0#0’"-",).2 ); "!’ 42 $ =879 2#-(’,-0’ ().(’#"6-**1 -./ ()5#6"-",).-**1 2"0-,+!";)0%-0/& 4!’2’
,./’3’2 -0’ :-2’/ ). 2-"’**,"’>/’0,<’/ ,.;)05-",). ).*1,-./ "!’ #)"’.",-* ;)0 -##*,(-",). ,2 "!’0’;)0’ *-0+’&
?2,.+ @A 4’00- $ BC89D ,5-+’2 ;0)5 EFF@ ); "!’ 2’5,-0,/ =)0"! D!--.3, G)’22 H*-"’-6,=879 (BC8II)-./ 42
(BC8I@)%’0’ ,.<’2",+-"’/ -2 ,./,(-")02 ); <’+’"-",). -:6./-.(’ -./ *-./ 260;-(’ "’5#’0-"60’& 4) ’3-5,.’ 2),* 5),2"60’
()./,",).2 -" "!’ 260;-(’,"!’ 2#-",-* #-""’0. -./ "’5#)0-* ’<)*6",). ,. 7J9,4J9 -./ 4789 %-2 -.-*12’/ %,"! 42 -./
=879& 4!’ 0’26*"2 %’0’:(I)G-./ 260;-(’ 5),2"60’ 0’#0’2’."’/ :1 7J9 )0 4J9 2’#-0-"’*1 ()6*/ .)" ’*,5,.-"’ "!’ ’;;’("2
(-62’/ ;0)5 2’-2).-* (!-.+’ ); 42 -./ =879& K)%’<’0,4789,/’0,<’/ ;0)5 -.-*12,2 ); "!’ "’5#’0-"60’ $ <’+’"-",). ,./’3
2#-(’ -./ (0’-"’/ :1 ,."’+0-"’/ -.-*12’2 ); "!’2’ "%) L,./2 ); /-"-,!’*#’/ ") /’"’05,.’ 5-,. #0,.(,#*’2 ); "!’ "’5#)0-* -./
2#-",-* <-0,-",). ); 260;-(’ 2),* 5),2"60’& (E)M0)5 "’5#)0-* ’<)*6",). ); 4789,," ()6*/ :’ ,.;’00’/ "!-" "!’ "0’./ ,. 4789
,2 !,+! <-*6’2 ,. "!’ /01 2’-2). (2#0,.+ )0 -6"65.)-./ *)% <-*6’2 ,. "!’ 0-,.1 2’-2). (2655’0 )0 %,."’0)& (@)D#-",-*
’<)*6",). ); 4789 ,./,(-"’2 "!-" 260;-(’ 2),* 5),2"60’ ); ’-(! 0’+,). !-2 ,"2 )%. 2’-2).-* 0!1"!5& B’-. 4789 ,. "!’
.)0"!’0. 0’+,). ,.(0’-2’/ 0-#,/*1 ;0)5 %,."’0 ") 2#0,.+ -./ /’(0’-2’/ ;0)5 2655’0 ") -6"65. ,. 2)6"!;,. "!’ 5,//*’ -./
%’2","!’ ,.(0’5’." ); 4789 ;0)5 %,."’0 ") 2#0,.+ %-2 -*5)2" ’N6-* ") /’(0’5’." ;0)5 2655’0 ") -6"65.& (O)D,5,*-0*1,
2#-",-* <-0,-",). ); 260;-(’ 2),* 5),2"60’ ,. %-"’02!’/2 %-2 2’-2).-*& 41#,(-**1,,. C("):’0 "!’ <-0,-",). ,. 4789 %-2
+0’-"’2" ,. 5)2" %-"’02!’/2,’2#’(,-**1 ,. .)0"!&
!"# $%&’(:4789;*-./ 260;-(’ 5),2"60’;G)’22 H*-"’-6
地表水分是监测土地退化的一个重要指标,是气候、水文、生态、农业等领域的主要参数,在地表与大气界
面的物质和能量交换中起重要作用[I P @]。基于小范围内的采样数据来预测地表水分的程度及分布范围对于
区域以及全球尺度的研究而言,可操作性不强,要实现大范围内的地表水分监测,还需要借助遥感技术[O P A]。
近年来,基于遥感技术的地表水分估测方法被广泛应用于森林火灾中对树木可燃性的评估[Q,R]、水文模型的
参数验证[S]以及沙漠化或干旱对植被的水分压力的估测等[IF]。
目前遥感监测地表水分的方法主要通过实际蒸散和潜在蒸散比值计算作物缺水指数(J0)# T-"’0 D"0’22
9./’3,简称 JTD9)评价土壤水分状况[II P I@],但是该比值只适用于裸土或稀疏植被覆盖[IO]。在植被覆盖条件
下,通常采用其他的方法代替,例如利用时间序列植被指数构建的植被状态指数(7J9)[IU],温度状态指数
(4J9)[IA],以及 B(7,(-0等[IQ]建立的归一化温度指数(=)05-*,V’/ 8,;;’0’.(’ 4’5#’0-"60’ 9./’3,简称 =849),
其中 =849能很好地描述土壤供水能力,对变化的环境反应比 =879灵敏[IR]。最近的研究结果表明结合植被
指数和地表温度研究区域地表水分状况会得到更加合理的结果[IS,EF]。其中具有代表性的是 D-./!)*"等[EI]基
于植被指数和地表温度的关系,提出的温度植被干旱指数(4789)。
利用 4789反映地表水分状况的原理主要是植被状态及其水分压力可以通过热红外波段和可见光 $近红
外比值之间的互补性得到很好地表现[IR]。大量相关研究表明样点在地表温度 $植被指数空间中的位置受很
多因素影响[IO],其中具有代表性的是对土壤>植被>大气传输模型(2),*><’+’"-",).>-"5)2#!’0’ "0-.2;’0,简称
D7W4)的模拟结果的研究[EI,EE]。利用 D7W4模型可以估算地表温度 $植被指数特征空间的形状[E@],进而结合
对该特征空间生态特征的解释可以构建对土壤和作物水分含量具有一定的指示意义的 4789 指数[IR,EO]。通
过与地面采样获取的地表水分比较可以发现 4789所表征的地表水分空间分布及细节变化更为显著[EU,EA]。
本文选择陕北黄土高原地区作为研究区域,首先采用基于植被覆盖特征的 7J9 指数和基于地表温度特
征的 4J9指数分别评价了区域地表水分状况的季节变化;其次分析了植被覆盖与地表温度特征线性关系的季
节变化,在此基础上采用能同时反映植被 $温度特征的 4789 指数,从而综合评价了植被 $温度特征对黄土高
原区地表水分状况季节变化的影响。
OAUO X 生X 态X 学X 报X X X EQ 卷X
!""#:$ $ %%%& ’()*)+,(-& (.
!" 研究区概况
研究区地处陕西省北部,东经 /01234’至 ///2/5’,北纬 6523/’至 67268’间,东隔黄河,与山西省相望;西
以子午岭为界,与宁夏、甘肃接壤;北邻内蒙;南接渭南、铜川、咸阳三市。全区总面积 40909:;3,分属延安和
榆林市,下辖 35 个县区、645 个乡镇、7//0 个行政村。全区属暖温带和温带半干旱大陆性季风气候,四季分
明。区域内不同地区的气温除随纬度的增加而降低外,地势高低也具有显著的影响。通常高温出现在黄土高
原的东南部和低平地区,低温出现在西北部和较高的山区。区内东南部,年均气温一般在 /3& 5<以上,局部
地区可达 /8& 6<,西北部和北部地区一般在 3& 5<以下。1 月份平均气温大多超过了 33<,/ 月份平均气温各
地变化在 = / > /9<之间[31,34]。
图 /? 研究区流域划分
@,+& /? A!’ *)(-",). )B C"DEF -G’- -.E ,E’.",B,(-",). )B %-"’GC!’E
本区植被类型多样,具有明显的过渡性特点。随着气候特征的带状更迭,植被变化从东南向西北,由森林
草原、干草原、荒漠草原依次出现。森林草原主要分布在洛河中游地区,植被为白羊草、铁杆蒿、艾蒿、长芒草
为优势组成的草甸草原或草原。沙棘、荆条、酸枣、狼牙刺等较耐旱的灌木比较发达,侧柏、油松等耐旱树种在
林地中分布十分普遍。干草原主要分布在窟野河、无定河中下游、洛河上游等地区。荒漠草原分布在孤山川、
无定河上游等地区。由于该区人类活动历史悠久,强度大,自然植被破坏比较严重。目前该区的森林覆盖率
不足 /0H,坡耕地面积所占比例较大,导致水土流失十分严重[37]。
由于流域具有相对明确的自然边界,同时也是完整的景观生态系统[39,31]。为更好地反映景观格局及过
程的动态变化,本文以流域作为基本研究单元。首先基于 /I5 万的数字地形图,利用 JG(KLM 平台的水文分析
功能划分了研究区的主要流域,依次为北部区域内的孤山川、窟野河、秃尾河、佳芦河、榆溪河流域,西部区域
的无定河上游、头道川、洛河上游及部分内陆河流域,中部区域的清涧河、大理河、无定河中、下游、延河上、下
游流域及南部区域的云岩河、仕望河、洛河中游、葫芦河、沮水流域等(图 /)。
#" 数据及处理方法
#& !" 数据预处理
在研究中,为了精确评价研究区植被 $温度特征的空间分布,本研究使用国家农业部资源遥感与数字农业
重点开放实验室提供的 A’GG- $ NOPLM数据产品,包括地表温度和 QPRL 植被指数两种产品格式,空间分辨率
分别为 /000;和 350;,影像为 3006 年逐旬合成数据。数据预处理主要利用 SQRL 遥感影像处理平台,对照
/I5万地形图选取控制点,进行几何精校正,误差在半个象元以内。采用双标准纬线等积圆锥投影
(JTUSVM),椭球体为 WVJMORMWX,坐标系为 U’,Y,.+/758。通过与陕北黄土高原区行政边界掩膜
(NJMWLQK)处理,得到研究区范围内各时相的遥感影像。
5958? //期 ? ? ? 李正国? 等:基于植被 $温度特征的黄土高原地表水分季节变化 ?
!""#:$ $ %%%& ’()*)+,(-& (.
!& !" 基于植被指数的地表水分估测模型
由于从植被指数反演出的土地表面绿度与植物的生长状态及其密度密切相关,因此,植被指数可用于监
测对作物生长的环境条件,尤其是对地表水分环境的监测[/]。影响植物生长的因素很多,主要有气候、土壤、
天气和人类活动等,在这些因素中,在一定的连续时间内,可以认为气候和土壤处于相对不变的状态,只有天
气变化对植物生长具有短期的效应。因此,植被指数可用于表示植物的水分状况[0/]。在应用植被指数测量
地表水分状况的研究中,较有代表性的为植被状态指数 !"#,应用 !"# 动态监测地表水分变化的范围比应用
其它方法如 $%!#和降水量的监测更有效,同时认为 $%!#适用于研究大尺度范围的气候变异,而 !"#适用于
估算区域级的干旱程度[1]。植被状态指数的定义为[02]:
!"# 3
$%!#& 4 $%!#5,.
$%!#5-6 4 $%!#5,.
7 088 (0)
式中,$%!#&为某一特定年第 , 个时期的 $%!# 值,$%!#5-6和 $%!#5,.分别代表所研究年限内第 & 个时期
$%!#的最大值和最小值。
式中的分母部分是在研究年限内第 , 个时期植被指数的最大值和最小值之差,它在一定意义上代表了
$%!#的最大变化范围,反映了当地植被的生境;分子部分在一定意义上表示了某一特定年第 &个时期的当地
气象信息,若 $%!#&和 $%!#5,.之间差值小,表示该时段植物长势很差,此时 !"#偏低,地表水分状况较差。
!& #" 基于地表温度的地表水分估测模型
地表温度是控制地球表面大多数物理、化学和生物过程的参数之一,对裸土来说地表温度指的是土壤表
面温度,浓密植被覆盖的地表温度可以认为是植物冠层的表面温度。植物冠层温度升高是植物受到水分胁迫
和干旱发生的指示器,这是因为植物叶片气孔的关闭可以降低由于蒸腾所造成的水分损失,进而造成地表潜
热通量的降低,从而将导致地表感热通量的增加,而感热通量的增加又进一步造成植被冠层温度的升
高[00,09]。因此,地表温度可用于对地表水分状况的监测,其中较有代表性的为 ’"# 与 $%’#[/,0:]。’"# 的定义
与 !"#的定义相似,但它强调了温度与植物生长的关系,即高温对植物生长不利。’"#的定义为[01]:
’"# 3
’(& 4 ’(5,.
’(5-6 4 ’(5,.
(9)
式中,’(&为某一特定年第 & 个时期的地表温度值,’(5-6和 ’(5,.分别表示所研究年限内第 & 个时期地表温
度的最大值和最小值。’"#愈大,表示愈干旱。’"#的缺点是未考虑白天的气象条件,如净辐射、风速、湿度等
对热红外遥感的影响及地表温度的季节性变化[9;]。
!& $" 基于 ’( ) $%!#特征空间的地表水分估测模型
如果单独以遥感获取的地表温度作为指标,在植被覆盖不完全条件下,较高的土壤背景温度会严重干扰
土壤湿度信息[/]。植被指数提供了绿色植被的生长状况和覆盖度信息,如果结合光谱植被指数和陆表温度
的综合信息监测土壤湿度,可消除土壤背景的影响[<<]。
国内外学者研究了各种空间尺度和时间分辨率的地表温度和植被指数的关系,发现 ’(和 $%!#之间存在
明显的负相关关系[壤水分(土壤水分增加可以加大土壤的热惯量)更为直接相关,而不是作为对潜在热能的限制性控制[<1]。相
关研究表明,以 $%!#和 ’(为横纵坐标的散点图呈三角形[!)*" 等基于上述关系,提出了
温度植被干旱指数(’!%#)估测土壤表层水分状况[0;]。相关研究认为,在相同大气和地表湿度状况下,不同
的地表类型有着不同的 ’( ) $%!#斜率和截距[98]。
图 9 展示了 ’( ) $%!#特征空间的概念框架,主要体现了 ’(与 $%!#的关系。左侧边代表了不同湿度的裸
土的温度变化,横轴代表随着植被绿度的上升,最大地表温度下降。斜边表示在干旱条件下,对于某一给定的
地表类型和气候条件,地表温度所能达到的极限,可通过在 ’( ) $%!#特征空间中定义等值线来代表不同的干
旱程度。例如,’!%#值为 0 是干边(?@A ’>+’),代表有限的水分供应;’!%#值为 8 则是湿边(B’" ’>+’),具有
最大的土壤蒸发蒸腾总量和无限的水分供应。’!%#的计算公式为:
112/ C 生C 态C 学C 报C C C 9D 卷C
!""#:$ $ %%%& ’()*)+,(-& (.
图 /0 1234原理示意图[56]
7,+& /0 3’8,.,",). )8 "!’ 1234[56]
!"#$ 9
!% : !%;,.
& < ’(#"$ : !%;,.
(=)
式中,!%;,.为三角形中最小的地表温度,定义了相
应的湿边,!% 为给定像元的观测温度,(#"$ 为观测的
归一化植被指数值,&和 ’分别为定义干边的线性拟合
方程(!%;-> 9 & < ’(#"$)中的参数,!%;->为给定 (#"$
值下的最大地表温度。通过在从湿到干、从裸土到全植
被覆盖的各种条件下,对大范围区域内的像元采样来估
算参数 &和 ’。
需要指出的是,像元在 !% ) (#"$特征空间位置受到
诸多因素的影响。首先,!"#$ 的计算主要取决于 !% 和
(#"$;而稀疏植被地区的 !% 和 (#"$ 的提取同时受到
土壤和植被的影响而变得更为复杂[5?,5@],更重要的是
!%;->与 (#"$的线性关系随着季节不同发生显著的变
化(图 =),直接影响对参数 &和 ’的估计。
图 =0 1A;->与 (#"$各月线性关系
7,+& =0 B,.’-C C’*-",).A!,#A D’"%’’. 1A;-> -.E (#"$ ,. ;)."!A
!" 结果与分析
!& #" "*$的时空分布
本文基于 FG34H数据的植被指数产品,计算了研究区逐旬的 2I4数据集(图 6)。统计结果表明(表 5),
"*$均值的季节变化并不显著,其中在植被生长季节(J K ? 月份)均处于 L& / K L& 6 之间,表明期间植被生长
均不同程度受到水分条件的约束;另一方面,"*$方差的季节变化也不甚明显,仅在植被生长季节略高于其他
月份,表明期间植被生长受约束程度空间差异较大。
MNJ60 55期 0 0 0 李正国0 等:基于植被 $温度特征的黄土高原地表水分季节变化 0
!""#:$ $ %%%& ’()*)+,(-& (.
图 /0 研究区地表土壤水分特征的空间分布(1223)
4,+& /0 5,6"7,89",). ): *-.; 697:-(’ %-"’7 6"-"96 ,. 6"9;< -7’- (1223)
为研究 !"#的时空分布规律,对各流域 !"#的均值和方差的统计结果表明(表 =),北部区域秃尾河和榆
溪河流域一月的 !"#均值高于其他流域,变化范围多在 2$ 12 > 2$ /? 间,各流域 !"#方差值均在 2$ = 左右;/ 月
份各流域均值维持在 2$ 1? > 2$ 3? 间,方差值下降到 2$ 2? 左右;@ 月份各流域均值基本在 2$ 1? 左右,方差值保
持不变;=2 月份除佳芦河流域有所上升外,均值保持在 2$ 1? 左右,方差值稍有下降。南部区域各流域 = 月份
的 !"#均值多在 2$ /? 左右,方差值也在 2$ =2 左右;/ 月份均值略有下降至 2$ /2 左右,方差值基本保持不变;@
月份均值则大幅上升至 2$ A2 > 2$ @2 间,方差值也有所上升;=2 月份各流域的均值下降至 2$ A2 左右,而方差
BA?/ 0 生0 态0 学0 报0 0 0 1@ 卷0
!""#:$ $ %%%& ’()*)+,(-& (.
值回复至 /! 0/ 左右。西部区域各流域 0 月份均值多在 /! 1/ 2 /! 3/ 间,方差值也在 /! 0 左右;3 月份均值下降
至 /! 45 2 /! 1/ 间,方差值亦降至 /! /5 以下;6 月份均值回升至 /! 1/ 左右,方差值略有上升;0/ 月份均值基本
维持不变,方差值维持在 /! /5 左右。中部区域各流域 0 月份均值多在 /! 4/ 2 /! 1/ 间,方差值也在 /! 0/ 以
下;3 月份均值略有下降,方差值亦降至 /! /5 以下;6 月份均值进一步上升至 /! 1/ 2 /! 3/ 左右,方差值在 /! 0/
以下;0/ 月份均值略有上升,方差值则小幅下降。
表 !" 研究区各流域 #$%的统计值
&’()* !" +,’,-.,-/. 01 #$% -2 3’,*4.5*6. -2 .,768 ’4*’
项目 7"’8
北部区域 9):"! #-:"
孤山川
90
窟野河
94
秃尾河
91
榆溪河
93
佳芦河
95
南部区域 ;)<"! =-:"
云岩河
;0
洛河中游
;4
葫芦河
;1
仕望河
;3
沮水
;5
0 月 均值 >’-. /& 46 /& 14 /& 35 /& 50 /& 44 /& 3/ /& 36 /& 36 /& 31 /& 51
?-.& 方差 @-:,-.(’ /& /A /& 00 /& 05 /& 00 /& /6 /& /B /& 0/ /& /6 /& 0/ /& 04
3 月 均值 /& 10 /& 14 /& 4B /& 46 /& 43 /& 13 /& 3/ /& 3/ /& 1B /& 35
C#:& 方差 /& /5 /& /D /& /5 /& /1 /& /D /& /D /& /B /& /D /& /6 /& /A
6 月 均值 /& 4D /& 46 /& 4D /& 4A /& 45 /& 50 /& D/ /& DA /& DD /& 63
?<*& 方差 /& /3 /& /3 /& /B /& 00 /& /3 /& 01 /& 04 /& /A /& 05 /& /A
0/ 月 均值 /& 45 /& 45 /& 4D /& 4D /& 1/ /& 3A /& 55 /& D/ /& D4 /& 6/
E("& 方差 /& /1 /& /1 /& /5 /& /1 /& /5 /& /A /& 0/ /& /6 /& 00 /& 0/
项目 7"’8
西部区域 F’G" #-:"
头道川
F0
洛河上游
F4
内陆河
F1
无定河
上游 F3
中部区域 >,HH*’ #-:"
清涧河
>0
大理河
>4
无定河中游
>1
无定河下游
>3
延河上游
>5
延河下游
>D
0 月 均值 >’-. /& 11 /& 11 /& 3/ /& 3/ /& 46 /& 46 /& 41 /& 41 /& 1/ /& 4A
?-.& 方差 @-:,-.(’ /& /D /& /B /& 0/ /& 01 /& /6 /& /D /& /5 /& /5 /& /6 /& /B
3 月 均值 /& 46 /& 4A /& 4D /& 4D /& 45 /& 41 /& 41 /& 41 /& 4D /& 4B
C#:& 方差 /& /4 /& /3 /& /1 /& /3 /& /1 /& /3 /& /5 /& /1 /& /3 /& /D
6 月 均值 /& 10 /& 13 /& 14 /& 4B /& 10 /& 41 /& 45 /& 43 /& 10 /& 1B
?<*& 方差 /& /D /& /6 /& /B /& /A /& /3 /& /3 /& /B /& /5 /& /6 /& /B
0/ 月 均值 /& 1/ /& 1D /& 4D /& 4B /& 10 /& 4B /& 4A /& 11 /& 13 /& 1B
E("& 方差 /& /3 /& /D /& /3 /& /3 /& /3 /& /1 /& /3 /& /3 /& /D /& /D
9& :" IJ7的时空分布
基于 >EK7;数据提取的地表温度,计算了研究区逐旬的 IJ7数据集(图 3)。"#$均值的季节变化则较为
显著,在植被生长季节均处在 /& 3 以上,其中 5 月份达 /& 55,表明期间地表蒸发强烈,导致地表水分条件较
差,而其他月份均在 /& 45 左右,此时地表水分相对较为湿润;另一方面,"#$ 方差在植被生长季节较高,普遍
在 /L 05 以上,表明期间地表水分状况的区域内部差异较大。
为研究 "#$的时空分布规律,对各流域 IJ7的均值和方差的统计结果表明(表 4),北部区域各流域 0 月
份的 "#$均值在 /& 0/ 2 /& 4/ 间,"#$方差值均在 /& /5 左右;3 月份各流域均值上升至 /& 5/ 左右,方差值上升
到 /& /5 2 /& 0/ 间;6 月份各流域均值基本维持不变,方差值则略有下降;0/ 月份回落至 /& 3/ 以下,方差值上
升至 /L 0/ 2 /& 45 间。南部区域各流域 0 月份均值为 /& 1/ 2 /& 15 间,方差值均在 /& /5 左右;3 月份各流域均
值上升至 /L 55 2 /& D/ 间,方差值略有上升;6 月份各流域均值明显回落至 /& 15 以下,方差值则基本不变;0/
月份进一步降至 /& 0/ 2 /& 4 间,方差值差别较大,取值在 /& /5 2 /& 05 间不等。西部区域各流域 0 月份均值多
在 /& 45 2 /& 15 间,方差值也在 /& 0 以下;3 月份均值上升至 /& 5/ 2 /& 55 间,方差值则基本不变;6 月份均值除
无定河上游流域外均下降至 /& 5/ 以下,方差值变化很小;0/ 月份均值则进一步降至 /& 1D 以下,方差值上升
至 /& 0/ 以上。中部区域各流域 0 月份均值多在 /& 4/ 2 /& 1/ 间,方差值均在 /& /5 左右;3 月份均值则大幅上
升至 /& 3/ 以上,方差值轻微上升至 /& /5 2 /& 0/ 间;6 月份清涧河与延河流域的 "#$ 均值有所下降,其他流域
AD53M 00期 M M M 李正国M 等:基于植被 $温度特征的黄土高原地表水分季节变化 M
!""#:$ $ %%%& ’()*)+,(-& (.
均继续上升,而方差值基本维持不变;/0 月份各流域均值明显下降,方差值保持在 0& /0 以下。
表 !" 研究区各流域 !"#的统计值
#$%&’ !" ()$)*+)*,+ -. !"# */ 0$)’1+2’3+ */ +)435 $1’$
项目 1"’2
北部区域 3)4"! #-4"
孤山川
3/
窟野河
35
秃尾河
36
榆溪河
37
佳芦河
38
南部区域 9):"! ;-4"
云岩河
9/
洛河中游
95
葫芦河
96
仕望河
97
沮水
98
/ 月 均值 <’-. 0& 0= 0& // 0& /6 0& /> 0& /? 0& 6/ 0& 66 0& 66 0& 68 0& 6@
A-.& 方差 B-4,-.(’ 0& 07 0& 07 0& 08 0& 08 0& 07 0& 08 0& 0> 0& 07 0& 0> 0& 07
7 月 均值 0& 8/ 0& 87 0& 86 0& 87 0& 7? 0& 8@ 0& 8> 0& 8? 0& 8@ 0& 8@
C#4& 方差 0& 08 0& 0= 0& 0? 0& 0= 0& 0= 0& 0= 0& 0> 0& 0@ 0& 0> 0& 0@
> 月 均值 0& 75 0& 80 0& 8@ 0& 8? 0& 8@ 0& 67 0& 58 0& /= 0& 5> 0& />
A:*& 方差 0& 0? 0& 0@ 0& 0@ 0& 0> 0& 0> 0& 0> 0& 08 0& 08 0& 0= 0& 07
/0 月 均值 0& 6/ 0& 6> 0& 6/ 0& 55 0& 58 0& /7 0& /> 0& /8 0& /5 0& /@
D("& 方差 0& 56 0& 5/ 0& /= 0& /5 0& /= 0& 07 0& /8 0& 0= 0& 0= 0& /5
项目 1"’2
西部区域 (E’F" #-4")
头道川
E/
洛河上游
E5
内陆河
E6
无定河
上游 E7
中部区域 (<,GG*’ #-4")
清涧河
大理河
<5
无定河中游
<6
无定河下游
<7
延河上游
<8
延河下游
<@
/ 月 均值 <’-. 0& 68 0& 6/ 0& 66 0& 57 0& 5? 0& 5> 0& /? 0& 56 0& 6/ 0& 5=
A-.& 方差 B-4,-.(’ 0& 0? 0& 0@ 0& 0@ 0& 0@ 0& 08 0& 08 0& 07 0& 08 0& 0@ 0& 0@
7 月 均值 0& 86 0& 87 0& 8/ 0& 8/ 0& 86 0& 7= 0& 80 0& 77 0& 87 0& 8@
C#4& 方差 0& 0? 0& 0> 0& 0> 0& 0? 0& /0 0& 0@ 0& 0@ 0& 0= 0& 0> 0& 0?
> 月 均值 0& 70 0& 68 0& 7> 0& 87 0& 77 0& 8/ 0& 86 0& 8? 0& 67 0& 6?
A:*& 方差 0& 0> 0& 0@ 0& 0> 0& 0= 0& 0@ 0& 0@ 0& 08 0& 0> 0& 0@ 0& 0>
/0 月 均值 0& 5> 0& 57 0& 68 0& 6@ 0& /7 0& 57 0& 50 0& /6 0& 56 0& /7
D("& 方差 0& // 0& 0= 0& /6 0& /@ 0& 0> 0& /0 0& 0= 0& 06 0& 0= 0& 0>
6& 6" !"2-H与 #$%&线性关系的季节变化
采用光谱植被指数和地表温度的综合信息监测地表水分,必须首先确定在不同的季节,某一植被覆盖条
件下,最大可能的地表温度,从而确立定义干边的线性拟合方程参数。
从对研究区像元随机采样的统计结果分析来看,在一月,#$%& 与 !"2-H呈现正相关关系(图 6-),当 #$%&
值每上升 /0(0 I 588),!"2-H升高 /& @J,两者相关系数高达 0& ?6。表明在冬季月份,植被覆盖较好的地区,地
表温度相应也较高。这主要是因为期间南部的平均气温要明显高于北部,而植被覆盖良好的南部地区其水分
条件也相应较好,从而起到了保温的作用,两者叠加使得 !"2-H与 #$%&呈现正相关关系。
在 7 月份,#$%&与 !"2-H呈现的相关程度较低(图 6K),当 #$%&值每上升 /0,!"2-H升高 0& @J,两者相关系
数仅 0& /=。主要是因为此时研究区北部的气温要略高于南部,一定程度上抵消了南部植被覆盖的保温效果,
两者共同作用导致 !"2-H与 #$%&几乎不相关。
在 > 月份,#$%&与 !"2-H呈现负相关关系(图 6(),当 #$%&值上升每上升 /0,!"2-H下降 5& 0J,两者相关系
数高达 0& =8。表明在夏季月份,植被覆盖较好的地区,地表温度相应也较低。这主要是因为此时研究区北部
的气温要显著高于南部,而植被覆盖良好的南部地区水分条件也相应较好,其蒸发效应起到了降温的作用,两
者共同作用使得 !"2-H与 #$%&呈现显著的负相关关系。
在 /0 月份,#$%&与 !"2-H的相关关系亦偏低(图 6G),当 #$%&值每上升 /0,!"2-H升高 0& >J,两者相关系
数仅 0& 68。表明在秋季月份,植被覆盖的程度与地表温度的线性关系并不明显。这主要是因为此时研究区
南北的气温差异并不显著,而植被的凋零使得植被覆盖的差异也不缩小,两者共同作用导致 !"2-H与 #$%& 几
乎不相关。
6& 7" !%$&的时空分布
基于
0>87 M 生M 态M 学M 报M M M 5> 卷M
!""#: %%%& ’()*)+,(-& (.
!"#数据集(图 /)。!"#均值的季节变化非常显著,在植被生长季节均在 0& 1 以上,其中 2 月份达 0& 30,表
明期间地表水分状况较差;另一方面,!"#方差也表现为在植被生长季节较高,普遍在 0& 41 以上,表明期间 地表水分状况的区域差异较大。 为研究地表水分特征的时空分布规律,对不同区域内各流域 !"# 的均值和方差加以统计分析(表 5)。
结果表明,北部区域内的各流域 !"#值变化的趋势主要表现为在冬季(如4月份)整体均值最低,同时空间差异性最小;到/月份时均值达到较高的水平(0的均值下降到 0% 51 以下,方差值上升到 0% 4 6 0% 8 间。南部区域各
流域 4 月份的 !"#均值多在0的均值上升到 0% 11 6 0% 20 间,方
差值基本保持不变;3 月份均值基本维持在 0% 20 左右,方差值均不足 0% 40;40 月份各流域的均值下降至 0% 80
6 0% 81 间,方差值基本维持不变。西部区域各流域 4 月份均值多在 0% 50 6 0% /0 间,方差值也在 0% 4 以下;/
月份均值上升到 0% 10 以上,方差值保持不变;3 月份均值进一步上升至 0% 30 左右,方差值维持不变;40 月份
均值下降到 0% /0 以下,方差值也上升到 0% 40 6 0% 41 间。中部区域各流域 4 月份均值多在 0% 80 6 0% 50 间,方
差值也在 0% 40 以下;/ 月份均值上升到 0% 10 以上,方差值维持不变;3 月份均值进一步上升至 0% 20 左右,方
差值不变;40 月份均值下降到 0% 80 6 0% 50 间,方差值不变。
表 !" 研究区各流域 !"#的统计值 #%&’ !" ())*+)*,+ -. !"# */ 0)′1+2′3+∗/+)4351’项目 9"’: 北部区域 ;)<"! #-<" 孤山川 ;4 窟野河 ;8 秃尾河 ;5 榆溪河 ;/ 佳芦河 ;1 南部区域 =)>"! ?-<" 云岩河 =4 洛河中游 =8 葫芦河 =5 仕望河 =/ 沮水 =1 4 月 均值 @’-. 0& 0A 0& 48 0& 4/ 0& 4B 0& 4A 0& 58 0& 5/ 0& 55 0& 51 0& 51 C-.& 方差 D-<,-.(’ 0& 0/ 0& 0/ 0& 02 0& 01 0& 0/ 0& 01 0& 03 0& 0/ 0& 03 0& 0/ / 月 均值 0& 15 0& 11 0& 1/ 0& 12 0& /A 0& 13 0& 1B 0& 1A 0& 12 0& 12 E#<& 方差 0& 01 0& 0A 0& 0A 0& 0A 0& 0A 0& 0A 0& 03 0& 02 0& 03 0& 02 3 月 均值 0& 25 0& 21 0& 2B 0& 3/ 0& 23 0& 24 0& 20 0& 20 0& 12 0& 1A C>*& 方差 0& 02 0& 02 0& 01 0& 02 0& 01 0& 03 0& 03 0& 02 0& 02 0& 01 40 月 均值 0& 5/ 0& 5A 0& 5/ 0& 83 0& 8A 0& 84 0& 81 0& 85 0& 80 0& 8/ F("& 方差 0& 80 0& 4B 0& 42 0& 40 0& 42 0& 05 0& 45 0& 0B 0& 0A 0& 40 项目 9"’: 西部区域 G’H" #-<" 头道川 G4 洛河上游 G8 内陆河 G5 无定河 上游 G/ 中部区域 @,II*’ #-<" 清涧河 @4 大理河 @8 无定河中游 @5 无定河下游 @/ 延河上游 @1 延河下游 @2 4 月 均值 @’-. 0& 53 0& 55 0& 5/ 0& 82 0& 50 0& 8B 0& 80 0& 8/ 0& 58 0& 54 C-.& 方差 D-<,-.(’ 0& 0B 0& 02 0& 02 0& 03 0& 02 0& 01 0& 0/ 0& 01 0& 02 0& 02 / 月 均值 0& 1/ 0& 12 0& 18 0& 18 0& 1/ 0& 14 0& 14 0& /2 0& 12 0& 1B E#<& 方差 0& 0A 0& 03 0& 03 0& 0B 0& 40 0& 02 0& 02 0& 0A 0& 0B 0& 03 3 月 均值 0& 30 0& 2B 0& 38 0& 30 0& 28 0& 2/ 0& 23 0& 28 0& 28 0& 28 C>*& 方差 0& 0/ 0& 01 0& 0/ 0& 0B 0& 01 0& 02 0& 03 0& 02 0& 01 0& 0/ 40 月 均值 0& 58 0& 8A 0& 53 0& 5B 0& 84 0& 8B 0& 82 0& 80 0& 8B 0& 84 F("& 方差 0& 0A 0& 0B 0& 44 0& 4/ 0& 02 0& 0B 0& 0B 0& 08 0& 0B 0& 02 !& 6" "&、!&与 !"#的相关关系
为比较各地表水分指标的相关性,对各月 "&、!&及 !"#的相关系数进行统计分析(表 /)。计算结果 表明,"&与 !&的相关系数有明显的季节变化,在春夏季(如462月份)及秋冬季(如44648月份),两者均表现为相关关系不明显,表明在上述季节采用植被覆盖与温度表征地表水分的差异较大;而在夏秋季(如36A月份),两者则表现为较强的负相关关系,其中A月份的相关系数值达J0温度特征的黄土高原地表水分季节变化 K
!""#: %%%& ’()*)+,(-& (.
!"#与"的相关系数总体偏低,并存在一定程度的季节变化。在冬春季(如 // 0 1 月份),两者表现为
微弱的负相关关系,表明在上述季节两者对地表水分的表征能力差异较大;在夏秋季(如 2 0 3 月份),两者相
关程度较高,其中 3 月份达 4 5& 65,表明在该季节两者对地表水分状况的表征能力大致相当。
!"#与!的相关系数全年均在 5& 65 以上,部分月份甚至达 5& 72 以上,表明两者对地表水分状况的表
征能力完全相当。总的来看,单独采用 !%与"表征地表水分会受到明显的季节影响,而 !"#虽然和!
关系更为密切,但能在不同季节综合体现植被覆盖和地表温度特征对地表水分的影响,从而能较好地反映区
域地表水分状况的时空变化特征。
表 !" !"#、"#与!%#各月相关系数
#+,)- .)’//,.,’-+0 )/ !"#,"#-1 !与 !%相关系数!! 5& // 5& /2 5& ;/ 5& ;5 5& /C 5& 52 4 5& ;@ 4 5& @7 4 5& 7/ 4 5& 1@ 5& // 5& 53 "%与 !"#相关系数" 4 5& 5C 4 5& 53 4 5& /; 4 5& ;6 4 5& 13 4 5& @6 4 5& @5 4 5& 11 4 5& 65 4 5& ;7 4 5& /5 4 5& 5C !%与 !"#相关系数# 5& 77 5& 63 5& 7C 5& 76 5& 63 5& 72 5& 7C 5& 76 5& 6C 5& 63 5& 67 5& 62 J J ! K)??’*-",). ()’LL,(,’."M =’"%’’. "% -.N !%;"K)??′∗−",).()′LL,(,′."M=′" -.N !"#;# K)??’*-",). ()’LL,(,’."M =’"%’’. !% -.N
!"#;!通过 5& 52 显著性水平检验 E,+.,L,(-." (5& 5/ O ’ O 5& 52) 4& 5" !"#的相关检验
上文的分析表明 !"#相对于其他指数能较好地体现地表水分的分布特征,为了进一步验证其与地表水 分的对应关系,本文通过利用相关的气候因子对 !"# 进行了空间验证。首先对各月 !"#进行累积并求取 其均值,进一步对该计算结果按密度分级,共划分为 /5 级(图 2-)。结合该地区年均降雨量、年均温以及相应 的湿润系数(图 2=)的空间插值结果,统计了不同分级内各指标的均值与方差(表 2)。结果表明,随着 !"#
值增加,地表干燥度随之升高,表现为年均温从 ;& 1CP持续上升到 /;& 51P,且空间变异有所加大;年均降雨
量从 22C& 3199下降至 1@3& 1299;而湿润系数则从 4 2& 2@ 降低至 4 ;2& 65。上述指标的变化规律表明,!"#能有效地反映研究区地表温度、降水以及植被覆盖对地表水分空间分布的综合影响。表6"!%#与验证因子的关系
#+,)-03,8 )/ !-1 ,+0 9&,1+,-: /.+)*0 分级 K*-MM QRI8 均值 >’-. 方差 R-?,-.(’ 年均温 AH’?-+’ -..D-* "’9#’?-"D?’ 均值(P) >’-. 方差 R-?,-.(’ 年均降雨量 AH’?-+’ -..D-* #?’(,#,"-",). 均值(99) >’-. 方差 R-?,-.(’ 湿润系数 >),M"D?’ 8.N’S 均值 >’-. 方差 R-?,-.(’ / 5& /7 5& 51 ;& 1C 5& 7/ 22C& 31 5& @2 4 2& 2@ 5& /@ ; 5& ;2 5& 5; 1& C1 5& 32 2;3& @7 5& 23 4 7& 7; 5& ;2 1 5& 15 5& 5/ @& 7C 5& 73 2/2& 2/ 5& 2C 4 //& /@ 5& ;7 @ 5& 12 5& 5; C& /; 5& 31 @3/& 25 5& C5 4 /1& 1C 5& ;C 2 5& @/ 5& 5; C& 66 5& 77 @26& C1 5& 25 4 /2& 76 5& /3 C 5& @6 5& 5; 6& @/ 5& 7@ @1@& 12 5& 17 4 /6& 2@ 5& /6 6 5& 21 5& 5; 7& /3 5& 3@ @56& 13 5& ;7 4 ;5& ;1& 5& /1 7 5& C5 5& 5; 3& 56 /& 5/ 135& ;@ 5& ;; 4 ;;& 57& 5& 53 3 5& CC 5& 5; /5& 22 5& 3C 162& 7C 5& ;C 4 ;1& ;C 5& 57 /5 5& 62 5& 52 /;& 51 /& /5 1@3& 12 5& 11 4 ;2& 65 5& /5 !" 结论与讨论 从利用植被指数测量地表水分状况的结果来看,各区域 "%均值的季节分布上存在较大差异,其中在北
部区域表现为 @ 月份较高,/5 月份较低的分布规律;中部区域与北部恰恰相反,为 @ 月份偏低,/5 月份较高;
南部区域与之类似,仅 6 月份偏高;西部区域则基本是 @ 月份较低而 / 月份偏高。同时各区地表水分状况的
;62@ J 生J 态J 学J 报J J J ;6 卷J
!""#: %%%& ’()*)+,(-& (.
图 /0 !"#分级及其相关验证因子 1,+& /0 2*-33,4,(-",). )4 !"# -.5 ,"3 6-*,5-",.+ 4-(")73
季节变化亦存在不同的特点,其中北部区域中秃尾河和榆溪河在 8 9 : 月份的 "%均值下降最为明显,反映期间上述流域地表水分状况有明显变干的趋势;南部区域的地表水分状况则是在:9;月份上升最为显著;西部及中部区域地表水分状况全年均处于偏低的水平。地表水分特征的空间变异性也存在季节差异,区内大部分流域在8月份的"方差值达 <& 8< 以上,其中北部区域尤为明显,反映了各流域内部地表水分状况差异在冬
季较为显著。
从利用地表温度测量地表水分状况的结果来看,各区域的 !%均值皆为:9;月份高于8<9翌年8月份,表现了区域地表水分状况的季节分布规律。但各区地表水分状况的季节变化存在不同的特点,其中北部区域89:月份的!均值上升最为明显,反映期间区域地表水分状况明显变干;南部及中部区域的地表水分状况
则在 ; 9 8< 月份间显著趋湿;而西部区域地表水分状况在 8 9 : 月份的变化幅度与 ; 9 8< 月份大致相当。
从植被指数与地表温度的线性关系来看,在 8 9 ; 月份,区域植被覆盖的分布格局会进一步增加南部与北
部气温的差异;而 : 月份和 8< 月份,则会缩减南北气温的差异。因此,采用 !"#评价地表水分状况时,必须 考虑地表气温与植被覆盖线性关系的季节变化。同时地表水分特征的空间变异性也存在季节差异,区内大部 分流域在 8< 9翌年 8 月份的 !"#方差值达 <& 8< 以上,其中北部区域尤为突出,表明期间各流域内部地表水
分状况差异较为显著。
从综合反映植被覆盖和地表温度特征的 !"#的时空分布来看,各区域的 !"#均值皆为 : 9 ; 月份高于
8< 9翌年 8 月份,表现了区域地表水分状况的季节分布规律。但各区地表水分状况的季节变化存在不同的特
点,其中北部区域 8 9 : 月份的 !"#均值上升最为明显,反映期间区域地表水分状况明显变干;南部区域的地 表水分状况则在 ; 9 8< 月份间显著趋湿;而西部及中部区域地表水分状况在 8 9 : 月份的上升幅度与 ; 9 8< 月份的下降幅度大致相当。同时地表水分特征的空间变异性也存在季节差异,区内大部分流域在 8< 月份的 !"#方差值达 <& 8< 以上,其中北部区域尤为突出,反映期间各流域内地表水分状况差异较显著的特征。
从分析结果来看,单独采用 !%与"表征地表水分会受到明显的季节影响,而 !"#虽然和!关系更
=;/:0 88期 0 0 0 李正国0 等:基于植被 温度特征的黄土高原地表水分季节变化 0 !""#: 能有效地反映研究区地表温度、降水及植被
覆盖对地表水分空间分布的综合影响。
!"#原理简单,应用性较强,但尚存在不足之处,主要表现在,在较高的 值下,!"#等值线间距变 小,导致 !"#不确定性增加;其次,!"#方法将 !& ’ %#"特征空间由梯形简化为三角形更增加了在 %#"高 值情况下,!"#的不确定性;再次,其将湿边看作水平直线,而不是在梯形特征空间中的斜线,也导致在 %#"低值区对 !"#的高估。另外,由于 !"#方法的内在限制,其提取的水分信息同时包含了地表土壤水分和植 被表层水分,因此,在未来研究中应通过对地表覆盖类型的细分、地面采样数据的时空匹配等途径来提升 !"#的应用水平。
!"#"$"%&"’:
[/ ]0 1-) 2 3,4!-.+ 4 5,6-.+ 5& 78-*9-",.+ :),* ;),:"9<’ :"-"9: ,. 5,.=,-.+ 9:,.+ "!’ >’;#’<-"9<’ ?’+’"-",). @?@B)& E’;)"’
3’.:,.+ >’(!.)*)+A -.C F##*,(-",).,GHHI,/J(K):ILM ILN&
[G ]0 E-. O,4!-.+ 4 5,4!-.+ P Q,() *+& @7R ()<<’(",). 9:,.+ >?@B ") ’8-*9-"’ :),* ;),:"9<’ :"-"9: ,. 2!,.-& 3(,’.(’ )S 3),* -.C 6-"’<
2).:’<8-",).,GHHT,M(G):MG MK&
[M ]0 U9 V W,6-.+ W,R- X R& 7SS’(" )S *-.C 9:’ ). :),* %-"’< ,. *)’:: !,** -<’-& 3(,’.(’ U)9.C-",). ,. 2!,.-,/JJJ,/M(I):GGT GGL&
[I ]0 6-.+ Q 5,6-. 4 R,P).+ W 1,() *+& FC8-.(’: ,. C<)9+!" ;).,")<,.+ YA 9:,.+ <’;)"’*A :’.:’C .)<;-*,Z’C C,SS’<’.(’ 8’+’"-",). ,.C’D -.C *-.C
:9[T ]0 [-;Y,. 7 U,7!<*,(! @& >!’ :93’.:,.+,/JJK,/L:IKM INL&
[K ]0 [,9 6,X)+-. U ]& R).,")<,.+ <’+,).-* C<)9+!" 9:,.+ "!’ 8’+’"-",). ().C,",). ,.C’D& B."’<.-",).-* W)9<.-* )S E’;)"’ 3’.:,.+,/JJK,/L:GLK/
GLNG&
[L ]0 ?’(’C&)& E)""’[N ]0 2!98,’() 7,F+9-C) B,2)(’<) @,() *+& @’:,+. )S -. ’;#,<,(-* ,.C’D ") ’:",;-"’ S9’* ;),:"9<’ ()."’." S<); ]_FF F?‘EE ,;-+’: ,. S)<’:" S,<’
C-.+’< :"9C,’:& B."’<.-",).-* W)9<.-* )S E’;)"’ 3’.:,.+,GHHM,GI:/KG/ /KML&
[J ]0 F.C’<:’. W,3-.C!)*" B,W’.:’. X ‘,() *+& Q’<:#’(",8’: ,. 9:,.+ - <’;)"’*A :’.:’C C:(-*’& ‘AC<)*)+,(-* Q<)(’::’:,GHHG,/K:GJLM GJNL&
[/H]0 X-<.,’*, F,@-** a_*;) P& E’;)"’ :’.:,.+ ;).,")<,.+ )S C’:’<",S,(-",).,#!’.)*)+A,-.C C<)9+!":& R-.-+’;’." )S 7.8,<).;’."-* O9-*,"A:F.
B."’<.-",).-* W)9<.-*& GHHM,/I:GG MN&
[//]0 W-(^:). E @,BC:) 3 V,E’+,.-") E W& 2-.)#A "’;#’<-"9<’ -: - (<)# %-"’< :"<’:: ,.C,(-")<& 6-"’< E’:)9<(’ E’:’-<(!,/JN/,/L://MM //MN&
[/G]0 4!-.+ E ‘& F .’% ;)C’* S)< ’:",;-",.+ (<)# %-"’< :"<’:: Y-:’C ). ,.S<-<’C <-C,-",). ,.S)<;-",). & 3(,’.(’ ,. 2!,.- V,/JNK,L:LLK LNI&
[/M]0 >,-. P9)*,-.+& 7:",;-",). )S ’8-#)"<-.:#,<-",). -.C :),* ;),:"9<’ -.C C<)9+!" ;).,")<,.+ 9:,.+ <’;)"’ :’.:,.+ ,. .)<"! 2!,.- #*-,.& 3#-(’ -.C
7.8,<).;’.",BFU,P<-Z F9:"<,-,/JJM,/GN /MT&
[/I]0 R)<-. R 3,2*-<^’ > E,B.)9’ 1,() *+& 7:",;-",.+ (<)# %-"’< C’S,(," 9:,.+ "!’ <’*-",). Y’"%’’. :9,.C’D& E’;)"’ 3’.:,.+ )S 7.8,<).;’.",/JJI,IJ:GIK GKM&
[/T]0 X)+-. U ]& E’;)"’ :’.:,.+ )S %’-"!’< ,;#-(": ). 8’+’"-",). ,. .)<\!);)+’.’)9: -<’-:& B."’<.-",).-* W)9<.-* )S E’;)"’ 3’.:,.+,/JJH,//:/IHT
/I/J&
[/K]0 X)+-. U ]& F##*,(-",). )S 8’+’"-",). ,.C’D -.C Y<,+!".’:: "’;#’<-"9<’ S)< C<)9+!" C’"’(",).& FC8-.(’: ,. 3#-(’ E’:’-<(!,/JJT,/T:J/ /HH&
[/L]0 R(?,(-< > E,W9## @ [ V,1-.+ 5,() *+& [,.^,.+ <’+,). %-"’< Y-*-.(’ ;)C’*: %,"! <’;)"’ :’.:,.+& B.:Q<)(’’C,.+ )S "!’ /M"! F:,-. 2).S’<’.(’
). E’;)"’ 3’.:,.+& b*--.Y--"-<,R).+)*,-:/JJG& VK& / VK& K&
[/N]0 3-.C!)*" B,E-:;9::’.- X,F.C’<:’.Y W& F :,;#*’ ,."’<#<’"-",). )S "!’ :9;),:"9<’ :"-"9:& E’;)"’ 3’.:,.+ )S 7.8,<).;’.",GHHG,LJ:G/M GGI&
[/J]0 3;,"! E 2 P,2!)9C!9W)9<.-* )S E’;)"’ 3’.:,.+,/JJ/,/G (/H):GHG/ GHII&
[GH]0 P)’"Z 3 W& R9*",:’.:)< -.-*A:,: )S ]@?B,:93’.:,.+,/JJL,/N(/T):L/ JI&
ILTI 0 生0 态0 学0 报0 0 0 GL 卷0
!""#: %%%& ’()*)+,(-& (.
[/0]1 2,**,’3 4 4,5-6*3). 7 8,29, :,!" #& ; <’6,=,(-",). )= "!’ "6,-.+*’>?’"!)@ =)6 )A"-,.,.+ 396=-(’ 3),* %-"’6 ()."’." -.@ ’.’6+B =*9C’3 =6)? 6’?)"’ ?’-396’?’."3 )= "!’ 8)6?-*,D’@ E,==’6’.(’ F’+’"-",). G.@’C (8EFG)-.@ 396=-(’ 6-@,-." "’?#’6-"96’& G."’6.-",).-* :)96.-* )= 4’?)"’ H’.3,.+, 0IIJ,0K (0L):M0NL M0OO& [//]1 P)6-. P H,5*-6Q’ 7,R93"-3 S T,!" #& ;& U<-*9-",). )= !B@6)*)+,( #-6-?’"’63 ,. - 3’?,-6,@ 6-.+’*-.@ 93,.+ 6’?)"’*B 3’.3’@ 3#’("6-* @-"-&
S-"’6 4’3)96(’3 4’3’-6(!,0IIN,MV (L):0/KJ 0/IJ&
[/M]1 5-6*3). 7 8,2,**,’3 4 4,H(!?9++’ 7 :& ;. ,."’6#6’"-",). )= ?’"!)@)*)+,’3 =)6 ,.@,6’(" ?’-396’?’." )= 3),* %-"’6 ()."’."& ;+6,(9*"96-* -.@ W)6’3"
P’"’)6)*)+B,0IIJ,JJ:0I0 /VL&
[/N]1 W6,’@* P ;,E-<,3 W S& H)96(’3 )= <-6,-",). ,. 6-@,)?’"6,( 396=-(’ "’?#’6-"96’ )<’6 - "-**+6-33 #6-,6,’& 4’?)"’ H’.3,.+ )= U.<,6).?’.",0IIN,NK:
0 0J&
[/L]1 S-.+ 2 H,X,- :,8,9 5 S& W*)% 6)9",.+ ?’"!)@ -.@ ,"3 -##*,(-",). ,. @,3"6,A9"’@ !B@6)*)+,(-* ?)@’*,.+& 2’)+6-#!,(-* 4’3’-6(!,/VVN,/M(/):
0JL 0K/&
[/O]1 Y!’.+ Z X,[,9 5 P,S-.+ Y 2& H,?9*-",). )= !B@6)*)+,(-* #6)(’33’3 ,. [93!, \-3,. A-3,.+ @,3"6,A9"’@ !B@6)*)+,(-* ?)@’*& 2’)+6-#!,(-*
4’3’-6(!,/VVN,/M(N):NNJ NLN&
[/J]1 [)’33 T*-"’-9 G."’+6-"’@ H96<’B 26)9# )= 5!,.’3’ ;(-@’?B )= H(,’.(’3& H),* 6’3)96(’3 -.@ ,"3 6’-3).-A*’ 93’ ). [)’33 T*-"’-9& \’,],.+:5!,.’3’
H(,’.(’ -.@ 7’(!.)*)+B T6’33,0II0&
[/K]1 [)’33 T*-"’-9 G."’+6-"’@ H96<’B 26)9# )= 5!,.’3’ ;(-@’?B )= H(,’.(’3& E-"- ;"*-3 )= 4’3)96(’3,U.<,6).?’.",H)(,’"B -.@ U().)?B ,. [)’33
T*-"’-9& \’,],.+:5!,.- U().)?B T6’33,0II/&
[/I]1 [, Y 2,S-.+ ^ [,5!-.+ Z W& [-.@3(-#’ =6-+?’."-",). -.@ 3),* A-6’.’33 3"9@B ,. .)6"! H!--.A’, [)’33 T*-"’-9,;("- U()*)+,(- H,.,(-,/VVL,
/L(M):N/0 N/J&
[MV]1 W9 \ :& 7!’ 3#-",-* #-""’6. -.-*B3,3 )= -+6,(9*"96-* *-.@3(-#’ ,. "!’ *)’33 -6’-& ;("- U()*)+,(- H,.,(-,0IIL,0L(/):00M 00I&
[M0]1 ^-. 8 [,^ 9 X 2& 2)-*3,#6,.(,#*’3,-.@ 3B3"’?3 )= ’()_=9.(",).-* 6’+,).-*,D-",). ,. 5!,.-& 4’3)96(’3 -.@ U.<,6).?’." ,. "!’ ^-.+"D’ \-3,.,
/VVM,0/(O):LJI LKL&
[M/]1 P(F,(-6 7 4,:9## E [ \& 7!’ (966’." -.@ #)"’.",-* )#’6-",).-* 93’ )= 6’?)"’ 3’.3,.+ ") -,@ @’(,3,).3 ). @6)9+!" ’C(’#",).-* (,6(9?3"-.(’3 ,.
;93"6-*,-:; 6’<,’%& ;+6,(9*"96-* HB3"’?,0IIK,LJ:MII NOK&
[MM]1 ‘, H Z,[, 2 5,S-.+ 5 ^,!" #& H"9@B ). ?).,")6,.+ @6)9+!" ,. 5!,.- %,"! PaEGH #6)@9("& ;@<-.(’3 ,. S-"’6 H(,’.(’,/VVL,0O(0):LO O0& [MN]1 T6,(’ : 5& b3,.+ 3#-",-* ()."’C" ,. 3-"’**,"’ @-"- ") ,.=’6 6’+,).-* 3(-*’ ’<-#)"6-.3#,6-",).& GUUU 76-.3-(",).3 ). 2’)3(,’.(’ -.@ 4’?)"’ H’.3,.+, 0IIV,/K:INV INK& [ML]1 5-6*3). 7 8,2,**,’3 4 4,T’66B U P& ; ?’"!)@ ") ?-Q’ 93’ )= "!’6?-* ,.=6-6’@ "’?#’6-"96’ -.@ 8EFG ?’-396’?’."3 ") ,.=’6 396=-(’ 3),* %-"’6 ()."’." -.@ =6-(",).-* <’+’"-",). ()<’6& 4’?)"’ H’.3,.+ 4’<,’%3,0IIN,I:0O0 0JM& [MO]1 W6,’@* P ;,E-<,3 W S& H)96(’3 )= <-6,-",). ,. 6-@,)?’"6,( 396=-(’ "’?#’6-"96’ )<’6 - "-**+6-33 #6-,6,’& 4’?)"’ H’.3,.+ )= U.<,6).?’.",0IIN,NK: 0 0J& 参考文献: [0 ]1 姚春生,张增祥,汪潇& 使用温度植被干旱指数法(7FEG)反演新疆土壤湿度&遥感技术与应用,/VVN,0I(O):NJM c NJK& [/ ]1 冉琼,张增祥,张国平,等& 温度植被干旱指数反演全国土壤湿度的 EUP订正& 中国水土保持科学,/VVL,M(/):M/ c MO& [M ]1 傅伯杰,王军,马克明&黄土丘陵区土地利用对土壤水分的影响&中国自然科学基金,0III,0M(N)://L c //J& [N ]1 王鹏新,S-. Y!’.+ P,.+,龚健雅,等& 基于植被指数和土地表面温度的干旱监测模型& 地球科学进展,/VVM,0K(K):L/J c LMM& [/0]1 王纲胜,夏军,牛存稳& 分布式水文模拟汇流方法及应用&地理研究,/VVN,/M(/):0JL c 0K/& [//]1 郑红星,刘昌明,王中根& 黄河典型流域分布式水文过程模拟& 地理研究,/VVN,/M(N):NNJ c NLN& [/M]1 中国科学院黄土高原综合科学考察队& 黄土高原地区土壤资源及其合理利用& 北京:中国科学技术出版社,0II0& [/N]1 中国科学院黄土高原综合科学考察队& 黄土高原地区资源环境社会经济数据集& 北京:中国经济出版社,0II/& [/L]1 李正国,王仰麟,张小飞& 陕北黄土高原景观破碎化及其土壤裸露效应研究,生态学报,/VVL,/L(M):N/0 c N/J& [/O]1 傅伯杰& 黄土区农业景观空间格局分析&生态学报,0IIL,0L(/):00M c 0/V& [/J]1 燕乃玲,虞孝感& 我国生态功能区划的目标、原则与体系,长江流域资源与环境,/VVM,0/(O):LJI c LKL& [/I]1 齐述华,李贵才,王长耀,等& 利用 PaEGH数据产品进行全国干旱监测的研究& 水科学进展,/VVL,0O(0):LO c O0& LJLN1 00期 1 1 1 李正国1 等:基于植被温度特征的黄土高原地表水分季节变化 1