基于METRIC模型的江漢平原蒸散量估算適應性分析

2021-04-02 03:51:30陳植華

安全與環境工程 2021年2期

吳凡，陳植華，胡成

(中國地質大學(武漢)環境學院，湖北武漢 430078)

江漢平原地處華中地區，西起宜昌，東至武漢，是國內重要的稻產水產基地之一。江漢平原屬亞熱帶季風氣候，年平均降雨量達1 000 mm以上。整個江漢平原河流縱橫交錯，湖泊星羅棋布，水是整個平原生態系統及資源環境的重要樞紐。水資源事關生態系統的循環及穩定，是區內重要的生態資源之一，正確地評價區域水資源量，可以更加明晰地揭示區域水文循環及能量循環過程，對合理配置水資源，保障區域生態用水提供更為科學的建議，也是實現水循環的可再生持續利用，保障生態、社會和諧發展的重要手段之一。蒸散發作為水文循環中的重要消耗途徑之一，濕潤及半濕潤地區可占降水總量的約60%，干旱地區可占降水總量的近90%，而完整的蒸散發過程更能夠將近50%的太陽輻射進行對流交換，因此蒸散發也是能量平衡以及水量平衡的重要影響因子，準確、快捷地評價區域蒸散量，對正確評價區域水資源量有著至關重要的作用。蒸散發表征了植物表面以及地表的水分蒸散以及植物體內的水分蒸騰，其物理過程十分復雜，與土壤熱交換、大氣湍流交換、植被冠層水汽交換等都有著緊密的聯系，也是整個水文循環中最難直接測量的部分。作為水資源主要消耗途徑之一的蒸散發也日益受到全球學者的關注，對其研究也隨著關注度日益加深。由于蒸散發物理過程復雜，涉及環境變量多，如何精確地估算蒸散量已經成為水文學、氣象學等多個學科的重要研究課題之一。傳統手段通過經驗模型估算的蒸散量結果，其精度已難以滿足現代水資源量評價的需求，氣象站通過蒸發皿測得的蒸發量數據，也只能表征該區域內的蒸發潛力，與實際蒸散發仍然有著較大的差異，而通過儀器測定如蒸滲儀法、渦度相關儀法等，其空間尺度仍然局限在點或小區域上，難以實現大面積的蒸散發測定。針對以上問題，本文采用遙感蒸散發模型反演江漢平原實際蒸散量，并探討METRIC模型對江漢平原蒸散量估算的適用性，對江漢平原的水資源量評價具有重要的意義。此次遙感反演的江漢平原屬半濕潤地區，且未進行地面蒸散測量設備如渦度儀、蒸滲儀等的建設，其基礎氣象數據暫無法滿足雙層及多層模型的使用條件，而由于冷熱像元的選取原則導致SEBAL模型在半濕潤地區的適用性較弱，此次引進了采用全新冷熱像元選取原則的METRIC模型，理論上能夠在半濕潤地區獲得更高的蒸散量遙感反演精度，同時能夠有效增強遙感反演在此類區域的可用性及實用性，為區域提供更為可靠的蒸散量數據，并提高區域水資源規劃精度，具有較高的研究意義。

1 研究數據與方法

1. 1 遙感數據源和氣象數據選擇

在本次研究區域蒸散量遙感反演中，主要使用的遙感數據為4期Landsat8遙感影像數據。Landsat8是美國陸地衛星計劃中的第8顆，于2013年2月11日發射升空，以填補自2003年Landsat7攜帶的SLC傳感器故障導致的遙感影像缺失。自1972年陸地衛星計劃首顆衛星發射以來，Landsat系列衛星已成為全球范圍內重要的環境監測及資源調查遙感衛星，從Landsat1/2/3只攜帶覆蓋了三波段的RBV傳感器及紅外波段的MSS傳感器，到如今Landasat8攜帶OLI光學傳感器及TIRS熱紅外傳感器，覆蓋波段高達11個波段，精度最高達15 m，Landsat系列衛星已在幾十年間取得了長足的進步。截止目前Landsat全系列衛星數據仍面向全球免費提供，至今仍然是進行區域蒸散量遙感反演的理想數據源之一。此次遙感反演選擇了研究區域2015年10月25日、2016年3月1日、2018年3月23日、2018年4月8日共4期的Landsat8遙感影像數據，此4期遙感影像數據為研究時間段內影像質量最佳的4期數據，并能夠一定程度反映研究區域蒸散量的季節性變化。本研究主要利用Landasat8攜帶的OLI光學傳感器的2～7波段進行植被指數(NDVI)、土地利用類型等的遙感反演以及TIRS熱紅外傳感器的10及11波段進行地表溫度等參數的遙感反演。氣象數據選取研究區域周邊孝感市氣象站和孝昌縣氣象站共兩個站點的氣溫、濕度、風速、降雨量等各類氣象數據輸入模型參與蒸散量遙感反演。

1.2 METRIC模型和SEBAL模型介紹

遙感蒸散發模型的發展距今已有30余年的歷史，其研究的開端為采用手持式紅外輻射儀測量植被冠層溫度并與蒸散發建立經驗公式，直至1983年，Seguin等首次利用衛星傳感器測得的地表溫度對蒸散發進行了估算，正式拉開了現代對遙感蒸散發模型的研究。

1988年，由Bastiaanssen提出了迄今為止使用最為廣泛的遙感蒸散發模型即SEBAL(Surface Energy Balance Algorithm for Land)模型。該模型基于地表能量平衡，巧妙避免了氣象數據在空間插值并訂正溫度時引起的誤差，具有十分重要的物理意義，目前在世界范圍內都得到了廣泛應用。美國的愛達荷州水利部門曾經對SEBAL模型的估算精度進行了大面積、高精度的地面驗證，結果表明其估算誤差可以控制在4%以內。在國內，由中國科學院蘭州大氣物理研究所在黑河地區的實驗研究表明：在估算單日蒸散量時，SEBAL模型的估算誤差保持在8%以內。上述研究說明，SEBAL模型在區域蒸散量估算上具有實際可用性。2005年，Allen等利用與SEBAL模型相同的原理，開發了METRIC模型，克服了SEBAL模型在山區反演精度不足的問題，并對模型中的干濕邊界采用了新的選取方式，同時在時間尺度上采取了全新的瞬時蒸散擴展方法，更好地考慮了風速、濕度等氣象因素帶來的影響。

METRIC模型和SEBAL模型的原理都是基于能量平衡方程，只需要遙感影像中的可見光、紅外光以及極少量氣象參數(氣溫、風速)等即可計算出地表各項輻射通量(包括凈輻射通量、土壤熱通量、顯熱通量)，同時根據這些輻射通量計算出表征蒸散發部分的輻射通量、潛熱通量，并利用比例不變法將瞬時蒸散量擴展到日蒸散量。

2 模型各項參數求解

遙感蒸散發模型的基礎原理為能量平衡，區域內的能量平衡方程如下：

λET

(1)

式中：

為凈輻射通量(W/m)，表征了太陽能夠有效輻射到地表的輻射能量;

為土壤熱通量(W/m)，表征了土壤與地面之間的熱量交換;

為顯熱通量(W/m)，表征了大氣與地面之間以湍流形式進行的感熱交換;

為水的汽化潛熱[W/(m·mm)]；

為蒸散量(mm)。

以上各項參數均可以通過遙感影像結合地表氣象數據反演得到，并具有較大的空間尺度。

在利用METRIC模型和SEBAL模型對區域內蒸散量進行遙感反演時，其主要需求的參數為太陽高度角、太陽天頂角、日地距離、歸一化植被指數、地表反照率、地表比輻射率，地表溫度等，以上各項地表參數需要結合遙感數據和氣象數據，并利用遙感軟件分別進行反演。此次遙感軟件平臺選用ENVI。METRIC模型的主要計算流程見圖1。

圖1 METRIC模型的計算流程圖Fig.1 Flow chart of METRIC model calculation

2. 1 模型參數獲取

2.1.1 歸一化植被指數(NDVI)

根據蒸散發的基礎原理，植被蒸騰作用是蒸散發的主要占比之一，區域內的植被覆蓋情況以及植被生長狀況等對蒸散發有著極大的影響。而在遙感反演中，體現植被覆蓋情況以及植被生長狀態的主要參數就是植被指數，又被稱作歸一化植被指數(NDVI)，其主要原理來自植被葉面對于光的吸收有著其獨特的特征，其中葉面對可見光中的紅光波段吸收能力很強，但與此同時，它對于不可見的近紅外波段有著極強的反射能力。根據以上兩個特征，對區域遙感影像的近紅外波段與紅光波段進行比值計算即可反映區域的植被覆蓋狀況，其計算公式如下：

NDVI=(NIR-R)/(NIR+R)

(2)

式中：NDVI為歸一化植被指數;NIR為近紅外波段的像元值(DN值);R為可見光紅波段DN值。

根據NDVI，可以計算得到區域的植被覆蓋度

，其計算公式如下：

(3)

式中:

為植被覆蓋度；NDVI為植被指數最小值；NDVI為植被指數最大值。2.1.2 地表反照率(

)地表反照率(

)也是整個蒸散量遙感反演中極為重要的地表參數之一，因為其表征的是區域內反射輻射通量總和與入射輻射通量總和的比值，該比值的大小決定了整個區域總入射能量的多少。

主要受到來自土地利用類型，即下墊面變化的影響，不同下墊面的

值差別較大，在Landsat8衛星影像的

計算中，采用如下計算公式：

(4)

式中：

為地表反照率;

為大氣外反照率;

為大氣程輻射(取值在0.025～0.04之間);

為大氣透過率。其中，大氣外反照率

的計算公式為

(5)

式中:

分別為Landsat8遙感影像2/4/5/6/7波段的DN值。大氣透過率

的計算公式為

=0.75+2×10×

(6)

式中:

為地面高程(m)，代指氣象觀測站的海拔高度。2.1.3 地表比輻射率(

)地表比輻射率(

)的物理意義表征在物體表面單位面積上所輻射出的輻射通量與在相同溫度下黑體輻射出的輻射通量的比值。

的值域受到很多因素的影響，通常有物質的界定常數、波長等不同條件的影響，通常

的取值范圍在0～1之間。正常情況下

難以進行量化計算，而多年來學者已經總結出地表輻射率與植被覆蓋率之間的經驗公式，并根據該經驗公式對

進行估算，經多位學者實踐證明其效果較好。但在對

進行估算時，需根據波段和土地利用類型分別進行計算，具體如下：水體：Landsat8 BAND10下，

取值為0.995；Landsat8 BAND11下，

取值為0.986。水體

取定值其主要原因為水體在熱紅外波段下，該

值遠遠高出其他土地利用類型，無限接近于黑體，因此根據不同紅外波段及不同衛星傳感器對其做定值處理。建筑用地：建筑用地地表比輻射率

的計算，公式為

(7)

式中：

為地表比輻射率;

為植被覆蓋度。

其他自然地物：其他自然地物地表比輻射率ε的計算公式為

(8)

2.1.4 地表溫度(

)

地表溫度的反演在遙感反演中起步較早，多年來已經開發出多種穩定算法用以計算地表溫度，均具有較好的效果,目前主流地表溫度反演算法有大氣校正法、分裂窗算法等。本次地表溫度的遙感反演根據傳感器類型以及區域氣候狀況采取了大氣校正法。

主要代表由Landsat8影像傳感器接收到的熱紅外輻射亮度，對比Landsat8中兩熱紅外波段的大氣透過情況可以發現，BAND10有更優異的性能表現，因此在日常使用中通常采取第10波段作為地表溫度反演的基礎波段。Landsat8影像傳感器接受到的輻射組成主要有三部分，分別是大氣向上輻射亮度、大氣向下輻射反射能量和地表熱輻射，其計算公式如下：

(

)+(1-

)

↓]·

↑

(9)

式中：

為熱紅外輻射亮度[W/(m·μm·sr]；

(

)為地表輻射亮度[W/(m·μm·sr]；

↑為大氣向上輻射亮度[W/(m·μm·sr]；

↓為大氣向下輻射亮度[W/(m·μm·sr]；

為地表溫度(K)。

對公式(9)進行數學轉換，可得到：

(10)

其中，

可以利用普朗克公式計算：

(11)

式中：

和

需根據傳感器及波段取定值，對于Landsat8的TRIS傳感器Band10波段，

為774.89W/(m·μm·sr)，

為1 321.08 K。

2. 2 能量平衡各分量反演

根據上述所得到的各地表特征參數，結合氣象參數(溫度、濕度、風速等)，即可對地表各項能量通量進行反演，并根據能量平衡計算潛熱通量。

2.2.1 凈輻射通量(

)

凈輻射通量(Net Radiant Flux，NRF)反映的是到達地表的太陽短波輻射和大氣長波輻射的凈收支，它是驅動大氣運動的主要能量。根據地表輻射平衡原理，地表單位面積上接收到的凈輻射通量的計算公式如下：

=(1-

)

εR

(12)

式中：

為凈輻射通量(W/m)；

為入射的太陽短波輻射(W/m);

為大氣下行長波輻射(W/m);

為大氣上行長波輻射(W/m)。

其中，入射的太陽短波輻射的計算公式如下：

×cos×

(13)

式中：

為太陽常數(W/m)，其取定值1 367 W/m;

為太陽天頂角(°)，可通過查詢遙感數據頭文件獲得；

為日地距離系數(無固定量綱),其計算公式為

(14)

式中：DOY為太陽歷日序。

大氣下行長波輻射的計算公式如下：

(15)

式中：

為玻爾茲曼常數(W/m·K)，其取定值5.67×10W/m·K；

為近地處氣溫(K)；

為大氣比輻射率，其計算公式為

=1.08(-ln

)

(16)

大氣上行長波輻射計算公式如下：

(17)

式中：ε為地表比輻射率。

2.2.2 土壤熱通量(

)通過儀器觀測可知，土壤在不同深度下其溫度也會隨深度發生改變，因此在不同深度下土壤發生了熱量交換，此過程中產生的通量叫做土壤熱通量(Soil Heat Flux，SOHF)。土壤熱通量與其他3個通量相比較小，而且它與凈輻射通量之間有一定的關系。Teixeira等提出考慮地表溫度、歸一化植被指數、地表反照率的土壤熱通量

計算公式為

(18)

式中：

為土壤熱通量(W/m)。上式僅針對區域內有植被覆蓋的區域，無植被覆蓋的區域中，水體：

=0.3

，建筑用地：

=0.41

-51。2.2.3 顯熱通量(

)

由地表流向大氣的熱量，一般被稱為顯熱通量，其主要計算公式及過程如下：

(19)

式中:

為顯熱通量(W/m)；

為空氣密度(kg/m)；

為空氣定壓比熱容[J/(kg·K)]，其取值為1 004 [J/(kg·K)];d

為不同高度處的溫差(K)，一般取0.01～2 m處的溫差；

為空氣動力學阻抗(s/m)。其中，空氣密度

的計算公式為

(20)

式中：

為近地表處氣溫(K);

為地面高程(m)。空氣動力學阻抗

的計算公式為

(21)

式中：

為卡曼常數，取值為0.41；

為摩擦速度(m/s)；

和

代表公式(19)中的兩個不同高度，

為2 m，

為0.01 m。其中，摩擦速度

的計算公式為

(22)

式中：

為

高度處的平均風速(m/s);

為動力粗糙度(m)，其暫無量化計算公式，可以由氣象站點的植被高度和歸一化植被指數按如下經驗公式擬合得到：

=exp(1.442 1×NDVI-3.916 7)

(23)

其中，相關系數

=0.898 8。

在此類單層遙感蒸散量估算中，使用了一個假設量，即在距離地面200 m高程的高空，區域內所有像元的風速均是相等的，根據此假設，即可方便地計算出整個區域內的摩擦速度，具體計算過程有以下幾個步驟：第一步，根據動力粗糙度計算公式計算出氣象站附近的動力粗糙度，再結合氣象數據，計算出氣象站點的摩擦速度；第二步，根據以上公式及氣象站的數據，計算出氣象站上空200 m處的風速，即為全研究區上空200 m處的風速，此風速值在整個反演過程中保持不變；第三步，根據以上數據及公式，即可以計算出整個研究區的摩擦速度。

基于能量平衡的單層模型中，假設溫度差d

與地表溫度

之間存在線性關系，并假設其關系式如下：d

(24)

式中：

和

均為線性系數，通過選擇干濕邊界，即冷熱像元可以求得。一般情況下，“冷點”即代表擁有大蒸散量的像元點，該類像元點一般擁有大量水分且生長旺盛的植被，且地表溫度較其他地區明顯偏低，該類像元一般假設其近地層氣溫與其地表溫度相等，對于METRIC模型，一般設定“冷點”像元的潛熱通量值為此像元點的參考蒸散量的1.05倍；“熱點”的選擇一般選擇裸土，因為此類像元點一般無植被覆蓋，無表面的水分來源，其蒸散量可以視作0，則

λET

≈0。本文中，將冷像元點選取在區域中部較其他地區溫度明顯較低的耕地內，而熱像元點則選取在西南地區遠離水源且溫度較周圍地區明顯偏高的建筑用地內。結合土壤熱通量計算公式(18)和顯熱通量計算公式(19)，根據能量平衡原理，針對冷熱像元的輻射特征，對線性系數

和

進行數學推導，可得到：

(25)

冷

(26)

式中：

熱和

分別為熱像元的凈輻射通量(W/m)和土壤熱通量(W/m);

熱為熱像元的空氣動力學阻抗(s/m);

和

分別為熱像元和冷像元的地表溫度(K)。在得到線性系數

和

之后，即可根據溫度差計算公式(24)計算得到全區域的溫度差。

由于低層大氣存在不穩定性，因此根據上述公式計算所得的顯熱通量也存在一定的誤差，METRIC模型等單層能量平衡模型引入了莫寧-奧布霍夫(Monin-Obukhov)的循環遞歸計算對顯熱通量的結果進行了校正，通過5次以上的循環計算至莫寧-奧布霍夫長度值大于0，即可得到正確的顯熱通量值。

2.2.4 潛熱通量

根據能量平衡公式(1)，即可計算潛熱通量

λET

。

2.2.5 日蒸散量

通過閱讀眾多學者前期的研究成果可以發現，估算所得蒸散量與參考作物蒸散量的比值在24 h內，其值的變化較小，因此在計算日蒸散量時，假設衛星影像過境瞬間的蒸發率在當天是保持不變的，此假設引起的誤差在天及更長的時間尺度上對最終結果的影響較小，其計算公式如下：

(27)

式中：

為瞬時蒸發率；

為日蒸發率。

則過境當天的日蒸散量計算公式如下：

(28)

式中：

為當日日蒸散量(mm/d)，86 400為時間系數；

24為日平均凈輻射通量(W/m);

為日平均土壤熱通量(W/m)。

3 研究結果與討論

本次以位于江漢平原的孝感市孝昌縣肖港鎮作為遙感反演的主要區域，其地理位置位于北緯30°00′～31°10′、東經113°45′～114°00′之間，包括了1∶5萬水文地質調查肖港鎮幅圖幅所在的整個區域，總面積為430 km。整個研究區域從氣候帶分帶劃分，屬亞熱帶季風性氣候，一年中四季雨水量均比較充足，日照時間長，氣候較溫暖。區域內風向主要以東南風為主，一年平均氣溫一般在15～16℃之間，冬季平均氣溫一般在1～4℃之間，夏季平均氣溫一般在28～29℃之間，其年內日平均氣溫一般不低于10℃，整個孝感市年平均降雨量在1 200 mm左右，一年之中約有70%的降水集中在屬于夏季的7月至9月，其多年年平均蒸散量約為1 435 mm。

3. 1 遙感反演結果

以上述模型公式及理論為基礎，利用ESRI旗下ENVI軟件平臺和ArcGIS空間分析平臺分別根據METRIC模型和SEBAL模型的原理，對研究區域進行日蒸散量遙感反演，并根據世界糧食組織(FAO)推薦的P-M模型對METRIC模型和SEBAL模型最終得到的研究區域蒸散量估算結果進行精度檢驗與對比，其中METRIC模型額外加入DEM數據對區域地形坡度加以修正，以降低高程差引起的冷卻效應。

本文利用METRIC模型反演得到的研究區域4日日蒸散量頻率分布圖和4日蒸散量空間分布圖，見圖2和圖3。

圖2 利用METRIC模型反演得到的研究區域4日日蒸散量頻率分布圖Fig.2 Four-day daily evapotranspiration frequency distribution diagram of the study area retrieved from METRIC model

圖3 利用METRIC模型反演得到的研究區域4日蒸散量空間分布圖Fig.3 Four-day daily evapotranspiration spatial distribution diagram of the study area retrieved from METRIC model

由圖2和圖3可見，研究區域內日蒸散量主要集中在0～7 mm/d，在空間分布上，城鎮、裸地等遠離水體及植被覆蓋度較低的地區，其日蒸散量較小，而河流及湖泊等因其水源直接暴露在太陽輻射下，無土地及植被等固定的地區，其日蒸散量較大；因研究區域內作物下墊面存在水稻田等含水量較大的下墊面，其蒸散規律與水體類似，導致耕地中也存在部分大蒸散量像元，其蒸散量達6 mm/d，總體上研究區域內日蒸散量符合以潛熱通量交換為主的客觀規律；在時間分布上，研究區域2018年4月8日的蒸散量要顯著大于其他3日的蒸散量,觀察當天平均氣溫達19.7℃，要顯著高于其他3日，且其時間節點剛好位于2018年清明節之后，氣溫顯著提升，日照條件良好，風速較大，其蒸散量為4日最大，符合客觀規律，而秋季研究區域2015年10月25日的蒸散量最小，符合其季節條件下的客觀規律。

SEBAL模型的區域蒸散量估算方法，與MEITRIC模型類似，均為利用能量平衡余項對潛熱通量進行估算，但未引入DEM數據進行精校正，且冷熱像元的選取原則以地表溫度大小為主，以干燥的裸地作為熱像元，以濕潤的低溫度水體作為冷像元。利用SEBAL模型反演得到的研究區域4日日蒸散量頻率分布圖，見圖4。

圖4 利用SEBAL模型反演得到的研究區域4日日蒸散量頻率分布圖Fig.4 Four-day daily evapotranspiration frequency distribution diagram of the study area retrieved from SEBAL model

利用FAO提供的改進P-M模型以及氣象站氣象數據估算研究區域4日遙感影像當天氣象站點參考蒸散量，并作為此次遙感反演的標準蒸散量，與遙感反演所得的氣象站點實際蒸散量進行了對比，同時進行了估算誤差分析。P-M模型4日日蒸散量的估算結果，見表1。

3. 2 討論

本文分別統計了各遙感蒸散發模型(即METRIC模型和SEBAL模型)反演得到的研究區域4日遙感影像當天氣象站點蒸散量，并與P-M模型日蒸散量估算結果進行了對比同時進行了估算誤差分析，其結果見表1。

表1 3種模型研究區域日蒸散量的估算結果及相對誤差

由表1可知，在日蒸散量較小的時間段，METRIC模型日蒸散量的估算誤差要低于SEBAL模型，其中在2016年3月1日，METRIC模型日蒸散量的估算誤差較SEBAL模型低約4%，在2018年3月23日，METRIC模型日蒸散量的估算誤差較SEBAL模型低約13%，在2015年10月25日，METRIC模型日蒸散量的估算誤差較SEBAL模型低約10%，而在2018年4月8日，其日蒸散量較大，METRIC模型日蒸散量的估算結果出現了較大誤差，其估算誤差較SEBAL模型高約33%，而SEBAL模型在蒸散量較大的時間段與P-M模型日蒸散量的估算結果保持了高度的一致，其相對誤差僅為0.28%，受限于Landsat8圖像分辨率原因，其在點位上的估算精度與實際測量結果仍存在一定的差距，導致其與P-M模型日蒸散量的估算誤差在大部分時間段維持在10%以上，此處點位估算僅用于誤差趨勢分析，其實際應用模式應為區域性估算。

通過以上數據可以發現，在蒸散量較小的時間段內，METRIC模型的估算精度要高于SEBAL模型，METRIC模型在蒸散量較小的3日內平均估算誤差較SEBAL模型低約9%，說明METRIC模型在蒸散量較小的季節，在反演區域所屬氣候環境下，相較于SEBAL模型具有更好的區域適用性；在蒸散量較大的季節，SEBAL模型表現出了更好的適用性，以P-M模型為基準，其估算誤差僅為0.28%，METRIC模型的估算精度較其有較大的差距。從METRIC模型和SEBAL模型的原理分析兩個模型在不同蒸散條件下估算精度表現不一致的原因，可以發現其主要原因集中在冷熱像元的選取原則上，即模型干濕邊界的選取原則上，兩模型均通過人工圈定方式來選擇區域內的干濕邊界，但選取原則具有一定的差別，SEBAL模型的濕邊界一般以溫度較低的水體為主，而METRIC模型則以溫度較低且水分較多的濕潤地帶，如林地、耕地等作為其濕邊界，不同土地利用類型的蒸散發特征表明，水體的蒸散發在區域內總蒸散量的占比較高；而METRIC模型的干濕邊界選取原則，雖避免了SEBAL模型在植被覆蓋度較高的區域估算時受水分脅迫像元的影響而帶來的誤差，但因未直接選擇水體作為濕邊界，降低了對水體蒸散部分權重的考慮，導致其在蒸散量較大的季節出現了較大的估算誤差，反之SEBAL模型的干濕邊界選取原則導致其在此條件下表現出了較高的估算精度；兩模型的估算誤差表現符合兩個模型不同的干濕邊界選取原則。此外，METRIC模型中通過加入DEM數據，考慮了高程、地形坡度的影響以及“高程差”帶來的冷卻效應，進一步減小了因下墊面不平整而帶來的估算誤差。

4 結論與建議

METRIC模型一般用于干旱及半干旱、地勢起伏較大地區的蒸散量遙感估算，此次將METRIC模型首次運用于江漢平原地區蒸散量估算，并采用全球內運用范圍較廣、效果較好的SEBAL模型進行蒸散量估算并進行估算結果對比，同時利用P-M模型進行估算誤差驗證，總體上取得了較好的實驗效果。結果表明：METRIC模型在蒸散量較小的時間段內，在研究區域范圍內具有較好的適用性，以P-M模型估算結果為基準，其平均估算精度較SEBAL模型提高約9%；而在蒸散量較大的時間段內，其估算結果出現了一定的偏差，與SEBAL模型的估算精度有較大的差距。因此，在研究區域所屬江漢平原地帶遙感蒸散量估算模型選擇中，在蒸散量較小的時間及季節，推薦以METRIC模型作為遙感反演模型，而在蒸散量較大的時間及季節，如夏季等，推薦以SEBAL模型作為遙感反演模型，兩種模型交替使用能夠有效提高本區域蒸散量的估算精度。

在本研究中，仍存在諸多不足之處，例如：雖同時采用了METRIC模型和SEBAL模型雙模型作為對比，但仍未解決干濕邊界選擇方式為人為圈定這個兩個模型都存在的重要限制問題，此次遙感反演因條件所限，未進行嚴格的地面蒸散量及溫度實測對比驗證，雖采用聯合國糧農組織提供的改進P-M模型估算參考蒸散量以驗模型估算精度，但仍與地面實際驗證數據存在差距，今后的研究中需提前對驗證方式提出方案，并對驗證所需儀器進行布置，以進一步提高驗證數據的精度及驗證方法的科學性；受限于國內衛星影像數據獲取途徑，此次僅獲得4期可用遙感影像數據，后期需通過更多途徑獲取更多可用遙感影像數據。

相較于傳統蒸散量測定方式，區域蒸散量遙感反演的便捷性及大尺度性賦予了其極高的研究價值，隨著近年來各類高分辨率傳感器反射、各種基于大數據統計以及AI神經網絡算法的涌現，賦予了蒸散量遙感估算長足的進步空間。因此，在未來的研究中還需加強對此類單層模型的優化，主要可集中在以下幾點：對干濕邊界選取流程化、數據化確定；進一步加強復雜下墊面光譜分析；對模型中仍使用經驗公式計算的部分參數提出量化計算方法；鋪設大孔徑閃爍通量儀，結合小流域范圍水文閉合，對遙感反演結果進行高精度驗證，以此為遙感蒸散發模型估算的精度提高提供更可信的數據。相信在解決以上難點之后，遙感蒸散發模型的估算精度及適用范圍將會得到進一步的提升。