劉星含，張佳華，2，許曉宏

（1.長江大學(xué)地球科學(xué)學(xué)院，武漢430100；2.中國科學(xué)院遙感與數(shù)字地球研究所數(shù)字地球?qū)嶒?yàn)室，北京100094）

MODIS－TVDI指數(shù)監(jiān)測新疆干旱動(dòng)態(tài)

劉星含1，張佳華1，2，許曉宏1

（1.長江大學(xué)地球科學(xué)學(xué)院，武漢430100；2.中國科學(xué)院遙感與數(shù)字地球研究所數(shù)字地球?qū)嶒?yàn)室，北京100094）

針對近年頻發(fā)的干旱情況不能準(zhǔn)確及時(shí)監(jiān)測評估的問題，該文以新疆為研究區(qū)域，基于溫度植被干旱指數(shù)方法，利用2007年到2012年3月～8月MODIS合成產(chǎn)品數(shù)據(jù)獲取歸一化植被指數(shù)和陸地地表溫度，構(gòu)建LST－NDVI特征空間，得到全區(qū)的溫度植被干旱指數(shù)和旱情等級空間分布圖，分析了新疆干旱變化趨勢，驗(yàn)證了溫度植被干旱指數(shù)和降水因子的關(guān)系。結(jié)果表明：2007年～2012年新疆的干旱面積逐年趨于平穩(wěn)，空間上表現(xiàn)為南疆旱情高于北疆，春季旱情高于夏季，降水量是影響溫度植被干旱指數(shù)的重要因子。該研究為政府部門對新疆旱情嚴(yán)重地區(qū)治理提供了有效數(shù)據(jù)保證。

歸一化植被指數(shù)；陸地表面溫度；溫度植被干旱指數(shù)；干旱監(jiān)測；干旱評估

1 引 言

干旱與其他自然災(zāi)害相比，出現(xiàn)的頻率最高，持續(xù)時(shí)間最長，影響范圍最大，除了危害作物生長，造成作物減產(chǎn)，同時(shí)還危害居民生活，影響工業(yè)生產(chǎn)和其他社會(huì)經(jīng)濟(jì)活動(dòng)。

基于地面觀測的干旱監(jiān)測方法主要是基于地面站點(diǎn)的土壤墑情數(shù)據(jù)，其準(zhǔn)確性、代表性和完整性都存在局限，遙感技術(shù)的優(yōu)勢在于能夠利用傳感器獲取及時(shí)、客觀的大范圍綜合地表信息，同時(shí)監(jiān)測陸地表面溫度和植被生長狀態(tài)，這使遙感技術(shù)成為了區(qū)域旱情監(jiān)測的重要手段。

國外學(xué)者利用遙感進(jìn)行土壤水分監(jiān)測研究開始于20世紀(jì)60年代末。20世紀(jì)70年代，Bowers和Hanks等［1］對土壤光譜反射特性進(jìn)行了研究；Watsom等［2］首次提出了用地表溫度日較差推算熱慣量的簡單模式。Carlson T N［3］在估算土壤有效水分和熱慣量研究中，較早開展了利用NOAA／AVHRR遙感資料進(jìn)行大面積作物旱情監(jiān)測的嘗試。1990年之后，國外在土壤水分遙感定量反演方面得到了新的發(fā)展，并逐步實(shí)現(xiàn)了實(shí)用化系統(tǒng)應(yīng)用。從Moran對歸一化植被指數(shù)（Normalized Difference Vegetation Index，NDVI）和陸地地表溫度（Land Surface Temperature，LST）的相關(guān)性研究開始，多項(xiàng)研究表明NDVI和LST之間存在著顯著的負(fù)相關(guān)性［4－5］。Price和Carlson等［3，6］研究發(fā)現(xiàn)，若研究區(qū)域包含地表覆蓋植被類型從裸土到密閉植被冠層，土壤濕度由干旱到濕潤，則該區(qū)域內(nèi)每個(gè)像元的植被指數(shù)和地表溫度組成的散點(diǎn)圖呈現(xiàn)為梯形，在此基礎(chǔ)上Sandholt等［7］提出了溫度植被干旱指數(shù)（Temperature Vegetation Dryness Index，TVDI）的概念。

國內(nèi)學(xué)者進(jìn)入90年代后，在利用NOAA／AVHRR資料進(jìn)行土壤水分或干旱的宏觀監(jiān)測研究工作取得了很大進(jìn)步。肖乾廣等［8］利用NOAA／AVHRR數(shù)據(jù)對我國干旱情況進(jìn)行了長期監(jiān)測；查勇［9］利用衛(wèi)星遙感圖像創(chuàng)建了統(tǒng)一的草地植被變化反演模型；柳欽火等［10］利用地表溫度和植被指數(shù)，將全國分為6個(gè)區(qū)進(jìn)行全國農(nóng)業(yè)旱情監(jiān)測；張順謙等［11］利用NOAA／AVHRR資料，結(jié)合四川省農(nóng)業(yè)氣象臺測定的土壤濕度數(shù)據(jù)，通過對溫度植被旱情指數(shù)計(jì)算方法的改進(jìn)和旱情分級標(biāo)準(zhǔn)的研究，建立了四川伏旱遙感監(jiān)測模型。田慶久和杜靈通［12］構(gòu)建了基于多遠(yuǎn)空間信息的干旱監(jiān)測模型。

新疆是干旱災(zāi)害頻發(fā)的地區(qū)，由于特殊的地理和氣候環(huán)境，水資源時(shí)空分布不均勻，季節(jié)性供需矛盾大，社會(huì)經(jīng)濟(jì)發(fā)展和不合理的人類運(yùn)動(dòng)等因素［13］，造成了新疆干旱常年發(fā)生。近年來，新疆地區(qū)頻繁的旱情讓研究者對該地區(qū)的干旱情況進(jìn)行了研究。姚春生等［14］使用溫度植被干旱指數(shù)法完成了新疆土壤濕度的反演。包安明等［15］利用MODIS數(shù)據(jù)對塔里木河下游的植被生長狀況進(jìn)行了監(jiān)測，完成了新疆生態(tài)環(huán)境建設(shè)的研究。

本文基于MODIS－TVDI方法開展2007年到2012年的新疆干旱監(jiān)測研究，構(gòu)建多年的LSTNDVI特征空間，得到新疆干旱情況分布圖。分析了2007年～2012年3月～8月新疆干旱變化情況，研究結(jié)果為政府部門對新疆旱情嚴(yán)重地區(qū)治理提供了有效的數(shù)據(jù)保證。

2 研究區(qū)域概況

新疆維吾爾自治區(qū)位于中國西北邊陲，總面積為1.6649×106km2。由于遠(yuǎn)離海洋、深居內(nèi)陸，四周有高山阻隔，海洋氣流不易到達(dá)，所以形成了典型的溫帶大陸性干旱氣候。區(qū)內(nèi)降水量山區(qū)多于盆地、西部多于東部、北部多于南部，氣候干燥，天山以南的南疆氣溫高于天山以北的北疆。

新疆的土地利用類型分為農(nóng)用地（37.85%）、建設(shè)用地（0.72%）和未利用地（61.43%）。農(nóng)用地分為耕地、園地、林地、牧草地和其他農(nóng)業(yè)用地，其中牧草地所占比例最大，達(dá)到30.82%；建筑用地分為居民點(diǎn)及工礦用地、交通用地和水利設(shè)施用地，所占比例依次為0.58%、0.03%、0.11%［16］。新疆的耕地主要分布在伊犁地區(qū)、昌吉自治州、塔城地區(qū)等，園地主要分布在巴州、阿克蘇、克州、喀什地區(qū)等，林地主要分布在伊犁地區(qū)、和田地區(qū)、巴州、阿勒泰地區(qū)，牧草地主要分布在伊利地區(qū)、阿勒泰地區(qū)。

通過新疆1961年～2005年88個(gè)氣象站多年的平均氣溫和降水量［17］，可以得到45年中的平均氣溫平均值為7.8℃，其中吐魯番最高（14.5℃），天山大西溝最低（－5.1℃）。平均降水量為154.6mm，其中天山小渠子最大（543.3mm），托克遜最?。?.8mm）。近年來，在新疆部分地區(qū)和小型灌溉區(qū)等易旱區(qū)域連年發(fā)生持續(xù)性干旱，嚴(yán)重影響了農(nóng)牧業(yè)生產(chǎn)，給新疆的社會(huì)經(jīng)濟(jì)持續(xù)發(fā)展造成了嚴(yán)重危害。

3 數(shù)據(jù)與方法

3.1 數(shù)據(jù)處理

本文收集了研究區(qū)2007年～2012年3月～8月的MODIS11A2和MODIS13A3。MODIS11A2是8天合成空間分辨率為1km的陸地表面溫度產(chǎn)品，包括白天地表溫度、夜間地表溫度、31和32波段通道發(fā)射率等資料。該產(chǎn)品中的地表溫度是通過建立31、32通道亮溫線性組合的劈窗算法計(jì)算得到。MODIS13A3是月合成的空間分辨率為1km的植被指數(shù)產(chǎn)品。該產(chǎn)品包括NDVI、EVI、紅光、近紅外、中紅外、藍(lán)光等波段反射率以及其他輔助信息。MODIS13A3使用新的合成算法減小了隨觀測角度的變化和太陽－目標(biāo)－傳感器幾何學(xué)因素引起的變化。在生成植被指數(shù)格點(diǎn)資料時(shí)將應(yīng)用分子散射、臭氧吸收、氣溶膠訂正算法、用BRDF模式將觀測量訂正到天頂角。將數(shù)據(jù)拼接、投影轉(zhuǎn)換，利用波段計(jì)算和裁剪完成預(yù)處理。本文遙感數(shù)據(jù)來源于美國NASA網(wǎng)站（http：／／ladsweb.nascom.nasa.gov）。

3.2 溫度植被干旱指數(shù)

Sandholt［7］等在研究中發(fā)現(xiàn)，LST－NDVI特征空間中有很多等值線，于是提出溫度植被干旱指數(shù)（TVDI）的概念。TVDI由植被指數(shù)和地表溫度計(jì)算得到，其定義為：

TVDI＝（T－Tmin）／（Tmax－Tmin） （1）

其中T為任意像元的地表溫度；Tmin為NDVI對應(yīng)的最低地表溫度，即濕邊，Tmin＝a＋b× NDVI；Tmax為NDVI對應(yīng)的最高地表溫度，即干邊，Tmax＝c＋d×NDVI；a、b、c、d是干濕邊擬合方程的系數(shù)。

TVDI的取值在0～1之間，TVDI值越大，表明該區(qū)域干旱越嚴(yán)重。在濕邊上TVDI＝0，在干邊上TVDI＝1，T越接近干邊，TVDI越大，土壤干旱情況越嚴(yán)重；相反，T越接近濕邊，TVDI越小，土壤含水量越高。對參數(shù)的計(jì)算要求研究區(qū)域范圍足夠大，地表覆蓋從裸土變化到稠密的植被覆蓋，土壤表層含水量從干燥到濕潤。

3.3 干旱等級劃分

由于不同地區(qū)的氣候條件不同，各種干旱監(jiān)測指數(shù)的適應(yīng)范圍差異大，在劃分干旱等級和監(jiān)測、評估干旱發(fā)生時(shí)，難以滿足不同地區(qū)的災(zāi)害需求。李國正［18］將TVDI按密度劃分為10個(gè)等級，然后又將其合并為濕潤區(qū)（1～3）、偏濕潤區(qū)（4～6）、偏干旱區(qū)（7～8）、干旱區(qū)（9～10）4個(gè)等級。吳孟泉等［19］將干旱劃分為5個(gè)等級，分別是水域（0～0.2）、濕潤（0.2～0.4）、正常（0.4～0.6）、缺水（0.6～0.8）和干旱（0.8～1.0）。新疆的沙漠面積較大，干旱常年發(fā)生，利用李國正和吳孟泉的干旱等級劃分方法，研究區(qū)域除了盆地邊緣和高山腳下，基本被劃分為缺水和干旱區(qū)，無法得到直觀的新疆干旱圖。張順謙［11］根據(jù)中國氣象局提供的由土壤濕度Y進(jìn)行旱情分級標(biāo)準(zhǔn)，結(jié)合土壤濕度和研究區(qū)TVDI指數(shù)的線性擬合，反演得到各旱情等級的TVDI分級標(biāo)準(zhǔn)。該標(biāo)準(zhǔn)利用在新疆地區(qū)的TVDI分級效果較好，因此本文采用張順謙TVDI對干旱等級的劃分（表1）。

表1 干旱等級劃分表

4 結(jié)果與分析

經(jīng)過數(shù)據(jù)預(yù)處理之后，應(yīng)用MODIS數(shù)據(jù)構(gòu)建2007年～2012年春季和夏季的LST－NDVI特征空間。在ENVI軟件下，使用IDL語言編程，經(jīng)過運(yùn)行輸出LST－NDVI散點(diǎn)圖，并得到所提取的NDVI下最大和最小陸地地表溫度Tmax和Tmin，并擬合回歸獲得干濕邊擬合方程（圖1）。特征空間中的干邊斜率小于0，說明隨著植被覆蓋度的增加，Tmax越??；濕邊的斜率大于0，說明隨著植被覆蓋度的增加，Tmin越大。

圖1 干濕邊擬合方程

在ArcGIS中將全區(qū)的TVDI圖像制成干旱等級圖，干旱分級標(biāo)準(zhǔn)如表1所示。從整體研究區(qū)域來看，新疆無旱和輕旱區(qū)域主要分布在盆地邊緣的高山腳下，夏季的高山融雪形成了豐富的水資源，中旱和重旱地區(qū)主要分布在人口較為密集的平原地區(qū)，特旱區(qū)域主要分布在塔克拉瑪干沙漠、古爾班通古特沙漠以及吐魯番盆地。

從每年7月份的TVDI干旱等級劃分圖像（圖2）可以看出，近年來，新疆自治區(qū)的無旱面積變化并不明顯，但是重旱地區(qū)面積呈明顯下降的趨勢。在2007年和2008年7月，新疆多處人口密集地區(qū)處于TVDI等級的重旱和特旱等級中，例如喀什、和田、阿勒泰等地區(qū)，而吐魯番地區(qū)常年處于特旱等級。

通過對2007年～2012年春季和夏季的TVDI求均值，以及等級劃分，從而獲得2007年～2012年3月～8月的各干旱等級面積比重圖（圖3）。2009年無旱面積占全區(qū)總面積的比重最高，為74.65%，2008年的無旱面積占全區(qū)總面積的比重最低，僅為37.16%。特旱和重旱面積比重在2012年為26.97%，而在2009年僅為0.61%，輕旱和中旱面積比重在6年里變化不明顯。從數(shù)據(jù)可以看出，新疆在2009年的干旱情況相較于2007年和2008年得到了緩解，2008年和2012年的干旱程度相比其他年份更加嚴(yán)重。

圖2 2007年～2012年7月新疆地區(qū)干旱情況比較

圖3 2007年～2012年各干旱等級面積比重圖

將2007年～2012年春、夏季36個(gè)月的干旱面積比重完成折線圖（圖4），從圖中顯示可以看出，每年的干旱高峰基本發(fā)生在春季，由此可以得出，新疆的春旱比夏旱情況嚴(yán)重。并且在2010年之后，干旱面積所占比重逐漸趨于平穩(wěn)，波動(dòng)范圍逐漸縮小。

圖4 2007年～2012年3月～8月干旱面積比重圖

圖5 降水量與TVDI的關(guān)系

降水量是影響干旱的直接因子，選取烏魯木齊、阿勒泰和喀什地區(qū)的歷史降水量數(shù)據(jù)，分析降水量和TVDI的關(guān)系。從降水量數(shù)據(jù)來看，烏魯木齊市和阿勒泰地區(qū)的降水量明顯高于喀什地區(qū)，喀什地區(qū)常年處在低降水量高TVDI的環(huán)境下，可見南疆地區(qū)的干旱情況比北疆更加嚴(yán)峻。從圖5看出，降水量與TVDI有很好的相關(guān)性，隨著降水量的增大，TVDI有下降的趨勢，但有一定的時(shí)間滯后性，當(dāng)降水量趨于平穩(wěn)時(shí)，TVDI的波動(dòng)變化也不明顯，可以認(rèn)為降水量是影響TVDI的重要因子之一。

5 結(jié)束語

利用TVDI方法，獲得了2007年～2012年3月～8月新疆地區(qū)旱情變化圖，并通過多年的干旱面積變化情況對比，以及與降水量的相關(guān)性分析，對新疆存在的干旱情況進(jìn)行了研究。結(jié)果表明，新疆地區(qū)的無旱區(qū)域面積穩(wěn)定，集中在盆地邊緣、山脈腳下，特旱區(qū)域集中在沙漠地區(qū)附近。結(jié)合降水量因子可以看出，南疆的旱情比北疆嚴(yán)重。新疆地區(qū)干旱面積波動(dòng)幅度趨于平穩(wěn)，近幾年來特旱地區(qū)面積比重卻在下降。

目前，對新疆的干旱情況監(jiān)測主要還在使用單一指數(shù)的方法，主要監(jiān)測數(shù)據(jù)都來自氣象站。而單一指數(shù)的監(jiān)測不能夠全面反映新疆的干旱情況，在人口較少區(qū)域，更是無法了解當(dāng)?shù)氐恼鎸?shí)旱情水平。MODIS數(shù)據(jù)獲取快捷，解決了對大面積區(qū)域宏觀監(jiān)測的問題，TVDI方法原理簡單，應(yīng)用型強(qiáng)，解決了僅使用單一數(shù)據(jù)進(jìn)行新疆地區(qū)干旱監(jiān)測的問題。

該方法也存在諸多不足之處，TVDI方法僅使用了地表溫度數(shù)據(jù)和歸一化植被指數(shù)數(shù)據(jù)，并未考慮其他干旱影響因子。要解決這一問題，必須用其他影響因子對TVDI方法進(jìn)行驗(yàn)證，同時(shí)可以結(jié)合其他干旱指數(shù)對新疆干旱情況進(jìn)行更全面的評估，以提升干旱模型的應(yīng)用水平。

［1］ BOWERS S A，HANKS R J.Reflection of radiant energy from soil［J］.Soil Science，1965，100（2）：130－138.

［2］ WATSON K，POHN H A.Thermal inertia mapping from satellites discrimination of geologic unitin oman［J］.Journal of Research Geology Surveying，1974，2（2）：147－158.

［3］ CARLSON T N，PERRY E M.Remote estimation of soil moisture availability and fractional vegetation cover for agriculture field［J］.Agriculture and Forest Meteorology，1990，（52）：45－49.

［4］ MORAN M S，CLARKE T R.Estimating crop water deficit using the relation between surface－air temperature and spectral vegetation index［J］.Remote Sensing of Environment，1994，（49）：246－263.

［5］ BOEGH E，SOEGAARD H，HANAN N，et al.A remote sensing study of the NDVI－Ts relationship and the transpiration from sparse vegetation in the sahel based on high resolution satellite data［J］.Remote Sensing of Environment，1998，（69）：224－240.

［6］ PRICE J C.Using spatial context in satellite data to infer regional scale evapotranspiration［J］.IEEE Transactions on Geoscience and Remote Sensing，1990，（28）：940－948.

［7］ SANDHOT L，RASMUSSEN K.A simple interpretation of the surface temperature／vegetation index space for assessment of surface moisture status［J］.Remote Sensing of Environment，2002，（79）：213－224.

［8］ 肖乾廣，陳維英，盛永偉，等.用氣象衛(wèi)星監(jiān)測土壤水分的試驗(yàn)研究［J］.應(yīng)用氣象學(xué)報(bào)，1994，5（3）：312－318.

［9］ 查勇.草地植被變化遙感監(jiān)測方法研究［D］.南京：南京師范大學(xué)地理科學(xué)學(xué)院，2003.

［10］ 柳欽火，辛景峰，辛?xí)灾?，?基于地表溫度和植被指數(shù)的農(nóng)業(yè)干旱監(jiān)測方法［J］.科技導(dǎo)報(bào)，2007，25（6）：12－18.

［11］ 張順謙，卿清濤，侯美亭，等.基于溫度植被干旱指數(shù)的四川伏旱遙感監(jiān)測與影響評估［J］.農(nóng)業(yè)工程學(xué)報(bào)，2007，23（9）：141－147.

［12］ 杜靈通.基于多遠(yuǎn)空間信息的干旱模型構(gòu)建及其應(yīng)用研究［D］.南京：南京大學(xué)地理學(xué)，2013.

［13］ 白云崗，木沙·如孜，雷曉云，等.新疆干旱災(zāi)害的特征及影響因素分析［J］.人民黃河，2012，34（7）：61－63.

［14］ 姚春生，張?jiān)鱿椋魹t.用溫度植被干旱指數(shù)法反演新疆土壤濕度［J］.遙感技術(shù)與應(yīng)用，2004，19（6）：473－478.

［15］ 包安明，張小雷，方暉，等.MODIS數(shù)據(jù)在新疆生態(tài)環(huán)境建設(shè)中的應(yīng)用［J］.干旱區(qū)地理，2004，27（2）：256－260.

［16］ 楊燕玲.新疆土地利用動(dòng)態(tài)變化研究［D］.新疆：新疆農(nóng)業(yè)大學(xué)土地資源管理，2005.

［17］ 毛煒嶧，南慶紅，史紅政.新疆氣候變化特征及氣候分區(qū)方法研究［J］.氣象，2008，34（10）：67－73.

［18］ 李國正，王仰麟，吳健生，等.基于TVDI方法的黃土高原地表干燥度與土地利用的關(guān)系研究［J］.地理研究，2006，23（9）：913－919.

［19］ 吳孟泉，崔偉宏，李景剛.溫度植被干旱指數(shù)（TVDI）在復(fù)雜山區(qū)干旱監(jiān)測是應(yīng)用研究［J］.干旱區(qū)地理，2007，30（1）：30－35.

Drought Monitoring in Xinjiang by Using MODIS－TVDI

LIU Xing－h(huán)an1，ZHANG Jia－h(huán)ua1，2，XU Xiao－h(huán)ong1
（1.College of Earth Sciences，Yangtze University，Wuhan 430100；2.Lab of Digital Earth Sciences，Institute of Remote Sensing and Digital Earth，Chinese Academy of Sciences，Beijing100094）

Considering the conditions of droughts occurred frequently in recent years，Xinjiang province was selected as the study area in this study，and by using the compound MODIS data from March to August during 2007and 2012，the NDVI（normalized difference vegetation index）and LST（land surface temperature）were obtained，and the LST－NDVI feature space was firstly established.Then，the condition of TVDI（temperature vegetation dryness index）and spatial distribution of drought grade in the study area were obtained，the variation tendency of the drought in Xinjiang was analyzed，and the relationship of precipitation and TVDI was verified.The result showed that the drought area of Xinjiang had leveled off with each passing year from 2007to 2012，drought severity on south Xinjiang was severer than that on north Xinjiang in the spatial performance，and drought severity on Spring was severer than that on summer；the precipitation factor showed an important factor affecting the TVDI.

NDVI；land surface temperature；TVDI；drought monitoring；drought evaluation

10.3969／j.issn.1000－3177.2015.02.020

P412.15

A

1000－3177（2015）138－0111－05

2014－03－14

2015－01－22

中科院院地合作項(xiàng)目（Y42301101A）；CAS－TWAS項(xiàng)目（Y3YI2701KB）；中科院百人計(jì)劃項(xiàng)目（Y24002101A）。

劉星含（1989～），女，碩士，研究方向?yàn)楦珊颠b感監(jiān)測。

E－mail：lxhan0820＠126.com

張佳華（1966～），男，博士，研究方向?yàn)樯鷳B(tài)與氣象災(zāi)害遙感研究。

E－mail：zhangjh＠radi.ac.cn