白 娟，楊勝天，董國(guó)濤，2，郝芳華，王鳴程，4

（1.北京師范大學(xué) 地理學(xué)與遙感科學(xué)學(xué)院 遙感科學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室，北京100875；2.黃河水利委員會(huì) 黃河水利科學(xué)研究院，鄭州450003；3.北京師范大學(xué) 環(huán)境學(xué)院 水環(huán)境模擬國(guó)家重點(diǎn)實(shí)驗(yàn)室，北京100875；4.中國(guó)科學(xué)院 昆明植物研究所 山地生態(tài)系統(tǒng)研究中心，昆明650204）

蒸散發(fā)是影響區(qū)域水熱平衡的主要因素，也是水文循環(huán)的重要環(huán)節(jié)之一。區(qū)域蒸散過程涉及土壤、植被和大氣等多個(gè)要素，時(shí)空尺度陸面蒸散量的精確評(píng)價(jià)對(duì)區(qū)域水循環(huán)研究和植被生產(chǎn)力評(píng)價(jià)具有重要意義。蒸散發(fā)的估算方法目前有很多，如點(diǎn)尺度的波文比能量平衡法、空氣動(dòng)力學(xué)方法和渦度相關(guān)法，區(qū)域尺度的水量平衡理論、SPAC理論和參考作物騰發(fā)量等［1］。聯(lián)合國(guó)糧農(nóng)組織（FAO）于1998年推出的修正Pennman－Monteith（P－M）模型是目前最常用的蒸散發(fā)估算方法，常用于檢驗(yàn)其他蒸散發(fā)模型的精度［2－3］。P－M 模型綜合了能量平衡（熱量平衡）方程和空氣動(dòng)力學(xué)方法，具有較明確的物理意義，用于估算較長(zhǎng)時(shí)期（月或年）的蒸散發(fā)時(shí)較為精確，被廣泛應(yīng)用于區(qū)域參考作物蒸散量的估算［4－8］。目前大部分研究主要是基于氣象觀測(cè)資料，通過估算氣象站點(diǎn)的蒸散量，然后對(duì)站點(diǎn)的蒸散量進(jìn)行空間插值推廣到區(qū)域上，這種方法在站點(diǎn)上具有很高的精度，但在區(qū)域尺度受下墊面幾何結(jié)構(gòu)與物理性質(zhì)的空間異質(zhì)性影響，較難取得準(zhǔn)確的結(jié)果。

隨著遙感技術(shù)的發(fā)展和應(yīng)用，衛(wèi)星遙感技術(shù)被廣泛應(yīng)用于區(qū)域蒸散發(fā)的研究中。雖然遙感不能直接測(cè)量蒸散發(fā)，但是利用遙感技術(shù)可以定量反演用于計(jì)算陸面蒸散的地表參數(shù)和地表通量。不少學(xué)者采用P－M模型結(jié)合遙感的方法對(duì)區(qū)域蒸散發(fā)量進(jìn)行估算［9－10］，但這些研究多利用遙感數(shù)據(jù)進(jìn)行通量反演，氣象數(shù)據(jù)仍采用站點(diǎn)數(shù)據(jù)的插值結(jié)果，在資料缺乏地區(qū)或者氣象觀測(cè)站點(diǎn)分布不均勻的地區(qū)的應(yīng)用仍存在局限。

本研究在FAO Penman－Monteith模型的基礎(chǔ)上，采用多源遙感數(shù)據(jù)，對(duì)三江平原2011年生長(zhǎng)季（5—9月）的日實(shí)際蒸散量進(jìn)行估算，并結(jié)合地表連續(xù)觀測(cè)數(shù)據(jù)對(duì)模型精度進(jìn)行評(píng)價(jià)?；诖?，進(jìn)一步分析研究區(qū)實(shí)際蒸散發(fā)的時(shí)空分布特征和氣象因素對(duì)實(shí)際蒸散發(fā)的影響，驗(yàn)證了遙感驅(qū)動(dòng)的Penman－Monteith模型在三江平原的適用性，同時(shí)為三江平原的生態(tài)水文過程模擬提供了重要的基礎(chǔ)數(shù)據(jù)，對(duì)研究三江平原的生態(tài)水文過程具有重要意義。

1 研究區(qū)與實(shí)驗(yàn)區(qū)概況

三江平原位于黑龍江省東部，是由黑龍江、烏蘇里江及松花江沖積而成的低平原。該區(qū)北起黑龍江，南抵興凱湖，西鄰小興安嶺，東至烏蘇里江，地理坐標(biāo)為43°49′55″—48°27′40″N，129°11′20″—135°05′26″E，總面積10.89萬km2。氣候?qū)儆跍貛駶?rùn)、半濕潤(rùn)大陸性季風(fēng)氣候區(qū)，全年日照時(shí)數(shù)2 400～2 500h，1月平均氣溫低于－18℃，7月平均氣溫21～22℃，無霜期120～140d，凍結(jié)期140～190d，季節(jié)性凍土深度為1.4～2.5m。年降水量500～650mm，降水時(shí)空分布不均勻，75%～85%的降水集中在6—10月，空間上為東部降水多于西部。土壤類型主要有暗棕壤、黑土、沼澤土、白漿土、草甸土。

為了獲取地面實(shí)際蒸散發(fā)數(shù)據(jù)，在三江平原859農(nóng)場(chǎng)內(nèi)建立實(shí)驗(yàn)區(qū)。實(shí)驗(yàn)區(qū)位置為47°24′20.42″—47°24′21.38″N，134°07′11.12″—134°07′12.84″E。試驗(yàn)田面積600m2，主要作物為小麥和大豆。區(qū)內(nèi)架設(shè)有RR－9310換位式波文比通量觀測(cè)系統(tǒng)。通量觀測(cè)系統(tǒng)采樣頻率設(shè)置為10min，觀測(cè)時(shí)間為2011年6—9月，觀測(cè)因子包括地表凈輻射、土壤熱通量、不同高度層的大氣溫度和相對(duì)濕度。

2 研究方法與數(shù)據(jù)處理

2.1 研究方法

2.1.1 實(shí)際蒸散量 FAO P－M模型通過參考作物蒸散量和作物特性對(duì)需水量的影響即作物系數(shù)來計(jì)算潛在蒸散發(fā)量，然后結(jié)合實(shí)際土壤水分狀況計(jì)算實(shí)際蒸散量。實(shí)際蒸散量和潛在蒸散量表示為［11］：

式中：ET0——參考作物蒸散量（mm／d）；ETp——潛在蒸散量（mm／d）；ETa——實(shí)際蒸散量（mm／d）；Ks——土壤水分脅迫系數(shù)；Kc——作物系數(shù)，與土壤、氣候、作物種類、作物生長(zhǎng)狀況等因素有關(guān)，在沒有試驗(yàn)資料的情況下，可采用聯(lián)合國(guó)糧農(nóng)組織（FAO）給出的不同作物各發(fā)育階段作物系數(shù)經(jīng)驗(yàn)值。

Jensen等［12］的研究結(jié)果表明，土壤水分脅迫系數(shù)Ks可用式（4）表示：

式中：Av——相對(duì)有效含水率；Wm——萎蔫含水量（m3／m3），通過全球土壤數(shù)據(jù)庫(kù)HWSD查詢表層土壤的機(jī)械組成和有機(jī)碳含量，然后帶入SPAW模型計(jì)算獲得；Wf——田間持水量（m3／m3），同樣采用SPAW模型進(jìn)行計(jì)算；W——根區(qū)實(shí)際土壤含水量（m3／m3），采用遙感驅(qū)動(dòng)的分布式時(shí)變?cè)鲆嫠哪Ｐ停≧S－DTVGM）進(jìn)行估算［13］。

2.1.2 參考蒸散量 根據(jù)FAO的推薦，參考下墊面為表面開闊、具有充足水分供應(yīng)、植被高度為0.12 m的草地，具有固定的表面阻抗70s／m，反照率為0.23。參考蒸散量的計(jì)算式為［11］：

式中：ET0——參考蒸散量（mm／d）；Δ——?dú)鉁豑 時(shí)的飽和水汽壓曲線斜率（kPa／℃）；Rn——地表凈輻射［MJ／（m2·d）］；G——土壤熱通量［MJ／（m2·d）］；γ——干濕表常數(shù)（kPa／℃）；T——日平均溫度（℃）；U2——2m 高度處風(fēng)速（m／s）；es——飽和水汽壓（kPa）；ed——實(shí)際水汽壓（kPa）。

飽和水汽壓曲線斜率Δ：

式中：γ——干濕表常數(shù)（kPa／℃）；Cp——空氣定壓比熱，指一定氣壓下，單位體積的空氣溫度升高1℃所需的能量為1.013×10－3MJ／（kg·℃）；λ——蒸發(fā)潛熱，取2.45MJ／kg；ε——水汽分子量與干空氣分子量之比，為0.622。

地表熱通量G采用Su提出的計(jì)算方法［14］，對(duì)于有植被覆蓋的地面，計(jì)算公式如下：

式中：全植被覆蓋下，土壤熱通量與凈輻射的比值τc＝0.05；裸地情況下，土壤熱通量與凈輻射比值τs＝0.315；fc為植被蓋度，采用Nilson提出的計(jì)算方法［15］，表達(dá)式為：

式中：LAI——葉面積指數(shù)；k——與植被幾何結(jié)構(gòu)有關(guān)的系數(shù)；Ω——聚集指數(shù)，取值參見不同IGBP類型的典型聚集指數(shù)［16］；K——冠層消光系數(shù)，對(duì)于2a生或多年生樹木，假定葉片在空間為球狀分布，則其折射光的消光系數(shù)只取決于太陽(yáng)高度角或太陽(yáng)天頂角，表達(dá)式為：

式中：z——太陽(yáng)天頂角（°），Monsi認(rèn)為草本植物的K 為0.3～0.5［17］。

2.1.3 遙感驅(qū)動(dòng)的FAO P－M模型 遙感驅(qū)動(dòng)的FAO P－M模型計(jì)算實(shí)際蒸散發(fā)的流程如圖1所示。首先利用遙感數(shù)據(jù)估算出凈輻射、土壤熱通量、平均飽和水汽壓和實(shí)際水汽壓等參數(shù)，然后利用以上參數(shù)結(jié)合作物系數(shù)計(jì)算潛在蒸散量，最后通過土壤含水量計(jì)算出實(shí)際蒸散量。

圖1 基于多源遙感數(shù)據(jù)的P－M模型計(jì)算流程

2.2 數(shù)據(jù)處理

2.2.1 地面實(shí)際蒸散實(shí)驗(yàn)數(shù)據(jù)處理 波文比方法采用地表能量平衡方程和空氣動(dòng)力學(xué)方程來計(jì)算實(shí)際蒸散發(fā)，原理簡(jiǎn)單，精度較高，適用于常規(guī)地面觀測(cè)。由波文比計(jì)算蒸散耗熱通量的公式為［18］：

式中：E——蒸散量（mm）；λ——水的汽化潛熱（MJ／kg）；Rn——地表凈輻射［MJ／（m2·d）］；G——土壤熱通量［MJ／（m2·d）］；β為波文比。對(duì)由波文比方法計(jì)算得到的蒸散量進(jìn)行嚴(yán)格篩選，剔除與實(shí)際方向相反的數(shù)據(jù)和嚴(yán)重偏離正態(tài)分布的數(shù)據(jù)，對(duì)剔除后的缺值數(shù)據(jù)采用滑動(dòng)平均的方法進(jìn)行補(bǔ)值，然后將補(bǔ)值后的10min間隔的蒸散量累加轉(zhuǎn)化為日實(shí)際蒸散量。

2.2.2 多源遙感數(shù)據(jù)處理 本研究主要采用的遙感數(shù)據(jù)包括：（1）2010年環(huán)境衛(wèi)星影像，對(duì)影像數(shù)據(jù)進(jìn)行預(yù)處理，根據(jù)不同土地覆蓋類型的影像色調(diào)、紋理等特征，建立相應(yīng)的解譯標(biāo)志，通過人機(jī)交互式判讀，得到研究區(qū)的土地利用數(shù)據(jù)。土地利用類型主要包括：林地、草地、濕地、旱田、水田、裸地和水體。（2）采用“MODIS／Terra Leaf Area Index／FPAR 8－Day L4Global 1km SIN Grid V005”產(chǎn)品計(jì)算植被覆蓋度。產(chǎn)品空間分辨率為1km，時(shí)間分辨率為8 d，需要進(jìn)行時(shí)間序列插值得到日尺度數(shù)據(jù)。（3）采用GLDAS＿NOAH025SUBP＿3H 產(chǎn)品，包括氣溫、風(fēng)速、比濕、凈短波輻射和凈長(zhǎng)波輻射數(shù)據(jù)，分別用于計(jì)算日平均氣溫、日平均風(fēng)速、空氣相對(duì)濕度和凈輻射。由于GLDAS提供的數(shù)據(jù)是間隔25km的數(shù)據(jù)集點(diǎn)，不是連續(xù)的柵格數(shù)據(jù)，需要將原始數(shù)據(jù)提取生成1 km×1km分辨率的柵格數(shù)據(jù)，然后結(jié)合DEM做反距離加權(quán)插值，最后對(duì)不同類型的數(shù)據(jù)進(jìn)行求平均或累加得到日尺度數(shù)據(jù)，獲得最終需要的氣象數(shù)據(jù)。（4）降水?dāng)?shù)據(jù)采用FY－2D標(biāo)稱格式24h降水估計(jì)產(chǎn)品，空間分辨率為5.154km。（5）DEM數(shù)據(jù)采用SRTM 3數(shù)據(jù)，空間分辨率為90m。將不同空間分辨率的數(shù)據(jù)重采樣為1km×1km分辨率。

2.2.3 遙感數(shù)據(jù)精度評(píng)價(jià) 氣象數(shù)據(jù)來自GLDAS產(chǎn)品，有必要結(jié)合氣象站點(diǎn)數(shù)據(jù)和野外定點(diǎn)觀測(cè)數(shù)據(jù)對(duì)遙感氣象數(shù)據(jù)的精度進(jìn)行驗(yàn)證。驗(yàn)證數(shù)據(jù)來自由中國(guó)氣象局提供的三江平原7個(gè)基準(zhǔn)地面氣象觀測(cè)站的地面氣候資料日值數(shù)據(jù)和監(jiān)測(cè)站的波文比觀測(cè)數(shù)據(jù)。提取氣象站和地面觀測(cè)站所在像元的GLDAS數(shù)據(jù)的氣溫、風(fēng)速和凈輻射數(shù)據(jù)，分別與相應(yīng)站點(diǎn)的實(shí)測(cè)數(shù)據(jù)進(jìn)行比較。GLDAS氣溫?cái)?shù)據(jù)通過反距離加權(quán)并結(jié)合DEM插值后，與實(shí)測(cè)數(shù)據(jù)的一致性更好，GLDAS日均氣溫?cái)?shù)據(jù)與研究區(qū)內(nèi)各氣象站點(diǎn)的實(shí)測(cè)日均氣溫都呈現(xiàn)很好的相關(guān)性，R2均大于0.97；對(duì)比GLDAS風(fēng)速數(shù)據(jù)和研究區(qū)內(nèi)各氣象站點(diǎn)的風(fēng)速數(shù)據(jù)，發(fā)現(xiàn)GLDAS數(shù)據(jù)總體上較觀測(cè)值偏小，GLDAS風(fēng)速數(shù)據(jù)與各氣象站實(shí)測(cè)風(fēng)速的R2為0.416～0.608，相關(guān)性較好；將基于 GLDAS數(shù)據(jù)獲取的凈輻射與波文比實(shí)測(cè)的凈輻射數(shù)據(jù)進(jìn)行對(duì)比，發(fā)現(xiàn)二者比較接近，誤差在2.4～6.2MJ／m2之間，R2為0.704。以上驗(yàn)證結(jié)果表明，GLDAS產(chǎn)品精度較高，由GLDAS產(chǎn)品提取的氣象數(shù)據(jù)適用于三江平原蒸散發(fā)量的計(jì)算。

3 結(jié)果與分析

3.1 實(shí)際蒸散量估算結(jié)果驗(yàn)證

將P－M模型估算的旱田蒸散量與波文比實(shí)測(cè)的旱田蒸散量進(jìn)行比較?？紤]到影像的幾何糾正會(huì)帶來取樣誤差，觀測(cè)點(diǎn)處的遙感估算值采用以該像元為中心的3×3窗口內(nèi)的平均值作為該像元的遙感估算值，與實(shí)測(cè)值進(jìn)行比較。圖2為生長(zhǎng)季內(nèi)日尺度上旱田蒸散發(fā)的模擬值與實(shí)測(cè)值。由圖2可以看出：模擬值和實(shí)測(cè)值比較接近，旱田日蒸散發(fā)的模擬值與實(shí)測(cè)值的R2為0.824，RMSE為0.493；旱田月平均蒸散量的模擬值與實(shí)測(cè)值的R2為0.938，RMSE為0.256，P－M模型對(duì)旱田蒸散量的模擬精度較高。生長(zhǎng)季內(nèi)月尺度上，在6—8月份，旱田實(shí)際蒸散發(fā)的P－M模型的模擬值比實(shí)測(cè)值偏高，相對(duì)誤差分別為11.7%，7.63%和19.4%；9月份的模擬值則低于實(shí)測(cè)值，相對(duì)誤差為－0.13%。在實(shí)際蒸散量較小時(shí)，模型模擬誤差較大，可能是由于模擬值為像元尺度，雖然旱田實(shí)驗(yàn)區(qū)地勢(shì)平坦開闊，實(shí)測(cè)數(shù)據(jù)可以代表較大范圍的蒸散發(fā)情況，但仍存在尺度差異性，此外由于驅(qū)動(dòng)P－M模型的遙感產(chǎn)品的精度問題，會(huì)造成誤差累積?？傮w而言，在整個(gè)觀測(cè)時(shí)期內(nèi)，P－M模型的模擬值與實(shí)測(cè)值較為一致，表明本文所采用的多源遙感數(shù)據(jù)驅(qū)動(dòng)P－M模型的方法適用于日蒸散發(fā)的估算。

圖2 三江平原旱田日蒸散發(fā)的模擬值與實(shí)測(cè)值

3.2 三江平原蒸散發(fā)計(jì)算結(jié)果

采用多源遙感數(shù)據(jù)，對(duì)三江平原2011年生長(zhǎng)季內(nèi)（5—9月）月蒸散量的變化情況進(jìn)行估算。從時(shí)間上來看，三江平原的月蒸散發(fā)量分布呈單峰型，從5月到7月蒸散量呈顯著增長(zhǎng)趨勢(shì)，其中7月蒸散量最大，普遍達(dá)到60～100mm，8月蒸散量開始減少，整體上與6月相比差異不大，基本上達(dá)到40～80mm，9月蒸散量顯著減少，大部分地區(qū)的蒸散量介于20～40mm，與杜嘉［19］在三江平原采用SEBS模型估算蒸散發(fā)的結(jié)論相同。分析其原因?yàn)椋?月氣溫已經(jīng)回升，空氣飽和差變大，導(dǎo)致耗水增加；6—8月氣溫達(dá)到最高，降水量大，供水充分，是植被生長(zhǎng)的旺季，植被蒸騰顯著提高，導(dǎo)致蒸散發(fā)量大幅增加；9月氣溫開始緩慢降低，降水量減少，植株衰老，葉片功能喪失，蒸散發(fā)降低。

從空間分布來看，三江平原生長(zhǎng)季的蒸散發(fā)量與土地利用方式密切相關(guān)。位于以黑龍江、松花江、烏蘇里江及其支流和興凱湖為主的水域，蒸發(fā)量最高，生長(zhǎng)季的平均蒸發(fā)總量為445.72mm；其次為位于東北部地區(qū)的濕地和位于小興安嶺青黑山的林地覆蓋區(qū)，林地具有較高的蒸發(fā)系數(shù)，再加上該區(qū)降水充足，平均蒸散發(fā)量大于340mm；位于老爺嶺、張廣才嶺和完達(dá)山的林地覆蓋區(qū)，由于降水量較少，平均蒸散發(fā)量不足290mm；濕地由于有大面積水生植被覆蓋，且降水充沛，可提供充足的水分用于蒸騰，平均蒸散發(fā)量接近330mm；水田由于供水充足，生長(zhǎng)發(fā)育旺盛，平均蒸散發(fā)量為310mm；位于三江平原中部的旱田和草地，平均蒸散發(fā)量不足280mm，這是由于三江平原旱地大多不具備灌溉條件，大面積的旱田只能雨養(yǎng)，很多時(shí)候處于缺水狀態(tài)，故蒸散發(fā)量較低。總體上，2011年三江平原生長(zhǎng)季的蒸散發(fā)量主要與植被蓋度和水分供給狀況有關(guān)。三江平原水體的日蒸散量高于沼澤濕地的日蒸散量，與Burba［20］在美國(guó)內(nèi)布拉斯加州對(duì)水體和濕地蒸散發(fā)研究的結(jié)論相同。

3.3 實(shí)際蒸散量影響因素分析

在三江平原區(qū)域范圍內(nèi)，以5km為間隔進(jìn)行均勻采樣，共得到44個(gè)樣點(diǎn)。通過計(jì)算整個(gè)生長(zhǎng)季內(nèi)各樣點(diǎn)的實(shí)際蒸散量與4個(gè)氣象因子（氣溫、凈輻射、風(fēng)速和比濕）的相關(guān)系數(shù)，并對(duì)各氣象因子取其相關(guān)系數(shù)的平均值來衡量三江平原實(shí)際蒸散量與氣象因素的相關(guān)性。實(shí)際蒸散量與氣溫、凈輻射的相關(guān)性較好，其平均相關(guān)系數(shù)分別為0.49和0.63；實(shí)際蒸散量與比濕和風(fēng)速的相關(guān)性較差，其平均相關(guān)系數(shù)分別為0.32和0.14；將比濕轉(zhuǎn)換為相對(duì)濕度，與實(shí)際蒸散量呈負(fù)相關(guān)，平均相關(guān)系數(shù)為0.15。各氣象因子與實(shí)際蒸散發(fā)的相關(guān)性分析表明，凈輻射和氣溫是影響三江平原實(shí)際蒸散發(fā)的兩個(gè)主要因子。各氣象因子按照對(duì)蒸散發(fā)量的貢獻(xiàn)大小排序?yàn)椋簝糨椛洌練鉁兀颈葷瘢撅L(fēng)速。

圖3是三江平原生長(zhǎng)季內(nèi)觀測(cè)站點(diǎn)的實(shí)際蒸散發(fā)與氣溫、凈輻射和降水隨時(shí)間的變化規(guī)律，可以看出在降水較多時(shí)，氣溫和凈輻射與實(shí)際蒸散發(fā)的變化趨勢(shì)大體相同，實(shí)際蒸散發(fā)的峰值相對(duì)降水量呈現(xiàn)出一定的滯后性；7月15日到7月28日缺少降水，雖然期間氣溫和凈輻射值仍然比較高，潛在蒸散的平均值也達(dá)到4.87mm，但是由于土壤水分含量較低，實(shí)際蒸散量銳減，平均只有2.22mm，說明在氣溫、凈輻射和比濕較為穩(wěn)定的情況下，次降雨可以明顯增加蒸散量。三江平原降水時(shí)空分布不均勻，降水集中在6—10月，東多西少，而實(shí)際蒸散發(fā)的時(shí)空分布與降水的時(shí)空分布密切相關(guān)，因此降雨是影響區(qū)域?qū)嶋H蒸散發(fā)的制約性因素。

圖3 三江平原旱田日實(shí)際蒸散發(fā)和氣象要素的變化趨勢(shì)

4 結(jié)論

（1）采用P－M模型，結(jié)合多源遙感數(shù)據(jù)，對(duì)三江平原生長(zhǎng)季內(nèi)（5—9月）的日實(shí)際蒸散量進(jìn)行估算，并用波文比通量觀測(cè)系統(tǒng)得到的實(shí)測(cè)數(shù)據(jù)對(duì)模擬結(jié)果進(jìn)行了驗(yàn)證。結(jié)果表明，在整個(gè)觀測(cè)時(shí)期內(nèi)，P－M模型的模擬值與實(shí)測(cè)值比較一致，說明本文所采用的方法適用于日實(shí)際蒸散發(fā)的估算。

（2）蒸散發(fā)受土地利用和覆被類型的影響，同時(shí)與植被的生長(zhǎng)季節(jié)有很好的相關(guān)關(guān)系。生長(zhǎng)季內(nèi)，三江平原蒸散發(fā)量呈明顯的季節(jié)變化，總體上表現(xiàn)為5月份最低，7月份最高；蒸散發(fā)量的空間分布與植被蓋度和水分供給狀況密切相關(guān)。

（3）三江平原生長(zhǎng)季內(nèi)各氣象因子與實(shí)際蒸散發(fā)的相關(guān)性分析表明，凈輻射和氣溫是影響三江平原實(shí)際蒸散發(fā)的兩個(gè)主要因子。各氣象因子對(duì)實(shí)際蒸散量的貢獻(xiàn)大小為：凈輻射＞氣溫＞比濕＞風(fēng)速。此外，分析觀測(cè)站點(diǎn)氣象因素對(duì)實(shí)際蒸散發(fā)的影響得出，在其他氣象因素較為穩(wěn)定的情況下，次降雨可以明顯增加實(shí)際蒸散發(fā)量。三江平原降水的時(shí)空分布不均勻，降水是影響區(qū)域?qū)嶋H蒸散發(fā)的關(guān)鍵因素。

［1］ 武夏寧，胡鐵松，王修貴，等.區(qū)域蒸散發(fā)估算測(cè)定方法綜述［J］.農(nóng)業(yè)工程學(xué)報(bào)，2006，22（10）：257－262.

［2］ 劉安花，李英年，薛曉娟，等.高寒草甸蒸散量及作物系數(shù)的研究［J］.中國(guó)農(nóng)業(yè)氣象，2010，31（1）：59－64.

［3］ 左德鵬，徐宗學(xué)，李景玉，等.氣候變化情景下渭河流域潛在蒸散量時(shí)空變化特征［J］.水科學(xué)進(jìn)展，2011，22（4）：455－461.

［4］ Tong Ling，Kang Shaozhong，Zhang Lu.Temporal and spatial variations of evapotranspiration for spring wheat in the Shiyang river basin in northwest China［J］.Agricultural Water Management，2007，87（3）：241－250.

［5］ Suleiman A A，Tojo S C M，Hoogenboom G.Evaluation of FAO－56crop coefficient procedures for deficit irrigation management of cotton in a humid climate［J］.Agricultural and Forest Meteorology，2007，91（1／3）：33－42.

［6］ 卓瑪蘭草，劉普幸，張亞寧，等.甘肅黃土高原區(qū)潛在蒸散量時(shí)空變化與成因研究［J］.水土保持研究，2012，19（1）：70－76.

［7］ 侯蘭功，肖洪浪，鄒松兵，等.黑河流域水循環(huán)特征研究［J］.水土保持研究，2010，17（3）：254－258.

［8］ 曾麗紅，宋開山，張柏，等.松嫩平原參考作物蒸散量變化及其影響因素分析［J］.農(nóng)業(yè)系統(tǒng)科學(xué)與綜合研究，2010，26（1）：41－48.

［9］ Cleugh H A，Leuning R，Mu Q，et al.Regional evaporation estimates from flux tower and MODIS satellite data［J］.Remote Sensing of Environment，2007，106（3）：285－304.

［10］ 田雷，楊勝天，王玉娟.應(yīng)用遙感技術(shù)研究貴州春季蒸散發(fā)空間分異規(guī)律［J］.水土保持研究，2008，15（1）：87－91.

［11］ Allen R G，Pereira L S，Raes D，et al.Crop Evapotranspiration－guidelines for Computing Crop Water Requirements［R］.FAO Irrigation and Drainage Paper 56.FAO，1998.

［12］ Kristensen K J，Jensen S E.A model for estimating actual evapotranspiration from potential evapotranspiration［J］.Nordic Hydrology，1975，6（3）：170－188.

［13］ 楊勝天，等.生態(tài)水文模型與應(yīng)用［M］.北京：科學(xué)出版社，2012.

［14］ Su Z.The Surface Energy Balance System （SEBS）for estimation of turbulent heat fluxes［J］.Hydrology and Earth System Sciences Discussions，2002，6（1）：85－100.

［15］ Nilson T.A theoretical analysis of the frequency of gaps in plant stands［J］.Agricultural Meteorology，1971，8：25－38.

［16］ 唐世浩，朱啟疆，孫睿.基于方向反射率的大尺度葉面積指數(shù)反演算法及其驗(yàn)證［J］.自然科學(xué)進(jìn)展，2006，16（3）：331－337.

［17］ Monsi M，Saeki T.The light factor in plant communities and its significance for dry matter production［J］.Japanese Journal of Botany，1953，14：22－52.

［18］ 孫曉敏，朱治林，張仁華.生態(tài)系統(tǒng)中蒸發(fā)過程的精確測(cè)定：換位式波文比觀測(cè)儀介紹［J］.資源生態(tài)環(huán)境網(wǎng)絡(luò)研究動(dòng)態(tài)，1995，6（4）：44－47.

［19］ 杜嘉，張柏，宋開山，等.基于NOAA／AV HRR數(shù)據(jù)估算三江平原蒸散量研究初探［J］.水土保持研究，2009，16（2）：56－62.

［20］ Burba G G，Verma S B，Kim J.A comparative study of surface energy fluxes of three communities（Phragmites australis，Scirpus acutus，and open water）in a prairie wetland ecosystem［J］.Wetlands，1999，19（2）：451－457.