免费文献传递   相关文献

Dynamic monitoring of drought based on Temperature Vegetation Dryness Index in Guangdong Province

基于温度植被干旱指数的广东省旱情动态监测



全 文 :基于温度植被干旱指数的广东省旱情动态监测
王莺1,王劲松1,姚玉璧1,2,赵福年1
(1.中国气象局兰州干旱气象研究所 甘肃省干旱气候变化与减灾重点实验室 中国气象局干旱气候变化与
减灾重点开放实验室,甘肃 兰州730020;2.甘肃省定西市气象局,甘肃 定西743003)
摘要:广东省饱受季节性干旱的困扰,开展遥感旱情监测工作具有重要的现实意义。利用 MODIS数据构建温度植
被指数特征空间,采用线性拟合的方法提取干湿边,计算广东省2011年温度植被干旱指数(TVDI),以该指数为依
据分析了广东省旱情的时空分布,得到了以下主要结论:1)随着EVI的增加,陆地表面最大温度减小,最大地表温
度和最小地表温度差值呈减小趋势,且地面温度的最大和最小值与EVI呈近似线性关系。EVITS 特征空间的季
节变化明显,随着温度降低,EVITS 特征空间明显萎缩。随着年内温度的变化,EVITS 特征空间的干湿边截距也
发生相应变化,即冬季截距较小,夏季截距较大。2)2011年广东省的冬旱非常严重,春季和秋季也有一定程度的干
旱。其冬旱和春旱的地域分布相似,均呈现自北向南逐渐加重的趋势,且沿海地区重于内陆。秋旱的地区分布特
点与冬春旱相反,大致呈自南向北逐渐加重的趋势。结合广东省历史气象资料对干旱监测结果进行评价,结果表
明遥感监测结果与实际旱情较吻合。3)比较了TVDI与土壤湿度的相关性,结果表明TVDI可以体现土壤湿度状
况,两者呈负相关关系。将TVDI和主要气象因子做了相关性分析,从结果可以看出,TVDI与降水量、温度和相对
湿度间都呈负相关关系。从相关系数来看,TVDI与降水量之间的相关性最高。从气象站尺度来说,TVDI对降水
的变化是敏感的,即连续降水可导致TVDI值下降,旱情得到缓解;持续无降水可使TVDI值增加,旱情加重。
关键词:广东省;温度植被干旱指数(TVDI);MODIS;EVITS 特征空间;旱情动态监测
中图分类号:S503.4;S165+.2  文献标识码:A  文章编号:10045759(2014)02009810
犇犗犐:10.11686/cyxb20140212  
  自有人类社会起,干旱就与其相伴相生,并阻碍着社会的进步与发展。作为一种复杂的自然灾害,干旱在每
一种气候区域内都有可能发生。据世界气象组织的报告,自1963年开始,干旱已经发展成为世界范围内影响最
为广泛的自然灾害[1]。中国是一个干旱灾害频发的国家。据统计,自20世纪90年代以来,我国因旱年均粮食损
失高达278亿kg,因旱年均工业损失超过2000亿元,因旱年均饮水困难人口达2746万人[2]。更需要注意的是,
在近年来全球气候变暖的大背景下,我国年均气温升高,年降水量时空变异增大,极端天气气候事件增多,干旱区
域不断扩大,并有从干旱区向湿润区发展的趋势。
广东省地处我国华南地区,地势北高南低,属于亚热带季风气候区,是中国多雨地区之一。但是在季风、地理
和地质条件的多重因素影响下也饱受季节性干旱的困扰。据《中国气象灾害大典·广东卷》记录,1950年以来,
全省每年受旱面积在十几万至几十万公顷,干旱严重的年份受旱面积可超过100万公顷,到90年代止,全省230
多万公顷的耕地中,约24%的耕地经常受旱。2004年,广东省就发生了半个世纪以来最严重的农业旱灾,受灾面
积7280km2,成灾面积3200km2,直接经济损失超过35亿元[3]。王静爱等[4]分析了1949-2000年中国旱灾动
态变化趋向的区域分异,发现广东省为旱灾增加区。随着社会经济的发展和城市化程度的不断提高,广东省对淡
水的需求量剧增,干旱在威胁着全省农业生产的同时,也越来越严重地影响到社会生活的各个方面。因此进行有
效的干旱监测可以帮助决策部门及时准确地了解旱情的时空分布及其变化特征,采取积极有效的防旱抗旱措施,
将干旱灾害带来的损失降到最低。
传统的干旱监测主要依靠位点尺度的土壤水分含量数据来表征研究区的干旱程度和范围,这种方法的单点
98-107
2014年4月
   草 业 学 报   
   ACTAPRATACULTURAESINICA   
第23卷 第2期
Vol.23,No.2
收稿日期:20130527;改回日期:20130609
基金项目:国家重点基础研究发展计划(2013CB430200,2013CB430206),国家重大科学研究计划(2012CB955903),中国清洁发展机制基金项目
“面向适应的气候灾害风险评估与管理机制研究”和兰州干旱气象研究所博士科研启动项目(KYS2012BSKYO2)资助。
作者简介:王莺(1984),女,甘肃兰州人,助理研究员,博士。Email:wangyn924@gmail.com
测量精度较高,但缺点是空间代表性差。近年来,空间遥感技术以其在时间和空间上可快速获得大面积地物光谱
信息的优势弥补了传统监测方法的不足,对实现干旱的动态监测具有十分重要的意义[5]。自20世纪70年代开
始,国内外对遥感监测土壤水分的方法开展了大量研究,应用较广泛的方法主要有热惯量法、热红外法、距平植被
指数法、作物缺水指数法、植被供水指数法等。热惯量法只适用于裸土或稀疏植被覆盖情况,若植被覆盖度大,热
惯量法对土壤水分的监测就会受到很大限制;归一化差值植被指数(normalizeddifferencevegetationindex,ND
VI)作为水分胁迫指标具有一定的滞后性[68];温度作为水分胁迫指标具有较高敏感性,但遥感获得的地表温度
数据受到土壤背景和植被覆盖的影响。2002年,Sandholt等[9]在简化的植被指数-地表温度特征空间的基础上
提出了温度植被干旱指数(temperaturevegetationdrynessindex,TVDI),该指数既考虑了区域内NDVI的变
化,又考虑了在NDVI相同情况下地表温度的变化,对区域农业旱情和地表土壤水分的监测效果得到了广大学
者的肯定[1016]。但是计算TVDI指数时用到的NDVI数据在植被覆盖度低的区域受到土壤背景的影响,在植被
覆盖度高的区域有饱和现象,且存在对大气影响的纠正不彻底、“最大值合成法”不能保证选择最佳像元等缺
点[17],而基于 MODIS数据的增强型植被指数(enhancedvegetationindex,EVI)在刻画高覆盖度植被方面比
NDVI更加敏感,且能有效抵御大气干扰,因此可以更好地描述区域植被在不同季节的差异[18]。
基于上述原因,尝试采用 MODIS-EVI与陆地表面温度相结合构建EVITS特征空间的方法开展广东省的
干旱动态监测,并结合野外实测土壤湿度数据对模型模拟结果做精度验证,评估该方法对广东省干旱监测工作的
适宜性,以期为广东省的防旱抗旱工作提供科技支撑。
1 材料与方法
1.1 研究区简介
广东省位于109°40′~117°20′E、20°08′~25°32′N之间,海岸线曲折绵长,境内山岭众多,其走向主要为东
北-西南方向,地势呈由北向南伸向海洋的斜坡,山地丘陵约占全境面积的2/3,台地平原占1/3(图1)。广东省
南北跨温带、亚热带和热带,年平均气温19~23℃,1月份平均气温12~16℃,热量资源十分丰富,不少经济作物
可以越冬。广东省濒临南海,年降水量1300~2500mm,居全国前列。但在季风、热带气旋和地形的作用下,该
地区热量和降水量在时空分配上不均匀,常有季节性的春旱、秋旱、冬春连旱和夏秋连旱。广东省的作物主要以
水稻(犗狉狔狕犪狊犪狋犻狏犪)为主,占粮食总产量的80%以上,其抗旱能力低,对水分要求高,极易受旱成灾。
1.2 遥感数据
在研究中使用的遥感数据为美国宇航局/中分辨率成像光谱辐射计(NationalAeronauticsandSpaceAd
ministration/ModerateResolutionImagingSpectroradiometer,NASA/MODIS)陆地产品组按照统一算法开发
的2011年16日最大合成产品 MOD13A2(1000m×1000m)中的EVI数据和8日最大合成产品 MOD11A2
(1000m×1000m)中的LST(LandSurfaceTemperature)数据。以上数据在全球正弦曲线投影SIN(SINusoidal
projection)系统中的编号为h28v06,资料版本为5.0,数据格式为EOS-HDF。该数据可以从http://reverb.
echo.nasa.gov/reverb/网站下载。MOD13A2数据包含 NDVI、EVI、红光、近红外、中红外、蓝光波段反射率及
其他辅助信息,使用最新合成算法减小随观测角度变化和太阳-目标-传感器几何学因素引起的变化,用双向反
射分布函数(BRDF)模式将观测量订正到天顶角。MOD11A2数据包含白天和夜间地表温度,其地表温度是用
31和32通道亮温线性组合的劈窗算法计算获取的,通道亮温值由辐射度与0.1K步长亮温的查找表来确定,发
射率是由 MODIS土地覆盖产品来确定[19]。
在 MRT(modisreprojectiontool)软件中将EOS-HDF格式的 MODIS图像转成TIFF格式,并定义为大
地坐标 WGS84;然后在ArcGIS9.3中用datamanagementtools将图像坐标系统转换为Albers地图投影,并用
celstatistics工具对 MOD11A2数据进行16日最大值合成;最后在spatialanalysttools中用标准的广东省行政
边界掩膜提取研究区EVI和LST数据。
1.3 土壤湿度数据
广东省气象数据和土壤湿度数据分别来源于中国气象科学数据共享服务网(http://cdc.cma.gov.cn/
99第23卷第2期 草业学报2014年
图1 各站点位置图
犉犻犵.1 犔狅犮犪狋犻狅狀狅犳狋犺犲狊狅犻犾犿狅犻狊狋狌狉犲狅犫狊犲狉狏犪狋犻狅狀狊犻狋犲狊犪狀犱犮犾犻犿犪狋犲狊狋犪狋犻狅狀狊
 
home.do)的《中国地面气候资料日值数据集》和《中国农作物生长发育和农田土壤湿度旬值数据集》。该数据集
包含了2011年广东省内36个气象站点的日降水量和温度,以及10个农气站观测的10cm土壤相对湿度(图
1)。
1.4 MODIS-EVI
MODIS-EVI数据是综合处理土壤、大气饱和问题的增强型植被指数。它采用SAVI(soiladjustedvegeta
tionindex)对土壤背景影响进行了校正[20],除了“云处理”和大气校正处理外,还用 ARVI(atmosphericalyre
sistantvegetationindex)对残留气溶胶做了处理。EVI就是通过参数构建的一个同时校正土壤和大气影响的植
被指数。MODIS-EVI的合成算法是根据合成期16d内数据质量状况分4步逐步合成的,即当16d内超过
30%的数据符合质量要求,则根据BRDF把不同视角换算为星下点反射值,分别计算植被指数,然后采用限定视
角内最大值原理合成;若16d内小于30%且多于2d的数据符合质量要求,选其中视角最小的2个植被指数,取
两者中最大值合成;若只有1d无云,则直接使用这天的数据计算植被指数;若观测期内均有云,则取合成期所有
植被指数的最大值[21]。其计算公式为:
犈犞犐=2.5(ρNIR-ρRed)/(ρNIR+犮1ρRed-犮2ρBlue+犔) (1)
式(1)中,ρ为经大气校正的各波段(NIR、Red和Blue)的反射率,犔=1为土壤调节参数,犮1 和犮2 是描述通过Blue
波段修正大气对Red波段影响的参数,分别是6.0和7.5。
该数据在算法方面的改进为遥感定量研究奠定了良好的基础。王正兴等[1718]的研究发现,在南亚热带和热
带地区,EVI与地表温度的线性关系明显高于NDVI,且对植被不同季相的变化反应更加灵敏。
1.5 TVDI计算方法
很多研究人员发现陆地表面温度(TS)和植被指数之间存在明显的负相关关系,这种关系与土壤湿度密切相
关[2223]。Carlson等[24]和Price[25]发现若研究区植被覆盖包含裸土到全覆盖,土壤湿度包含极干旱到极湿润,则
以遥感资料获得的NDVI和TS数据的散点图呈三角形,而 Moran等[26]发现其表现为梯形。这就是 NDVITS
特征空间(图2)。Sandholt等[9]在研究土壤湿度时发现NDVITS 特征空间中有很多等值线,于是提出了TVDI
指数,用以估算土壤表层水分状况,从而进行遥感干旱监测。
001 ACTAPRATACULTURAESINICA(2014) Vol.23,No.2
  由于1.4中提到的EVI的优势,所以应用EVI来
图2 犖犇犞犐犜犛特征空间
犉犻犵.2 犐犾狌狊狋狉犪狋犻狅狀狅犳狊犻犿狆犾犻犳犻犲犱犖犇犞犐犜犛狊狆犪犮犲
 
替代NDVI建立EVITS 特征空间,进而构建 TVDI
指数。TVDI的计算见公式(2):
犜犞犇犐=(犜犛-犜犛min)/(犜犛max-犜犛min) (2)
式中,犜犞犇犐的取值在0~1之间,湿边的犜犞犇犐最小,
为0,土壤含水量接近田间持水量,干边的犜犞犇犐最
大,为1,土壤含水量接近萎蔫点;犜犛 是任意像元的地
表温度;犜犛min是某一EVI对应的最低地表温度,对应
EVITS空间的湿边;犜犛max是某一EVI值对应的最高
地表温度,对应EVITS 空间的干边。犜犛min和犜犛max的
拟合方程如下:
犜犛min=犪1+犫1×犈犞犐 (3)
犜犛max=犪2+犫2×犈犞犐 (4)
式中,犪1、犫1 是湿边拟合方程的系数,犪2、犫2 是干边拟合方程的系数。
2 结果与分析
2.1 特征空间参数
提取广东省相同EVI下不同像元对应的所有温度中的最大和最小陆地表面温度,然后通过EVI和LST二
维散点图建立每16d1次的EVITS特征空间图(图3)。2011年共获得23张EVITS 特征空间图,图3为其中
的4张。从图中可以看出,EVITS散点图有相似的分布特征,即随着EVI的增加,陆地表面最大温度减小,最大
图3 犈犞犐犜犛特征空间
犉犻犵.3 犈犞犐犜犛狊狆犪犮犲 
101第23卷第2期 草业学报2014年
地表温度和最小地表温度差值呈减小趋势,且地面温
度的最大和最小值与EVI呈近似线性关系。另外,
EVITS特征空间的季节变化也非常明显,随着温度降
低,EVITS 特征空间明显萎缩。和图2相比,若将湿
边描述成与EVI轴平行的直线会给结果带来较大误
差,所以有必要对湿边进行线性拟合。
按照TVDI原理,当植被指数增加时,陆地表面最
大温度逐渐降低。这个假设基于植被指数与地表植被
覆盖度呈线性关系这一前提。但实际情况并非如
此[27]。设EVI的精度为0.01,提取相应EVI在特征
空间中的最大和最小陆地表面温度,并对特征空间进
行最优目视判读,回归拟合获得每个特征空间的干湿
边方程(表1),线性拟合结果经狋检验发现均达到0.05
的显著性水平。从表中可以看出,EVITS特征空间干
边拟合的整体效果很好,平均拟合系数犚2>0.87;湿
边的拟合效果也比较好,平均拟合系数犚2>0.55。但
是在某些时段干湿边拟合精度仍比较低,造成这一现
象的主要原因是 MODIS产品存在的噪声会影响线性
拟合精度。
犪1 和犪2 是干湿边拟合方程的常数项,在特征空
间中是干湿边在犢 轴(地表温度)的截距。它代表了
裸土像元在水分充足(犪1)和水分匮乏(犪2)时的地表温
度值。从表1中可以看到,随着年内温度的变化,
EVITS特征空间的干湿边截距也发生相应变化,即冬
季截距较小,夏季截距较大。犫1 和犫2 为干湿边拟合方
程的斜率,其变化可以归因于蒸散[28]、冠层传导度[29]
和土壤湿度[30]。在多重因素的影响下,犫1 和犫2 随时
间变化的规律性不明显。
表1 犈犞犐犜犛特征空间干湿边方程
犜犪犫犾犲1 犜犺犲犲狇狌犪狋犻狅狀狊狅犳狋犺犲犱狉狔犲犱犵犲犪狀犱
狑犲狋犲犱犵犲犳狅狉犈犞犐犜犛狊狆犪犮犲
日序数
Julianday
干边方程
Dryedge
犚2 湿边方程
Wetedge
犚2
1 狔=-9.33狓+24.68 0.92 狔=22.78狓+0.36 0.94
17 狔=-12.46狓+24.88 0.87 狔=5.31狓+4.67 0.40
33 狔=-21.03狓+32.99 0.90 狔=3.97狓+14.79 0.51
49 狔=-21.13狓+34.95 0.91 狔=17.41狓+11.29 0.74
65 狔=-16.86狓+36.16 0.86 狔=-4.36狓+17.98 0.35
81 狔=-16.36狓+34.55 0.95 狔=24.28狓+9.75 0.86
97 狔=-30.24狓+43.54 0.96 狔=-7.41狓+25.31 0.72
113 狔=-19.51狓+40.35 0.92 狔=2.95狓+20.91 0.43
129 狔=-14.22狓+38.96 0.92 狔=3.11狓+20.95 0.34
145 狔=-15.54狓+41.20 0.93 狔=-4.64狓+27.04 0.51
161 狔=-4.02狓+34.21 0.50 狔=12.05狓+15.43 0.76
177 狔=-18.66狓+42.76 0.91 狔=3.63狓+22.24 0.40
193 狔=-13.78狓+39.24 0.93 狔=2.34狓+24.19 0.42
209 狔=-15.67狓+40.72 0.95 狔=3.29狓+23.28 0.34
225 狔=-16.10狓+41.19 0.95 狔=3.23狓+22.78 0.44
241 狔=-11.75狓+38.17 0.97 狔=7.77狓+20.74 0.61
257 狔=-18.24狓+39.85 0.97 狔=4.22狓+20.71 0.38
273 狔=-13.54狓+35.26 0.93 狔=6.74狓+17.38 0.58
289 狔=-13.81狓+34.45 0.96 狔=5.28狓+19.97 0.42
305 狔=-11.26狓+33.53 0.92 狔=6.05狓+18.26 0.82
321 狔=-5.65狓+28.80 0.77 狔=4.96狓+17.96 0.42
337 狔=-2.73狓+23.74 0.43 狔=3.06狓+11.50 0.58
353 狔=-4.18狓+23.75 0.87 狔=4.32狓+10.88 0.60
2.2 旱情等级分布
根据表1获得的干湿边方程和公式(2),分别计算各时段的TVDI值。以TVDI值作为旱情分级指标,将旱
情分为5级,分别是湿润(0~0.2)、正常(0.2~0.4)、轻旱(0.4~0.6)、中旱(0.6~0.8)和重旱(0.8~1.0)。由此
可得2011年23个时段16d的广东省旱情分布图(图4)。从图中可以看出,2011年广东省的冬旱非常严重,春
季和秋季也有一定程度的干旱。广东省冬旱和春旱的地区分布相似,均呈现自北向南逐渐加重的趋势,且沿海地
区重于内陆。秋旱的地区分布特点与冬春旱相反,大致呈自南向北逐渐加重的趋势。具体来看,1月,重旱区主
要位于广东省西南部的徐闻县、雷州市、湛江市、遂溪县、廉江市、吴川市、化州市、茂名市、电白县、高州市和阳西
县,还有阳东县和信宜市的大部分地区,以及广东省东部的海丰县、汕尾市、陆丰市、陆河县、惠来县、普宁市、潮阳
区、揭西县、揭阳市、揭东县、汕头市和澄海区。2月,广东西南部的重旱区仍然存在,同时广东北部的龙川县、兴
宁市和梅州市也发展成了重旱区。3月,重旱区主要分布在广东的东南沿海。4月,重旱区主要位于广东的西南
部和东部。5月中上旬,广东西南部的重旱区基本消失,而到了5月下旬,广东西南部化州市、廉江市、茂名市、电
白县、吴川市、湛江市、遂溪县、雷州市、徐闻县又出现重旱。6月,广东西南部的重旱范围向南推移。7月,广东西
南部的重旱范围进一步缩小。8月,广东西南部的重旱基本得到缓解,而广东中部和北部出现局部的重旱。9月
201 ACTAPRATACULTURAESINICA(2014) Vol.23,No.2
和10月,广东西南部、中偏东部和东部出现局部重旱。11月,广东西南部又出现大面积重旱。12月,重旱分布范
围在11月的基础上有了明显扩大。
《中国气象灾害大典·广东卷》中记载,广东省的重春旱区位于雷州半岛、周江口一带及东南部沿海地区,而
雷州半岛南端春旱是固有的灾害现象;重秋旱区主要位于北部偏北及潮汕一带沿海地区;重冬春连旱区主要位于
西南部偏西地区和珠江口附近沿海地区。本文的监测结果与上述规律基本吻合。而《中国气象灾害年鉴(2012)》
中记载,广东省自2010年10月1日至2011年4月28日平均降水量为225.6mm,比常年同期偏少61%,为历史
同期最少,持续偏少的降水导致广东出现了历史罕见的秋冬春连旱,而旱情从3月开始从粤东到粤西蔓延加剧。
本文的研究结果较准确的反映了上述记录。因此,可以应用 MODIS数据进行广东省旱情的动态监测,为防灾减
灾和灾后救援提供及时准确的基础资料。
图4 2011年广东省旱情遥感监测分布
犉犻犵.4 犜犺犲犱狉狅狌犵犺狋犱犻狊狋狉犻犫狌狋犻狅狀犿狅狀犻狋狅狉犲犱犫狔狉犲犿狅狋犲狊犲狀狊犻狀犵犻狀犌狌犪狀犵犱狅狀犵犻狀2011
 
2.3 TVDI与土壤湿度
图5 犜犞犇犐与土壤湿度
犉犻犵.5 犜犺犲狉犲犾犪狋犻狅狀狊犺犻狆犫犲狋狑犲犲狀犜犞犇犐犪狀犱狊狅犻犾犿狅犻狊狋狌狉犲
 
提取和土壤湿度数据采集时间相近的相应像元
的TVDI数据,分析TVDI值与站点0~10cm土壤
湿度间的相关性(图5)。从图中可以看出,TVDI和
土壤湿度间呈负相关关系,犚2 为0.16,经过狋检验
发现线性回归方程达到0.01的显著性水平。即
TVDI值越高,土壤湿度越低。这一结论满足TV
DI值越接近1,土壤水分越少的原理。这进一步说
明TVDI可以反映广东省地表土壤水分状况,用该
指标做广东省旱情动态监测是可行的,具有应用价
值。
301第23卷第2期 草业学报2014年
2.4 TVDI与气象因子
图6 犜犞犇犐与气象因子
犉犻犵.6 犜犺犲狉犲犾犪狋犻狅狀狊犺犻狆犫犲狋狑犲犲狀犜犞犇犐犪狀犱
犿犲狋犲狅狉狅犾狅犵犻犮犪犾犳犪犮狋狅狉狊 
在本研究中主要分析了降水量、平均温度和平均
相对湿度与相应像元的 TVDI值之间的相关性(图
6)。从图中可以看出,TVDI与降水量、温度和相对湿
度都呈负相关关系。说明随着降水量减少、温度降
低和平均相对湿度减小,TVDI值增加,广东省干旱程
度增加。经狋检验表明,方程达到0.05的显著性水
平。从相关系数来看,TVDI与降水量之间的相关性
最高,犚2 为0.35,与平均温度间的相关性最低,犚2 为
0.21。因此,广东省旱情发展的主要影响因子是降水
量。
在气象站尺度上进一步分析了TVDI和降水量的
时间变化趋势。选择了 2 个站点,分 别 是 罗 定
(111.57°E,22.77°N)和惠阳(114.42°E,23.08°N)。
这两个站所在像元在23幅图像中均有有效数据,降水
量有较明显差异,且分别位于广东省的西南和东北区,
具有地域代表性。从图7中可以看出,TVDI对降水
的变化是敏感的,在连续降水之后,土壤湿度增加,植
物受胁迫的情况得到缓解,TVDI值降低;当降水减
少,土壤湿度减小,植被生长受到胁迫,TVDI值增加。
从相关系数来看,罗定站的 TVDI和降水量间的犚2
为0.0396,惠阳站的 TVDI和降水量间的 犚2 为
0.0118,均未通过显著性检验。出现这种情况的主要
原因是空间尺度的不同。TVDI的空间分辨率为1
km,代表的是1km2 的地表干湿状况,而降水量数据
只代表气象站一个点的降水状况。但从图7中仍可以
看出,TVDI指数对降水是敏感的。
图7 降水量与犜犞犇犐在站点尺度随时间变化关系图
犉犻犵.7 犜狉犲狀犱狊狅犳狉犪犻狀犳犪犾犪狀犱犜犞犇犐犳狉狅犿狋犺犲狊狋犪狋犻狅狀狊犮犪犾犲 
401 ACTAPRATACULTURAESINICA(2014) Vol.23,No.2
3 结论
通过建立EVITS特征空间算法,获得近似三角形的温度植被指数特征空间,然后采用线性拟合的方法提取
了干湿边,构建了广东省2011年的温度植被干旱指数TVDI,得到了以下主要结论。
1)随着EVI的增加,陆地表面最大温度减小,最大地表温度和最小地表温度差值呈减小趋势,且地面温度的
最大和最小值与EVI呈近似线性关系。EVITS 特征空间的季节变化也非常明显,随着温度降低,EVITS 特征
空间明显萎缩。随着年内温度的变化,EVITS特征空间的干湿边截距也发生相应变化,即冬季截距较小,夏季截
距较大。
2)以TVDI值作为旱情分级指标,建立广东省2011年旱情遥感监测图,从图中可以看出2011年广东省的冬
旱非常严重,春季和秋季也有一定程度的干旱。冬旱和春旱的地区分布相似,均呈现自北向南逐渐加重的趋势,
且沿海地区重于内陆。秋旱的地区分布特点与冬春旱相反,大致呈自南向北逐渐加重的趋势。结合广东省历史
气象资料对干旱监测结果进行评价,结果表明遥感监测结果与实际旱情较吻合,因此应用 MODIS数据构建
TVDI指数的方法可以用于广东省旱情动态监测。
3)比较TVDI与土壤湿度的相关性的结果表明,TVDI可以体现土壤湿度状况,两者呈负相关关系。从相关
分析结果可以看出,TVDI与降水量、温度和相对湿度间都呈负相关关系。从相关系数来看,TVDI与降水量之
间的相关性最高。在气象站尺度上TVDI对降水的变化最敏感。
从以上结论可以看出,TVDI是一种适用于大尺度的简单有效且物理意义明确的旱情动态监测方法,具有快
速、及时、宏观等优势。而利用EVI代替NDVI建立EVITS 特征空间,既改善了NDVI的植被饱和、对大气纠
正不彻底等问题,也降低了土壤背景对植被指数的影响,使其更适宜于广东省的旱情监测与评估。因为定性反映
土壤湿度情况的干边和湿边拟合方程是由EVITS特征空间数据得到的,不需要其他辅助数据,所以在不同时间
和不同区域不具有可比性,且对边界的拟合具有一定的主观性。另外,在复杂地表和多变的气候环境条件下,遥
感影像的部分像元数值质量较差,从而影响特征空间干湿边线性拟合的精度,在今后的工作中也要注意加强对无
效像元分析剔除等方法的研究。
参考文献:
[1] 李华朋,张树清,高自强,等.MODIS植被指数监测农业干旱的适宜性评价[J].光谱学与光谱分析,2013,33(3):756
761.
[2] 亚行技援中国干旱管理战略研究课题组.中国干旱灾害风险管理战略研究[M].北京:中国水利水电出版社,2011:12.
[3] 高懋芳,张虹鸥,秦晓敏,等.广东省农业旱灾遥感监测[J].国土资源遥感,2008,(3):9499.
[4] 王静爱,孙恒,徐伟,等.近50年中国旱灾的时空变化[J].自然灾害学报,2002,11(2):16.
[5] 李慧,王鹏新.基于TerraMODIS和NOAAAVHRR数据的条件植被温度指数干旱监测及其对比分析[J].干旱区资源与
环境,2013,27(3):6166.
[6] 花立民.玛曲草原植被NDVI与气候和载畜量变化的关系分析[J].草业学报,2012,21(4):224235.
[7] 周伟,王倩,章超越,等.黑河上游草地NDVI时空变化规律及其对气候因子的响应分析[J].草业学报,2013,22(1):138
147.
[8] 王浩,李文龙,杜国祯,等.基于3S技术的甘南草地覆盖度动态变化研究[J].草业学报,2012,21(3):2637.
[9] SandholtI,RasmussenK,AndersenJ.Asimpleinterpretationofthesurfacetemperature/vegetationindexspaceforassess
mentofsurfacemoisturestatus[J].RemoteSensingofEnvironment,2002,79(2):213224.
[10] 冯蜀青,殷青军,肖建设,等.基于温度植被旱情指数的青海高寒区干旱遥感动态监测研究[J].干旱地区农业研究,2006,
24(5):141145.
[11] PatelNR,AnapashshaR,KumarS,犲狋犪犾.AssessingpotentialofMODISderivedtemperature/vegetationconditionindex
(TVDI)toinfersoilmoisturestatus[J].InternationalJournalofRemoteSensing,2009,30(1):2339.
[12] 齐述华,王长耀,牛铮.利用温度植被旱情指数 (TVDI)进行全国旱情监测研究[J].遥感学报,2003,7(5):420427.
[13] GaoZ,GaoW,ChangNB.Integratingtemperaturevegetationdrynessindex(TVDI)andregionalwaterstressindex(RW
501第23卷第2期 草业学报2014年
SI)fordroughtassessmentwiththeaidofLANDSATTM/ETM+images[J].InternationalJournalofAppliedEarthObser
vationandGeoinformation,2011,13(3):495503.
[14] 吴孟泉,崔伟宏,李景刚.温度植被干旱指数 (TVDI)在复杂山区干旱监测的应用研究[J].干旱区地理,2007,30(1):30
35.
[15] 赵莉荣,武伟,刘洪斌,等.基于温度植被干旱指数法的农业干旱研究[J].西南师范大学学报(自然科学版),2009,34(2):
8084.
[16] WanZ,WangP,LiX.UsingMODISlandsurfacetemperatureandnormalizeddifferencevegetationindexproductsformoni
toringdroughtinthesouthernGreatPlains[J].USA,InternationalJournalofRemoteSensing,2004,25(1):6172.
[17] 王正兴,刘闯,HueteA.植被指数研究进展:从AVHRRNDVI到 MODISEVI[J].生态学报,2003,23(5):979987.
[18] 王正兴,刘闯,陈文波,等.MODIS增强型植被指数EVI与 NDVI初步比较[J].武汉大学学报(信息科学版),2006,
31(5):407410.
[19] WanZ.MODISlandsurfacetemperatureproductsusers’guide[Z].InstituteforComputationalEarthSystemScience,Uni
versityofCalifornia,SantaBarbara,CA.http://www.icess.ucsb.edu/modis/LstUsrGuide/usrguide.html,2006.
[20] HueteAR.Asoiladjustedvegetationindex(SAVI)[J].RemoteSensingofEnvironment,1988,(25):295309.
[21] 李文梅,覃志豪,李文娟,等.MODISNDVI与 MODISEVI的比较分析[J].遥感信息,2010,6:7378.
[22] 李云鹏,司瑶冰,刘朋涛,等.基于空间信息的内蒙古农业干旱监测研究[J].干旱区资源与环境,2011,25(11):125131.
[23] 于敏,程明虎.基于NDVITS 特征空间的黑龙江省干旱监测[J].应用气象学报,2010,21(2):221228.
[24] CarlsonTN,GiliesRR,PerryEM.AmethodtomakeuseofthermalinfraredtemperatureandNDVImeasurementstoin
fersurfacesoilwatercontentandfractionalvegetationcover[J].RemoteSensingReviews,1994,9(12):4559.
[25] PriceJC.Usingspatialcontextinsatelitedatatoinferregionalscaleevapotranspiration[J].GeoscienceandRemoteSens
ing,IEEETransactionson,1990,28(5):940948.
[26] MoranMS,ClarkeTR,InoueY,犲狋犪犾.Estimatingcropwaterdeficitusingtherelationbetweensurfaceairtemperatureand
spectralvegetationindex[J].RemoteSensingofEnvironment,1994,49(3):246263.
[27] 叶琦,赵萍,孙静.基于 MODIS/NDVI与EVI的皖江流域植被覆盖比较分析[J].长江流域资源与环境,2012,21(3):
361368.
[28] HopeAS,PetzoldDE,GowardSN,犲狋犪犾.Simulatingcanopyreflectanceandthermalinfraredemissionsforestimatinge
vapotranspiration[J].WaterResourcesBuletin,1987,22(6):10111019.
[29] NemaniRR,RunningS W.Estimationofregionalsurfaceresistancetoevapotranspirationfrom NDVIandthermalIR
AVHRRdata[J].JournalofAppliedMeteorology,1989,28(4):276284.
[30] CarlsonT.Anoverviewofthe“trianglemethod”forestimatingsurfaceevapotranspirationandsoilmoisturefromsatelite
imagery[J].Sensors,2007,7(8):16121629.
601 ACTAPRATACULTURAESINICA(2014) Vol.23,No.2
犇狔狀犪犿犻犮犿狅狀犻狋狅狉犻狀犵狅犳犱狉狅狌犵犺狋犫犪狊犲犱狅狀犜犲犿狆犲狉犪狋狌狉犲犞犲犵犲狋犪狋犻狅狀犇狉狔狀犲狊狊犐狀犱犲狓犻狀犌狌犪狀犵犱狅狀犵犘狉狅狏犻狀犮犲
WANGYing1,WANGJingsong1,YAOYubi1,2,ZHAOFunian1
(1.KeyLaboratoryofAridClimaticChangeandReducingDisasterofGansuProvince,KeyOpen
LaboratoryofAridChangeandDisasterReductionofChinaMeteorologicalAdministration,
InstituteofAridMeteorology,ChinaMeteorologicalAdministration,Lanzhou
730020,China;2.MeteorologicalBureauofDingxiCity,
Dingxi743003,China)
犃犫狊狋狉犪犮狋:GuangdongProvincehasbeensufferingfromseasonaldroughtandthedynamicmonitoringofdrought
byremotesensingdatacouldbeanimportanttool.TheEVITSspacewasconstructedbyModerateResolution
ImagingSpectroradiometer(MODIS)data,thedryandwetedgeswereextractedwithalinearfittingmethod
andthentheTemperatureVegetationDrynessIndex(TVDI)wascalculatedforGuangdongProvincein2011
andtheTemporalspatialdistributionofdroughtwasanalyzed.1)WithanincreaseoftheEnhancedVegetation
Index(EVI),themaximumlandsurfacetemperaturedecreasedandthedifferencevaluebetweenthemaximum
andminimumlandsurfacetemperaturealsodecreased.Themaximumandminimumlandsurfacetemperature
bothhadanapproximatelylinearrelationtotheEVI.TheseasonalvariationofEVITSspacewasobvious.
Whenthetemperaturedecreased,theEVITSspaceshrankandwithachangeoftemperature,theinterceptof
thedryandwetedgeschangedcorrespondingly,i.e.theinterceptwassmalerinwinterandbiggerinsummer.
2)Winterdroughtwasveryserious,whilespringandautumndroughtwerealsoserioustosomedegreein
GuangdongProvincein2011.Thedroughtinwinterandspringhadsimilardistributionareas,andthegrade
weakenedgradualyfromsouthtonorth.Thedroughtincoastalareaswasmoreseriousthanininlandareas.
Thegeographicaldistributionofautumndroughtapproximatelyincreasedfromsouthtonorth,whichwasthe
oppositetowinterandspringdroughts.AccordingtotheweatherdatainGuangdongProvince,theresultsof
thedroughtmonitoringbyremotesensingwereinconformitywithreality.3)Therewasasignificantnegative
correlationbetweenTVDIandsoilmoistureandtheresultsshowedthatthesoilmoisturecanbesuccessfuly
reflectedbytheTVDIdata.TheTVDIwasnegativelycorrelatedwithprecipitation,temperatureandmeanrel
ativehumidity.TVDIhadahighcorrelationcoefficientwithprecipitation.TVDIwassensitivetoprecipitation
onthestationscale,whichisasimilartrendtothechangeofTVDIandprecipitationinthetimeseries.Contin
uousprecipitationforacertaintimewasfolowedbyadeclineintheTVDIwhichindicatedslightdroughtonly
andviceversa.
犓犲狔狑狅狉犱狊:GuangdongProvince;TemperatureVegetationDrynessIndex(TVDI);ModerateResolutionIma
gingSpectroradiometer(MODIS);EVITSspace;dynamicmonitoringofdrought
701第23卷第2期 草业学报2014年