基于同化數(shù)據(jù)的標準化土壤濕度指數(shù)監(jiān)測農(nóng)業(yè)干旱的適宜性研究

周洪奎,武建軍,*,李小涵,劉雷震,楊建華,韓忻憶

1 北京師范大學環(huán)境演變與自然災害教育部重點實驗室,北京 1008752 北京師范大學地理科學學部,北京 100875

干旱是一種多發(fā)的極端氣候事件,同時也是最具破壞性的自然災害之一。頻繁發(fā)生的干旱災害給農(nóng)業(yè)生產(chǎn)帶來嚴重影響,威脅糧食安全[1- 2]。黃淮海平原是我國最重要的糧食生產(chǎn)基地之一,準確、有效地監(jiān)測黃淮海平原的農(nóng)業(yè)干旱狀況對于及時制定區(qū)域的減災策略,降低災害損失具有重要意義。

農(nóng)業(yè)干旱是作物生長過程中因土壤水分不足,阻礙作物正常生長,影響糧食產(chǎn)量的水量供需不平衡現(xiàn)象[3- 4]。在農(nóng)業(yè)干旱監(jiān)測中,土壤濕度扮演著十分重要的角色[5- 7]。當土壤水分低于土壤田間持水量時,會產(chǎn)生水分虧缺,植被開始不處于最佳生長狀態(tài)；當水分進一步降低,即會產(chǎn)生水分脅迫,影響植被生長,進而影響植被生產(chǎn)力。因此,利用土壤濕度進行農(nóng)業(yè)干旱監(jiān)測也具有農(nóng)學意義。

近年來,已經(jīng)有學者提出了一些基于土壤濕度的農(nóng)業(yè)干旱監(jiān)測指數(shù),主要分為兩類：(1)基于歷史時間序列數(shù)據(jù),確定當前狀態(tài)相對于歷史正常范圍的偏離程度建立的干旱指數(shù),如標準化土壤濕度指數(shù)(SSMI,standardized soil moisture index、土壤濕度距平(SMA,soil moisture anomaly)、土壤濕度百分位數(shù)(SMP,soil moisture percentile)等。該類指數(shù)僅利用土壤濕度數(shù)據(jù),無需獲取其他數(shù)據(jù)源,并且不同地區(qū)的監(jiān)測結果具有可比性,能夠用于區(qū)域干旱監(jiān)測與評估。Mishra等[8]基于歷史時間序列數(shù)據(jù),利用標準化土壤濕度指數(shù)詳細剖析了局地尺度的農(nóng)業(yè)干旱。(2)從土壤可利用水的角度,基于土壤濕度和土壤特性參數(shù)構建的干旱指數(shù),如土壤濕度指數(shù)(SMI,soil moisture index)[9- 10]、土壤水分虧缺指數(shù)(SWDI,soil water deficit index)[5]等。該類指數(shù)的優(yōu)點在于從土壤可利用水的角度出發(fā),考慮了不同地區(qū)土壤性質(zhì)對于水分虧缺的影響差異,但對于區(qū)域農(nóng)業(yè)干旱監(jiān)測而言,準確獲取不同地區(qū)土壤屬性參數(shù)是十分困難的,因而在區(qū)域或更大尺度上的應用受到一定的限制。

標準化土壤濕度指數(shù)(SSMI)是基于歷史土壤濕度時間序列構建的一種農(nóng)業(yè)干旱指數(shù),具有計算簡單易行,考慮數(shù)據(jù)分布特征等優(yōu)點,探討其在區(qū)域農(nóng)業(yè)干旱監(jiān)測中的適宜性能夠為區(qū)域業(yè)務化的農(nóng)業(yè)干旱監(jiān)測以及干旱影響評估提供基礎。然而,由于根區(qū)土壤濕度獲取的不易性以及土壤濕度的估算精度不高等問題,目前該指數(shù)在農(nóng)業(yè)干旱監(jiān)測中的適宜性研究還十分缺乏。此外,現(xiàn)有的農(nóng)業(yè)干旱指數(shù)評價研究中,主要通過與其他常用干旱指數(shù)或者氣象要素的對比來評價[5,11],考慮干旱災害記錄以及干旱影響的評價研究較少。實際上,我國的主要農(nóng)業(yè)區(qū)擁有比較完善的農(nóng)業(yè)氣象觀測網(wǎng)絡,記錄了作物生長發(fā)育狀況和農(nóng)業(yè)氣象災害發(fā)生情況,這對于評價農(nóng)業(yè)干旱指數(shù)的適宜性是十分有用的數(shù)據(jù)源。因此,本文采用數(shù)據(jù)同化后的根區(qū)土壤濕度數(shù)據(jù)構建標準化土壤濕度指數(shù)來監(jiān)測黃淮海平原的農(nóng)業(yè)干旱狀況,通過與標準化降水蒸散指數(shù)(SPEI,standardized precipitation-evapotranspiration index)、農(nóng)業(yè)干旱災害數(shù)據(jù)的對比以及與冬小麥產(chǎn)量的關系分析,綜合評價其農(nóng)業(yè)干旱監(jiān)測的適宜性,以期為黃淮海平原農(nóng)業(yè)干旱監(jiān)測業(yè)務化運行以及防災減災策略的制定提供理論依據(jù)和科學參考。

1 研究區(qū)概況

黃淮海平原是我國重要的糧食主產(chǎn)區(qū)且干旱災害頻發(fā),已有許多學者以此為研究區(qū)開展干旱相關研究[12- 14]。黃淮海平原地處我國北方,位于32°—41°N和112°—123°E之間,面積約為39萬km2。在行政區(qū)劃上包括天津、山東以及北京南部、河北大部、河南大部、安徽和江蘇北部地區(qū)(圖1)。根據(jù)我國農(nóng)業(yè)綜合區(qū)劃圖,黃淮海平原可分為燕山太行山山麓平原區(qū)、冀魯豫低洼平原區(qū)、黃淮平原區(qū)和山東丘陵農(nóng)林區(qū)。

圖1 研究區(qū)位置示意圖Fig.1 Location of the study area

冬小麥和夏玉米的輪作是該地區(qū)主要種植模式,冬小麥生長季在區(qū)域南北部稍有差異,主要集中在每年9月下旬/10月中上旬—次年5月下旬/6月中上旬；夏玉米生長季則為6月中上旬—9月中下旬。年降水量差異較大,在500—950 mm之間,且60%—70%降水集中在夏季。因此,相對而言,冬小麥更容易受到干旱威脅[14]。根據(jù)中國氣象局農(nóng)氣站點記錄的農(nóng)業(yè)氣象災害中,黃淮海平原經(jīng)常發(fā)生不同程度和持續(xù)時間的干旱事件,在本文研究時段(2002—2010年)內(nèi),2002、2004、2006年發(fā)生了較為嚴重的干旱事件。

2 資料與方法

2.1 數(shù)據(jù)來源與預處理

本文所用到的數(shù)據(jù)包括根區(qū)土壤濕度同化數(shù)據(jù)、標準化降水蒸散指數(shù)數(shù)據(jù)、農(nóng)業(yè)氣象站點災害數(shù)據(jù)和冬小麥產(chǎn)量數(shù)據(jù)。

(1)根區(qū)土壤濕度同化數(shù)據(jù)

本文所使用的土壤濕度數(shù)據(jù)是將LPRM模型反演的微波土壤濕度同化到雙層帕默爾水量平衡模型中得到的根區(qū)土壤濕度數(shù)據(jù)[15- 16]。已有研究表明,根區(qū)土壤濕度相比于某一固定深度的土壤濕度數(shù)據(jù)更適合于農(nóng)業(yè)干旱監(jiān)測[8]。該根區(qū)土壤濕度數(shù)據(jù)代表的是土壤深度為1 m內(nèi)的土壤水分狀況,該深度能夠滿足大部分作物根部水分、營養(yǎng)吸收的需求[17]。該數(shù)據(jù)覆蓋時段為2002年6月至2010年12月,分辨率為0.25度的逐日數(shù)據(jù),數(shù)據(jù)格式為NetCDF格式,可從Reverb平臺下載(http://reverb.echo.nasa.gov/reverb/)。數(shù)據(jù)預處理過程為：首先利用計算機程序?qū)etCDF格式數(shù)據(jù)轉(zhuǎn)換為TIFF格式,然后對數(shù)據(jù)無效值(負值)進行剔除,最后利用GIS軟件裁剪獲得黃淮海平原根區(qū)土壤濕度數(shù)據(jù)。

(2)標準化降水蒸散指數(shù)數(shù)據(jù)

標準化降水蒸散指數(shù)(SPEI)是Vicente-Serrano等[18]提出的一種應用非常廣泛的干旱指數(shù)。SPEI是在考慮降水的基礎上,加入了溫度信息,利用降水與潛在蒸散發(fā)之間的水量平衡關系來表征干旱。SPEI具有多尺度特征,研究表明,3個月尺度SPEI(SPEI- 3)適合監(jiān)測農(nóng)業(yè)干旱[12,19],因此,文中將SPEI- 3與SSMI進行對比分析,作為評價SSMI監(jiān)測農(nóng)業(yè)干旱適宜性的一個方面。SPEI的計算方法：(1)獲得某一時間尺度(如3個月)內(nèi)降水累計值與潛在蒸散發(fā)累計值的差值的多年時間序列,一般要求時間序列的長度為30年以上[18]。根據(jù)氣象數(shù)據(jù)的可獲得性,潛在蒸散發(fā)的計算可以采用Penman-Monteith、Thornthwaite、Hargreaves等方法[20- 21]。(2)對該時間序列進行Log-Logistic概率分布擬合,獲得概率分布函數(shù)。(3)最后將概率分布函數(shù)進行正態(tài)標準化得到SPEI值,具體計算公式參考Vicente-Serrano等[18]的文獻。本研究使用的SPEI數(shù)據(jù)集是由提出者Vicente-Serrano等生產(chǎn)的全球SPEI數(shù)據(jù)集[22],可以在以下網(wǎng)址(http://sac.csic.es/spei/index.html)下載。該數(shù)據(jù)格式為NetCDF格式,空間分辨率為0.5度,時間分辨率為每月,數(shù)據(jù)時段為1901年1月—2015年12月,包含1—48月尺度的SPEI數(shù)據(jù)。為了與土壤濕度數(shù)據(jù)時段保持一致,本研究選取黃淮海平原2002年6月至2010年12月的SPEI格網(wǎng)數(shù)據(jù)集。

(3)農(nóng)業(yè)氣象站點災害數(shù)據(jù)

由于大范圍地獲取完整、詳細的災害記錄數(shù)據(jù)是較為困難的,因此,干旱觀測數(shù)據(jù)對于評價干旱指數(shù)的監(jiān)測效果是十分寶貴的數(shù)據(jù)源。我國黃淮海平原擁有比較密集的農(nóng)業(yè)氣象站點,記錄了農(nóng)業(yè)干旱的相關信息。中國氣象局根據(jù)農(nóng)業(yè)氣象臺站上報的農(nóng)業(yè)氣象旬月報報文資料整理形成農(nóng)業(yè)氣象站點災害旬值數(shù)據(jù)集,該數(shù)據(jù)集包括災害名稱、災害發(fā)生日期、災害強度、災害面積以及受害百分比等主要字段,該數(shù)據(jù)可以從中國氣象數(shù)據(jù)網(wǎng)下載(http://data.cma.cn/)。根據(jù)農(nóng)業(yè)干旱國家標準中作物形態(tài)指標,災害記錄中干旱強度分為輕度、中度和重度干旱[23]。本研究在全國農(nóng)氣站點的災害記錄中篩選出了2002—2010年間黃淮海平原共40個農(nóng)氣站點的1244條干旱災害記錄。

(4)冬小麥產(chǎn)量數(shù)據(jù)

本文還收集了黃淮海平原主要省市天津、河北、山東、河南、安徽35個縣級行政單位的1994—2013年冬小麥產(chǎn)量數(shù)據(jù)來評價干旱指數(shù)與作物產(chǎn)量之間的關系。數(shù)據(jù)來源是從各省市統(tǒng)計年鑒資料中獲取。在預處理過程中,對產(chǎn)量數(shù)據(jù)進行了多項式去趨勢處理,從實際產(chǎn)量中去除科技進步等人為因素對冬小麥產(chǎn)量的影響,得到氣候要素為主要影響因素的產(chǎn)量(即氣象產(chǎn)量)。文中用減產(chǎn)率作為農(nóng)業(yè)干旱對作物造成的影響的表達指標,減產(chǎn)率(YLR,yield loss ratio)可通過公式(1)計算：

(1)

2.2 標準化土壤濕度指數(shù)(SSMI)

本文選用標準化土壤濕度指數(shù)作為評價農(nóng)業(yè)干旱的指標,對根區(qū)土壤濕度進行標準化。對數(shù)據(jù)標準化的過程中,需確定土壤濕度數(shù)據(jù)的概率分布。常用于干旱指數(shù)構建的概率分布函數(shù)包括伽瑪分布、皮爾遜III型分布、經(jīng)驗累積概率分布、正態(tài)分布、Log-Logistic分布等[24]。利用Kolmogorov-Smirnov方法檢驗,結果顯示土壤濕度數(shù)據(jù)符合正態(tài)分布,與Mishra等的結果相同[8]。因此,SSMI的構建方法如下：

(2)

表1 標準化土壤濕度指數(shù)(SSMI)干旱等級劃分

3 結果與分析

3.1 基于SSMI的黃淮海平原干旱監(jiān)測結果

根據(jù)2.2節(jié)描述的SSMI的計算方法,利用根區(qū)土壤濕度數(shù)據(jù)計算得到了2002—2010年間黃淮海平原SSMI的干旱監(jiān)測結果。圖2展示了黃淮海平原內(nèi)所有柵格數(shù)據(jù)SSMI月均值數(shù)據(jù)序列。根據(jù)SSMI的干旱等級劃分,黃淮海平原在2002—2010年間共發(fā)生12次干旱事件,其中達到重度干旱事件1次(2002.6—2003.8),中度干旱事件3次(2004.1—2004.7、2005.3—2005.8、2006.8—2007.1),其余為輕度干旱事件。

圖2 黃淮海平原SSMI月均值時間序列Fig.2 Time series of monthly SSMI in the Huang-Huai-Hai Plain

3.2 SSMI與SPEI的對比分析

由于SPEI和SSMI數(shù)據(jù)在時間尺度和空間分辨率上不同,為了對二者進行對比分析,需要將二者整合成統(tǒng)一的時空分辨率上。在空間分辨率上,SPEI數(shù)據(jù)的空間分辨率為0.5度,SSMI數(shù)據(jù)的空間分辨率為0.25度,文中將SSMI數(shù)據(jù)按照雙線性內(nèi)插法重采樣為0.5度數(shù)據(jù)用于對比分析。在時間分辨率上,由于SPEI是逐月數(shù)據(jù),為與SPEI數(shù)據(jù)序列保持一致,采用SSMI月值數(shù)據(jù)用于分析二者之間的關系,SSMI月值數(shù)據(jù)的計算方法見2.2節(jié)。

本研究選用3個月尺度的SPEI(SPEI- 3)來與SSMI進行對比分析。SPEI- 3即綜合考慮當前月份前3個月的降水和蒸散發(fā)量之差,按照上述SPEI的計算方法得到當前月份SPEI值,按照月份順序依次滑動計算得到不同月份的SPEI- 3值,最終形成逐月的SPEI- 3時間序列。首先從整個黃淮海平原尺度上對比分析SSMI和SPEI之間的關系,將黃淮海平原內(nèi)所有柵格數(shù)據(jù)進行取均值處理,得到2002年6月至2010年12月共103個月的SSMI和SPEI時間序列數(shù)據(jù),結果如圖3左圖所示。結果顯示,從整個黃淮海平原尺度來看,SPEI與SSMI具有較好的一致性,相關系數(shù)達到0.63,呈現(xiàn)極顯著的相關關系。在2002、2004、2006年發(fā)生的農(nóng)業(yè)干旱中,二者均具有良好地識別效果。在柵格尺度上,進一步分析SSMI與SPEI的相關程度,2002—2010年月尺度SSMI與SPEI的相關性如圖3右圖所示?？梢园l(fā)現(xiàn),SSMI與SPEI之間的相關系數(shù)在0.3—0.8之間,平均相關系數(shù)為0.52,每個柵格參與相關性計算的樣本數(shù)量為103個,所有柵格均通過0.001顯著性水平檢驗。以上分析結果表明,SSMI與SPEI之間具有極顯著的相關性和較好的一致性,以SPEI作為參照指標,SSMI能夠反映黃淮海平原的干濕狀況。

圖3 黃淮海平原月尺度SSMI與SPEI的對比分析Fig.3 Comparisons of the monthly SSMI and SPEI in the Huang-Huai-Hai Plain

3.3 SSMI與干旱災害觀測記錄的比較

本節(jié)分別從黃淮海平原整體情況和典型農(nóng)氣站點上對干旱災害觀測記錄與SSMI進行對比分析。在分析黃淮海平原整體情況時,由于災害觀測記錄為旬值數(shù)據(jù)集,時間序列較長,為了便于空間上展示,按照逐月方式對干旱災害記錄進行了整理歸納；在典型農(nóng)氣站點尺度上,采用旬尺度數(shù)據(jù)進行對比分析,從相對較小的時間尺度上進一步檢驗SSMI的干旱監(jiān)測效果。表2是根據(jù)農(nóng)氣站點的災害記錄整理出的2002—2010年間的主要干旱事件。從表中可以發(fā)現(xiàn),黃淮海平原每年都會有干旱發(fā)生,僅是在持續(xù)時間和發(fā)生范圍上有所差別,文中重點分析了2004、2006年較長持續(xù)時間的干旱過程。圖4展示了2004和2006年黃淮海平原從干旱發(fā)生到逐步加重直至干旱減輕的時空演變過程。按照農(nóng)氣站點災害記錄,從2004年3月開始,黃淮海平原就已經(jīng)有零星分布的干旱發(fā)生,至4—6月,干旱范圍已經(jīng)擴展為黃淮海平原大部分地區(qū)。在2004年3—6月間,黃淮海平原共有18個站點記錄發(fā)生干旱,平均受害百分比為40%—50%,平均受旱面積超過2.5萬hm2,山東菏澤、安徽蒙城、宿縣、亳州等地受害百分比甚至達到90%—100%,平均受旱面積超過6.7萬hm2。在2004年7月—11月,干旱范圍有所減小,主要集中在安徽、江蘇北部,以及山東、河南部分地區(qū)。在2006年4—12月間,共有30個站點顯示發(fā)生干旱災害,平均受害百分比為47%—57%,平均受旱害面積超過2.8萬hm2。在干旱空間分布來看,在7—9月干旱得到部分緩解；在10—11月,由于降水不足,干旱又逐漸加重；直至12月,黃淮海平原大部分地區(qū)干旱逐漸解除,只在山東膠東半島和河北部分地區(qū)干旱持續(xù)。綜合2004、2006年兩次農(nóng)業(yè)干旱事件,從圖4可以看出,在空間分布上,SSMI的干旱監(jiān)測結果與農(nóng)氣站點災害記錄是基本符合的。

表2 農(nóng)氣站點記錄的黃淮海平原主要農(nóng)業(yè)干旱事件

圖4 2004、2006年黃淮海平原干旱空間演變過程Fig.4 The evolutions of drought in 2004 and 2006 in the Huang-Huai-Hai Plain

為了進一步評估SSMI在農(nóng)業(yè)干旱監(jiān)測中的準確性,從站點尺度上分析干旱指數(shù)監(jiān)測結果與災害記錄的一致性。根據(jù)農(nóng)氣站點災害數(shù)據(jù)的完整性,選取德州、菏澤、泰安、商丘、蒙城、亳州6個站點用于對比分析。圖5中分別用數(shù)值-1、-2、-3來表示農(nóng)氣站點記錄的輕度、中度和重度干旱事件。另外,農(nóng)氣站點還記錄了洪澇災害,其輕、中、重強度分別用1、2、3表示。從圖5中可以發(fā)現(xiàn),各站點SSMI所代表的干旱強度與干旱災害記錄是基本一致的,總體上SSMI能夠較為準確地反映區(qū)域農(nóng)業(yè)干旱強度特征。從干旱監(jiān)測的時效性方面,旬值尺度的SSMI與干旱災害記錄也基本吻合,說明利用SSMI在旬尺度上監(jiān)測農(nóng)業(yè)干旱也具有較大的應用潛力。

圖5 SSMI旬值與農(nóng)氣站點干旱記錄的對比Fig.5 Comparison of the 10-day SSMI and drought records from agro-meteorological sites

3.4 農(nóng)業(yè)干旱對冬小麥產(chǎn)量的影響分析

在進行農(nóng)業(yè)干旱監(jiān)測時,分析干旱指數(shù)與作物產(chǎn)量的關系是檢驗干旱指數(shù)監(jiān)測效果的重要方面。本文搜集到黃淮海平原35個縣市級1994—2013共20年的冬小麥產(chǎn)量數(shù)據(jù)。為了使不同縣市之間的產(chǎn)量數(shù)據(jù)能夠進行類比,本文采用減產(chǎn)率來分析農(nóng)業(yè)干旱監(jiān)測結果與冬小麥產(chǎn)量的關系。首先,計算得到了所有縣市2002—2010年8個生長季SSMI與冬小麥減產(chǎn)率的相關性(圖6左圖)。結果顯示,SSMI與冬小麥減產(chǎn)率呈現(xiàn)極顯著的正相關關系(R=0.68,P<0.001)。圖6右圖展示了各縣市生長季SSMI與減產(chǎn)率之間的相關性,在35個縣市中,有24個可以通過0.05水平顯著性檢驗,11個未通過顯著性檢驗。根據(jù)農(nóng)氣站點災害記錄顯示,河北保定、石家莊、邢臺等地的冬小麥在2002—2003生長季遭遇病蟲害、大風的災害影響,山東菏澤在2003—2004生長季遭受漬害影響,并非都因干旱引起冬小麥減產(chǎn)。因此,總體上,生長季SSMI與冬小麥減產(chǎn)率是顯著相關的,SSMI能夠?qū)r(nóng)業(yè)干旱引起的冬小麥減產(chǎn)起到一定的指示作用。

圖6 冬小麥生長季SSMI與減產(chǎn)率的相關關系Fig.6 Correlations of the SSMI during the growing season and winter wheat yield loss ratio

本文選取2002—2003年發(fā)生在冬小麥生長季內(nèi)的干旱事件,進一步分析SSMI干旱監(jiān)測結果與冬小麥減產(chǎn)率的關系。從圖7中可以發(fā)現(xiàn),黃淮海平原2002—2003生長季冬小麥產(chǎn)量基本都呈現(xiàn)減產(chǎn)狀態(tài),減產(chǎn)幅度在15%以內(nèi)。根據(jù)SSMI識別結果,除了河北保定、石家莊、邢臺(圖中綠色顯示區(qū)域)外,其余區(qū)域均處于干旱狀態(tài)?？傮w上,基于統(tǒng)計數(shù)據(jù)得到的冬小麥發(fā)生減產(chǎn)的區(qū)域與基于SSMI識別的干旱發(fā)生區(qū)域具有較好地一致性。在數(shù)值上,有些區(qū)域減產(chǎn)率和SSMI可能無法很好地一一對應,即減產(chǎn)程度并不完全隨著SSMI的減小而增加。主要原因在于糧食減產(chǎn)可能有多種原因引起,在農(nóng)業(yè)氣象災害范疇內(nèi),除干旱外,還可能由大風、漬害、干熱風、冰雹等引起減產(chǎn)。如圖7中河北保定、石家莊、邢臺(綠色區(qū)域)SSMI顯示為正?；驖駶櫊顟B(tài),而產(chǎn)量卻呈現(xiàn)減產(chǎn)狀況,其主要原因在于受到病蟲害、大風的災害影響。2002—2003生長季內(nèi)保定、石家莊、邢臺的降水量分別為101.7 mm、247.4 mm和177.6 mm,占平均降水量的94%、185%和134%,顯示該生長季降水接近或者多于平均降水量,冬小麥產(chǎn)量減少可能并非受到干旱影響。另外,不同生育期內(nèi)即使出現(xiàn)相同等級的干旱發(fā)生,對作物產(chǎn)量的影響也不盡相同,這也是導致SSMI與減產(chǎn)率無法很好地一一對應的原因之一。

圖7 2002—2003生長季SSMI區(qū)域統(tǒng)計值與冬小麥減產(chǎn)率Fig.7 The zonal SSMI statistics and winter wheat yield loss ratio during the growing season of 2002—2003

4 結論和討論

本文利用數(shù)據(jù)同化后的根區(qū)土壤濕度數(shù)據(jù)構建了標準化土壤濕度指數(shù)(SSMI),通過與常用的干旱指數(shù)SPEI、農(nóng)業(yè)干旱災害數(shù)據(jù)以及冬小麥減產(chǎn)率的對比分析,綜合評價了利用SSMI監(jiān)測農(nóng)業(yè)干旱的適宜性。研究表明,SSMI能夠有效反映黃淮海平原的農(nóng)業(yè)干旱狀況,利用SSMI監(jiān)測農(nóng)業(yè)干旱是適宜的。具體研究結論如下：

(1)黃淮海平原格網(wǎng)尺度SSMI與SPEI之間的平均相關系數(shù)達到0.52,具有極顯著相關關系(P<0.001)。從整個黃淮海平原來看,SSMI與SPEI也具有良好的一致性,能夠準確識別該區(qū)大范圍農(nóng)業(yè)干旱。

(2)從區(qū)域尺度上看,SSMI能夠有效反映干旱發(fā)生、發(fā)展直至減輕的演變過程。通過與德州、菏澤、泰安、商丘、蒙城、亳州6個農(nóng)氣站點干旱災害記錄的對比分析表明：在站點尺度上,利用SSMI識別的農(nóng)業(yè)干旱與農(nóng)氣站點干旱災害記錄是基本一致的,SSMI能夠較為準確地監(jiān)測農(nóng)業(yè)干旱的強度。

(3)冬小麥生長季SSMI與減產(chǎn)率具有顯著的相關性和良好的對應關系,基于SSMI識別的農(nóng)業(yè)干旱發(fā)生區(qū)域與基于統(tǒng)計數(shù)據(jù)得到的冬小麥減產(chǎn)區(qū)域是基本一致的,SSMI能夠?qū)r(nóng)業(yè)干旱引起的冬小麥減產(chǎn)起到一定的指示作用。

相比于氣象干旱,農(nóng)業(yè)干旱是一種更為復雜的干旱類型,與氣象條件和作物生長狀況密切相關[5,8],因而,在評價農(nóng)業(yè)干旱指數(shù)的適宜性時,需要從不同方面綜合評價[27]。在以往的農(nóng)業(yè)干旱指數(shù)評價中,往往通過與其他指數(shù)的對比來評價干旱指數(shù)的監(jiān)測效果[5,11]。由于不同干旱指數(shù)的構建原理不同,僅利用對比分析指數(shù)間的一致性不足以說明干旱指數(shù)的適宜性。本文從干旱指數(shù)、災害觀測記錄對比以及干旱對產(chǎn)量影響的角度綜合評價了標準化土壤濕度指數(shù)(SSMI)監(jiān)測農(nóng)業(yè)干旱的適宜性,在干旱指數(shù)的適宜性評價方法上更加全面,研究結果可以為黃淮海平原農(nóng)業(yè)干旱監(jiān)測業(yè)務化運行以及防災減災策略的制定提供理論依據(jù)和科學參考。

本研究尚存在一些不足之處。由于采用的根區(qū)土壤濕度數(shù)據(jù)是同化模型模擬結果和遙感數(shù)據(jù)所得,土壤濕度數(shù)據(jù)時間序列相對較短。隨著SMOS(soil moisture ocean salinity)、SMAP(soil moisture active passive)等衛(wèi)星土壤濕度產(chǎn)品的增多,空間分辨率也會更高,后續(xù)研究還需采用更長時間序列數(shù)據(jù)用于農(nóng)業(yè)干旱監(jiān)測及影響評價。在探討干旱對作物影響時,還需進一步考慮灌溉等人為因素的影響。本文研究區(qū)黃淮海平原的大部分地區(qū)為灌溉農(nóng)業(yè)區(qū),為了減輕干旱對產(chǎn)量的影響,大部分區(qū)域會進行人為灌溉,所以在發(fā)生干旱時導致產(chǎn)量可能并未降低,在分析干旱對作物生長及產(chǎn)量影響時,人為灌溉則會增加分析結果的不確定性。本文所利用的土壤濕度為同化后的根區(qū)土壤濕度數(shù)據(jù),數(shù)據(jù)同化過程中用到了微波遙感土壤濕度數(shù)據(jù),微波遙感信號可以捕捉到地面灌溉的信息[28],因而該數(shù)據(jù)在一定程度上考慮了人為灌溉的信息。此外,在利用作物產(chǎn)量統(tǒng)計數(shù)據(jù)時,先進行了去趨勢處理,去除了人為因素導致的科技進步對產(chǎn)量的影響,其中也包含人為灌溉因素,因此,本研究在分析干旱對作物產(chǎn)量影響時部分考慮了灌溉因素的影響。即便如此,在分析SSMI與冬小麥減產(chǎn)率之間的關系時,有些區(qū)域仍然出現(xiàn)減產(chǎn)率和呈現(xiàn)的干旱狀況不完全一致的情況。因而,未來還需獲取農(nóng)田灌溉數(shù)據(jù)深入考慮灌溉因素可能帶來的影響。目前研究工作僅在黃淮海平原開展,進一步研究還需在更多區(qū)域(如雨養(yǎng)農(nóng)業(yè)區(qū)、不同氣候區(qū))進行SSMI的驗證評價工作。