吳 黎, 解文歡, 張有智, 宋麗娟, 李 巖, 呂志群

(黑龍江省農(nóng)業(yè)科學(xué)院 農(nóng)業(yè)遙感與信息研究所, 哈爾濱 150086)

農(nóng)業(yè)干旱是作物生長(zhǎng)過程中因供水不足，阻礙作物正常生長(zhǎng)而發(fā)生的水量供應(yīng)不平衡的現(xiàn)象[1]。黑龍江省是典型的旱作農(nóng)業(yè)區(qū)，近20 a來黑龍江省旱災(zāi)占總農(nóng)業(yè)氣象災(zāi)害總面積的63%，占據(jù)省內(nèi)三大氣象災(zāi)害首位。遙感監(jiān)測(cè)技術(shù)以其宏觀、快速、動(dòng)態(tài)、經(jīng)濟(jì)的優(yōu)勢(shì)，成為旱災(zāi)監(jiān)測(cè)研究的有效重要手段。眾多學(xué)者研究表明利用地面溫度和植被指數(shù)特征空間耦合而成的溫度植被干旱指數(shù)(TVDI)監(jiān)測(cè)旱情方法優(yōu)越、可行性高[2]。利用TVDI模型監(jiān)測(cè)地表旱情時(shí)，其受傳感器類型和數(shù)據(jù)空間分辨率的影響較小，因此被廣泛應(yīng)用到干旱監(jiān)測(cè)中[3-8]。

眾多學(xué)者對(duì)黑龍江省的干旱研究主要?dú)w納為兩類：一是基于長(zhǎng)時(shí)間的氣象站點(diǎn)數(shù)據(jù)采用不同干旱方法進(jìn)行時(shí)空特征、演變規(guī)律、變化趨勢(shì)、災(zāi)害評(píng)價(jià)等研究，二是基于遙感數(shù)據(jù)進(jìn)行干旱監(jiān)測(cè)及預(yù)測(cè)模型分析。分析干旱時(shí)空特征研究基本都是以氣象數(shù)據(jù)為基礎(chǔ)的，其原因是氣象數(shù)據(jù)歷史時(shí)間更久，從1951年開始我國便已有氣象站觀測(cè)數(shù)據(jù)，時(shí)間越長(zhǎng)對(duì)規(guī)律和趨勢(shì)分析就越有利，這是氣象數(shù)據(jù)的優(yōu)點(diǎn)。但氣象站屬于點(diǎn)數(shù)據(jù)，不能代表面上干旱情況，學(xué)者需采用一系列方法進(jìn)行面數(shù)據(jù)的轉(zhuǎn)化，這樣在監(jiān)測(cè)精度上有所損失。遙感數(shù)據(jù)屬于面數(shù)據(jù)，以像元為單元，可監(jiān)測(cè)每一個(gè)像元的干旱程度彌補(bǔ)氣象數(shù)據(jù)的不足，且MODIS數(shù)據(jù)具有時(shí)效性好、幅寬大、波段范圍廣，數(shù)據(jù)開放的優(yōu)勢(shì)，可快速實(shí)現(xiàn)全省范圍內(nèi)土壤干旱監(jiān)測(cè)，在防治干旱，指導(dǎo)精準(zhǔn)灌溉，節(jié)約水資源，降低農(nóng)民損失方面均具有重要意義。目前對(duì)黑龍江省利用遙感數(shù)據(jù)進(jìn)行20 a的時(shí)空變化特征研究特別缺少，學(xué)者主要集中在對(duì)單一時(shí)刻或某一時(shí)段的干旱監(jiān)測(cè)上。針對(duì)如上問題本文基于MODIS時(shí)序數(shù)據(jù)計(jì)算逐月TVDI，累計(jì)20 a間黑龍江省干旱發(fā)生頻率及干旱強(qiáng)度，分析黑龍江省干旱時(shí)間特征，找到規(guī)律并同時(shí)開展空間變化特征研究，為后續(xù)作者進(jìn)行黑龍江省干旱預(yù)警研究奠定基礎(chǔ)，為政府抗旱救災(zāi)、制定農(nóng)業(yè)政策和管控農(nóng)業(yè)項(xiàng)目提供科學(xué)數(shù)據(jù)支撐。

1 研究區(qū)概況

黑龍江省地理位置介于121°11′—135°05′E和43°25′—53°33′N。年平均氣溫在-6～4℃，無霜期多在100～140 d，年降雨量介于400～650 mm，中部最多，東部次之，北部和西部最少，黑龍江省地勢(shì)西北高東南低，農(nóng)業(yè)地域特征是東西兩大平原、南北兩大山地。

2 數(shù)據(jù)來源與研究方法

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

中國科學(xué)院計(jì)算機(jī)網(wǎng)絡(luò)信息中心(http:∥www.nsdata.cn/)獲取MODIS標(biāo)準(zhǔn)產(chǎn)品中16 d合成的植被指數(shù)MOD13A2，8 d合成的地表溫度MOD11A2數(shù)據(jù)。構(gòu)建溫度植被干旱指數(shù)(TVDI)需將植被指數(shù)(MOD13A2)及地表溫度(MOD11A2)數(shù)據(jù)統(tǒng)一到相同的時(shí)間分辨率和相同的空間分辨率上，將臨近的兩個(gè)8 d的MOD11A2以最大值法(同一點(diǎn)兩景影像中取最大值作為新影像的相元值·)合成16 d的數(shù)據(jù)。兩種數(shù)據(jù)的空間分辨率均為1 km。數(shù)據(jù)時(shí)間從2000—2019年，每年的5月初至9月末，涵蓋農(nóng)作物的整個(gè)生育期。從(https:∥ladsweb.nascom.nasa.gov/data/search.html)網(wǎng)站獲取土地覆蓋類型產(chǎn)品數(shù)據(jù)(MCD12Q1)，獲取黑龍江省耕地區(qū)域界限范圍。根據(jù)國際地圈生物圈計(jì)劃(IGBP)的全球植被分類方案，選取DN值為12的像元為耕地。統(tǒng)一全部圖幅投影，后進(jìn)行圖幅拼接，再利用黑龍江省行政區(qū)劃圖進(jìn)行黑龍江省遙感影像的裁剪，得到黑龍江省歸一化植被指數(shù)和地表溫度數(shù)據(jù)；最后利用MCD12Q1數(shù)據(jù)提取黑龍江省耕地，對(duì)MODIS所有時(shí)序數(shù)據(jù)進(jìn)行掩膜處理，獲得最終的農(nóng)作物研究區(qū)的歸一化植被指數(shù)和地表溫度時(shí)序數(shù)據(jù)集。

從中國氣象數(shù)據(jù)共享網(wǎng)(http:∥cdc.nmic.cn/)獲取黑龍江省40個(gè)氣象站2000—2019年氣象觀測(cè)降雨量數(shù)據(jù)用于計(jì)算標(biāo)準(zhǔn)化降水指數(shù)(SPI)[9]。

2.2 研究方法

(1) TVDI模型法。國內(nèi)外眾多學(xué)者研究不同時(shí)間分辨率和空間尺度的植被指數(shù)(NDVI)和地表溫度(Ts)的關(guān)系表明，二者存在明顯的負(fù)相關(guān)關(guān)系。Sandholt等[10]在研究NDVI-Ts特征空間時(shí)發(fā)現(xiàn)了很多條直線，據(jù)此提出了溫度植被干旱指數(shù)(TVDI)監(jiān)測(cè)地表相對(duì)含水量狀況，計(jì)算公式為：

(1)

Tsmax=a1+b1×NDVI

(2)

Tsmin=a2+b2×NDVI

(3)

將公式(2—3)代入公式(1)

(4)

式中:Ts為地表溫度;Tsmin為相同NDVI條件下的最小地表溫度;Tsmax為相同NDVI條件下的最大地表溫度;a1，a2，b1，b2為擬合方程系數(shù)。

TVDI值為0～1，引用作者較前TVDI等級(jí)劃分結(jié)果[11]，共分5個(gè)強(qiáng)度其中無旱對(duì)應(yīng)的強(qiáng)度等級(jí)I(Intensity)為1，輕旱為2，中旱為3，重旱為4，特旱為5，相對(duì)應(yīng)的TVDI值范圍見表1；20 a總干旱強(qiáng)度為I20，定義I20<60為干旱頻率低強(qiáng)度弱，60≤I20<80為干旱發(fā)生頻率中強(qiáng)度中，80≤I20為干旱頻率高強(qiáng)度大。

表1 溫度植被干旱指數(shù)干旱強(qiáng)度

為了研究黑龍江省干旱特征，分析年平均干旱強(qiáng)度和年內(nèi)極端干旱(年內(nèi)月平均TVDI最大值)情況，本文將2000—2019年每年5—9月月溫度植被干旱指數(shù)數(shù)據(jù)求平均獲得年平均溫度植被干旱指數(shù)(TVDImean)，其代表的是作物生育期內(nèi)平均受旱程度；一年中5—9月的每月TVDI值中最大溫度植被干旱指數(shù)為TVDImax，其代表的是作物生育期內(nèi)極端受旱情況。采用一元線性回歸分析方法研究時(shí)序TVDI數(shù)據(jù)的趨勢(shì)特征。

(2) 標(biāo)準(zhǔn)化降水指數(shù)。降雨分布是一種偏態(tài)分布，因此直接用降雨量很難在不同時(shí)間、不同地區(qū)上進(jìn)行相互比較，標(biāo)準(zhǔn)化降水指數(shù)(SPI)是表示某時(shí)段降水量出現(xiàn)概率的指標(biāo)之一，因其具有概率屬性，其能很好的比較不同區(qū)域的干旱狀況[12]。通過計(jì)數(shù)出累積降水的概率密度函數(shù)Γ，進(jìn)行正態(tài)標(biāo)準(zhǔn)化處理，其計(jì)算方法是在特定時(shí)間段內(nèi)產(chǎn)生的降水量與平均值之差除以標(biāo)準(zhǔn)差，其中降水量平均值和標(biāo)準(zhǔn)差是根據(jù)過去的歷年氣象觀測(cè)記錄來確定的[13]。公式如下：

式中:x為降水樣本值;H(x)為與Γ函數(shù)相關(guān)的降水分布累積概率;c0=2.515 517;c1=0.802 853;c2=0.010 328;d1=1.432 788;d2=0.189 269;d3=0.001 308。

本文計(jì)算的是1個(gè)月時(shí)間尺度SPI下黑龍江各站點(diǎn)的SPI指數(shù)值。

3 結(jié)果與分析

3.1 精度驗(yàn)證

標(biāo)準(zhǔn)化降水指數(shù)(SPI)是氣象干旱監(jiān)測(cè)中被廣泛應(yīng)用的指標(biāo)，本研究利用SPI與TVDI進(jìn)行相互精度驗(yàn)證，可明確遙感TVDI干旱指標(biāo)的指導(dǎo)意義[14-16]。SPI選取2000—2019年黑龍江省40個(gè)氣象站的5—9月的月降水量數(shù)據(jù)，計(jì)算1個(gè)月尺度SPI值，求40個(gè)站點(diǎn)平均值代表全省1個(gè)月尺度SPI情況，最后將1個(gè)月SPI值累積求平均獲得年SPI值。同時(shí)獲取與SPI時(shí)間相匹配的月TVDI值，將月TVDI值累積求平均獲得年平均TVDI(TVDImean)。對(duì)年SPI與TVDImean進(jìn)行關(guān)系分析，可見年平均SPI與TVDImean有負(fù)相關(guān)關(guān)系，相關(guān)系數(shù)為-0.67，公式為：SPI=-3.2737TVDI+1.9136，p值為0.002，表明年平均SPI與TVDImean存在極顯著負(fù)相關(guān)性。說明TVDImean可以很好的體現(xiàn)年干旱情況，即TVDImean越高，年平均標(biāo)準(zhǔn)化降水指數(shù)越小，越干旱，且氣象干旱也是造成農(nóng)業(yè)干旱的主要原因。

3.2 黑龍江省干旱時(shí)間特征分析

研究全省20 a來干旱時(shí)間特征變化，TVDImean時(shí)序數(shù)據(jù)(圖1)顯示，黑龍江省2006年是省內(nèi)干旱趨勢(shì)變化轉(zhuǎn)折年，干旱發(fā)生較重的在2000—2006年，此時(shí)間段內(nèi)TVDImean均值為0.63；到2013年干旱水平降到最低，2007—2013年此時(shí)間段內(nèi)TVDImean均值為0.57；2015年TVDImean值雖有所升高，但2014—2019年時(shí)間段內(nèi)TVDImean均值依然有所降低，為0.56。3個(gè)時(shí)段內(nèi)TVDImean均值呈下降趨勢(shì)。對(duì)2000—2019年TVDImean值與時(shí)間進(jìn)行回歸分析，表明TVDImean值20 a間總體呈減弱的趨勢(shì)，經(jīng)過顯著性分析得出p<0.05，表明黑龍江省年平均干旱有顯著減弱趨勢(shì)。

TVDImax監(jiān)測(cè)(圖1)出黑龍江省2006年、2009年和2015年極端干旱情況最重(TVDImax最大)，且分布在3個(gè)不同干旱時(shí)間內(nèi)，TVDImean與TVDImax關(guān)系不大，從變化趨勢(shì)來看黑龍江省近20 a的TVDImax有減弱的趨勢(shì)，但未通過顯著性分析。

圖1 2000-2019年TVDI變化特征

分析2000—2019年TVDImean中干旱面積占比情況，由圖2可以看出，黑龍江省年干旱類型中，無旱與特旱面積占比較少，主要集中在輕旱、中旱和重上，其中中旱面積占比最大，表明黑龍江省年干旱主要以中旱為主。輕旱和無旱面積占比呈增加趨勢(shì)，其他干旱等級(jí)面積占比呈減少趨勢(shì)，該趨勢(shì)與TVDImean呈減弱趨勢(shì)相一致。

圖2 2000-2019年黑龍江省不同等級(jí)干旱面積占比

分析5—9月20 a平均干旱面積占比的月動(dòng)態(tài)情況，由圖3可知，黑龍江省無旱面積占比呈先增加后降低現(xiàn)象，在7月無旱面積占比達(dá)到峰值，8月無旱面積占比劇降其比例小于5月和6月，9月無旱面積占比降到最低；輕旱面積占比也呈先增加后降低現(xiàn)象，7月輕旱面積占比達(dá)到峰值，8月輕旱面積占比較7月減少但高于6月，9月輕旱面積占比降到最低；中旱面積占比呈先降低后升高現(xiàn)象，7月中旱面積占比降到最低，8月中旱面積占比逆轉(zhuǎn)升高比例高于5月和6月，9月中旱面積占比達(dá)到峰值；重旱面積占比6月份最高，5月重旱面積占比僅次于6月居第二，9月重旱面積占比小于6月和5月，7月和8月重旱面積占比特小；特旱面積占比呈降低現(xiàn)象5月最大，7月和8月特旱面積占比最小，幾乎接近0值，9月特旱面積占比有所增加但仍小于6月。中旱、重旱、特旱3個(gè)等級(jí)干旱面積占比總和由大到小：9月>5月>6月>8月>7月，其中雖然9月3個(gè)等級(jí)干旱面積占比總量最大，但以中旱為主，特旱特少，重旱占比也比5月和6月少，中旱本身級(jí)別不高，9月又屬秋收季節(jié)，因此對(duì)作物影響不是很大；5月雖不是3個(gè)等級(jí)干旱面積占比總量最大，但其重旱、特旱面積占比均較大。因此分析表明5月是黑龍江省干旱最嚴(yán)重期，6月是干旱期，7月是黑龍江省干旱最弱時(shí)期，8月是干旱較弱時(shí)期，9月是干旱期以中旱為主。

圖3 2000-2019年黑龍江省5-9月干旱面積占比

3.3 黑龍江省干旱空間特征分析

根據(jù)定義TVDImean值中各種干旱類型的干旱強(qiáng)度：無旱干旱強(qiáng)度為1，輕旱干旱強(qiáng)度為2，中旱干旱強(qiáng)度為3，重旱干旱強(qiáng)度為4，特旱干旱強(qiáng)度為5，利用公式(5)計(jì)算I20，本文假設(shè)20 a同一像元干旱強(qiáng)度都為4，那80就是重旱和中旱的臨界值，因此認(rèn)為80≤I20的像元點(diǎn)為干旱頻率高且干旱強(qiáng)度大的區(qū)域稱為重風(fēng)險(xiǎn)區(qū)，主要分布在訥河、甘南縣、市齊齊哈爾市、龍江縣、林甸縣、泰來縣、大慶市、肇源縣和杜爾伯特縣。假設(shè)20 a同一點(diǎn)干旱強(qiáng)度都為3，那60就是中旱和輕旱的臨界值，因此認(rèn)為60≤I20<80的像元點(diǎn)為干旱頻率中干旱強(qiáng)度中區(qū)域稱為中風(fēng)險(xiǎn)區(qū)，主要分布在西部地區(qū)、寶清縣、富錦市、友誼縣、集賢縣、勃利縣和寧安市干旱頻率中強(qiáng)度中。I20<60的像元點(diǎn)為干旱頻率輕干旱強(qiáng)度輕區(qū)域稱為低風(fēng)險(xiǎn)區(qū)，主要分布在北部、中部和東北部地區(qū)(圖4)。

分析5—9月逐月干旱空間特征情況，5月重風(fēng)險(xiǎn)區(qū)內(nèi)含縣市較多，主要分布在西部地區(qū)、寧安市、勃利縣、雞東縣、友誼縣、寶清縣和集賢縣，表明該地區(qū)每年5月份發(fā)生重旱、特旱的概率大；中風(fēng)險(xiǎn)區(qū)主要分布在北部、中部；干旱頻率輕干旱強(qiáng)度輕區(qū)域主要分布在東部地區(qū)。6月重風(fēng)險(xiǎn)區(qū)明顯減少，主要分布在西部部分地區(qū)；中風(fēng)險(xiǎn)區(qū)域增加，增加主要為5月份的重風(fēng)險(xiǎn)區(qū)，也說明6月份較5月份干旱有明顯減弱；低風(fēng)險(xiǎn)區(qū)分布在北部和東部地區(qū)。7月、8月份重風(fēng)險(xiǎn)區(qū)幾乎沒有，此結(jié)果與上文時(shí)間特征結(jié)果相吻合；中風(fēng)險(xiǎn)區(qū)零星分布在西部地區(qū)，該區(qū)面積8月份大于7月份；低風(fēng)險(xiǎn)區(qū)7月份幾乎覆蓋全省(除齊齊哈爾市和大慶市附近外)，8月份除西南部外的全部區(qū)域，表明7月、8月較少發(fā)生干旱。9月重風(fēng)險(xiǎn)區(qū)主要分布在大慶市周邊及泰來縣、龍江縣、甘南縣、齊齊哈爾市、訥河市等地；中風(fēng)險(xiǎn)區(qū)分布在西部地區(qū)和東南部地區(qū)；低風(fēng)險(xiǎn)區(qū)主要分布在北部和東北部(圖4)。

圖4 黑龍江省20 a及逐月干旱空間特征分布

4 討 論

作者較前研究[11]已對(duì)溫度植被干旱指數(shù)干旱監(jiān)測(cè)指標(biāo)的干旱強(qiáng)度進(jìn)行了劃分，且該劃分結(jié)果也在黑龍江省進(jìn)行應(yīng)用與驗(yàn)證[2,17]證明了劃分的TVDI干旱強(qiáng)度對(duì)黑龍江省土壤干旱情況的反演具有較好的適用性。本文就是基于此對(duì)黑龍江省開展2000—2019年的干旱時(shí)空變化特征研究，張劍俠等[18]以氣象數(shù)據(jù)和統(tǒng)計(jì)數(shù)據(jù)計(jì)算月干旱指數(shù)，統(tǒng)計(jì)30 a(1975—2005)內(nèi)月干旱發(fā)生的次數(shù)，區(qū)域劃分干旱多發(fā)區(qū)、干旱次發(fā)區(qū)和較少發(fā)生干旱區(qū)，區(qū)域劃分結(jié)果與本研究結(jié)果具有一致性?；赥VDI的干旱空間特征分析是以像元為基礎(chǔ)單元，屬于面狀數(shù)據(jù)，對(duì)于全省研究更具有代表性，而傳統(tǒng)的氣象數(shù)據(jù)和統(tǒng)計(jì)數(shù)據(jù)均是以氣象監(jiān)測(cè)站和縣為基礎(chǔ)單元，屬于點(diǎn)狀數(shù)據(jù)，對(duì)全省研究屬以點(diǎn)概面，因此TVDI監(jiān)測(cè)更具優(yōu)勢(shì)。

TVDI的構(gòu)建是基于溫度—植被的理想假設(shè)，即在相同植被覆蓋度的前提下，地表溫度越高，土壤水分越低，這一區(qū)域就越干旱，反之就越濕潤。本文研究的黑龍江省耕地范圍生長(zhǎng)季內(nèi)的干旱情況，僅從區(qū)域角度和時(shí)間跨度上進(jìn)行分析，并未考慮到地面附著植被類型，由于植被類型的不同，植株高矮存在一定的差異，對(duì)TVDI的計(jì)算會(huì)有一定的影響，因此后續(xù)作者會(huì)從植被覆蓋類型出發(fā)，研究不同植被類型、不同時(shí)期、不同區(qū)域?qū)Ω珊档捻憫?yīng)程度的研究將更有意義。

5 結(jié) 論

(1) 從時(shí)間尺度上來分析，2000—2019年黑龍江省20 a間2000—2006年干旱嚴(yán)重、2007—2013年干旱減弱、2014—2019年干旱震蕩減弱；TVDImax監(jiān)測(cè)出黑龍江省2006年、2009年、2015年三年中極端干旱情況最重；20 a間TVDImean有顯著減弱趨勢(shì)而TVDImax未通過顯著性檢驗(yàn)。分析20 a間TVDImean中干旱面積占比可知，全省中旱現(xiàn)象較為普遍。分析月內(nèi)干旱情況可知，5月是黑龍江省干旱發(fā)生的最嚴(yán)重時(shí)期，6月是干旱期，7月是黑龍江省干旱最弱時(shí)期，8月是干旱較弱時(shí)期，9月也是干旱期以中旱為主。

(2) 從空間尺度上來分析，20 a重風(fēng)險(xiǎn)區(qū)主要分布在訥河市、甘南縣、齊齊哈爾市、龍江縣、林甸縣、泰來縣、大慶市、肇源縣和杜爾伯特縣。中風(fēng)險(xiǎn)區(qū)主要分布在西部地區(qū)、寶清縣、富錦市、友誼縣、集賢縣、勃利縣和寧安市。低風(fēng)險(xiǎn)區(qū)主要分布在北部、中部和東北部地區(qū)。分析逐月干旱空間特征：西部地區(qū)、寧安市、勃利縣、雞東縣、友誼縣、寶清縣和集賢縣等地區(qū)5月干旱發(fā)生頻率高強(qiáng)度大，西部部分地區(qū)6月持續(xù)高頻率高強(qiáng)度干旱，大慶市周邊及泰來縣、龍江縣、甘南縣、齊齊哈爾市、訥河市等地區(qū)9月干旱發(fā)生頻率高強(qiáng)度大。