靖娟利,蔡江濤,2a,耿仁方,2b,王永鋒

(1.桂林理工大學(xué) a.測(cè)繪地理信息學(xué)院;b.廣西空間信息與測(cè)繪重點(diǎn)實(shí)驗(yàn)室,廣西 桂林 541006;2.南京信息工程大學(xué)a.地理科學(xué)學(xué)院；b.遙感與測(cè)繪工程學(xué)院，南京 210044)

0 引 言

地表蒸散(evapotranspiration, ET)是指在不同下墊面的氣候要素影響下,特定區(qū)域內(nèi)水分蒸發(fā)和散發(fā)的總和,包括地表水分蒸發(fā)和植物體內(nèi)水分的蒸騰。地表蒸散是全球水和能量循環(huán)的主要成分,轉(zhuǎn)移了陸地上約2/3的降水量,維持著地表能量的平衡,對(duì)監(jiān)測(cè)氣候變化和改善水資源管理具有重要意義[1],受到國(guó)內(nèi)外學(xué)者的廣泛關(guān)注。

傳統(tǒng)的蒸散研究方法大都基于“點(diǎn)” 尺度[2-4], 觀測(cè)成本高、 代表性差, 難以直接應(yīng)用到較大范圍的區(qū)域尺度。 遙感技術(shù)提供了大范圍地表特征信息, 具有較好的時(shí)效性和區(qū)域性特點(diǎn), 為大尺度非均勻下墊面的蒸散發(fā)監(jiān)測(cè)提供了新途徑[5]。 2011年美國(guó)NASA 團(tuán)隊(duì)基于Penman-Monteith 遙感模型和MODIS 數(shù)據(jù)研發(fā), 發(fā)布了全球陸地蒸散數(shù)據(jù)集產(chǎn)品MOD16, 該數(shù)據(jù)集已由全球通量塔數(shù)據(jù)驗(yàn)證, 模擬精度達(dá)到86%[6-7], 在全球范圍內(nèi)得到廣泛應(yīng)用。

珠江流域處于我國(guó)南方亞熱帶濕潤(rùn)地區(qū),其上游地處滇、黔、桂3省交匯處,是我國(guó)喀斯特環(huán)境的集中分布區(qū),自然生態(tài)環(huán)境脆弱,其人口、資源、環(huán)境與可持續(xù)發(fā)展矛盾突出[8]。近年來(lái),珠江流域旱澇災(zāi)害頻繁,造成了巨大經(jīng)濟(jì)損失[9]。因此,開展珠江流域地表蒸散發(fā)時(shí)空分布研究,在水資源管理、旱澇災(zāi)害預(yù)警等方面具有重要的應(yīng)用價(jià)值。

近年來(lái),國(guó)內(nèi)研究人員基于MOD16產(chǎn)品對(duì)鄱陽(yáng)湖流域[10]、中國(guó)[11]、陜西省[12]、三江平原[13]、渭河流域[14]、淮河流域[5]、漢江流域[15]等開展了不同尺度的地表蒸散發(fā)時(shí)空分布特征評(píng)估。國(guó)內(nèi)學(xué)者也相繼對(duì)珠江流域地表蒸散量時(shí)空特征進(jìn)行了研究：劉昌明等[16]采用優(yōu)化太陽(yáng)輻射計(jì)算的Penman-Monteith潛在ET計(jì)算方法, 分析發(fā)現(xiàn)珠江流域潛在ET呈下降趨勢(shì), 潛在ET對(duì)最高氣溫的變化最為敏感; 李修倉(cāng)等[17]采用水量平衡模型和Penman公式分別計(jì)算了珠江流域7個(gè)子流域1961—2000年實(shí)際ET和潛在ET。 但這些研究都是基于站點(diǎn)氣象數(shù)據(jù)進(jìn)行空間內(nèi)插, 進(jìn)而研究珠江流域的潛在ET或?qū)嶋HET的趨勢(shì), 其結(jié)果空間代表性差, 不能精準(zhǔn)反映區(qū)域尺度ET趨勢(shì)走向, 同時(shí)缺乏對(duì)流域內(nèi)ET的空間分布特征的探討。

本文以空間分辨率為1 km的MOD16A2月數(shù)據(jù)和MOD16A3年尺度數(shù)據(jù)為基礎(chǔ), 分析珠江流域蒸散發(fā)時(shí)空變化特征。 此外, 由于不同土地利用類型本身的生理生態(tài)特性及其所處區(qū)域的水熱條件差異, 其平均蒸散量分布特征會(huì)表現(xiàn)出不同的變化特點(diǎn)[18]。 分別提取了珠江流域各種土地利用類型年均ET值, 并分析其變化特征。 通過(guò)以上研究, 可以為認(rèn)識(shí)流域旱澇災(zāi)害成因、 加強(qiáng)生態(tài)保護(hù)、 合理配置水資源等方面提供科學(xué)的理論指導(dǎo)依據(jù)。

1 研究區(qū)概況

本文所研究的珠江流域范圍在北緯21°31′—26°49′ 、東經(jīng)102°14′—115°53′,流域覆蓋滇、桂、粵、黔、湘、贛等省區(qū),總面積達(dá)44.78萬(wàn)km2,是中國(guó)流域面積第三、流量第二的南方大河流域,主要的支流包括西江、北江和東江。珠江流域地勢(shì)大體上西高東低、北高南低(圖1),從西到東橫跨云貴高原、兩廣丘陵和珠江三角洲平原,流域內(nèi)山地、丘陵面積占94.4%,平原面積僅占5.6%。流域地處濕熱多雨的熱帶、亞熱帶氣候區(qū),年平均氣溫在14～22 ℃,年際變化不大,但地區(qū)差異大；流域內(nèi)雨量豐沛,多年平均降水量在1 200～2 200 mm,降水量由東向西遞減,降雨季節(jié)分配不均,地區(qū)差異和年際變化較大,致使流域洪、澇、旱等自然災(zāi)害頻繁。

2 研究方法

2.1 數(shù)據(jù)來(lái)源及預(yù)處理

數(shù)據(jù)包括MOD16地表蒸散數(shù)據(jù)、土地利用數(shù)據(jù)和流域邊界數(shù)據(jù),詳細(xì)數(shù)據(jù)來(lái)源見(jiàn)表1。MOD16產(chǎn)品是基于Penman-Monteith方程,結(jié)合MODIS的日常氣象再分析資料和8 d遙感植被屬性動(dòng)態(tài)信息得到的,其空間分辨率為1 km。本文選取2000年1月—2014年12月連續(xù)15年的MOD16A3年產(chǎn)品和MOD16A2月產(chǎn)品,利用MRT(MODIS Reprojection Tool)工具將原始的HDF數(shù)據(jù)進(jìn)行投影轉(zhuǎn)換和拼接, 將原有的正弦曲線(sinusoidal)投影轉(zhuǎn)換為等積圓錐(albers equal area conic)投影，最后在ArcGIS 10.2 平臺(tái)中去除數(shù)據(jù)中的無(wú)效值，并進(jìn)行裁剪操作得到研究區(qū)1 km分辨率的ET。結(jié)合珠江流域?qū)嶋H情況,從中國(guó)1∶100萬(wàn)植被類型空間分布數(shù)據(jù)中選擇林地、灌木、草地、荒地和耕地5類土地利用類型(圖2)。從全球流域數(shù)據(jù)庫(kù)GDBD 中提取珠江流域邊界,并分為北盤江流域、南盤江流域、紅水河流域等9個(gè)子流域。

表1 研究所用數(shù)據(jù)來(lái)源Table 1 Data sources for research

圖2 研究區(qū)土地利用現(xiàn)狀Fig.2 Land use in study area

2.2 Theil-Sen 中值趨勢(shì)分析與 Mann-Kendall顯著性檢驗(yàn)分析

本文采用Theil-Sen中值趨勢(shì)分析法研究珠江流域年均ET時(shí)間序列的變化特征,并結(jié)合Mann-Kendall檢驗(yàn)法對(duì)變化趨勢(shì)的顯著性進(jìn)行檢驗(yàn)。

Theil-Sen 趨勢(shì)度ρ計(jì)算公式為[19]

(1)

嚴(yán)格來(lái)說(shuō),ρ=0的像元非常稀少, 故釆用±0.000 5為閾值, 劃分兩閾值之間的像元為穩(wěn)定區(qū)，ρ>0.000 5的像元為ET增加區(qū),ρ值<-0.000 5值的像元為ET減少區(qū)。

Mann-Kendall 檢驗(yàn)常廣泛應(yīng)用于降雨、氣溫和徑流等水文氣象要素的趨勢(shì)分析及變異檢驗(yàn),其優(yōu)點(diǎn)是無(wú)論序列是否屬于正態(tài)分布,或受少數(shù)異常值的干擾,都可用于變量的趨勢(shì)檢驗(yàn),計(jì)算方便,可用于變量的突變分析[20]。

Mann-Kendall檢驗(yàn)的計(jì)算公式為[21]

(2)

(3)

(4)

(5)

其中:Z為標(biāo)準(zhǔn)化后的檢驗(yàn)統(tǒng)計(jì)量;S為檢驗(yàn)統(tǒng)計(jì)量; sgn是符號(hào)函數(shù);θ=ETj-ETi;n為時(shí)間序列長(zhǎng)度(本研究中n=15)。 當(dāng)Z為正時(shí), 表示序列呈增加趨勢(shì);Z為負(fù)時(shí), 表示序列呈減少趨勢(shì);Z的絕對(duì)值越大, 表示序列變化趨勢(shì)越顯著。 在給定顯著性水平α下, 在正態(tài)分布表中查臨界值Z1-α/2,當(dāng)|Z|>|Z1-α/2|時(shí),認(rèn)為時(shí)間序列變化趨勢(shì)顯著。本文選擇置信水平α=0.05,Z1-α/2=1.96。

為了研究珠江流域2000—2014年間ET動(dòng)態(tài)變化趨勢(shì),本文基于 Theil-Sen 中值趨勢(shì)分析和 Mann-Kendall顯著性檢驗(yàn)方法,利用ArcGIS 10.2中的Arcpy編程工具對(duì)年均ET進(jìn)行逐像元編程,得到Theil-Sen 中值趨勢(shì)度ρ值,并進(jìn)行Mann-Kendal檢驗(yàn)得到檢驗(yàn)統(tǒng)計(jì)量Z值。綜合ρ和Z值,將ET變化趨勢(shì)劃分為顯著增加、輕微增加、基本穩(wěn)定、輕微減少、顯著減少5類。

3 結(jié)果與分析

3.1 蒸散量時(shí)間變化特征

由圖3可看出, 2000—2014年珠江流域ET年際波動(dòng)較大, 年均值變化范圍在941.00～1 002.93 mm/a, 變化率為0.32 mm/a, 整體呈微弱增加趨勢(shì)。 多年ET均值為959.33 mm/a, 明顯超出多年平均值的年份為2003和2007年, 其中2003年ET最高, 超出平均值43.60 mm, 相對(duì)變化率達(dá)到了4.54%; ET最小值出現(xiàn)在2010年, 為941.00 mm/a。 相關(guān)資料表明, 2003年珠江流域出現(xiàn)了罕見(jiàn)的持續(xù)高溫天氣, 降水嚴(yán)重偏少, 各地先后呈現(xiàn)夏秋連續(xù)干旱現(xiàn)象; 2007年珠江流域降雨量偏少嚴(yán)重, 流域內(nèi)旱情發(fā)展蔓延,局部區(qū)域旱情較嚴(yán)重[22],這可能是ET波動(dòng)大的重要影響因素。

圖3 珠江流域多年平均ET的年際變化Fig.3 Inter-annual variation of average annual ET in the pearl River Basin

圖4顯示了珠江流域ET的年內(nèi)變化特征,總體上呈先增大后減小的單峰型分布。ET主要集中在 5—10月份, 4—5月快速增長(zhǎng), 7月份達(dá)到峰值121.60 mm, 10—11月迅速下降, 1月降至年內(nèi)最低值 39.29 mm。 這主要是因?yàn)?、 2和12月份流域內(nèi)的氣溫相對(duì)較低, 不利于地表蒸散發(fā); 5—9月份是該流域的雨季, 水熱條件充分，再加上風(fēng)速大、 日照充足, 提供了有利于ET的條件; 10月后氣溫逐漸回落, 又向著不利于蒸散發(fā)的條件轉(zhuǎn)變, ET逐漸下降到最小值。 受流域內(nèi)海拔、 氣候因子、 太陽(yáng)輻射等多種因素的綜合影響[23], 珠江流域四季ET差異顯著,整體表現(xiàn)為:夏季(349.75 mm)>秋季(248.06 mm)>春季 (237.05 mm)>冬季(123.92 mm), 多年均值分別占年蒸散量的36.48%、 25.87%、 24.72%、 12.92%。 ET在春季(3—5月) 增長(zhǎng)迅速, 夏季(6—8月)達(dá)到峰值, 秋季(9—11月)急劇減小, 冬季(12—2月)變動(dòng)較平緩且全年最低。

圖4 珠江流域多年月平均ET年內(nèi)分布Fig.4 Distribution of average annual monthly ET in the Pearl River Basin

3.2 蒸散量空間分布特征

3.2.1 蒸散量多年平均及季節(jié)分布特征 珠江流域多年ET均值在401.2～1 665.1 mm/a, 表現(xiàn)出顯著的空間差異性, 整體上呈現(xiàn)東高西低、 南高北低,流域上游低、 中下游高的分布格局。

由圖5可知,流域內(nèi)多年平均ET峰值出現(xiàn)在北江流域, 達(dá)到1 030.3 mm/a,ET高值區(qū)集中分布在北江流域中南部、桂江流域與黔洵江及西江流域接壤地區(qū)、右江流域北部及左江流域南部等地區(qū),而低值區(qū)主要分布在北盤江流域和南盤江流域的西部,其ET值在400～600 mm/a。

圖5 珠江流域多年ET平均值空間分布Fig.5 Spatial distribution of average annual ET in the Pearl River Basin

ET空間分布格局主要受地形地貌、地表植被類型、大氣濕度等因子的影響。ET高值區(qū)主要位于珠江三角洲、廣西盆地,這些地區(qū)海拔一般低于200 m,地表植被主要為常綠闊葉林,而且南臨南海,水汽供應(yīng)充足。而ET低值區(qū)集中分布在流域源頭,位于云貴高原南緣,海拔大都在2 000 m以上,地處喀斯特地區(qū),石灰?guī)r廣泛分布,土層薄且土壤侵蝕退化嚴(yán)重,地表植被主要以亞熱帶天然次生植被和人工植被為主,蓄水能力較差,而且由于地形抬升,水汽供應(yīng)不足,造成多年平均ET較小。

為更好地反映ET的分布,將四季平均ET值每隔100 mm重分類成不同等級(jí)。珠江流域ET存在明顯的季節(jié)性差異。春季流域內(nèi)ET值在50～350 mm波動(dòng)(圖6a), 隨著溫度回升, 地表蒸散量逐漸增加。 蒸散量介于250～350 mm的高值區(qū)主要分布在兩廣丘陵地區(qū)海拔較高的山區(qū),這些地區(qū)主要分布林地,蒸散量普遍較高;而<150 mm的低值區(qū)集中分布在珠江源頭,這些地區(qū)處于云貴高原,海拔相對(duì)較高,遠(yuǎn)離海洋,空氣濕度相對(duì)流域其他部分較低,同時(shí)該地區(qū)耕地廣泛分布,春季尚未種植農(nóng)作物,因此ET低于流域其他部分。夏季降水豐沛,溫度隨之升高,流域內(nèi)ET達(dá)到全年最大值,低值區(qū)縮小至流域西北一角,紅水河流域的中部和右江流域北部形成高值區(qū)中心,高值區(qū)域的植被類型大都為林地,林地的葉面積指數(shù)遠(yuǎn)高于其他植被類型,夏季水熱條件優(yōu)異的情況下,植被蒸騰作用較強(qiáng),而植物蒸騰在蒸散中占據(jù)最大比重,因此產(chǎn)生高值集中區(qū)域。此外,右江流域高值區(qū)旁邊有一條明顯的西北—東南走向的狹長(zhǎng)低值帶,原因是該地區(qū)處于百色盆地,因此ET較周邊區(qū)域低,這與多年ET均值的空間分布也相符合(圖6b)。秋季,流域內(nèi)氣溫逐漸降低,地表蒸散量逐漸回落,ET的高值中心縮小至百色中部,ET的整體分布與春季相似,低值區(qū)零星分布在流域西部(圖6c)。冬季流域內(nèi)ET達(dá)到全年最小值,50～150 mm的低值區(qū)在流域內(nèi)廣泛分布，僅在東南沿海地區(qū)地表蒸散量相對(duì)較高(圖6d),原因是這些區(qū)域冬季水熱、太陽(yáng)輻射等相對(duì)于流域其他區(qū)域都更優(yōu)異。

圖6 珠江流域四季地表蒸散量空間分布Fig.6 Spatial distribution of land-surface evapotranspiration in different seasons of the Pearl River Basin

3.2.2 蒸散量空間變化趨勢(shì) 由表2、圖7和圖8可知,珠江流域近15年蒸散量在空間上總體呈增加趨勢(shì)。ET增加區(qū)域占54.13%,顯著增加區(qū)域占7.15%;ET減少的區(qū)域占43.55%,顯著減少的區(qū)域占4.68%;ET基本穩(wěn)定的區(qū)域僅占2.32%。ET_Sen趨勢(shì)度介于-57.3～44.2,經(jīng)統(tǒng)計(jì)分析發(fā)現(xiàn),低值區(qū)集中分布在南盤江流域、右江流域和北江流域,高值區(qū)散布于各流域。

云貴高原東南部和廣西盆地等地區(qū)ET變化趨勢(shì)復(fù)雜,這些區(qū)域是典型的喀斯特地區(qū),生態(tài)環(huán)境脆弱,植被覆蓋度較低,在人類工農(nóng)業(yè)生產(chǎn)活動(dòng)持續(xù)干擾下,改變下墊面的性質(zhì),ET波動(dòng)較大。ET顯著增加的區(qū)域主要分布于貴州境內(nèi),集中在安順、黔南和黔西南三者交界處,桂林和柳州的東南部地區(qū)也存在小范圍帶狀的顯著增長(zhǎng)區(qū),這些區(qū)域地表植被覆蓋良好,不同植被類型分布較為均衡。ET顯著減少的區(qū)域主要集中在云南曲靖西部、紅河哈尼族彝族自治州東部,廣西百色中部的西北東南走向的帶狀區(qū)域以及桂林、南寧等城市中心周圍的環(huán)狀區(qū)域,這些地區(qū)受人類活動(dòng)持續(xù)干擾, 地表植被覆蓋度低, 且植被類型較為單一。ET基本穩(wěn)定的像元零星分散在研究區(qū),圍繞在變化趨勢(shì)不顯著的像元。

圖7 ET_Sen變化趨勢(shì)Fig.7 Trends of ET_Sen

圖8 Mann-Kendall顯著性檢驗(yàn)Fig.8 Significance test of Mann-Kendall

3.2.3 不同土地利用類型蒸散量變化特征 為了揭示珠江流域不同土地利用類型蒸散特征,基于GIS中的區(qū)域統(tǒng)計(jì)功能,分別提取不同土地利用類型的年均蒸散量變化特征。

從圖9可知,不同土地利用類型年均ET值存在一定差異性,林地的年均ET均值遠(yuǎn)高于其他土地利用類型,達(dá)到1 011.22 mm/a,其次為草地、灌木、耕地,荒地年均ET僅為837.32 mm/a。究其原因,是研究區(qū)林地資源豐富、土壤水分比較充沛、氣候濕潤(rùn),因此林地蒸散量比較高;而荒地由于植被覆蓋度低,土壤保水能力差,在相同的水熱條件下,其蒸散量相對(duì)較低。

從圖10可以進(jìn)一步看出,不同土地利用類型15年來(lái)年均ET變化趨勢(shì)大致類似?；牡氐腅T值年際波動(dòng)最為顯著,林地次之,兩者均呈微弱增長(zhǎng)趨勢(shì),與全流域年均ET的變化趨勢(shì)相同,而草地和灌木的年際波動(dòng)則不明顯。

圖9 不同土地利用類型的年平均蒸散量Fig.9 Annual average evapotranspiration of different land use types

圖10 不同土地利用類型蒸散量年際變化Fig.10 Interannual variation of evapotranspiration of different land use types

4 結(jié)論與討論

本文利用MOD16產(chǎn)品研究了2000—2014年珠江流域蒸散量時(shí)空分布特征和變化趨勢(shì),結(jié)論如下:

(1)時(shí)間格局上,2000—2014年珠江流域ET年際波動(dòng)較大,總體呈微弱增加趨勢(shì),變化率為0.321 1 mm/a。多年ET在941.0～1 002.9 mm/a,均值為959.33 mm/a。ET的年內(nèi)分布呈現(xiàn)先增大后減小的單峰型分布,各季節(jié)ET均值的大小關(guān)系為:夏季>秋季>春季>冬季。

(2)空間格局上,珠江流域多年平均ET空間差異顯著,變化范圍在 401.2～1 665.1 mm/a,整體上呈現(xiàn)東高西低、南高北低的空間分布格局。ET峰值出現(xiàn)在北江流域,而低值區(qū)集中分布在珠江源地區(qū)。四季平均ET夏季最大,高值區(qū)成片出現(xiàn)在紅水河流域和右江流域;春秋接近且在流域內(nèi)均衡分布;冬季最小,高低值的分布格局較其他季節(jié)有較大差異且空間分布差異不明顯。

(3)由Theil-Sen 趨勢(shì)分析結(jié)果可知:ET_Sen趨勢(shì)度介于-57.3～44.2,ET變化趨勢(shì)以輕微增加為主。其中ET增加區(qū)域(54.13%)大于減少區(qū)域(43.55%), 顯著增加區(qū)域(7.15%)高于顯著減少的區(qū)域(4.68%),而ET基本穩(wěn)定的區(qū)域僅占2.32%。顯著增加的區(qū)域主要分布于紅水河流域、柳江流域、南盤江流域;而顯著減少的區(qū)域主要分布在南盤江流域、右江流域、北江流域。

(4)受地表類型的影響,不同土地利用類型ET的分布特征具有差異性。林地的ET均值最高,荒地最低。不同土地利用類型年均ET波動(dòng)基本類似,荒地的ET值年際波動(dòng)最為顯著。

2000—2014年珠江流域ET年際呈微弱增加趨勢(shì),年內(nèi)變化在春季增長(zhǎng)迅速,夏季達(dá)到峰值,秋季急劇減小,冬季變動(dòng)較平緩且全年最低。這與鐘昊哲等[24]采用最小二乘法對(duì)Penman-Monteith-Leuning模型進(jìn)行參數(shù)優(yōu)化,并結(jié)合MOD15A2葉面積指數(shù)進(jìn)行空間外推,得出西南喀斯特區(qū)域蒸散發(fā)與季節(jié)的明顯相關(guān)性:夏季蒸散發(fā)最高,冬季蒸散發(fā)最低的研究結(jié)果一致?；赥heil-Sen中值趨勢(shì)分析和Mann-Kendall顯著性檢驗(yàn)方法分析研究區(qū)ET變化趨勢(shì),ET明顯增加的區(qū)域主要分布于貴州境內(nèi),集中在安順、黔南和黔西南三者交界處。這與戴明宏等[25]運(yùn)用反距離權(quán)重插值法得出多年平均ET最高值點(diǎn)分布于安順地區(qū),畢節(jié)南部、六盤水東部和黔西南州北部次之的結(jié)論相似。由于動(dòng)力和熱力性質(zhì)差異,不同土地覆蓋的蒸散發(fā)量差異顯著,林地>草地>灌木>耕地>荒地,且多年月平均蒸散量最大值均出現(xiàn)在7月。邴龍飛等[26]對(duì) NOAH 陸面模式模擬的近30年中國(guó)陸地蒸散量和土壤含水量進(jìn)行統(tǒng)計(jì)分析,認(rèn)為林地在所有的生態(tài)系統(tǒng)類型中平均蒸散量最大,其次為草地、農(nóng)田,除荒漠外,各生態(tài)系統(tǒng)類型蒸散量皆在7月達(dá)到最大,與本研究的結(jié)論相似。因此,本研究的結(jié)論對(duì)認(rèn)識(shí)珠江流域不同時(shí)空維度的蒸散發(fā)格局增添了新的證據(jù)。

MOD16數(shù)據(jù)集基于Penman-Monteith方程,較好地反映了區(qū)域水和能量平衡,可用于揭示區(qū)域蒸散發(fā)的時(shí)空動(dòng)態(tài)特征。隨著區(qū)域氣候系統(tǒng)的變化和人類活動(dòng)的影響,珠江流域內(nèi)ET顯著減少的區(qū)域是否會(huì)持續(xù)發(fā)展,尤其是城市擴(kuò)張?jiān)斐傻淖匀恢脖黄茐?還需要全面、長(zhǎng)序列的數(shù)據(jù)進(jìn)行深入研究。區(qū)域地表蒸散發(fā)過(guò)程十分復(fù)雜,影響因素眾多,主要包括水量供應(yīng)因素(土壤含水量)和能量因素(主要有太陽(yáng)輻射、 氣溫、 水汽壓差以及風(fēng)速等),本文尚未進(jìn)一步定量分析這些因素對(duì)ET變化趨勢(shì)的影響,有待進(jìn)一步建立模型進(jìn)行分析。