吳金亮 單麗佳
摘要:干旱監(jiān)測(cè)的主要方法包括土壤含水量的中子儀測(cè)量法、TDR測(cè)量法、濕度計(jì)法、稱重烘干發(fā)法以及幾種衛(wèi)星遙感監(jiān)測(cè)方法等。對(duì)各類干旱監(jiān)測(cè)方法存在的問題和發(fā)展趨勢(shì)進(jìn)行探討,并提出加強(qiáng)我國干旱遙感監(jiān)測(cè)技術(shù)研究的建議,為我國農(nóng)業(yè)旱情監(jiān)測(cè)提供一定的幫助。
關(guān)鍵詞:土壤墑情;干旱監(jiān)測(cè);遙感
中圖分類號(hào):X171 文獻(xiàn)標(biāo)識(shí)碼:A 文章編號(hào):1674-1161(2014)06-0065-04
干旱是指因天氣久晴無雨或少量降雨、土壤中水分缺失、空氣干燥等而引發(fā)的農(nóng)作物枯死、人畜飲水不足等現(xiàn)象。從天氣影響角度考慮,干旱還應(yīng)包括高溫、熱浪天氣和干熱風(fēng)等。干旱是全世界面臨的重大自然災(zāi)害之一,而中國也是受旱情災(zāi)害影響較大的國家,主要體現(xiàn)在危害中國農(nóng)業(yè)上。我國每年農(nóng)作物平均受旱耕地面積達(dá)0.20億hm2,成災(zāi)耕地面積達(dá)0.08億hm2,作物減產(chǎn)達(dá)100~150億kg,經(jīng)濟(jì)損失達(dá)2 000多億元。農(nóng)業(yè)干旱災(zāi)害面積受年際及季節(jié)變化影響較大。通常情況下,將年總降水量小于200 mm的地區(qū)定義為干旱區(qū)域,年總降水量在200~500 mm之間的地區(qū)定義為半干旱區(qū)域。
目前用于土壤干旱和水分監(jiān)測(cè)的方法主要有常規(guī)地面監(jiān)測(cè)和衛(wèi)星遙感監(jiān)測(cè)兩種。常規(guī)地面監(jiān)測(cè)是利用地面監(jiān)測(cè)站將各監(jiān)測(cè)點(diǎn)形成網(wǎng)絡(luò),實(shí)時(shí)對(duì)地面土壤濕度等情況進(jìn)行監(jiān)測(cè)記錄,常見測(cè)量方法包括中子儀測(cè)量法、TDR測(cè)量法、濕度計(jì)法、稱重烘干發(fā)法等。衛(wèi)星遙感監(jiān)測(cè)則借助于衛(wèi)星傳感器與遙感信息相結(jié)合,從而達(dá)到大面積監(jiān)測(cè)。常用測(cè)量方法有作物缺水指數(shù)和主動(dòng)、被動(dòng)微波法,特征波段空間法,基于多光譜數(shù)據(jù)信息的熱慣量方法,以及垂直干旱指數(shù)法等。
2 干旱監(jiān)測(cè)方法
2.1 傳統(tǒng)的土壤含水量測(cè)量方法
傳統(tǒng)土壤含水量測(cè)量方法主要包括野外烘干稱重法、光學(xué)測(cè)量法、張力計(jì)土壤法、駐波比法、中子法、電阻法、γ-射線法、時(shí)域反射法、TDR法等。我國大部分氣象站測(cè)定土壤水分采用的是烘干法,少部分地區(qū)采用中子法、TDR法等自動(dòng)測(cè)站。
野外烘干稱重法是直接測(cè)量土壤水分含量的有效方法,在數(shù)據(jù)精度上優(yōu)于其他方法,但費(fèi)時(shí)費(fèi)力,時(shí)效性較差,且連續(xù)取樣會(huì)損壞當(dāng)?shù)赝寥澜Y(jié)構(gòu)?;谳椛湓淼闹凶臃ê挺?射線法因其危害人體健康而被國外棄用,在我國也很少使用。張力計(jì)土壤法的取樣對(duì)象為基質(zhì)勢(shì),與其他測(cè)定土壤含水量的方法相結(jié)合,可測(cè)定土壤水分特征曲線。電阻法的標(biāo)定過程較為繁瑣,且隨著時(shí)間的推移其結(jié)果很快失效,加上測(cè)定精度較低,測(cè)量范圍也存在一定的限制,因此基本上已被淘汰。光學(xué)測(cè)量法具有非接觸的優(yōu)勢(shì),但適應(yīng)性低,易受土壤變異性影響,應(yīng)用前景并不樂觀。TDR的優(yōu)點(diǎn)是測(cè)量速度快、操作簡(jiǎn)便、精確度高(能達(dá)到0.5%)、可連續(xù)測(cè)量,既可測(cè)量土壤表層水分、也可用于測(cè)量剖面水分,既可用于手持式的實(shí)時(shí)測(cè)量、也可用于遠(yuǎn)距離多點(diǎn)自動(dòng)監(jiān)測(cè),測(cè)量數(shù)據(jù)易處理。近20 a來,基于測(cè)量土壤介電常數(shù)的方法與其他方法相比,在測(cè)量精度和實(shí)時(shí)性上更具優(yōu)勢(shì),在操作上更加方便靈活,已成為國內(nèi)外廣泛使用的一種土壤水分測(cè)量方法。
2.2 遙感干旱監(jiān)測(cè)原理及方法
遙感影像的植被相關(guān)信息主要通過植被冠層與葉片光譜特性及差異變化來獲得,不同光譜特性波段所讀取的植被信息與植被不同波段信息或某些特征狀態(tài)存在不同的相關(guān)性。研究表明,在可見光波段與近紅外波段,土壤含水量與其反照率密切關(guān)聯(lián),兩者協(xié)同變化,水量較高的濕土比干土的反照率高。為此,可通過遙感影像在可見光和近紅外波段的反照率差異來計(jì)算土壤水分。此外,地表反照率還與土質(zhì)種類、地表紋理、地表起伏、遙感影像太陽入射角以及植被、作物覆蓋程度等因素相關(guān)。通常采用幾個(gè)波段值組合的形式來反演土壤水分,這樣可以去除土質(zhì)類型不同等因素的干擾。
微波遙感是指借助傳感器來接收微波信號(hào),這些微波信號(hào)主要來自地面上各類物體的發(fā)射或反射作用,并從信號(hào)中研究出有價(jià)值的內(nèi)容。利用遙感進(jìn)行土壤水分監(jiān)測(cè)時(shí),微波對(duì)水成分存在明顯的感應(yīng)能力,這就使遙感監(jiān)測(cè)成為測(cè)定土壤含水量的較為合理的方式(包含主動(dòng)方式和被動(dòng)方式)。主動(dòng)微波遙感主要測(cè)量土壤的后散射參量,而土壤粗糙度和介電常數(shù)可以決定土壤后散射參量,土壤濕度決定介電常數(shù),因此土壤墑情可以較好地反演出來。被動(dòng)式微波遙感的傳感器接收來自地面物體自身輻射的微波,而這些微波土壤濕度由計(jì)算出的土壤亮溫來反演得到,介電常量和熱量共同評(píng)測(cè)出土壤亮溫,介電常量與溫度和土壤含水量關(guān)聯(lián),因此含水量可借助于土壤亮溫計(jì)算得出[1-3]。
2.2.1 物缺水指數(shù)法 1981年Jackson提出了作物缺水指數(shù)(CWSI),其基本原理是通過計(jì)算出來的臨近空氣溫度與植物冠層溫度間的差值來估算監(jiān)測(cè)土壤含水量。作物缺水指數(shù)介于0~1之間,指數(shù)的增加可以反映出水分的脅迫程度逐漸增加,指數(shù)等于1時(shí)為最大。
作物缺水指數(shù)法的最初假設(shè)是冠層溫度值可以測(cè)量得到,地表作物完全覆蓋。但這種假設(shè)被現(xiàn)實(shí)否定了,原因在于地表裸露會(huì)導(dǎo)致虛擬干旱。Moran又提出植被指數(shù)溫度梯形模型理論,即在二維特征空間中,冠氣溫差與植被覆蓋是以一種梯形分布而存在的,4種極端干旱情況可用梯形的4個(gè)頂點(diǎn)來表示:裸土在干旱狀態(tài)下、裸土在飽和水分狀態(tài)下、遭水分脅迫時(shí)植被完全覆蓋、水分供給充足時(shí)的植被完全覆蓋。在這個(gè)定理基礎(chǔ)上有人又提出水分虧缺指數(shù)(WDI),其定義公式為:
WDI=1-Ed/Ep
式中:Ed為白天實(shí)際蒸散發(fā)量;Ep為潛在蒸散發(fā)量。為了得到WDI值,首先要通過遙感熱紅外數(shù)據(jù)來反演地表作物的植被輻射率、冠層溫度、臨近冠層的氣溫等因子,其次再推算實(shí)際的蒸騰蒸發(fā)量,利用彭曼公式計(jì)算出潛在蒸散發(fā)能力的強(qiáng)弱,這樣,作物缺水指數(shù)就可以計(jì)算得出。但在實(shí)際應(yīng)用中,需要較多的氣象實(shí)測(cè)數(shù)據(jù)才能計(jì)算得到潛在蒸騰蒸發(fā)量,因此,水分虧缺指數(shù)的計(jì)算存在較大難度。
2.2.2 微波遙感法 微波遙感監(jiān)測(cè)水分的理論基礎(chǔ)是干燥土壤與水的介電常數(shù)具有強(qiáng)烈反差。土壤含水量直接影響土壤的介電特性,使雷達(dá)回波對(duì)土壤濕度極為敏感。微波遙感監(jiān)測(cè)土壤水分精度較高,且具備全天候和全天時(shí)觀測(cè)的特點(diǎn),不受光照條件限制,可穿透云層,長(zhǎng)波段微波還能夠穿透植被并對(duì)土壤具有一定的穿透能力。微波遙感監(jiān)測(cè)干旱主要有被動(dòng)微波法和主動(dòng)微波法兩種。對(duì)于被動(dòng)微波法,在植被覆蓋區(qū),微波輻射測(cè)量土壤濕度的有效采樣深度約為5~10 cm,且選擇較長(zhǎng)的波躍更具優(yōu)越性;主動(dòng)微波法以應(yīng)用x波段側(cè)視雷達(dá)為主,方法主要是后向散射系數(shù)法(Back Scattering Coefficient)。用微波遙感研究土壤水分始于20世紀(jì)70年代,現(xiàn)雷達(dá)遙感技術(shù)在土壤水分方面的研究已較為成熟。李杏朝根據(jù)微波后向散射系數(shù)法,用微波遙感監(jiān)測(cè)土壤水分的相對(duì)誤差率僅12%。Rajat Bindlish等(2001)在積分模型基礎(chǔ)上通過改進(jìn),將實(shí)際測(cè)得的土壤水分與雷達(dá)獲取的數(shù)據(jù)相關(guān)系數(shù)由0.84提高到0.95。Moeremans等(2000)利用衛(wèi)星和雷達(dá)遙感監(jiān)測(cè)田間和區(qū)域兩個(gè)不同尺度的土壤含水量,認(rèn)為在裸地或植被稀疏的地區(qū)近地表土壤含水量與后向散射系數(shù)有很高的相關(guān)性。Tansey等(2001)也得到了相同的結(jié)論,同時(shí)認(rèn)為地表的粗糙度對(duì)土壤水分的監(jiān)測(cè)有很大影響。施建成等(2002)利用目標(biāo)分解技術(shù)方法分離出植被散射部分,在一定程度上克服了植被層和表面粗糙度的變化對(duì)雷達(dá)后向散射對(duì)土壤水分的影響[4]。微波遙感監(jiān)測(cè)土壤水分盡管受地表參數(shù)影響較大,但估算土壤水分精度較高,可全天候使用,不受云的干擾,若綜合其他可見光與近紅外圖像,將是監(jiān)測(cè)土壤水分最有希望的方法。但要進(jìn)入實(shí)用階段還有許多問題需要解決,如怎樣選擇有關(guān)參數(shù),如何更好地消除植被、壟向、地表粗糙度、土壤質(zhì)地等因素的影響?,F(xiàn)有大多數(shù)模型都是經(jīng)驗(yàn)的線性統(tǒng)計(jì)模型,不具有普遍適應(yīng)性。另外,雷達(dá)遙感數(shù)據(jù)的獲取途徑較少,價(jià)格相對(duì)較高,這在一定程度上限制了其在監(jiān)測(cè)土壤水分上的應(yīng)用。
2.2.3 Ts-NDVI特征空間法 Carlson,Price和Moran等得出在二維特征空間中Ts-NDVI分布是三角形或梯形,不同學(xué)者據(jù)此提出了很多監(jiān)測(cè)指數(shù)。Inge Sanholt等提出了條件溫度植被指數(shù)(TVDI),計(jì)算公式如下:
濕邊方程是指歸一化植被指數(shù)等于某一特定值時(shí)的最小值溫度,干邊方程是指當(dāng)歸一化植被指數(shù)值等于某一特定值時(shí)地表的最大溫度。a,b,a',b'是4個(gè)待定系數(shù)。通常根據(jù)研究區(qū)域的地表溫度和歸一化植被指數(shù)二維散點(diǎn)圖來近似提取干、濕邊方程。
該指數(shù)對(duì)研究區(qū)域有如下要求:區(qū)域內(nèi)有植被處于凋萎狀態(tài)時(shí)的象元,此時(shí)土壤水分含量是比較低的,同時(shí)也存在促使植被良好生長(zhǎng)發(fā)育的象元,此時(shí)土壤水量較大,通常用于表示年內(nèi)某地區(qū)的干旱等級(jí)[5-6]。
2.2.4 土壤熱慣量法 當(dāng)土壤含水量較多時(shí),會(huì)產(chǎn)生更大的熱慣量,兩者具有一定的關(guān)聯(lián)性。熱慣量值的計(jì)算式如下:
通常計(jì)算時(shí)無法直接利用遙感數(shù)據(jù)計(jì)算以上參數(shù)。為此,Price等提出了表觀熱慣量(ATI)定義,可以利用遙感影像數(shù)據(jù)的熱紅外輻射溫度差、近紅外反射率和可見光數(shù)據(jù),計(jì)算熱慣量值。計(jì)算公式如下:
通過遙感可獲取地表的有關(guān)信息,例如MODIS遙感影像采用的是近紅外波段(第1通道、第2通道、第3通道、第4通道)和第31熱紅外通道值,來反演推算大氣透過率、地表溫度、反射率等參數(shù),進(jìn)而計(jì)算土壤表觀熱慣量(ATI),再通過表觀熱慣量與土壤水分量之間的擬合關(guān)系,建立干旱模型,反演土壤墑情變化程度和干旱程度。該原理僅適用于植被覆蓋度不是很高的情況,對(duì)于植被覆蓋度較高的區(qū)域,其監(jiān)測(cè)結(jié)果會(huì)明顯降低[7-8]。
2.2.5 垂直干旱指數(shù)法 植物葉表面對(duì)紅光和藍(lán)紫光的吸收能力強(qiáng),對(duì)近紅外波段具有較高的反射表現(xiàn),因此,從紅波段(Red)到近紅外波段(Nir),裸地反射率值較高但其增幅卻不明顯。植被對(duì)紅光波段的吸收很快就能達(dá)到飽和狀態(tài),所以植被的增加可通過近紅外反射的增加而得到反映。秦明和詹志明等在分析前人研究以及紅光和近紅外特征空間后,提出了基于土壤線斜率的垂直干旱指數(shù),計(jì)算公式如下:
該指數(shù)假想了一條通過原點(diǎn)的直線L與土壤線垂直,其二維特征空間中任意一個(gè)點(diǎn)到直線L的距離都可表示為水分的分布狀況。水分含量越高代表距離越近,水分含量越少則代表越遠(yuǎn)。該方法的優(yōu)點(diǎn)是用光譜曲線特征直接代替了地表溫度和反射率的反演,簡(jiǎn)便實(shí)效。
秦明等研究指出,用垂直干旱指數(shù)監(jiān)測(cè)裸地地表水分效果良好,但遇到植被覆蓋度較高的田間,垂直干旱指數(shù)監(jiān)測(cè)會(huì)受到很大影響。不同土壤類型下地形特征會(huì)分布成非平面,為了改進(jìn)模型,阿布都瓦斯提等提出了修正的垂直千旱指數(shù)法 (MPDI),定義為:
植被地面和部分植被覆蓋程度是植被冠層的兩個(gè)重要參量因子,前兩者與葉片面積指數(shù)相關(guān)性很大。植被光譜指數(shù)可以計(jì)算出地表的植被覆蓋度。例如,歸一化植被指數(shù)(NDVI)和植被覆蓋度關(guān)聯(lián)性非常高,國內(nèi)部分學(xué)者提出對(duì)植被覆蓋度與歸一化植被指數(shù)進(jìn)行冪次方變換和均方根可提升植被覆蓋度的精度,還有研究指出使用歸一化植被指數(shù)(NDVI)計(jì)算出的地表覆蓋度要低于土壤調(diào)節(jié)的植被指數(shù)(SAVI)。Jiang等在研究歸一化植被指數(shù)的空間屬性后,在多數(shù)情況下發(fā)現(xiàn)歸一化植被指數(shù)尺度方法會(huì)使植被覆蓋度升高。因此,不提倡用歸一化植被指數(shù)尺度模型來推算,而應(yīng)采用植被指數(shù)尺度差分方法(DVI)。Ghulam等采用Baret提出NDVI的x次冪來推算植被覆蓋度,公式如下:
減小大氣擾動(dòng)可以采用歸一化植被指數(shù)冪函數(shù)法,要最大限度降低大氣的影響,基于比值指數(shù)就可以實(shí)現(xiàn)。
假想一個(gè)圖像象元是混合的,即由完全植被覆蓋和裸露土地兩個(gè)部分組成,則該混合象元波段的反射率可通過植被與土壤線性函數(shù)推算,即R=fvRv,i+(1-fv)RS,i,則MPDI可以表示為:
假設(shè)在同種土壤類型下保持含水量的相對(duì)穩(wěn)定,干旱指數(shù)的影響主要發(fā)生在植被的fvPDIv和(1-fv)-1。相關(guān)研究指出, 隨著植被覆蓋度線性增加,在植被覆蓋度達(dá)到0.7之前, (1-fv)-1不會(huì)顯著變化,直到fv>0.7之后才會(huì)曲線增加。
3 結(jié)語
目前,利用遙感技術(shù)進(jìn)行干旱監(jiān)測(cè)逐步趨于成熟。為了提高國內(nèi)借助遙感技術(shù)進(jìn)行干旱實(shí)時(shí)監(jiān)測(cè)的研究水平,降低干旱對(duì)我國農(nóng)業(yè)發(fā)展的危害程度,提出如下建議: 一是由于下載的衛(wèi)星遙感數(shù)據(jù)往往是二級(jí)數(shù)據(jù)(系統(tǒng)幾何校正產(chǎn)品),因此需要進(jìn)行輻射定標(biāo)、幾何精校正、大氣校正、裁剪、云掩膜等多個(gè)前期處理工作[9]。二是由于不同研究區(qū)域的土壤下墊面土質(zhì)種類不同,因此選取的模型方法會(huì)嚴(yán)重影響最后反演旱情的精度和可信度。三是每個(gè)植被指數(shù)均涉及農(nóng)作物生長(zhǎng)狀況,對(duì)其進(jìn)行旱情監(jiān)測(cè),可拓寬遙感數(shù)據(jù)源在干旱監(jiān)測(cè)中的研究。四是微波遙感方法對(duì)反演土壤水量的精度較好,可連續(xù)使用且不受云層干擾,但相關(guān)系數(shù)較難獲取,應(yīng)與其他遙感數(shù)據(jù)結(jié)合使用。五是在特征空間法的應(yīng)用過程中,要針對(duì)不同時(shí)相、不同分辨率的遙感數(shù)據(jù),分析不同植被覆蓋類型和水分盈虧狀況下干、濕邊的變化規(guī)律及其在干旱監(jiān)測(cè)中的應(yīng)用,以提高模型精度。
參考文獻(xiàn)
[1] 李付琴,隋洪智.農(nóng)作物遙感估產(chǎn)研究進(jìn)展[J].地球科學(xué)進(jìn)展,1992,7(3):30-36.
[2] 陳述彭,童慶禧,郭華東.遙感信息機(jī)理研究[M].科學(xué)出版社,1998.
[3] 郭廣猛,趙冰茹.使用MODIS數(shù)據(jù)監(jiān)測(cè)土壤濕度[J].土壤,2004,36(2):219-221.
[4] 黃楊,耿淮濱,王松,等.潮濕裸露土壤的微波散射特性[J].遙感信息,1986(1):24-27.
[5] 劉良云,張兵,鄭蘭芬.利用溫度和植被指數(shù)進(jìn)行地物分類和土壤水分反演[J].紅外與毫米波學(xué)報(bào),2002,21(4):269-273.
[6] 張紅衛(wèi),陳懷亮,申雙和,等. NDVI-ST特征空間及干濕邊變化特征[J].氣象科技,2010,38(1):45-52
[7] 李星敏,劉安麟,張樹譽(yù),等.熱慣量法在干旱遙感監(jiān)測(cè)中的應(yīng)用研究[J].干旱地區(qū)農(nóng)業(yè)研究,2005,23(1):54-59.
[8] 張仁華.改進(jìn)的熱慣量模式及遙感土壤水分[J].地理研究,1990,9(2):101-112.
[9] 蔣耿明.MODIS數(shù)據(jù)基礎(chǔ)處理方法研究和軟件實(shí)現(xiàn)[D].北京:中國科學(xué)院遙感應(yīng)用研究所,2003.
2.2.3 Ts-NDVI特征空間法 Carlson,Price和Moran等得出在二維特征空間中Ts-NDVI分布是三角形或梯形,不同學(xué)者據(jù)此提出了很多監(jiān)測(cè)指數(shù)。Inge Sanholt等提出了條件溫度植被指數(shù)(TVDI),計(jì)算公式如下:
濕邊方程是指歸一化植被指數(shù)等于某一特定值時(shí)的最小值溫度,干邊方程是指當(dāng)歸一化植被指數(shù)值等于某一特定值時(shí)地表的最大溫度。a,b,a',b'是4個(gè)待定系數(shù)。通常根據(jù)研究區(qū)域的地表溫度和歸一化植被指數(shù)二維散點(diǎn)圖來近似提取干、濕邊方程。
該指數(shù)對(duì)研究區(qū)域有如下要求:區(qū)域內(nèi)有植被處于凋萎狀態(tài)時(shí)的象元,此時(shí)土壤水分含量是比較低的,同時(shí)也存在促使植被良好生長(zhǎng)發(fā)育的象元,此時(shí)土壤水量較大,通常用于表示年內(nèi)某地區(qū)的干旱等級(jí)[5-6]。
2.2.4 土壤熱慣量法 當(dāng)土壤含水量較多時(shí),會(huì)產(chǎn)生更大的熱慣量,兩者具有一定的關(guān)聯(lián)性。熱慣量值的計(jì)算式如下:
通常計(jì)算時(shí)無法直接利用遙感數(shù)據(jù)計(jì)算以上參數(shù)。為此,Price等提出了表觀熱慣量(ATI)定義,可以利用遙感影像數(shù)據(jù)的熱紅外輻射溫度差、近紅外反射率和可見光數(shù)據(jù),計(jì)算熱慣量值。計(jì)算公式如下:
通過遙感可獲取地表的有關(guān)信息,例如MODIS遙感影像采用的是近紅外波段(第1通道、第2通道、第3通道、第4通道)和第31熱紅外通道值,來反演推算大氣透過率、地表溫度、反射率等參數(shù),進(jìn)而計(jì)算土壤表觀熱慣量(ATI),再通過表觀熱慣量與土壤水分量之間的擬合關(guān)系,建立干旱模型,反演土壤墑情變化程度和干旱程度。該原理僅適用于植被覆蓋度不是很高的情況,對(duì)于植被覆蓋度較高的區(qū)域,其監(jiān)測(cè)結(jié)果會(huì)明顯降低[7-8]。
2.2.5 垂直干旱指數(shù)法 植物葉表面對(duì)紅光和藍(lán)紫光的吸收能力強(qiáng),對(duì)近紅外波段具有較高的反射表現(xiàn),因此,從紅波段(Red)到近紅外波段(Nir),裸地反射率值較高但其增幅卻不明顯。植被對(duì)紅光波段的吸收很快就能達(dá)到飽和狀態(tài),所以植被的增加可通過近紅外反射的增加而得到反映。秦明和詹志明等在分析前人研究以及紅光和近紅外特征空間后,提出了基于土壤線斜率的垂直干旱指數(shù),計(jì)算公式如下:
該指數(shù)假想了一條通過原點(diǎn)的直線L與土壤線垂直,其二維特征空間中任意一個(gè)點(diǎn)到直線L的距離都可表示為水分的分布狀況。水分含量越高代表距離越近,水分含量越少則代表越遠(yuǎn)。該方法的優(yōu)點(diǎn)是用光譜曲線特征直接代替了地表溫度和反射率的反演,簡(jiǎn)便實(shí)效。
秦明等研究指出,用垂直干旱指數(shù)監(jiān)測(cè)裸地地表水分效果良好,但遇到植被覆蓋度較高的田間,垂直干旱指數(shù)監(jiān)測(cè)會(huì)受到很大影響。不同土壤類型下地形特征會(huì)分布成非平面,為了改進(jìn)模型,阿布都瓦斯提等提出了修正的垂直千旱指數(shù)法 (MPDI),定義為:
植被地面和部分植被覆蓋程度是植被冠層的兩個(gè)重要參量因子,前兩者與葉片面積指數(shù)相關(guān)性很大。植被光譜指數(shù)可以計(jì)算出地表的植被覆蓋度。例如,歸一化植被指數(shù)(NDVI)和植被覆蓋度關(guān)聯(lián)性非常高,國內(nèi)部分學(xué)者提出對(duì)植被覆蓋度與歸一化植被指數(shù)進(jìn)行冪次方變換和均方根可提升植被覆蓋度的精度,還有研究指出使用歸一化植被指數(shù)(NDVI)計(jì)算出的地表覆蓋度要低于土壤調(diào)節(jié)的植被指數(shù)(SAVI)。Jiang等在研究歸一化植被指數(shù)的空間屬性后,在多數(shù)情況下發(fā)現(xiàn)歸一化植被指數(shù)尺度方法會(huì)使植被覆蓋度升高。因此,不提倡用歸一化植被指數(shù)尺度模型來推算,而應(yīng)采用植被指數(shù)尺度差分方法(DVI)。Ghulam等采用Baret提出NDVI的x次冪來推算植被覆蓋度,公式如下:
減小大氣擾動(dòng)可以采用歸一化植被指數(shù)冪函數(shù)法,要最大限度降低大氣的影響,基于比值指數(shù)就可以實(shí)現(xiàn)。
假想一個(gè)圖像象元是混合的,即由完全植被覆蓋和裸露土地兩個(gè)部分組成,則該混合象元波段的反射率可通過植被與土壤線性函數(shù)推算,即R=fvRv,i+(1-fv)RS,i,則MPDI可以表示為:
假設(shè)在同種土壤類型下保持含水量的相對(duì)穩(wěn)定,干旱指數(shù)的影響主要發(fā)生在植被的fvPDIv和(1-fv)-1。相關(guān)研究指出, 隨著植被覆蓋度線性增加,在植被覆蓋度達(dá)到0.7之前, (1-fv)-1不會(huì)顯著變化,直到fv>0.7之后才會(huì)曲線增加。
3 結(jié)語
目前,利用遙感技術(shù)進(jìn)行干旱監(jiān)測(cè)逐步趨于成熟。為了提高國內(nèi)借助遙感技術(shù)進(jìn)行干旱實(shí)時(shí)監(jiān)測(cè)的研究水平,降低干旱對(duì)我國農(nóng)業(yè)發(fā)展的危害程度,提出如下建議: 一是由于下載的衛(wèi)星遙感數(shù)據(jù)往往是二級(jí)數(shù)據(jù)(系統(tǒng)幾何校正產(chǎn)品),因此需要進(jìn)行輻射定標(biāo)、幾何精校正、大氣校正、裁剪、云掩膜等多個(gè)前期處理工作[9]。二是由于不同研究區(qū)域的土壤下墊面土質(zhì)種類不同,因此選取的模型方法會(huì)嚴(yán)重影響最后反演旱情的精度和可信度。三是每個(gè)植被指數(shù)均涉及農(nóng)作物生長(zhǎng)狀況,對(duì)其進(jìn)行旱情監(jiān)測(cè),可拓寬遙感數(shù)據(jù)源在干旱監(jiān)測(cè)中的研究。四是微波遙感方法對(duì)反演土壤水量的精度較好,可連續(xù)使用且不受云層干擾,但相關(guān)系數(shù)較難獲取,應(yīng)與其他遙感數(shù)據(jù)結(jié)合使用。五是在特征空間法的應(yīng)用過程中,要針對(duì)不同時(shí)相、不同分辨率的遙感數(shù)據(jù),分析不同植被覆蓋類型和水分盈虧狀況下干、濕邊的變化規(guī)律及其在干旱監(jiān)測(cè)中的應(yīng)用,以提高模型精度。
參考文獻(xiàn)
[1] 李付琴,隋洪智.農(nóng)作物遙感估產(chǎn)研究進(jìn)展[J].地球科學(xué)進(jìn)展,1992,7(3):30-36.
[2] 陳述彭,童慶禧,郭華東.遙感信息機(jī)理研究[M].科學(xué)出版社,1998.
[3] 郭廣猛,趙冰茹.使用MODIS數(shù)據(jù)監(jiān)測(cè)土壤濕度[J].土壤,2004,36(2):219-221.
[4] 黃楊,耿淮濱,王松,等.潮濕裸露土壤的微波散射特性[J].遙感信息,1986(1):24-27.
[5] 劉良云,張兵,鄭蘭芬.利用溫度和植被指數(shù)進(jìn)行地物分類和土壤水分反演[J].紅外與毫米波學(xué)報(bào),2002,21(4):269-273.
[6] 張紅衛(wèi),陳懷亮,申雙和,等. NDVI-ST特征空間及干濕邊變化特征[J].氣象科技,2010,38(1):45-52
[7] 李星敏,劉安麟,張樹譽(yù),等.熱慣量法在干旱遙感監(jiān)測(cè)中的應(yīng)用研究[J].干旱地區(qū)農(nóng)業(yè)研究,2005,23(1):54-59.
[8] 張仁華.改進(jìn)的熱慣量模式及遙感土壤水分[J].地理研究,1990,9(2):101-112.
[9] 蔣耿明.MODIS數(shù)據(jù)基礎(chǔ)處理方法研究和軟件實(shí)現(xiàn)[D].北京:中國科學(xué)院遙感應(yīng)用研究所,2003.
2.2.3 Ts-NDVI特征空間法 Carlson,Price和Moran等得出在二維特征空間中Ts-NDVI分布是三角形或梯形,不同學(xué)者據(jù)此提出了很多監(jiān)測(cè)指數(shù)。Inge Sanholt等提出了條件溫度植被指數(shù)(TVDI),計(jì)算公式如下:
濕邊方程是指歸一化植被指數(shù)等于某一特定值時(shí)的最小值溫度,干邊方程是指當(dāng)歸一化植被指數(shù)值等于某一特定值時(shí)地表的最大溫度。a,b,a',b'是4個(gè)待定系數(shù)。通常根據(jù)研究區(qū)域的地表溫度和歸一化植被指數(shù)二維散點(diǎn)圖來近似提取干、濕邊方程。
該指數(shù)對(duì)研究區(qū)域有如下要求:區(qū)域內(nèi)有植被處于凋萎狀態(tài)時(shí)的象元,此時(shí)土壤水分含量是比較低的,同時(shí)也存在促使植被良好生長(zhǎng)發(fā)育的象元,此時(shí)土壤水量較大,通常用于表示年內(nèi)某地區(qū)的干旱等級(jí)[5-6]。
2.2.4 土壤熱慣量法 當(dāng)土壤含水量較多時(shí),會(huì)產(chǎn)生更大的熱慣量,兩者具有一定的關(guān)聯(lián)性。熱慣量值的計(jì)算式如下:
通常計(jì)算時(shí)無法直接利用遙感數(shù)據(jù)計(jì)算以上參數(shù)。為此,Price等提出了表觀熱慣量(ATI)定義,可以利用遙感影像數(shù)據(jù)的熱紅外輻射溫度差、近紅外反射率和可見光數(shù)據(jù),計(jì)算熱慣量值。計(jì)算公式如下:
通過遙感可獲取地表的有關(guān)信息,例如MODIS遙感影像采用的是近紅外波段(第1通道、第2通道、第3通道、第4通道)和第31熱紅外通道值,來反演推算大氣透過率、地表溫度、反射率等參數(shù),進(jìn)而計(jì)算土壤表觀熱慣量(ATI),再通過表觀熱慣量與土壤水分量之間的擬合關(guān)系,建立干旱模型,反演土壤墑情變化程度和干旱程度。該原理僅適用于植被覆蓋度不是很高的情況,對(duì)于植被覆蓋度較高的區(qū)域,其監(jiān)測(cè)結(jié)果會(huì)明顯降低[7-8]。
2.2.5 垂直干旱指數(shù)法 植物葉表面對(duì)紅光和藍(lán)紫光的吸收能力強(qiáng),對(duì)近紅外波段具有較高的反射表現(xiàn),因此,從紅波段(Red)到近紅外波段(Nir),裸地反射率值較高但其增幅卻不明顯。植被對(duì)紅光波段的吸收很快就能達(dá)到飽和狀態(tài),所以植被的增加可通過近紅外反射的增加而得到反映。秦明和詹志明等在分析前人研究以及紅光和近紅外特征空間后,提出了基于土壤線斜率的垂直干旱指數(shù),計(jì)算公式如下:
該指數(shù)假想了一條通過原點(diǎn)的直線L與土壤線垂直,其二維特征空間中任意一個(gè)點(diǎn)到直線L的距離都可表示為水分的分布狀況。水分含量越高代表距離越近,水分含量越少則代表越遠(yuǎn)。該方法的優(yōu)點(diǎn)是用光譜曲線特征直接代替了地表溫度和反射率的反演,簡(jiǎn)便實(shí)效。
秦明等研究指出,用垂直干旱指數(shù)監(jiān)測(cè)裸地地表水分效果良好,但遇到植被覆蓋度較高的田間,垂直干旱指數(shù)監(jiān)測(cè)會(huì)受到很大影響。不同土壤類型下地形特征會(huì)分布成非平面,為了改進(jìn)模型,阿布都瓦斯提等提出了修正的垂直千旱指數(shù)法 (MPDI),定義為:
植被地面和部分植被覆蓋程度是植被冠層的兩個(gè)重要參量因子,前兩者與葉片面積指數(shù)相關(guān)性很大。植被光譜指數(shù)可以計(jì)算出地表的植被覆蓋度。例如,歸一化植被指數(shù)(NDVI)和植被覆蓋度關(guān)聯(lián)性非常高,國內(nèi)部分學(xué)者提出對(duì)植被覆蓋度與歸一化植被指數(shù)進(jìn)行冪次方變換和均方根可提升植被覆蓋度的精度,還有研究指出使用歸一化植被指數(shù)(NDVI)計(jì)算出的地表覆蓋度要低于土壤調(diào)節(jié)的植被指數(shù)(SAVI)。Jiang等在研究歸一化植被指數(shù)的空間屬性后,在多數(shù)情況下發(fā)現(xiàn)歸一化植被指數(shù)尺度方法會(huì)使植被覆蓋度升高。因此,不提倡用歸一化植被指數(shù)尺度模型來推算,而應(yīng)采用植被指數(shù)尺度差分方法(DVI)。Ghulam等采用Baret提出NDVI的x次冪來推算植被覆蓋度,公式如下:
減小大氣擾動(dòng)可以采用歸一化植被指數(shù)冪函數(shù)法,要最大限度降低大氣的影響,基于比值指數(shù)就可以實(shí)現(xiàn)。
假想一個(gè)圖像象元是混合的,即由完全植被覆蓋和裸露土地兩個(gè)部分組成,則該混合象元波段的反射率可通過植被與土壤線性函數(shù)推算,即R=fvRv,i+(1-fv)RS,i,則MPDI可以表示為:
假設(shè)在同種土壤類型下保持含水量的相對(duì)穩(wěn)定,干旱指數(shù)的影響主要發(fā)生在植被的fvPDIv和(1-fv)-1。相關(guān)研究指出, 隨著植被覆蓋度線性增加,在植被覆蓋度達(dá)到0.7之前, (1-fv)-1不會(huì)顯著變化,直到fv>0.7之后才會(huì)曲線增加。
3 結(jié)語
目前,利用遙感技術(shù)進(jìn)行干旱監(jiān)測(cè)逐步趨于成熟。為了提高國內(nèi)借助遙感技術(shù)進(jìn)行干旱實(shí)時(shí)監(jiān)測(cè)的研究水平,降低干旱對(duì)我國農(nóng)業(yè)發(fā)展的危害程度,提出如下建議: 一是由于下載的衛(wèi)星遙感數(shù)據(jù)往往是二級(jí)數(shù)據(jù)(系統(tǒng)幾何校正產(chǎn)品),因此需要進(jìn)行輻射定標(biāo)、幾何精校正、大氣校正、裁剪、云掩膜等多個(gè)前期處理工作[9]。二是由于不同研究區(qū)域的土壤下墊面土質(zhì)種類不同,因此選取的模型方法會(huì)嚴(yán)重影響最后反演旱情的精度和可信度。三是每個(gè)植被指數(shù)均涉及農(nóng)作物生長(zhǎng)狀況,對(duì)其進(jìn)行旱情監(jiān)測(cè),可拓寬遙感數(shù)據(jù)源在干旱監(jiān)測(cè)中的研究。四是微波遙感方法對(duì)反演土壤水量的精度較好,可連續(xù)使用且不受云層干擾,但相關(guān)系數(shù)較難獲取,應(yīng)與其他遙感數(shù)據(jù)結(jié)合使用。五是在特征空間法的應(yīng)用過程中,要針對(duì)不同時(shí)相、不同分辨率的遙感數(shù)據(jù),分析不同植被覆蓋類型和水分盈虧狀況下干、濕邊的變化規(guī)律及其在干旱監(jiān)測(cè)中的應(yīng)用,以提高模型精度。
參考文獻(xiàn)
[1] 李付琴,隋洪智.農(nóng)作物遙感估產(chǎn)研究進(jìn)展[J].地球科學(xué)進(jìn)展,1992,7(3):30-36.
[2] 陳述彭,童慶禧,郭華東.遙感信息機(jī)理研究[M].科學(xué)出版社,1998.
[3] 郭廣猛,趙冰茹.使用MODIS數(shù)據(jù)監(jiān)測(cè)土壤濕度[J].土壤,2004,36(2):219-221.
[4] 黃楊,耿淮濱,王松,等.潮濕裸露土壤的微波散射特性[J].遙感信息,1986(1):24-27.
[5] 劉良云,張兵,鄭蘭芬.利用溫度和植被指數(shù)進(jìn)行地物分類和土壤水分反演[J].紅外與毫米波學(xué)報(bào),2002,21(4):269-273.
[6] 張紅衛(wèi),陳懷亮,申雙和,等. NDVI-ST特征空間及干濕邊變化特征[J].氣象科技,2010,38(1):45-52
[7] 李星敏,劉安麟,張樹譽(yù),等.熱慣量法在干旱遙感監(jiān)測(cè)中的應(yīng)用研究[J].干旱地區(qū)農(nóng)業(yè)研究,2005,23(1):54-59.
[8] 張仁華.改進(jìn)的熱慣量模式及遙感土壤水分[J].地理研究,1990,9(2):101-112.
[9] 蔣耿明.MODIS數(shù)據(jù)基礎(chǔ)處理方法研究和軟件實(shí)現(xiàn)[D].北京:中國科學(xué)院遙感應(yīng)用研究所,2003.