楊秀芹 孫恒 王燕

摘要 基于2000—2014年MOD16遙感蒸散發(fā)數(shù)據(jù)，利用經(jīng)驗(yàn)正交函數(shù)（EOF），分析淮河流域年蒸散發(fā)時空分布特征。結(jié)果表明，淮河流域年蒸散發(fā)空間分布可由4個模態(tài)特征場來表現(xiàn)，其空間分布結(jié)構(gòu)為全區(qū)一致型、東南至西北差異型、東南至西北相反型、南北相反型。各模態(tài)2000—2010年蒸散發(fā)均呈上升趨勢，之后蒸散發(fā)呈下降趨勢。

關(guān)鍵詞 地表蒸散發(fā)；經(jīng)驗(yàn)正交函數(shù)法；時空分布；淮河流域

中圖分類號 S161.4 文獻(xiàn)標(biāo)識碼 A 文章編號 0517-6611（2017）02-0197-03

Abstract Annually spatialtemporal characteristics of land evapotranspiration （ET） over Huaihe River basin were analyzed based on MOD16 ET product during 2000-2014 using Empirical Orthogonal Function method （EOF）.Research results showed that annually spatial distribution of evapotranspiration over Huaihe River basin can be expressed by four model characteristics， including spatial distribution structure in the same type， the southeast and the northwest opposite type， the southeast and the northwest opposite type， the north and the south opposite type.Major modes presented that the land evapotranspiration had upward trend in the first decade of the 21st century， and then decreased.

Key words Land evapotranspiration；Empirical orthogonal function method；Spatialtemporal distribution；Huaihe River basin

在以變暖為主要特征的全球變化背景下，溫度和降水等要素變化會影響地表水分收支平衡，從而改變地表干濕狀況。降落到地球表面的降水有60%是通過蒸發(fā)或蒸散的作用返回大氣[1]，在部分干旱地區(qū)可高達(dá)90%[2]，蒸散發(fā)過程同時也消耗了大約3/5的地表凈輻射能量[3]。陸面蒸散發(fā)可以影響降水，同時伴隨的潛熱具有降溫作用，因此蒸散發(fā)作為地表能量平衡的重要組成部分和陸面水分平衡的重要組成部分，對地表干濕狀況起到了不可忽視的作用[4]?；春恿饔蚴菄抑匾纳唐芳Z棉油基地，地表干濕狀況及其變化對農(nóng)業(yè)生產(chǎn)、區(qū)域生態(tài)環(huán)境和國家經(jīng)濟(jì)有重要影響。因此，研究淮河流域蒸散發(fā)時空格局對該地區(qū)乃至全國都有重要意義。

不同遙感蒸散發(fā)模型的假設(shè)和主要控制因子的不同會造成陸面蒸散發(fā)估計的差異，經(jīng)過反演得到時空格局合理的數(shù)據(jù)集是研究蒸散發(fā)的前提。MOD16產(chǎn)品不僅提供了高時空分辨率的蒸散發(fā)數(shù)據(jù)集，且具有可以免費(fèi)在線獲取的特點(diǎn)[5]。目前已有很多學(xué)者使用MOD16產(chǎn)品對我國或某一區(qū)域地表蒸散量時空分布特征進(jìn)行分析[6-9]，但對淮河流域蒸散發(fā)時空分布特征的研究較少。筆者利用時間尺度為月的MOD16產(chǎn)品2000—2014年蒸散發(fā)數(shù)據(jù)，采用經(jīng)驗(yàn)正交函數(shù)（EOF），探討淮河流域年蒸散發(fā)時空變化特征。

1 資料與方法

1.1 研究區(qū)概況

淮河流域（112.125°～121.375° E、31.125°～36.375° N）跨河南、湖北、安徽、江蘇、山東共5個省，西高東低，在黃河與長江之間，水系分布密集。東北部及東部為平原，西部有桐柏山、伏牛山，南部有大別山、江淮丘陵，臨近長江，北部緊鄰泰山和黃河。淮河流域處于秦嶺—淮河一線南北分界線處，屬大陸性季風(fēng)氣候，全年降水充分，流域年徑流量大，土壤濕度大，年平均氣溫11～16 ℃。

1.2 研究方法及數(shù)據(jù)來源

1.2.1 MOD16產(chǎn)品及數(shù)據(jù)。MOD16產(chǎn)品數(shù)據(jù)集包括全球植被覆蓋區(qū)域的8 d、月、年時間尺度的地表蒸散發(fā)（ET），空間分辨率為1 km×1 km和0.05 °×0.05 °。該研究主要使用2000—2014年空間分辨率為0.05°×0.05°的Tiff格式的MOD16A2月尺度產(chǎn)品，通過該產(chǎn)品可獲得高時空分辨率的地表蒸散發(fā)，MOD16產(chǎn)品的蒸散發(fā)算法詳見文獻(xiàn)[10-11]。

1.2.2 EOF分析方法。

EOF（Empirical Orthogonal Function）方法，即經(jīng)驗(yàn)正交函數(shù)展開法，通過使用EOF分析方法，可以提取源數(shù)據(jù)矩陣的特征向量和時間系數(shù)。將其特征向量稱為空間模態(tài)，通過分析主要的空間模態(tài)和其相對應(yīng)的時間系數(shù)，對原始場做出估計和解釋[12]?？杀硎緸椋?/p>

X=ZV（1）

式中，X、Z、V分別稱為原始場矩陣、空間函數(shù)矩陣、時間函數(shù)矩陣[13]。由（1）式可知，若空間函數(shù)矩陣場為正時，時間系數(shù)為正，原始場矩陣場為正，表示該地區(qū)年蒸散發(fā)有增大的趨勢；時間系數(shù)為負(fù)，原始場矩陣場為負(fù)，表示該地區(qū)蒸散發(fā)有減小的趨勢。同理可得空間函數(shù)矩陣場為負(fù)時的情況。該研究采用距平化后2000—2014年逐年蒸散發(fā)數(shù)據(jù)，進(jìn)行EOF分析，得到2000—2014年淮河流域年蒸散發(fā)時空特征。

2 結(jié)果與分析

2.1 蒸散發(fā)空間變化特征

由表1可見，前4個模態(tài)的累積方差貢獻(xiàn)率為82.7%，因此利用前4個模態(tài)的時間序列及其特征場可較好地表現(xiàn)2000—2014年淮河流域年蒸散發(fā)的時空變化特征。

從圖1可看出，淮河流域第1特征向量基本為負(fù)值，表明2000—2014年淮河流域年蒸散發(fā)在空間上具有很好的一致性，定義為淮河流域年蒸散發(fā)分布全區(qū)一致型。年蒸散發(fā)第2模態(tài)特征向量，在淮河流域東南部出現(xiàn)正值中心，西北部出現(xiàn)負(fù)值中心，說明淮河流域東南部和西北部呈反位相關(guān)系，但東南部正值區(qū)較弱，定義為淮河流域年蒸散發(fā)分布東南至西北差異型。年蒸散發(fā)第3模態(tài)特征向量，在淮河流域西北部出現(xiàn)負(fù)值中心，東南部出現(xiàn)正值中心，說明淮河流域東南部與西北部呈反位相關(guān)系，東南部正值區(qū)較強(qiáng)且范圍大，定義為淮河流域年蒸散發(fā)分布東南至西北相反型。年蒸散發(fā)第4模態(tài)特征向量，在淮河流域北部出現(xiàn)正值中心，南部出現(xiàn)負(fù)值中心，說明淮河流域北部和南部呈反位相關(guān)系，定義為淮河流域年蒸散發(fā)分布南北相反型。

2.2 蒸散發(fā)時間變化特征

年第1模態(tài)中，年蒸散發(fā)同一變化趨勢持續(xù)時間長，年代際區(qū)分明顯。由圖2a可知，淮河流域地區(qū)2000、2001、2009、2011—2014年年蒸散發(fā)有增加趨勢，2002—2008、2010年年蒸散發(fā)有減少趨勢。2001—2010年淮河流域地區(qū)表現(xiàn)為由干轉(zhuǎn)濕的過程，2002—2008年年蒸散發(fā)持續(xù)下降，在2011年后，淮河流域地區(qū)表現(xiàn)為由濕轉(zhuǎn)干，蒸散發(fā)持續(xù)上升。

第2模態(tài)中，年蒸散發(fā)增加、減少趨勢呈周期變化。由圖2b可知，淮河流域西北地區(qū)2000—2003、2006、2009—2012年年蒸散發(fā)有減少趨勢，2004、2005、2007、2008、2013、2014年年蒸散發(fā)有增加趨勢。

2000—2014年淮河流域年蒸散發(fā)增加趨勢和減少趨勢的年份交替出現(xiàn)，年蒸散發(fā)呈現(xiàn)出周期性轉(zhuǎn)變，蒸散發(fā)減少周期為3—4年，蒸散發(fā)增加周期為1—2年。

第3模態(tài)中，年蒸散發(fā)增加、減少趨勢持續(xù)時間長，年代際區(qū)分明顯。由圖2c可知，2002、2007—2010、2012、2013年淮河流域西北部地區(qū)年蒸散發(fā)有減少趨勢，2000、2001、2003—2006、2011、2014年年蒸散發(fā)有增加趨勢。2000—2010年淮河流域西北地區(qū)表現(xiàn)為由干轉(zhuǎn)濕的過程，2003—2006年年蒸散發(fā)持續(xù)上升，2006年到達(dá)上升節(jié)點(diǎn)，2007—2010年年蒸散發(fā)持續(xù)下降。

第4模態(tài)中，淮河流域南部地區(qū)2000、2001、2006—2008、2011、2012年年蒸散發(fā)有減少趨勢，2002—2005、2009、2010、2013、2014年年蒸散發(fā)有增加趨勢（圖2d）。2000—2010年淮河流域南部地區(qū)表現(xiàn)為由干轉(zhuǎn)濕的過程，2000—2005年年蒸散發(fā)呈上升趨勢，2006—2010年年蒸散發(fā)持續(xù)下降；2011—2014年淮河流域南部地區(qū)表現(xiàn)為干濕周期性轉(zhuǎn)換，濕潤年持續(xù)時間為1年，干旱年持續(xù)時間為1年。

綜上所述，年蒸散發(fā)第1和第3模態(tài)中，2000—2010年蒸散發(fā)逐漸減少，淮河流域地區(qū)表現(xiàn)為由干轉(zhuǎn)濕的過程，在2010年左右，蒸散發(fā)開始逐漸增加，淮河流域地區(qū)表現(xiàn)為由濕轉(zhuǎn)干的過程。第2和第4模態(tài)中，年蒸散發(fā)的增加趨勢和減少趨勢的年份交替出現(xiàn)，并呈現(xiàn)出一定的周期性，即淮河流域地區(qū)干濕變化呈現(xiàn)出周期性，表明淮河流域干濕變化受到某特定規(guī)律的影響。

3 結(jié)論

（1）2000—2014年淮河流域年蒸散發(fā)空間場主要表現(xiàn)為4種形式，分別是全區(qū)一致型、東南至西北差異型、東南至西北相反型、南北相反型。

（2）通過經(jīng)驗(yàn)正交函數(shù)分析可得，2000—2010年淮河流域地區(qū)表現(xiàn)為由干轉(zhuǎn)濕的過程，2002—2008年年蒸散發(fā)持續(xù)下降，2010年后，淮河流域地區(qū)表現(xiàn)為由濕轉(zhuǎn)干。

參考文獻(xiàn)

[1] WANG K C，DICKINSON R E.A review of global terrestrial evapotranspiration：Observation，modeling， climatology，and climatic variability[J].Reviews of geophysics，2012，50（2）：1-54.

[2] 達(dá)日瑪.基于MODIS數(shù)據(jù)的區(qū)域蒸散發(fā)量遙感估算及其應(yīng)用：以巴彥淖爾市為例[D].呼和浩特：內(nèi)蒙古師范大學(xué)，2010.

[3] TRENBERTH K E，F(xiàn)ASULLO J T，KIEHI J.Earth′s global energy budget[J].Bulletin of the american meteorological society，2009，90（3）：311-323.

[4] 郭曉寅，程國棟.遙感技術(shù)應(yīng)用于地表面蒸散發(fā)的研究進(jìn)展[J].地球科學(xué)進(jìn)展，2004，19（1）：107-114.

[5] 喻元，白建軍，王建博，等.基于MOD16的關(guān)中地區(qū)實(shí)際蒸散發(fā)時空特征分析[J].干旱地區(qū)農(nóng)業(yè)研究，2015，33（3）：245-253.

[6] 范建忠，李登科，高茂盛.基于MOD16的陜西省蒸散量時空分布特征[J].生態(tài)環(huán)境學(xué)報，2014，23（9）：1536-1543.

[7] 馮飛，姚云軍，張彥彬，等.基于MOD16產(chǎn)品的三江平原蒸散量時空分布特征分析[J].生態(tài)環(huán)境學(xué)報，2015，24（11）：1858-1864.

[8] 吳桂平，劉元波，趙曉松，等.基于MOD16產(chǎn)品的鄱陽湖流域地表蒸散量時空分布特征[J].地理研究，2013，32（4）：617-627.

[9] 何慧娟，卓靜，李紅梅，等.基于MOD16產(chǎn)品的陜西關(guān)中地區(qū)干旱時空分布特征[J].干旱地區(qū)農(nóng)業(yè)研究，2016，34（1）：236-241.

[10] MU Q，HEINSCH F A，ZHAO M S，et al.Development of a global evapotranspiration algorithm based on MODIS and global meteorology data[J].Remote sensing of environment，2007，111（4）：519-536.

[11] MU Q Z，ZHAO M S，RUNNING S W.Improvements to a MODIS global terrestrial evapotranspiration algorithm[J].Remote sensing of environment，2011，115（8）：1781-1800.

[12] 黎偉標(biāo).廣西汛期降水的EOF分析[J].廣西氣象，1991，12（4）：12 -16.

[13] 鄧愛軍，陶詩言，陳烈庭.我國汛期降水的EOF分析[J].大氣科學(xué)，1989，13（3）：289-295.