申官正　王龍　余航　張素凡　高瑞

摘要：參考作物蒸散量是表征氣候干旱程度以及水資源供需平衡的重要指標(biāo)。本研究利用橫斷山區(qū)31個(gè)氣象站1960—2013年逐月氣象資料，采用聯(lián)合國糧農(nóng)組織推薦的彭曼公式計(jì)算研究區(qū)各站點(diǎn)逐月參考作物蒸散量，建立31個(gè)站點(diǎn)蒸散量序列，利用氣候傾向率、累積距平法、多元逐步回歸分析以及基于GIS中的反距離插值技術(shù)等方法，分析橫斷山區(qū)蒸散量的空間分布及時(shí)間變化趨勢。結(jié)果表明：橫斷山區(qū)參考作物蒸散量由南向北遞減。受氣溫上升、日照時(shí)數(shù)減少、風(fēng)速及相對(duì)濕度降低的影響，研究區(qū)蒸散量變化呈上升趨勢。1960—1968年研究區(qū)參考作物蒸散量逐步降低，1968—1984年研究區(qū)參考作物蒸散量逐步增加，1984—2004年研究區(qū)參考作物蒸散量逐步降低，2004—2013年研究區(qū)參考作物蒸散量逐步增加。參考作物蒸散量變化趨勢區(qū)域差異性明顯。就貢獻(xiàn)率而言，研究區(qū)北部參考作物蒸散量減少與平均風(fēng)速相關(guān)性大；研究區(qū)南部參考作物蒸散量增加與平均日照時(shí)數(shù)相關(guān)性大；研究區(qū)中部參考作物蒸散量增加與平均氣溫相關(guān)性大。

關(guān)鍵詞：橫斷山區(qū)；參考作物蒸散量；貢獻(xiàn)率；傾向率；時(shí)空變化

中圖分類號(hào)：S161.4 文獻(xiàn)標(biāo)識(shí)號(hào)：A 文章編號(hào)：1001-4942（2017）04-0104-08

Temporal and Spatial Variation Characteristics of Reference Crop

Evapotranspiration in the Hengduan Mountains

Shen Guanzheng， Wang Long， Yu Hang， Zhang Sufan， Gao Rui

（Department of Water Conservancy， Yunnan Agricultural University， Kunming 650201， China）

Abstract Reference evapotranspiration （ETo） is an important index for evaluating climate drought， vegetation water-consumption. Based on the meteorological data from 31 meteorological stations during 1960-2013 in the Hengduan Mountainous， the reference crop evapotranspiration was calculated with the Penman formula recommended by FAO. And then， the evapotranspiration sequence of 31 stations was set up to analyze the temporal and spacial changing trend of evapotranspiration in the Hengduan Mountainous by using linear regression method， cumulative anomaly method， multiple stepwise regression analysis and inverse distance interpolation technique based on GIS. The results showed that the reference crop evapotranspiration in the Hengduan Mountainous decreased from south to north. Influenced by temperature rising and decrease of sunshine duration， wind speed and relative humidity， the evapotranspiration of this area increased. The reference crop evapotranspiration gradually reduced from 1960 to 1968， gradually increased from 1968 to 1984， gradually reduced from 1984 to 2004， and gradually increased from 2004 to 2013. The spatial variation tendency of reference crop evapotranspiration was obvious. As for the contribution rate analysis， the decrease of reference crop evapotranspiration in the north research area was largely correlated with the average wind speed； the increase of reference crop evapotranspiration in the south of research area was largely correlated with average sunshine time； and the increase of reference crop evapotranspiration had great relevance with average temperature in the middle of research area.

Keywords The Hengduan Mountains； Reference crop evapotranspiration； Contribution rate； Tendency rate； Temporal and spatial variation

隨著對(duì)全球各地差異性氣候的研究，溫度遞增已成為全球范圍內(nèi)氣候變化的主誘因[1]。而氣候的變化很大程度上決定了水資源的時(shí)空分布。隨著人們對(duì)環(huán)境重視程度的不斷提升，氣候的變化也愈發(fā)得以重視[2，3]。

參考作物蒸散量又稱最大可能蒸散量，是表征氣候干旱程度以及水資源供需平衡的重要指標(biāo)，對(duì)水資源利用與規(guī)劃以及節(jié)水農(nóng)業(yè)的推廣有著深遠(yuǎn)的指導(dǎo)意義。聯(lián)合國糧農(nóng)組織（FAO）于1998年就參考作物蒸散量作出解釋，假設(shè)作物高度為0.12 m，并有固定的表面阻力為70 s/m，反照率為0.23的參考冠層的蒸散量，相當(dāng)于高度一致、生長旺盛、完全覆蓋地面而不缺水的開闊草地的蒸散量[4，5]。

近百年來，在全球變暖的大背景下，國內(nèi)外學(xué)者利用聯(lián)合國糧農(nóng)組織（FAO）推薦的公式分析世界各地參考作物蒸散量的變化規(guī)律。研究表明，全球范圍內(nèi)，大部分地區(qū)的參考作物蒸散量呈下降趨勢。就中國而言，參考作物蒸散量呈減小趨勢，但各地差異較大，區(qū)域性明顯[6-11]。

橫斷山區(qū)幅員遼闊，其范圍包括青藏高原東南緣、川滇藏境內(nèi)山川東西并列、南北縱橫的廣大地區(qū)。位于我國第一地形階梯與第二地形階梯交界處，是氣候變化反應(yīng)較為敏感的區(qū)域之一[12]。近年來，許多學(xué)者針對(duì)橫斷山區(qū)溫度及降水的變化趨勢及規(guī)律進(jìn)行了很多研究，但對(duì)于參考作物蒸散量的研究還較少[13]。本研究利用橫斷山區(qū)31個(gè)地面氣象觀測站點(diǎn)1960―2013年的逐月氣象資料，在采用公式計(jì)算出近54年橫斷山區(qū)各站點(diǎn)逐月參考作物蒸散量的基礎(chǔ)上，利用多元逐步回歸分析、氣候傾向率、累積距平法以及地理信息系統(tǒng)（GIS）等技術(shù)方法對(duì)橫斷山區(qū)參考作物蒸散量進(jìn)行分析，以期為保護(hù)橫斷山區(qū)脆弱生態(tài)環(huán)境及水資源開發(fā)利用提供科學(xué)依據(jù)，進(jìn)而促進(jìn)橫斷山區(qū)農(nóng)業(yè)及生態(tài)環(huán)境的可持續(xù)發(fā)展。

1 數(shù)據(jù)來源與研究方法

1.1 數(shù)據(jù)來源

利用橫斷山區(qū)31個(gè)地面氣象觀測站點(diǎn)1960―2013年的逐月平均氣溫、平均最高氣溫、平均最低氣溫、日照時(shí)數(shù)、平均相對(duì)濕度、平均風(fēng)速等氣象資料及各地面氣象觀測站點(diǎn)的經(jīng)緯度及海拔等地理信息，對(duì)橫斷山區(qū)參考作物蒸散量的計(jì)算及時(shí)空分布進(jìn)行研究?？紤]山脈及河流走向等因素的影響，確定橫斷山區(qū)范圍及所選站點(diǎn)分布情況如圖1。

1.2 研究方法

1.2.1 參考作物蒸散量的計(jì)算 利用聯(lián)合國糧農(nóng)組織推薦的彭曼公式[4，5]進(jìn)行計(jì)算，見式（1）。

1.2.2 反距離加權(quán)插值法 反距離加權(quán)插值法多應(yīng)用于山區(qū)或地面氣象觀測站點(diǎn)較少的區(qū)域，可以提升預(yù)測值的精確度，通過調(diào)整權(quán)重改變等值線分布狀況。

1.2.3 參考作物蒸散量變化趨勢分析方法 采用累積距平法、氣候傾向率、多元逐步回歸和貢獻(xiàn)率的方法分析1960―2013年橫斷山區(qū)參考作物蒸散量的變化趨勢[14，15]。

2 結(jié)果與分析

2.1 參考作物蒸散量的空間分布特征

通過對(duì)橫斷山區(qū)31個(gè)氣象站點(diǎn)的氣象要素統(tǒng)計(jì)結(jié)果匯總分析得（圖2），橫斷山區(qū)1960—2013年參考作物蒸散量均值在727～1 275 mm范圍內(nèi)波動(dòng)，呈階梯狀分布。就區(qū)域變化而言，橫斷山區(qū)年蒸散量呈現(xiàn)出由北向南帶狀遞增的趨勢，低值區(qū)位于研究區(qū)東北部若爾蓋、松潘、都江堰一帶，最低值出現(xiàn)在若爾蓋（660 mm，1962年）；高值區(qū)位于研究區(qū)南部元江一帶，最高值出現(xiàn)在元江（1 633 mm，1980年）。參考作物蒸散量均值隨緯度的降低而增大，但受經(jīng)度變化的影響較小。

相較于年平均氣溫、年平均風(fēng)速、年平均日照時(shí)數(shù)及年平均相對(duì)濕度的空間分布可知，氣溫較高、日照充足、風(fēng)速較大且空氣濕度較小的區(qū)域，其對(duì)應(yīng)的參考作物蒸散量較大。

利用多元逐步回歸分析的方法[15]對(duì)橫斷山區(qū)31個(gè)站點(diǎn)年參考作物蒸散量與各研究站點(diǎn)的年平均氣溫、年平均日照時(shí)數(shù)、年平均相對(duì)濕度和年平均風(fēng)速等氣象要素進(jìn)行統(tǒng)計(jì)分析，可以得到：

照時(shí)數(shù)，h；V為年平均風(fēng)速，m/s。上式通過α=0.05的顯著性檢驗(yàn)，其具體意義表現(xiàn)為，在橫斷山區(qū)范圍內(nèi)，其它氣候條件保持不變，年平均氣溫每升高1℃，年參考作物蒸散量將增加15.39 mm；年平均相對(duì)濕度每增大1%，10年參考作物蒸散量將降低8.49 mm；年平均日照時(shí)數(shù)每增加1 h，年參考作物蒸散量將增加38.24 mm；年平均風(fēng)速每增加1 m/s，年參考作物蒸散量將增加43.29 mm。

2.2 參考作物蒸散量時(shí)間分布趨勢及成因分析

2.2.1 年參考作物蒸散量變化趨勢 1960―2013年橫斷山區(qū)參考作物蒸散量平均值以4.5 mm/10a的傾向率呈顯著上升趨勢（圖3）。 橫斷山區(qū)參考作物蒸散量距平是對(duì)其1960―2013年54年均值的偏差值，從累積距平曲線（圖4）可得，1968、1984和2004年為參考作物蒸散量的轉(zhuǎn)折點(diǎn)。1960―1968年是參考作物蒸散量偏少的階段，以負(fù)距平為主，呈現(xiàn)較小的減少趨勢，其傾向率為-1.74 mm/10a；1968―1984年是參考作物蒸散量偏多的階段，以正距平為主，呈現(xiàn)微弱的增長趨勢，其傾向率為0.28 mm/10a；1984―2004年為參考作物蒸散量偏少的階段，以負(fù)距平為主，呈現(xiàn)微弱的增長趨勢，其傾向率為0.41 mm/10a；2004―2013年為參考作物蒸散量偏多的階段，以正距平為主，呈現(xiàn)較大的增長趨勢，其傾向率為35 mm/10a。近54年來，最小距平值（-43.88 mm）出現(xiàn)在2000年，最大距平值（52.27 mm）出現(xiàn)在2013年。

2.2.2 參考作物蒸散量變化的氣候成因 為確定近54年影響橫斷山區(qū)參考作物蒸散量變化的氣候成因，本研究利用1960―2013年橫斷山區(qū)31個(gè)氣象觀測點(diǎn)54年的氣溫、風(fēng)速、日照時(shí)數(shù)及相對(duì)濕度等實(shí)測氣象數(shù)據(jù)，統(tǒng)計(jì)分析各因子與參考作物蒸散量的相關(guān)關(guān)系。結(jié)果表明，橫斷山區(qū)參考作物蒸散量與年平均風(fēng)速相關(guān)不顯著，與年平均氣溫和年平均日照時(shí)數(shù)呈顯著正相關(guān)（P<0.05），與年平均相對(duì)濕度呈顯著負(fù)相關(guān)（見表1）。

統(tǒng)計(jì)分析1960―2013年各主要?dú)庀笠蜃拥淖兓贸觯▓D6），近54年來，橫斷山區(qū)年平均氣溫以0.239℃/10a的傾向率呈顯著增加趨勢（P<0.05）；年平均風(fēng)速以0.038 m/（s·10a）的傾向率呈顯著減少趨勢；年平均相對(duì)濕度以0.462%/10a的傾向率呈顯著減少趨勢；年平均日照時(shí)數(shù)以0.034 h/10a的傾向率呈顯著減少趨勢。結(jié)合參考作物蒸散量與各氣象要素的相關(guān)關(guān)系可知，年平均氣溫的上升趨勢與橫斷山區(qū)參考作物蒸散量上升趨勢保持一致，表明年平均氣溫是橫斷山區(qū)參考作物蒸散量變化的重要誘因。

2.2.3 各氣象因子貢獻(xiàn)率分析 統(tǒng)計(jì)分析各氣象因子對(duì)參考作物蒸散量變化的貢獻(xiàn)率得到（圖7），平均氣溫對(duì)參考作物蒸散量的貢獻(xiàn)率高值分布在研究區(qū)東北部及中部，結(jié)合橫斷山區(qū)參考作物蒸散量傾向率的變化得到影響研究區(qū)中部蒸散量增加的主要因素為平均氣溫；平均日照時(shí)數(shù)對(duì)參考作物蒸散量的貢獻(xiàn)率高值分布在研究區(qū)的南部及東北部的都江堰一帶，結(jié)合橫斷山區(qū)參考作物蒸散量傾向率的變化得到影響研究區(qū)南部蒸散量增加的主要因素為平均日照時(shí)數(shù)；平均相對(duì)濕度對(duì)參考作物蒸散量貢獻(xiàn)率高值分布在研究區(qū)中部及東北部若爾蓋一帶；平均風(fēng)速對(duì)參考作物蒸散量貢獻(xiàn)率高值主要分布在研究區(qū)中北部及南部元江一帶，結(jié)合橫斷山區(qū)參考作物蒸散量傾向率的變化得到影響研究區(qū)中北部蒸散量減少的主要因素是平均風(fēng)速。

3 結(jié)論

（1）參考作物蒸散量是一個(gè)衡量氣候變化的重要指標(biāo)，可以很清晰地顯示出氣候變化的進(jìn)程和類別。橫斷山區(qū)年平均蒸散量分布在727～1 275 mm范圍內(nèi)，從南到北隨緯度增大而減小。

就整個(gè)橫斷山區(qū)年際變化而言，1960―1968年屬于逐步降低階段，1968―1984年屬于穩(wěn)步上升階段，1984―2004年屬于穩(wěn)步降低階段，2004―2013年屬于穩(wěn)步上升階段。但是各個(gè)站點(diǎn)差異性明顯，對(duì)各個(gè)站點(diǎn)還有進(jìn)一步研究的必要。

（2）橫斷山區(qū)31個(gè)站點(diǎn)的平均傾向率為4.5 mm/10a，表明研究區(qū)內(nèi)蒸散量呈遞增趨勢。蒸散量呈現(xiàn)增加趨勢的區(qū)域位于研究區(qū)的中部以及南部，蒸散量呈現(xiàn)減少趨勢的區(qū)域集中于研究區(qū)北部。

（3）通過對(duì)研究區(qū)氣象因子與參考作物蒸散量相關(guān)性分析，參考作物蒸散量與平均相對(duì)濕度呈顯著負(fù)相關(guān)，與平均氣溫、平均風(fēng)速、平均日照時(shí)數(shù)呈正相關(guān)，其中與平均氣溫和平均日照時(shí)數(shù)的相關(guān)性達(dá)顯著水平（P<0.05）。研究區(qū)中北部參考作物蒸散量減少主要受平均風(fēng)速影響，研究區(qū)南部參考作物蒸散量的增加主要受平均日照時(shí)數(shù)影響，研究區(qū)中部參考作物蒸散量增加主要受平均氣溫影響。

參 考 文 獻(xiàn)：

[1] IPCC. Summary for policymakers of climate change 2007： the physical science basis[C]//Contribution of Working Group Ⅰ to the fourth assessment report of the intergovernmental panel on climate change. Cambridge： Cambridge University Press，2007.

[2] 施雅風(fēng)，沈永平，胡汝驥. 西北氣候由暖干向暖濕轉(zhuǎn)型的信號(hào)、影響和前景初步探討[J]. 冰川凍土，2002，24（3）：219-226.

[3] 張山清，普宗朝，伏曉慧，等. 氣候變化對(duì)新疆自然植被凈第一性生產(chǎn)力的影響[J].干旱區(qū)研究，2010，27（6）：905-914.

[4] Smith M，Allen R G，Monteith J L，et al. Report on the expert consultation on revision of FAO methodologies for crop water requirements [R]. Rome： FAO，1992.

[5] Allen R G，Luis S P，Durk R，et al. FAO Irrigation and Drainage Paper No.56—Crop evapotranspiration： guidelines for computing crop water requirements [R]. Food and Agriculture Organization of the United Nations， Rome， 1998：300.

[6] 牛紀(jì)蘋，粟曉玲. 石羊河流域參考作物蒸發(fā)蒸騰量對(duì)氣候變化的響應(yīng)模擬及預(yù)測[J]. 水利學(xué)報(bào)，2014，45（3）：286-295.

[7] 趙璐，梁川，崔寧博，等. 川中丘陵區(qū)參考作物蒸發(fā)蒸騰量近60年變化成因研究[J]. 水利學(xué)報(bào)，2013，44（2）：183-190.

[8] 陳超，龐艷梅，潘學(xué)標(biāo)，等. 四川地區(qū)參考作物蒸散量的變化特征及氣候影響因素分析[J]. 中國農(nóng)業(yè)氣象，2011，32（1）：35-41.

[9] 杜加強(qiáng)，舒儉民，劉成程，等. 黃河上游參考作物蒸散量變化特征及其對(duì)氣候變化的響應(yīng)[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2012，28（12）：92-100.

[10]劉定輝，劉永紅，熊洪，等. 西南地區(qū)農(nóng)業(yè)重大氣象災(zāi)害危害及監(jiān)測防控研究[J]. 中國農(nóng)業(yè)氣象，2011，32（增1）：208-212.

[11]詹存，梁川，趙璐. 川中丘陵區(qū)季節(jié)性干旱時(shí)空分布特征及成因分析[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2013，29（21）：82-90.

[12]中國科學(xué)院青藏高原綜合科學(xué)考察隊(duì). 橫斷山區(qū)自然地理[M]. 北京：科學(xué)出版社，1997.

[13]李宗省，何元慶，辛惠娟，等. 我國橫斷山區(qū)1960―2008年氣溫和降水時(shí)空變化特征[J]. 地理學(xué)報(bào)，2010，65（5）：563-579.

[14]謝云峰，張樹文. 基于數(shù)字高程模型的復(fù)雜地形下的黑龍江平均氣溫空間插值[J].中國農(nóng)業(yè)氣象，2007，28（2）：205-211.

[15]普宗朝，張山清.近48年新疆夏半年參考作物蒸散量時(shí)空變化[J].中國農(nóng)業(yè)氣象，2011，32（1）：67-72.