朱非林，王衛(wèi)光，孫一萌，鄭 強

(1.河海大學(xué)水文水資源與水利工程科學(xué)國家重點實驗室，江蘇 南京 210098;2.河海大學(xué)水文水資源學(xué)院，江蘇 南京 210098;3.河海大學(xué)水利水電學(xué)院，江蘇 南京 210098)

蒸散發(fā)在地球大氣系統(tǒng)的能量收支和水平衡變化中扮演著關(guān)鍵角色。作為能量平衡和水量平衡唯一的共同項，蒸散發(fā)是衡量水文循環(huán)和氣候變化的重要指標(biāo)。因此，確定蒸散發(fā)的長期變化趨勢已成為研究區(qū)域尺度氣候變化的水文響應(yīng)所關(guān)注的焦點［1-2］。表述蒸散發(fā)的變量主要包括水面蒸發(fā)、參考蒸散發(fā)和實際蒸散發(fā)。水面蒸發(fā)是指從蒸發(fā)皿觀察到的蒸發(fā)量，參考蒸散發(fā)表征充分供水下墊面條件下假想作物蒸散到空氣中的水量，而實際蒸散發(fā)則是直接參與水文循環(huán)的變量。

近幾年來，關(guān)于中國各大流域參考蒸散發(fā)和水面蒸發(fā)時空變化規(guī)律的研究已經(jīng)得到了廣泛開展。蔡輝藝等［3］結(jié)合敏感系數(shù)法，利用26個氣象站點的資料分析了淮河流域1960—2008年參考蒸散發(fā)的時空變化趨勢;Wang等［4］分析了海河流域1957—2008年參考蒸散發(fā)的變化規(guī)律及其影響因素;劉宏偉等［5］基于Penman-Monteith公式，評價了Hargreaves方法在太湖流域?qū)⒖颊羯l(fā)的模擬精度;Xu等［6］對長江流域1960—2000年蒸發(fā)皿蒸發(fā)和參考蒸散發(fā)的時空變化特征進(jìn)行研究，發(fā)現(xiàn)顯著下降的太陽輻射和風(fēng)速是造成長江流域參考蒸散發(fā)呈顯著下降趨勢的主要原因;雷慧閩等［7］采用生態(tài)水文模型分析了黃河下游大型引黃灌區(qū)1984—2006年蒸散發(fā)的變化規(guī)律;杜軍等［8］分析了藏北高原蒸發(fā)皿蒸發(fā)的變化特征，研究表明日照時數(shù)和平均氣溫日較差的顯著減小和相對濕度的明顯增加造成了蒸發(fā)皿蒸發(fā)呈現(xiàn)顯著的減少趨勢。

事實上，相對于參考蒸散發(fā)和水面蒸發(fā)而言，實際蒸散發(fā)是水平衡過程中的直接參與變量，因而更加受到關(guān)注。但由于受到氣候條件和地表狀況等眾多因素的共同作用，實際蒸散發(fā)被認(rèn)為是水文循環(huán)中最為復(fù)雜的環(huán)節(jié)，很難直接測出，因此通常由間接估算的方法得到。傳統(tǒng)的實際蒸散發(fā)估算方法(比如:波溫比和渦度相關(guān)法、SPAC水分傳輸綜合模擬法、定量遙感估算法等)由于所需的資料難以獲取和估算精度有限等，在實際應(yīng)用中受到了較大限制。Bouchet［9］于1963年提出了蒸散發(fā)互補相關(guān)原理，為實際蒸散發(fā)的估算開辟了一條新途徑。基于互補相關(guān)原理的概念，3種實際蒸散發(fā)計算模型(AA模型、GG模型、CRAE模型)被相繼提出。這3種模型僅需要常規(guī)的氣象數(shù)據(jù)，因此近年來被許多學(xué)者用于估算實際蒸散發(fā)［10-11］。然而，互補相關(guān)原理的成立需要一定的假設(shè)條件，加上蒸散發(fā)的各種影響因子隨氣候條件的不同會發(fā)生改變，因此3種模型的估算能力和適用性在不同的氣候類型區(qū)存在較大差別。Xu等［12］和劉紹民等［13］的研究均表明，基于3種模型原始參數(shù)的模擬效果較差，往往需要通過參數(shù)率定來提高估算精度，并且各模型的最優(yōu)參數(shù)值在不同的氣候類型區(qū)不一致。Gao等［14］利用參數(shù)率定后的AA模型和GG模型估算了海河流域的實際蒸散發(fā)，并基于Thornthwaite水量平衡法的計算結(jié)果對2種模型進(jìn)行評價，發(fā)現(xiàn)AA模型在海河流域的估算效果相對更優(yōu)。此外，劉波等［15］和劉健等［16］分別采用AA模型和GG模型估算長江流域和鄱陽湖流域的實際蒸散發(fā)，均取得了較好的模擬效果。

漢江是長江流域最長的支流，沿江分布有較大規(guī)模的工農(nóng)業(yè)，同時也是我國重要的農(nóng)業(yè)產(chǎn)區(qū)和商品糧基地。然而，漢江干流歷史上洪水災(zāi)害嚴(yán)重，是容易發(fā)生洪澇災(zāi)害的地區(qū)之一。研究漢江流域?qū)嶋H蒸散發(fā)的時空演變規(guī)律有助于進(jìn)一步理解氣候變化對水量平衡的影響，為流域的水資源配置和開發(fā)利用提供理論依據(jù)。目前，關(guān)于漢江流域?qū)嶋H蒸散發(fā)方面的深入研究還未見報道?；诖?，筆者分別利用參數(shù)率定后的AA模型和GG模型，根據(jù)漢江流域內(nèi)14個氣象站點1961—2011年的逐日氣象資料進(jìn)行實際蒸散發(fā)的模擬，基于水量平衡法的計算結(jié)果優(yōu)選出GG模型，分析了過去51a間實際蒸散發(fā)的時空演變規(guī)律和主要氣象因子的變化趨勢，并探討實際蒸散發(fā)演變的成因。

圖1 漢江流域及氣象站點分布示意圖Fig.1 Sketch map of location of Hanjiang River Basin and distribution of meteorological stations

1 研究區(qū)域和資料

漢江流域發(fā)源于陜西省漢中市嶓冢山，沿途流經(jīng)鄂、陜、豫、川、渝、甘6省市，位于東經(jīng)106°12'～114°14'、北緯30°08'～34°11'，干流全長1577km，流域面積15.9萬km2。漢江流域?qū)儆趤啛釒Ъ撅L(fēng)區(qū)，氣候比較溫和，年平均氣溫在15～17℃之間，水量較充沛。

研究資料包括向家坪水文站的徑流量資料和漢江流域內(nèi)14個氣象站點的逐日常規(guī)氣象資料。水文資料摘自《中華人民共和國水文年鑒》，時間序列為1961年1月至1986年12月;氣象資料由中國氣象局氣象信息中心氣象資料室提供，時間序列為1961年1月至2011年12月。資料包括日降水量、日最高氣溫、日最低氣溫、日平均氣溫、氣壓、相對濕度、日照時數(shù)。

2 研究方法

2.1 實際蒸散發(fā)模型

2.1.1 AA模型

Brutsaert等［17］根據(jù)Bouchet［9］的互補相關(guān)原理，提出了AA模型。AA模型中的參考蒸散發(fā)EP與濕潤環(huán)境蒸散發(fā)EW分別利用Penman公式和Priestly-Taylor公式計算。實際蒸散發(fā)Ea的計算公式為

式中:α——經(jīng)驗系數(shù)，其推薦取值為1.26;Δ——飽和水汽壓曲線斜率，kPa/℃;γ——干濕表常數(shù)，kPa/℃;Rn——地表凈輻射，MJ/(m2·d);M——土壤熱通量，MJ/(m2·d);ea——空氣干燥力，mm/d。

2.1.2 GG模型

Granger等［18］引進(jìn)相對蒸散發(fā)的概念，利用道爾頓蒸散定律推導(dǎo)出實際蒸散發(fā)和參考蒸散發(fā)的定量互補關(guān)系:

式中G為相對蒸散發(fā)參數(shù)，G=Ea/EP。

定義相對干燥力D=Ea/(Ea+Rn)，G與D之間存在以下關(guān)系:

式中a和b為經(jīng)驗系數(shù)，推薦取值分別為0.793和0.20。

2.2 Mann-Kendall檢驗和Sen’s坡度檢測法

Mann-Kendall檢驗［19-20］和Sen’s坡度檢測法［21］是基于秩和變化斜率的檢驗方法，其優(yōu)點在于不需要預(yù)先假定樣本的分布，不受異常值的干擾，能很好地處理非正態(tài)分布的數(shù)據(jù)。因此，這2種方法被廣泛地用于水文、氣象數(shù)據(jù)序列的變化趨勢研究［22-23］。同時，Mann-Kendall檢驗也是世界氣象組織積極推薦的非參數(shù)檢驗方法。

2.3 Pettitt突變檢測法

Pettitt于1979年提出的Pettitt突變檢測法［24］是一種非參數(shù)檢驗方法，其物理意義清晰，可以明確突變時間，已經(jīng)被廣泛地應(yīng)用于水文、氣象序列的突變研究［25］。該檢驗選用Mann-Whitney的統(tǒng)計函數(shù)Ut，T，認(rèn)為樣本x1，x2，…，xt，xt+1，…，xT獨立同分布(T為樣本容量)，對于連續(xù)序列，Ut，T和Vt，T的計算公式為

3 結(jié)果和討論

3.1 模型參數(shù)率定和實際蒸散發(fā)模擬

為了得到AA模型和GG模型在漢江流域的最優(yōu)參數(shù)值，選取漢江流域的子流域——洵河流域(面積為6448 km2)作為參數(shù)率定區(qū)域。

由洵河流域內(nèi)氣象站點1961—1986年的降雨量資料和向家坪水文站的同期徑流量資料求得多年平均降雨量P0=820.4mm，多年平均徑流深R0=312 mm，運用多年水量平衡方程E0=P0-R0得到洵河流域多年平均實際蒸散發(fā)E0=508.4 mm。采用模型的原始參數(shù)值進(jìn)行計算，AA模型和GG模型求得的E0分別為860.5 mm和597.3 mm，所得結(jié)果與多年水量平衡方程的計算結(jié)果均存在一定誤差，但GG模型的誤差明顯比AA模型小。因此，本文通過對2種模型參數(shù)的不斷調(diào)整，使模型估算出的E0與水量平衡法計算結(jié)果的閉合差盡可能小，最終將AA模型的參數(shù)α調(diào)整為0.965，將GG模型的參數(shù)a和b分別調(diào)整為1.30和0.21。利用參數(shù)率定后的AA模型和GG模型對洵河流域模擬的E0分別為514.2 mm和509.2 mm，與水量平衡方程計算結(jié)果的相對誤差僅為1.14%和0.16%，兩者的估算精度都有大幅度提高，并且表現(xiàn)出較好的適用性，但相對而言GG模型的模擬效果更優(yōu)。

從洵河流域?qū)嶋H蒸散發(fā)的年內(nèi)分布和年際變化(圖2)看，參數(shù)率定后AA模型和GG模型的模擬結(jié)果表現(xiàn)出高度的一致性，各月份的實際蒸散發(fā)最大僅相差1.9 mm，年實際蒸散發(fā)的變化趨勢基本一致。該結(jié)果進(jìn)一步驗證了利用參數(shù)率定后的兩種模型對漢江流域?qū)嶋H蒸散發(fā)進(jìn)行模擬的可靠性。

圖2 參數(shù)率定后AA模型和GG模型計算出的洵河流域?qū)嶋H蒸散發(fā)的年內(nèi)分布和年際變化對比Fig.2 Comparison of monthly and annual average actual evapotranspiration calculated by AA model and GG model using calibrated parameter values for Xunhe River Basin

根據(jù)漢江流域14個站點1961—2011年的逐日實測氣象資料，分別利用參數(shù)率定后的AA模型和GG模型對整個漢江流域的實際蒸散發(fā)進(jìn)行模擬。從年際變化(圖3)看，AA模型和GG模型估算出的年實際蒸散發(fā)的變化趨勢基本一致，各年份的實際蒸散發(fā)也大致相等。為了便于研究，本文根據(jù)GG模型的計算結(jié)果對漢江流域過去51a的實際蒸散發(fā)時空演變規(guī)律進(jìn)行分析。

圖3 參數(shù)率定后AA模型和GG模型計算出的漢江流域?qū)嶋H蒸散發(fā)的年際變化對比Fig.3 Comparison of annual average actual evapotranspiration calculated by AA model and GG model using calibrated parameter values for Hanjiang River Basin

3.2 漢江流域?qū)嶋H蒸散發(fā)的空間變化特征

本文利用反距離權(quán)重插值法分析漢江流域1961—2011年平均實際蒸散發(fā)的空間分布特征(圖4)。從總體上看，漢江流域多年平均實際蒸散發(fā)呈現(xiàn)從東南向西北逐漸減少的趨勢，其中位于漢江下游的鐘祥站和天門站的多年平均實際蒸散發(fā)在14個氣象站點中較大，分別為607.9 mm和607.5 mm;漢江上游各氣象站點的多年平均實際蒸散發(fā)相對較小，最小值出現(xiàn)在佛坪站，其多年平均實際蒸散發(fā)為484.5 mm。

圖4 漢江流域多年平均實際蒸散發(fā)的空間分布Fig.4 Spatial distribution of annual average actual evapotranspiration over Hanjiang River Basin

分別利用Mann-Kendall檢驗和Sen’s坡度檢測法對漢江流域14個氣象站點1961—2011年的逐年實際蒸散發(fā)進(jìn)行趨勢分析，圖5給出了Mann-Kendell檢驗的統(tǒng)計量Z值和Sen’s坡度值的空間分布。在過去的51 a間，除鎮(zhèn)安氣象站的年實際蒸散發(fā)呈顯著上升趨勢(置信度超過95%)外，其余13個站點均呈現(xiàn)下降趨勢，其中8個站點的下降趨勢較顯著(置信度超過95%)，主要分布在流域東部。Sen’s坡度檢測法的結(jié)果表明:流域東部的Sen’s坡度值較小，坡度范圍集中在-1.0～-1.5 mm/a;流域中部的Sen’s坡度值相對較大，局部地區(qū)甚至出現(xiàn)了正的上升趨勢，其中鎮(zhèn)安站的坡度值為0.812 mm/a，是所有站點中唯一出現(xiàn)正Sen’s坡度值的站點，這與Mann-Kendall檢驗的結(jié)果一致。

圖5 Z值和Sen’s坡度值的空間分布Fig.5 Spatial distribution of Z values and Sen’s slope values

3.3 漢江流域?qū)嶋H蒸散發(fā)的時間變化特征

圖6 漢江流域1961—2011年實際蒸散發(fā)的年際變化Fig.6 Annual variation of actual evapotranspiration over Hanjiang River Basin from 1961 to 2011

分別采用簡單線性回歸模型和二元回歸模型得出漢江流域?qū)嶋H蒸散發(fā)的回歸曲線(圖6)?？傮w上看，過去51a間漢江流域的實際蒸散發(fā)呈較明顯的下降趨勢，下降斜率為-0.671 4 mm/a。二元回歸曲線顯示漢江流域?qū)嶋H蒸散發(fā)的年際變化可大致分為2個階段:1961—1991年下降趨勢較顯著，曲線斜率較小;1991年之后下降幅度明顯減小，2001—2011年間呈現(xiàn)微弱的上升趨勢。

為了進(jìn)一步分析漢江流域?qū)嶋H蒸散發(fā)的年際變化規(guī)律(包括趨勢、坡度和突變)，分別采用Mann-Kendall檢驗法、Sen’s坡度檢測法和Pettitt突變檢測法對該流域1961—2011年的實際蒸散發(fā)序列進(jìn)行檢測，結(jié)果表明，過去51 a漢江流域年實際蒸散發(fā)呈現(xiàn)顯著的下降趨勢(置信度為99%)，其Sen’s坡度值為-0.6913 mm/a，與一元線性回歸曲線的斜率(-0.671 4 mm/a)基本一致。圖7(a)顯示年實際蒸散發(fā)在1979年發(fā)生顯著突變，P值為3.44×10-4(遠(yuǎn)小于99%置信度的閾值0.01)，Kt值為424。圖7(b)表明突變前后的Ea均值由575.5 mm降至547.3 mm，降幅達(dá)5%。

圖7 漢江流域1961—2011年實際蒸散發(fā)的Pettitt檢測結(jié)果Fig.7 Results of Pettitt test of actual evapotranspiration over Hanjiang River Basin from 1961 to 2011

從漢江流域的逐站點Pettitt檢測結(jié)果(圖8)看，流域內(nèi)所有14個站點的年實際蒸散發(fā)均呈顯著突變趨勢(置信度為95%)，其中有8個站點的突變達(dá)到了99%置信度，主要分布在流域的東南部和西北部地區(qū)。此外，共有11個站點的突變年份為1979年，占總數(shù)的78.6%。

3.4 實際蒸散發(fā)變化的成因分析

受供水條件和能量因素共同作用，流域?qū)嶋H蒸散發(fā)的機理一般比較復(fù)雜，由于土壤含水量資料往往難以搜集，因此以降水量代替［26］。研究表明，氣象因子變化是造成濕潤地區(qū)實際蒸散發(fā)時空波動的主要原因［27］?；诖?，為了探討漢江流域1961—2011年實際蒸散發(fā)時空演變的成因，筆者分別采用Mann-Kendall檢驗法、Sen’s坡度檢測法和Pettitt突變檢測法對風(fēng)速、氣溫、日照時數(shù)、氣壓、相對濕度、降水量這6個氣象因子進(jìn)行檢測。研究結(jié)果表明:(a)過去51 a間漢江流域的年平均風(fēng)速、年平均日照時數(shù)和年平均氣壓均呈現(xiàn)顯著的下降趨勢(置信度為99%)，并且分別于1981年、1979年和1993年發(fā)生顯著突變(置信度為99%);(b)在全球氣候變暖的背景下，漢江流域的年平均氣溫也顯著升高(置信度為99%)，其Sen’s坡度值為0.016℃/a，在1993年發(fā)生顯著突變(置信度為99%);(c)年平均相對濕度和年降水量分別呈微弱的上升趨勢和微弱的下降趨勢，兩者未表現(xiàn)出顯著的突變特征。

圖8 各氣象站點實際蒸散發(fā)的突變年份Fig.8 Years of mutation of actual evapotranspiration at 14 meteorological stations

風(fēng)速、氣溫、日照時數(shù)和降水量對實際蒸散發(fā)均產(chǎn)生正作用，而相對濕度和氣壓對實際蒸散發(fā)則產(chǎn)生負(fù)作用［24］。但是過去51 a間氣溫的顯著上升和氣壓的顯著下降并沒有導(dǎo)致漢江流域?qū)嶋H蒸散發(fā)上升，究其原因，主要是風(fēng)速和日照時數(shù)顯著下降對實際蒸散發(fā)的減弱作用補償了氣溫和氣壓引起的實際蒸散發(fā)增量，從而造成實際蒸散發(fā)呈顯著的下降趨勢。各氣象因子的Pettitt突變檢測結(jié)果也表明風(fēng)速和日照時數(shù)的突變年份(1981年、1979年)與實際蒸散發(fā)的突變年份(1979年)幾乎一致，兩者的顯著下降是導(dǎo)致漢江流域1961—2011年實際蒸散發(fā)顯著下降的主要原因。

4 結(jié) 語

對漢江流域1961—2011年實際蒸散發(fā)的時空演變規(guī)律和主要氣象因子的變化趨勢進(jìn)行分析，并探討了實際蒸散發(fā)變化的成因，得到以下主要結(jié)論:

a.參數(shù)率定后的AA模型和GG模型均表現(xiàn)出較好的適用性，兩種模型對整個漢江流域?qū)嶋H蒸散發(fā)的模擬結(jié)果基本一致。相對而言，GG模型的模擬效果更優(yōu)。

b.漢江流域1961—2011年平均實際蒸散發(fā)總體上呈現(xiàn)自東南向西北逐次減少的空間分布特征;除鎮(zhèn)安站的年實際蒸散發(fā)呈顯著上升趨勢外，其余站點均呈下降趨勢，其中8個站點的下降趨勢較顯著，主要分布在流域東部。

c.漢江流域1961—2011年平均實際蒸散發(fā)呈現(xiàn)顯著的下降趨勢，其Sen’s坡度值為-0.691 3 mm/a，并于1979年發(fā)生顯著突變，突變前后的均值由575.5 mm降至547.3 mm，降幅達(dá)5%。

d.過去的51 a間，風(fēng)速和日照時數(shù)顯著下降對實際蒸散發(fā)的減弱作用補償了氣溫和氣壓引起的實際蒸散發(fā)的增量，是導(dǎo)致漢江流域年實際蒸散發(fā)顯著下降的主要原因。

［1］WANGWeiguang，PENGShizhang，YANGTao，et al.Spatial and temporal characteristics of reference evapotranspiration trends in the Haihe River Basin，China［J］.Journal of Hydrologic Engineering，2011，16:239-252.

［2］JHAJHARIA D，DINPASHOH Y，KAHYA E，et al.Trends in reference evapotranspiration in the humid region of northeast India［J］.Hydrol Processes，2012，26:421-435.

［3］蔡輝藝，余鐘波，楊傳國，等.淮河流域參考蒸散發(fā)量變化分析［J］.河海大學(xué)學(xué)報:自然科學(xué)版，2012，40(1):76-82.(CAI Huiyi，YU Zhongbo，YANG Chuanguo，et al.Analysis of variation of reference evapotranspiration in Huaihe Basin［J］.Journal of Hohai University:Natural Sciences，2012，40(1):76-82.(in Chinese))

［4］WANGWeiguang，SHAOQuanxi，PENGShizhang，et al.Reference evapotranspiration change and the causes across the Yellow River Basin during 1957—2008 and their spatial and seasonal differences［J］.Water Resource Research，2012，48:1-27.

［5］劉宏偉，余鐘波.用Hargreaves法與Penman-Monteith法計算ET0:以太湖流域的應(yīng)用為例［J］.水資源保護(hù)，2010，26(1):6-8.(LIU Hongwei，YU Zhongbo.Application of Hargreaves and Penman-Monteith Equation to estimating ET0:a case study in Taihu Basin［J］.Water Resources Protection，2010，26(1):6-8.(in Chinese))

［6］XU CY，GONG L，JIANG T，et al.Analysis of spatial distribution and temporal trend of reference evapotranspiration and pan evaporation in Changjiang(Yangtze River)catchment［J］.Journal of Hydrology，2006，327:81-93.

［7］雷慧閩，蔡建峰，楊大文，等.黃河下游大型引黃灌區(qū)蒸散發(fā)長期變化特性［J］.水利水電科技進(jìn)展，2012，32(1):13-17.，YANG Dawen，et al.Long-term variability of evapotranspiration in a large irrigated area in lower reach of Yellow River［J］.Advances in Science and Technology of Water Resources，2012，32(1):13-17.(in Chinese))

［8］杜軍，邊多，鮑建華，等.藏北高原蒸發(fā)皿蒸發(fā)量及其影響因素的變化特征［J］.水科學(xué)進(jìn)展，2008，19(6):786-791.(DU Jun，BIAN Duo，BAOJianhua，et al.Changes of pan evaporations and its impact factors over northern Tibet in 1971-2006［J］.Advances in Water Science，2008，19(6):786-791.(in Chinese))

［9］BOUCHET R J.Evapotranspiration reele et potentielle，signification climatique［J］.General Assembly Berkeley，Int Ass Sci Hydrol，1963，62:134-142.

［10］趙玲玲，陳喜，夏軍，等.氣候變化下烏江流域蒸散發(fā)互補關(guān)系變化及成因辨識［J］.河海大學(xué)學(xué)報:自然科學(xué)版，2011，39(6):629-634.(ZHAO Lingling，CHEN Xi，XIA Jun，et al.Variation of complementary relationship of evapotranspiration with climate changes and its cause analysis in Wujiang Basin［J］.Journal of Hohai University:Natural Sciences，2011，39(6):629-634.(in Chinese))

［11］王艷君，姜彤，劉波.長江流域?qū)嶋H蒸散發(fā)量的變化趨勢［J］.地理學(xué)報，2010，65(9):1079-1088.(WANG Yanjun，JIANG Tong，LIU Bo.Trends of estimated and simulated actual evapotranspiration in the Yangtze River Basin［J］.Acta Geographica Sinica，2010，65(9):1079-1088.(in Chinese))

［12］XU C Y，SINGH V P.Evaluation of three complementary relationship evapotranspiration models by water balance approach to estimate actual regional evapotranspiration in different climatic regions［J］.Journal of Hydrology，2005，308:105-121.

［13］劉紹民，孫睿，孫中平，等.基于互補相關(guān)原理的區(qū)域蒸散量估算模型比較［J］.地理學(xué)報，2004，59(3):331-340.(LIU Shaomin，SUN Rui，SUN Zhongping，et al.Comparison of different complementary relationship models for regional evapotranspiration estimation［J］.Acta Geographica Sinca，2004，59(3):331-340.(in Chinese))

［14］GAOGe，XU Chongyu，CHEN Deliang.Spatial and temporal characteristics of actual evapotranspiration over Haihe River Basin in China［J］.Stochastic Environmental Research and Risk Assessment，2012，26(5):655-669.

［15］劉波，翟建青，高超，等.基于實測資料對日蒸散發(fā)估算模型的比較［J］.地球科學(xué)進(jìn)展，2010，25(9):974-980.(LIU Bo，ZHAI Jianqing，GAO Chao，et al.A comparison of daily actual evapotranspiration evaluation models based on field observational data［J］.Advances in Earth Science，2010，25(9):974-980.(in Chinese))

［16］劉健，張奇，許崇育，等.近50年鄱陽湖流域?qū)嶋H蒸散發(fā)量的變化及影響因素［J］.長江流域資源與環(huán)境，2010，19(2):139-145.(LIU Jian，ZHANGQi，XU Chongyu，et al.Change of actual evapotranspiration of Poyang Lake watershed and associated influencing factors in the past 50 years［J］.Resources and Environment in the Yangtze Basin，2010，19(2):139-145.(in Chinese))

［17］BRUTSAERT W，STRICKER H.An advection-aridity approach to estimate actual regional evapotranspiration［J］.Water Resource Research，1979，15(2):443-450.

［18］GRANGER R J，GRAY D M.Evaporation from natural nonsaturated surfaces［J］.Journal of Hydrology，1989，111:9-19.

［19］MANN H B.Nonparametric tests against trend［J］.Journal of the Econometric Society，1945，133(3):245-259.

［20］KENDALL M G.Rank correlation methods［M］.London:Griffin Press，1970.

［21］SEN PK.Estimates of the regression coefficient based on Kendall’s Tau［J］.Journal of the American Statistical Association，1968，63:1379-1389.

［22］楊肖麗，任立良，江善虎，等.西遼河源頭流域徑流變化趨勢及影響因素分析［J］.河海大學(xué)學(xué)報:自然科學(xué)版，2012，40(1):37-41.(YANG Xiaoli，REN Liliang，JIANG Shanhu，et al.Analysis of trend and influencing factors of runoff variation in headwater basin of West Liaohe River［J］.Journal of Hohai University:Natural Sciences，2012，40(1):37-41.(in Chinese))

［23］陸志華，夏自強，于嵐嵐，等.近51年松花江流域氣溫時空變化特征［J］.河海大學(xué)學(xué)報:自然科學(xué)版，2012，40(6):629-635.(LU Zhihua，XIA Ziqiang，YU Lanlan，et al.Characteristics of spatio-temporal variation of temperature in Songhua River Basin from 1960 to 2010［J］.Journal of Hohai University:Natural Sciences，2012，40(6):629-635.(in Chinese))

［24］PETTITT A N.A non-parametric approach to the change point problem［J］.Applied Statistics，1979，28(2):126-135.

［25］ZHANG Shurong，LU X X.Hydrological responses to precipitation variation and diverse human activities in a mountainous tributary of the lower Xijiang，China［J］.Catena，2009，77(2):130-142.

［26］GAO Ge，CHEN Deliang，XU Chongyu，et al.Trend of estimated actual evapotranspiration over China during 1960-2002［J］.Journal of Geophysical Research，2007，112(D11):1-8.

［27］COHEN S，STANHILL G.Evaporative climate changes at Bet-Dagan，Israel，1964-1998［J］.Agricultural and Forest Meteorology，2002，111(2):83-91.(LEI Huimin ，CAI Jianfeng