基于MOD16產(chǎn)品的陜西關(guān)中地區(qū)干旱時(shí)空分布特征

何慧娟，卓　靜，李紅梅，權(quán)文婷

(陜西省農(nóng)業(yè)遙感信息中心， 陜西 西安 710014)

?

基于MOD16產(chǎn)品的陜西關(guān)中地區(qū)干旱時(shí)空分布特征

何慧娟，卓靜，李紅梅，權(quán)文婷

(陜西省農(nóng)業(yè)遙感信息中心， 陜西 西安 710014)

摘要：基于MOD16產(chǎn)品，利用作物缺水指數(shù)法，結(jié)合氣象站點(diǎn)觀測資料及MOD17植被指數(shù)信息，分析了2000—2013年陜西關(guān)中地區(qū)干旱的時(shí)空分布特征。結(jié)果表明：基于MOD16產(chǎn)品的作物缺水指數(shù)的計(jì)算結(jié)果與站點(diǎn)土壤相對濕度的變化規(guī)律一致，并且兩者的相關(guān)系數(shù)通過了P<0.001的顯著性檢驗(yàn)，說明MOD16數(shù)據(jù)可以用于關(guān)中地區(qū)干旱的時(shí)空變化分析； 作物缺水指數(shù)的空間分布不均勻，在關(guān)中地區(qū)的東部和北部值較高，較為干旱，特別是在關(guān)中東北部和西安市周邊干旱特別嚴(yán)重，作物缺水指數(shù)多年平均值在0.795～0.832，達(dá)到重旱等級(jí)，并且區(qū)域內(nèi)干旱程度呈顯著上升趨勢； 作物缺水指數(shù)有明顯的年內(nèi)變化波動(dòng)，在3—6月及10—11月值最高，作物缺水指數(shù)達(dá)0.7左右，其中3—6月是重旱和特旱集中發(fā)生的時(shí)期，而10—11月以輕旱和中旱為主； 2000—2013年作物缺水指數(shù)在2003年出現(xiàn)最低值(0.65左右)，后逐年升高，在2013年達(dá)最高值(0.8左右)； 2000—2013年輕旱面積呈減小的趨勢，中、特旱呈增加趨勢； 作物缺水指數(shù)與溫度、降水及植被的關(guān)系密切，高溫少雨時(shí)，作物缺水指數(shù)偏高，干旱容易發(fā)生，同時(shí)歸一化植被指數(shù)較低。

關(guān)鍵詞：MOD16；干旱；時(shí)空分布；作物缺水指數(shù)；關(guān)中地區(qū)

干旱是陜西省主要的自然災(zāi)害，有“十年九旱”之說。據(jù)史料記載，陜西省干旱的危害程度占所有氣象災(zāi)害的50%[1]。陜西關(guān)中地區(qū)是陜西省主要的糧食產(chǎn)區(qū)，也是陜西省干旱發(fā)生最嚴(yán)重的區(qū)域[2]，在歷史上干旱發(fā)生的次數(shù)占總年數(shù)的42.25%[3]，因此對干旱準(zhǔn)確監(jiān)測對農(nóng)業(yè)生產(chǎn)至關(guān)重要。

傳統(tǒng)干旱監(jiān)測的方法多采用氣象和水文參數(shù)[2,4-5]，其數(shù)據(jù)是基于站點(diǎn)的，而觀測站點(diǎn)較少且分布不均勻，因此不能全面細(xì)致地反映區(qū)域干旱的時(shí)空變化特征?，F(xiàn)如今遙感技術(shù)的迅速發(fā)展，提供了大范圍高分辨率多時(shí)象的地表特征信息，為干旱監(jiān)測提供了新途徑。大量的研究[6-11]表明遙感數(shù)據(jù)已經(jīng)能夠較為準(zhǔn)確地反映土壤干旱信息，但是由于遙感數(shù)據(jù)信息量大，處理相對復(fù)雜，已有的遙感監(jiān)測研究多基于干旱發(fā)生的一個(gè)時(shí)段，沒有長期連續(xù)系統(tǒng)的干旱遙感監(jiān)測數(shù)據(jù)，因此對分析區(qū)域干旱的時(shí)空分布特征有一定的局限性。2011年美國NASA研究團(tuán)隊(duì)在MODIS遙感數(shù)據(jù)蒸散反演算法上取得新進(jìn)展[12-13]，并且發(fā)布了2000—2013年全球MODIS陸地蒸散參數(shù)產(chǎn)品數(shù)據(jù)(MOD16)，具有較高的時(shí)間分辨率和空間分辨率，以及免費(fèi)獲取等特點(diǎn)，在研究水分循環(huán)等方面具有一定的優(yōu)勢[14]。另外數(shù)據(jù)集中包含蒸散和潛在蒸散兩個(gè)參數(shù)，正好可以用于作物缺水模型的計(jì)算，為干旱監(jiān)測研究提供了方便可行的數(shù)據(jù)。因此本文將基于MOD16蒸散數(shù)據(jù)產(chǎn)品，利用作物缺水指數(shù)法，探討2000—2013年陜西關(guān)中地區(qū)干旱的時(shí)空變化特征。

1研究區(qū)概況與數(shù)據(jù)來源

1.1研究區(qū)域概況

研究區(qū)域沒有按照行政邊界進(jìn)行劃分，而是利用2014年國家基礎(chǔ)地理信息中心制作的2010年土地利用數(shù)據(jù)(http://www.globallandcover.com/GLC30Download/index.aspx)提取了關(guān)中平原作為研究區(qū)域。從圖1可以看到，關(guān)中地區(qū)位于陜西省中部，西起寶雞、東至潼關(guān)，南接秦嶺，北到陜北高原，西窄東寬，土地覆蓋類型為明顯的大片耕地區(qū)域，農(nóng)作物以冬小麥和夏玉米為主，為一年兩熟制，是陜西省主要的糧食產(chǎn)區(qū)。研究區(qū)域是斷層陷落區(qū)即地塹，后經(jīng)渭河及其支流涇河、洛河等沖積而成。區(qū)域內(nèi)地形平坦，土層深厚，蓄水性好，但是區(qū)域內(nèi)水資源相對不足，干旱是造成糧食減產(chǎn)的主要原因[15]。因此這樣的劃分更符合自然環(huán)境規(guī)律，利于區(qū)域干旱的研究。

圖1陜西關(guān)中地區(qū)土地利用和氣象站點(diǎn)分布

Fig.1Land use and distribution of meteorological stations in Guanzhong region of Shaanxi Province

1.2數(shù)據(jù)來源

MOD16數(shù)據(jù)是利用MODIS衛(wèi)星傳感器數(shù)據(jù)結(jié)合氣象數(shù)據(jù)反演得到的地面蒸散數(shù)據(jù)集。文章所用的蒸散數(shù)據(jù)是2011年Mu等在2007年算法的基礎(chǔ)上改進(jìn)的蒸散算法[12-13]。蒸散算法基于Penman-Monteith公式。蒸散包括地表蒸發(fā)和植物蒸騰。陸地蒸散水分來源于土壤、植被攔截未達(dá)地面的雨水、植被莖葉氣孔的呼吸。計(jì)算需要反照率、葉面積指數(shù)、地表覆蓋類型等遙感數(shù)據(jù)和氣壓、氣溫、濕度、輻射等氣象數(shù)據(jù)。本文使用的是2000—2013年月、年合成1 km分辨率蒸散(Evapotranspiration，ET)、潛在蒸散(Potential evapotranspiration，PET)數(shù)據(jù)，軌道號(hào)是h26v05和h27v05，下載地址為：ftp://ftp.ntsg.umt.edu。數(shù)據(jù)在MODLAND提供的MRT(Modis Reprojection Tool)投影轉(zhuǎn)換工具中對h26v05和h27v05兩區(qū)圖像進(jìn)行拼接和投影轉(zhuǎn)換，轉(zhuǎn)換為等經(jīng)緯度坐標(biāo)投影，基準(zhǔn)面為WGS-84坐標(biāo)系。

使用歸一化植被指數(shù)(Normailized difference vegetation index，NDVI)來反應(yīng)植被長勢信息，文章使用的是2000—2013年MODIS植被指數(shù)數(shù)據(jù)產(chǎn)品(MOD13Q1)，為16天合成250 m分辨率數(shù)據(jù)，軌道號(hào)是h26v05和h27v05，下載地址為：http://e4ftl01.cr.usgs.gov。同樣使用MRT工具提取了NDVI數(shù)據(jù)并進(jìn)行了拼接和投影轉(zhuǎn)換。為了與MOD16數(shù)據(jù)對應(yīng)，將16天合成數(shù)據(jù)內(nèi)插成月值。年NDVI數(shù)據(jù)處理是將一年NDVI值進(jìn)行了最大值合成，用來反映一年中植被生長最好的狀況。

氣象數(shù)據(jù)包括2000—2013年33個(gè)氣象站點(diǎn)觀測的降水量、平均氣溫的日值數(shù)據(jù)及5個(gè)農(nóng)業(yè)氣象站點(diǎn)2013年土壤相對濕度時(shí)值數(shù)據(jù)。數(shù)據(jù)來源于陜西省氣象局，站點(diǎn)分布情況見圖1。為了與蒸散數(shù)據(jù)進(jìn)行對比分析，將數(shù)據(jù)處理為月值和年值。

1.3作物缺水指數(shù)法

用于干旱監(jiān)測的方法很多，有植被指數(shù)法、熱慣量法、作物缺水指數(shù)法等。其中作物缺水指數(shù)法由于物理意義明確、適用范圍廣等特點(diǎn)在干旱研究方面得到廣泛應(yīng)用[9-11]。由于水分供應(yīng)是蒸散過程的基本條件之一，土壤含水量對蒸散速率有一定的影響[10]，因此作物缺水指數(shù)(Crop water stress index，CWSI)的定義為：

CWSI=1-ET/PET

式中：ET為實(shí)際蒸散量(mm)；PET為潛在蒸散量(mm)，即在區(qū)域供水充足條件下的蒸散量。實(shí)際蒸散量與潛在蒸散量的比值可以直觀地反應(yīng)出區(qū)域水分的盈虧狀況，從而可以用于土壤干旱的監(jiān)測，從式中可見，CWSI值在0～1之間，值越大，表示區(qū)域越干旱缺水，值越小表示越濕潤。在MOD16數(shù)據(jù)集中包含ET和PET數(shù)據(jù)，正好能夠方便用于作物缺水指數(shù)的計(jì)算。

2數(shù)據(jù)分析

2.1干旱檢測結(jié)果檢驗(yàn)

為了保證遙感監(jiān)測數(shù)據(jù)的有效性，首先需要對區(qū)域MOD16產(chǎn)品計(jì)算的作物缺水指數(shù)進(jìn)行檢驗(yàn)。對于干旱的監(jiān)測最可靠的方法是對土壤水分的直接測量，國家氣象局以20 cm土壤相對濕度作為旱情的分級(jí)標(biāo)準(zhǔn)[16],因此，基于遙感數(shù)據(jù)的CWSI將與測站測量的20 cm土壤相對濕度數(shù)據(jù)進(jìn)行對比。圖2給出了2013年區(qū)域內(nèi)五個(gè)測站CWSI與土壤相對濕度數(shù)據(jù)的對比結(jié)果。

圖22013年作物缺水指數(shù)與土壤相對濕度的比較

Fig.2Comparison between CWSI and soil relative humidity in 2013

從圖2中可以看出兩者的結(jié)果對于干旱的反映比較吻合。在較干旱的時(shí)期，土壤相對濕度較低，CWSI較高；在較濕潤的時(shí)期，土壤相對濕度較高，CWSI較低。并且，測站之間的對比也符合這一規(guī)律，例如在關(guān)中最東側(cè)的韓城站與最西側(cè)的隴縣站比較：韓城站較為干旱，20 cm土壤相對濕度整體偏低，CWSI偏高；而隴縣站較為濕潤，20 cm土壤相對濕度整體偏高，CWSI偏低。 這一變化也與2013年關(guān)中地區(qū)春季干旱情況較為一致，2013年4月1日陜西省氣象局啟動(dòng)了重大氣象災(zāi)害(干旱)Ⅲ級(jí)應(yīng)急響應(yīng)，至4月24日解除了干旱Ⅲ級(jí)應(yīng)急響應(yīng)狀態(tài)[8]。從測站CWSI的變化情況來看，3月測站CWSI處于高位，到5月各測站CWSI都有不同程度的降低，這與實(shí)際的干旱發(fā)展情況比較一致。

再對CWSI與20 cm土壤相對濕度的散點(diǎn)圖進(jìn)行分析(圖3)，發(fā)現(xiàn)呈顯著的負(fù)相關(guān)性，相關(guān)系數(shù)通過了0.001的顯著性檢驗(yàn)。說明用MOD16蒸散產(chǎn)品計(jì)算的作物缺水指數(shù)用于關(guān)中地區(qū)干旱的監(jiān)測是可行的。另外通過對CWSI與20 cm土壤相對濕度建立回歸方程，再參照旱情的劃分標(biāo)準(zhǔn)[16]，可以將20 cm土壤相對濕度的干旱劃分標(biāo)準(zhǔn)轉(zhuǎn)換為CWSI干旱等級(jí)的劃分標(biāo)準(zhǔn)，詳細(xì)劃分見表1。

圖3　作物缺水指數(shù)與土壤相對濕度的關(guān)系

2.2區(qū)域干旱的時(shí)空變化特征

2.2.1多年平均干旱的空間分布狀況從2000—2013年CWSI多年平均值反映的干旱空間分布來看(圖4a)，關(guān)中地區(qū)干旱呈現(xiàn)明顯的空間分異性規(guī)律，關(guān)中地區(qū)的東部和北部較為干旱，特別是在關(guān)中東北部干旱特別嚴(yán)重，CWSI多年平均值在0.795～0.832，達(dá)到重旱等級(jí)。這一結(jié)果與譚學(xué)志等[4]和蔡新玲等[5]的研究結(jié)果較為一致。另外西安市周邊，也較為干旱，特別是城市北部地區(qū)。干旱的空間分布規(guī)律與氣候條件密切相關(guān)，結(jié)合溫度、降水的空間分布規(guī)律(圖4c、圖4d)可以看到，降水從西南向東北呈減小的趨勢，因此CWSI從西南向東北增加，而溫度則在關(guān)中盆地中央最高，向南向北逐漸降低，尤其在城市周邊，特別是西安城市周邊溫度最高，因此在西安周邊干旱也較為嚴(yán)重。結(jié)合NDVI的空間分布(圖4b)也可以看到，CWSI高值區(qū)干旱較重，NDVI也較低，而在CWSI的低值區(qū)水分較充足，NDVI也較高。

2.2.2干旱的月變化圖5為根據(jù)MOD16數(shù)據(jù)計(jì)算的2000—2013年關(guān)中地區(qū)CWSI干旱月平均變化?？梢钥吹紺WSI有明顯的年內(nèi)變化波動(dòng)，在3—6月及10—11月CWSI值較高，CWSI達(dá)0.7左右，6月CWSI最高，超過了0.8。CWSI的低值出現(xiàn)在夏季，8月最低，在0.55左右。從干旱情況來看，干旱發(fā)生在春季、夏初及秋季。3—6月是重旱和特旱集中發(fā)生的時(shí)期，而10—11月以輕旱和中旱為主。這與周丹等[2]的陜西省春季和秋季干旱較嚴(yán)重的研究結(jié)果較為一致。

再結(jié)合溫度、降水及植被的變化情況來分析，溫度和降水呈單峰變化，溫度在1月最低，7月達(dá)到最高值，降水集中在7、8、9月，9月降水最多。植被呈雙峰變化，有兩個(gè)峰值(5月和8月)，這是因?yàn)殛P(guān)中地區(qū)作物以冬小麥和夏玉米為主，6月冬小麥?zhǔn)斋@，NDVI值出現(xiàn)明顯的下降。ET的變化與降水的變化相似，8月ET值最高，而ET在6月出現(xiàn)一個(gè)低谷與NDVI的急劇減小有關(guān)。PET與溫度的變化較為相似，月變化較均勻，在6月達(dá)到最高值。從CWSI、ET、PET與溫度、降水及NDVI的月變化相關(guān)系數(shù)計(jì)算結(jié)果(表2)也可以看到，雖然CWSI與溫度、降水及NDVI的相關(guān)系數(shù)不高，但是CWSI計(jì)算的兩個(gè)分量ET、PET卻與溫度、降水及NDVI密切相關(guān)。年內(nèi)變化ET與降水、NDVI相關(guān)性強(qiáng)，而PET與溫度相關(guān)性最強(qiáng)。

表2　干旱月變化下的相關(guān)系數(shù)

注：**代表P<0.05; ***代表P<0.01。

Note：** meansP<0.05; *** meansP<0.01.

圖4　干旱(a)、歸一化植被指數(shù)(b)、降水量(c)、平均溫度(d)多年平均的空間分布

圖5　陜西關(guān)中地區(qū)干旱月變化

圖6陜西關(guān)中地區(qū)干旱年際變化

Fig.6Annual variations of drought in Guanzhong

region of Shaanxi Province

2.2.3干旱的年際變化特征從區(qū)域干旱的年際變化情況來看(圖6)，CWSI在2003年出現(xiàn)了一個(gè)最低值0.65左右，后逐年升高，在2013年最高，達(dá)0.8左右。從干旱程度來看，輕旱面積呈減小的趨勢，中旱、重旱和特旱形勢從2003年開始呈加劇的趨勢。這一研究結(jié)果與周丹等[2]近20年陜西省干旱發(fā)生頻率呈明顯增加趨勢的研究結(jié)果較為一致。將CWSI與溫度降水及NDVI的年際變化結(jié)合分析來看(表3)，ET與NDVI、降水的相關(guān)系數(shù)較高，而PET與溫度相關(guān)性最為密切，因此CWSI與降水、NDVI呈負(fù)相關(guān)，與溫度呈正相關(guān)。也就是在降水充足、溫度較低的年份，CWSI低干旱程度輕，NDVI值也較高；而在溫度較高、降水缺乏的年份，CWSI高干旱程度重，NDVI值也較低。這一變化規(guī)律也可以說明利用 MOD16數(shù)據(jù)計(jì)算的作物缺水指數(shù)符合關(guān)中地區(qū)水熱條件及植被的變化規(guī)律，可以反映關(guān)中地區(qū)干旱的變化情況。

圖7作物缺水指數(shù)年際變化趨勢的空間分布

Fig.7Spatial distribution of CWSI annual variation tendency

最后計(jì)算了關(guān)中地區(qū)每個(gè)像元CWSI的變化趨勢(圖7)，可以看到關(guān)中地區(qū)CWSI主要呈明顯上升趨勢，在區(qū)域東部和西安市周邊上升趨勢顯著。這一變化可能與全球氣候變暖有關(guān)，據(jù)IPCC最新報(bào)告[17]顯示全球溫度升高是顯著的，因此在全球變暖的背景下使得干旱加劇。再有在城市周邊干旱的加劇與城市熱島效應(yīng)有關(guān)，城市周邊屬于干旱脆弱性地區(qū)[18]，城市的擴(kuò)張使得周邊溫度上升，干旱也隨之加劇。

表3　干旱年際變化下的相關(guān)系數(shù)

注：*代表P<0.1; **代表P<0.05。

Note：* meansP<0.1; ** meansP<0.05.

3結(jié)論與討論

利用MOD16產(chǎn)品計(jì)算的作物缺水指數(shù)CWSI，結(jié)合氣象站點(diǎn)觀測資料及MOD17植被指數(shù)信息，分析了陜西關(guān)中地區(qū)干旱的時(shí)空分布特征，主要結(jié)論為：

1) CWSI計(jì)算結(jié)果與站點(diǎn)土壤相對濕度的變化規(guī)律一致，并且兩者的相關(guān)系數(shù)通過了P<0.001的顯著性檢驗(yàn)，說明基于MOD16產(chǎn)品的CWSI的計(jì)算結(jié)果可以用于關(guān)中地區(qū)干旱時(shí)空變化的分析。

2) 從空間變化來看，CWSI多年平均的空間分布不均勻，在關(guān)中地區(qū)的東部和北部較為干旱，特別是在關(guān)中東北部干旱特別嚴(yán)重，CWSI多年平均值在0.795～0.832，達(dá)到重旱等級(jí)。另外西安市周邊，也較為干旱，特別是在西安城區(qū)北部地區(qū)，達(dá)到重旱等級(jí)。關(guān)中地區(qū)CWSI變化趨勢主要呈上升趨勢，在區(qū)域東部和西安市周邊上升趨勢顯著。

3) 從時(shí)間變化來看，CWSI有明顯的年內(nèi)變化波動(dòng)，在3—6月及10—11月CWSI值較高，CWSI達(dá)0.7左右，3—6月是重旱和特旱集中發(fā)生的時(shí)期，而10—11月以輕旱和中旱為主。CWSI在2003年出現(xiàn)最低值0.65左右，后逐年升高，在2013年達(dá)最高值0.8左右。2000—2013年輕旱面積呈減小的趨勢，中旱、重旱和特旱從2003年開始呈加劇的趨勢。

4) CWSI與溫度、降水及植被的關(guān)系密切。溫度較高(低)、降水較少(多)時(shí)，CWSI較高(低)，干旱較嚴(yán)重(輕)，NDVI值也較低(高)。

從研究結(jié)果可以看到，基于MOD16遙感數(shù)據(jù)反映的干旱的時(shí)空變換規(guī)律與基于氣象站點(diǎn)的干旱的時(shí)空變化規(guī)律的研究較一致，而基于遙感數(shù)據(jù)反映的干旱時(shí)空分布更為細(xì)致，更有利于我們認(rèn)識(shí)區(qū)域干旱的分布和發(fā)展。例如，通過分析干旱的年內(nèi)變化規(guī)律可用于不同作物在不同時(shí)期的節(jié)水灌溉，在冬小麥的生長期(10月—次年6月)較為缺水，尤其冬前分蘗期(11月前后)和拔節(jié)抽穗灌漿期(3—6月)干旱最為嚴(yán)重，此時(shí)也是小麥耗水的高峰期[19]，應(yīng)注意隨時(shí)查看墑情及麥苗長勢，做到及時(shí)灌溉。而在夏玉米種植期(6—10月)，雖然降水較為充分，CWSI值較低，但是夏季多發(fā)生局地強(qiáng)降雨，因此田間會(huì)出現(xiàn)時(shí)澇時(shí)旱的現(xiàn)象，應(yīng)該隨時(shí)注意玉米對水分的要求。再有，通過CWSI的空間分布特征，我們可以看到關(guān)中東部和城市周邊是干旱發(fā)生的重點(diǎn)區(qū)域，應(yīng)該加強(qiáng)對水分利用的管理。另外人為因素是干旱加劇的主要原因，全球變暖和城市擴(kuò)張使得干旱進(jìn)一步惡化，這警示我們需要注意溫度上升對作物生長帶來的影響，加強(qiáng)對水分的利用管理，及時(shí)應(yīng)對干旱變化。

參 考 文 獻(xiàn)：

[1]《陜西歷史自然災(zāi)害簡要紀(jì)實(shí)》編委會(huì)．陜西歷史自然災(zāi)害簡要紀(jì)實(shí)[M]．北京:氣象出版社，2004:6-39．

[2]周丹，張勃，任培貴，等．基于標(biāo)準(zhǔn)化降水蒸散指數(shù)的陜西省近50年干旱特征分析[J]．自然資源學(xué)報(bào)，2014，29(4)：676-688．

[3]文星，廖杰，段翰晨，等．基于歷史文獻(xiàn)記載的關(guān)中地區(qū)中世紀(jì)干旱過程[J]．干旱區(qū)資源與環(huán)境，2012，26(3)：83-87．

[4]譚學(xué)志，粟曉玲，邵東國．基于SPI的陜西關(guān)中地區(qū)氣象干旱時(shí)空特征分析[J]．干旱地區(qū)農(nóng)業(yè)研究，2011，29(3)：224-229．

[5]蔡新玲，葉殿秀，李茜，等．基于CI指數(shù)的陜西干旱時(shí)空變化特征分析[J]．干旱地區(qū)農(nóng)業(yè)研究，2013，31(5)：1-8．

[6]劉歡，劉榮高，劉世陽．干旱遙感監(jiān)測方法及其應(yīng)用發(fā)展[J]．地球信息科學(xué)學(xué)報(bào)，2012，14(2)：232-239．

[7]宋榮杰．基于遙感的旱區(qū)土壤濕度反演方法研究[D]．楊凌：西北農(nóng)林科技大學(xué)，2013．

[8]王鵬新，林巧，張樹譽(yù)，等．2013年關(guān)中地區(qū)春旱的實(shí)時(shí)監(jiān)測與剖析[J]．干旱地區(qū)農(nóng)業(yè)研究，2013，31(3)：234-259．

[9]王玲玲, 張友靜, 佘遠(yuǎn)見，等．遙感旱情監(jiān)測方法的比較與分析[J]．遙感應(yīng)用，2010，(5)：49-53．

[10]王玉娟，王樹東，曾紅娟，等．基于作物缺水指數(shù)法的渭河流域干旱特征[J]．干旱區(qū)研究，2014，31(1)：118-124．

[11]劉安麟，李星敏，何延波，等．作物缺水指數(shù)法的簡化及在干旱遙感監(jiān)測中的應(yīng)用[J]．應(yīng)用生態(tài)學(xué)報(bào)，2004，15(2)：210-214．

[12]Mu Q Z, Heinsch F A, Zhao M S, et al. Development of a global evapotranspiration algorithm based on MODIS and global mETeorology data[J]. Remote sensing of Environment,2007,111:519-536.

[13]Mu Q Z, Zhao M S, Running S W, Improvements to a MODIS global terrestrial evapotranspiration algorithm[J]. Remote Sensing of Environment,2011,115:1781-1800.

[14]吳桂平，劉元波，趙曉松，等．基于MOD16產(chǎn)品的鄱陽湖流域地表蒸散量時(shí)空分布特征[J]．地理研究，2013，32(4)：617-627．

[15]龐文保，魯淵平．渭河平原的水資源與農(nóng)作物需水量的研究[J]．陜西氣象，1996，(4)：19-22．

[16]張強(qiáng)，周旭凱，肖風(fēng)勁，等.GB/T20481-2006.氣象干旱等級(jí)[S].北京：中國標(biāo)準(zhǔn)出版社，2006.

[17]IPCC. Working Group Ⅰ Contribution to the IPCC Fifth Assessment Report (AR5). Climate Change 2013: The Physical Science Basis[R]. Cambridge: Cambridge University Press, 2013.

[18]喻忠磊，楊新軍，石育中．關(guān)中地區(qū)城市干旱脆弱性評(píng)價(jià)[J]．資源科學(xué)，2012，34(3)：581-588．

[19]李驥，皇甫庭，王志紅．小麥-玉米復(fù)種連作關(guān)鍵技術(shù)[M]．北京：中國農(nóng)業(yè)科學(xué)技術(shù)出版社，2012:15-18．

Spatial-temporal distribution characteristics of drought in

Guanzhong region of Shaanxi Province based on MOD16 products

HE Hui-juan, ZHUO Jing, LI Hong-mei, QUAN Wen-ting

(ShaanxiRemote-SensingInformationCenterforAgriculture，Xi'an,Shaanxi710014,China)

Abstract：Based on MOD16 products and the method of crop water stress index, combined with meteorological station observations and MOD17 vegetation index information, this article analyzed the Spatial-temporal variation characteristics of drought in the Guanzhong region of Shaanxi Province during 2000—2013. The results showed that crop water stress index based on MOD16 products was consistent with the variation of soil relative humidity, and the correlation coefficient between them was significant (P<0.001). This relationship indicated that MOD16 products could be used to analyze the variation features in Guanzhong region. In addition, the crop water stress index had not been evenly distributed. In eastern and northern of Guanzhong regions, its value was high, suggesting relatively dry condition. Particularly, in the northeast of Guanzhong region and around Xi’an city, the crop water stress index was 0.795～0.832, reaching a heavy drought level. In the meantime, drought severity was becoming elevated. Furthermore, the crop water stress index had obvious monthly fluctuation variations, reaching its highest value of 0.7 in March to June and October to November. From March to June, heavy and severe drought events mainly occurred. Rather from October to November, light and medium drought events mostly took place. In Guanzhong region, the crop water stress index exhibited an increasing trend, after minimum and maximum values of 0.65 in 2003 and 0.8 in 2013, respectively. The area suffered light drought showed a trend of becoming reduced. However, areas with medium to severe drought were increased during 2000—2013. Moreover, temperature, precipitation and vegetation were closely related to crop water stress index. During the period with higher temperature and less precipitation, the crop water stress index was high, drought was likely to happen, and the normalized difference vegetation index was low.

Keywords:MOD16; drought; Spatial-temporal distribution; CWSI; Guanzhong region

中圖分類號(hào)：S127；S165+.2

文獻(xiàn)標(biāo)志碼：A

作者簡介：何慧娟(1983—)，女，陜西岐山人，碩士，工程師，主要從事生態(tài)環(huán)境遙感監(jiān)測方面的研究。E-mail：393621703@qq.com。

基金項(xiàng)目：陜西省氣象局科技創(chuàng)新基金計(jì)劃項(xiàng)目(2014M-14)

收稿日期：2014-12-17

doi:10.7606/j.issn.1000-7601.2016.01.36

文章編號(hào)：1000-7601(2016)01-0236-06