滇中地區參考作物蒸散量時空變化特征

王 杰, 吳 灝, 黃 英, 王樹鵬

(云南省水利水電科學研究院, 昆明 650228)

?

滇中地區參考作物蒸散量時空變化特征

王 杰, 吳 灝, 黃 英, 王樹鵬

(云南省水利水電科學研究院, 昆明 650228)

參考作物蒸散量是灌溉設計、灌溉計劃等的基礎數據，利用滇中地區19個氣象臺站的觀測數據，計算了滇中地區的參考作物蒸散量(ET0)，分析了ET0時間和空間的變化特征及氣象要素對其的影響。結果表明：研究區的ET0于1982年發生突變，1960—1982年變化趨勢不明顯，1982—2002年呈現下降的趨勢，2003—2012年ET0呈現增加的趨勢，多年平均ET0約為1 223.7 mm。ET0的空間特征表現為中部高，東西低，春季最大，夏季高于秋季，冬季最小，高值區出現在元謀地區。ET0與風速、氣溫和日照時數呈現顯著的正相關關系，與相對濕度呈現極顯著的負相關關系。偏相關分析和逐步回歸分析顯示在年尺度上，風速、相對濕度和日照時數的組合可以預測ET0的年際變化。

參考作物蒸散量; 滇中地區; Penman-Monteith公式

在全球的降水量中，約有60%通過陸地生態系統蒸散返回到大氣中，被認為是水循環過程中最難的一個環節之一[1]，控制著陸地生態系統和大氣之間的物質和能量交換，成為水資源研究的重點領域之一。在蒸散量估算中，參考作物蒸散量是基礎參數，表征大氣蒸發能力的一個量度[2]，同時也是制定作物灌溉制度、灌區灌溉需水量的基礎，影響到區域水資源的優化配置[3]。對于農業用水、生態系統模型和徑流估算等具有重要的影響。國內外對于參考作物蒸散量進行計算的方法很多，其中彭曼—蒙特斯公式被認為是目前最好的方法，具有很強的理論基礎，所需數據能夠很容易得到，并得到廣泛的應用。

在氣候變暖大背景下，世界上大多數地區ET0呈下降趨勢[4]。在中國不同區域ET0空間差異十分明顯，總體呈下降趨勢[5-6]。云南水資源總量豐富，但時空分布不均，與人口、土地和生產力布局等不匹配。近年來隨著社會經濟的發展，農業用水量受到擠占，2009—2012年發生4 a連續干旱，使得農業缺水愈發嚴重。為此，本文利用滇中地區氣象臺站的日觀測資料，計算滇中地區的ET0，在此基礎上分析滇中地區ET0的時空分布規律及氣象要素對其的影響，以期為當地灌溉設計、灌溉計劃和水資源管理等提供理論依據當地和實踐指導。

1　研究區概況

滇中地區位于東經99°33′—103°40′，北緯24°23′—27°05′，包括昆明、玉溪、紅河、大理、曲靖、楚雄和麗江等，國土面積約為1.5×105km2。研究區屬于亞熱帶氣候，日照充足，四季如春，氣候宜人，干濕季分明，最熱月平均溫度約為19～22℃，最冷月平均溫度約為6～8℃。年降水量約為955.0 mm，其中85%以上的降雨量集中在5—10月，11月至次年4月的降水量僅占全年的15%左右。土壤類型以山原紅壤、棕壤和水稻土等為主。地形以山地和山間盆地為主，地勢起伏緩和。植被類型多樣，多為次生植被和人工植被。

2　數據及方法

2.1數據

本文采用滇中地區及臨近地區的19個氣象臺站，時間序列為1960—2012年，主要包括逐日最低氣溫和最高氣溫、相對濕度、風速和日照時數，各站點的經緯度和海拔等數據。

2.2方法

2.2.1ET0計算ET0的計算采用FAO Penman-Monteith公式[7]。通過matlab軟件編程，讀入計算所需要的日最低氣溫和最高氣溫、相對濕度、風速、日照時數、經緯度和海拔等數據，實現對日尺度上ET0的計算；通過使用ArcGIS 9.3軟件，利用kriging空間插值方法實現ET0的空間插值，生成年尺度和季節尺度上ET0的空間分布圖。

2.2.2Mann-Kendall突變檢驗本文采用Mann-Kendal檢驗對ET0的序列進行趨勢分析和突變分析。Mann-kendall檢驗不需要樣本遵從一定的分布，也不受少數異常值的干擾[8]，適用于分析實際中的非正態分布，在進行突變分析時還可以明確突變開始的時間，并指出突變區域[9]。

3　結果與分析

3.1ET0時間變化特征

由各氣象站1960—2012年ET0年均值，求平均值得到研究區近53 a來ET0的年均值。研究區多年平均ET0約為1 238.0 mm，最大值約為1 348.4 mm，出現在2012年；最小值約為1 115.9 mm，出現在1990年(圖1)。通過M-K檢驗得到統計量U的順序、逆序變化曲線UBk和UFk(圖2)，可知：取顯著水平0.05得到兩條臨界線y=±1.96，在臨界值之間UB和UF曲線有一個明顯的交點，位于1982年左右，即研究區平均ET0突變時間約為1982年，從1982年以后出現下降趨勢，從1990年以后出現顯著下降趨勢(圖2)。

從時間序列上可以看出1960—1982年ET0變化趨勢不明顯，MK檢驗同樣顯示不顯著的變化趨勢；1982—2002年ET0表現為下降趨勢，MK檢驗顯示1982—2002年表現出下降的趨勢；2003—2012年ET0表現為增加的趨勢，MK檢驗顯示相對于整個時間序列而言，其顯著性呈現減小的趨勢(圖1,2)。

圖1滇中地區參考蒸散量年際變化

圖2研究區平均ET0的Mann-Kendall檢驗

3.2ET0空間變化特征

利用Kriging空間插值方法將ET0空間柵格化，得到ET0在空間上的連續分布。滇中地區ET0在空間上表現出明顯的空間差異。多年平均ET0空間分布特征總體上表現為：中部高，東西低，大多數等值線呈現南北走向的半環分布，研究區多年平均ET0均超過了1 000 mm(圖3)。高值區出現在元謀地區，范圍為1 310.1～1 863.7 mm，這是由于元謀地區位于干熱河谷，相對于研究區的其他地區，氣溫高、日照時間長，接受的太陽輻射能多，加之降水稀少導致相對濕度較低，因而ET0較大。

春、夏、秋和冬的ET0空間分布與年的ET0比較，變化特征基本一致，在空間上同樣表現出明顯的空間差異，分布特征為中間高，東西低，高值區同樣出現在元謀地區。春季(3—5月)變化范圍334.6～474.3 mm，約占全年的34.0%；夏季(6—8月)變化范圍306.6～354.2 mm，約占全年的27.7%；秋季(9—11月)變化范圍225.9～263.5 mm，約占全年的20.3%；冬季(12月—次年2月)變化范圍180.7～245.5 mm，約占全年的18.0%；總體上表現為春季最大，夏季高于秋季，冬季最小(圖4)。

圖3滇中地區年ET0空間分布

圖4滇中地區ET0季節空間分布

3.3氣象要素對ET0的影響

ET0受當地氣候條件的影響，反映作物在水肥條件充足的條件的蒸散潛力，本文討論風速、氣溫、相對濕度和日照時數在年尺度上對ET0的影響。對研究區而言，風速和日照時數與ET0呈現極顯著的正相關關系決定系數分別為0.43，0.58，通過0.01的顯著性檢驗；氣溫與ET0呈現顯著的正相關關系，決定系數均為0.20，通過0.05的顯著性檢驗；而相對濕度與ET0呈現極顯著的負相關關系，決定系數均為0.58，通過0.01的顯著性檢驗(圖5)。通過ET0與氣象要素偏相關分析可知，與相對濕度和日照時數的偏相關系數最大，其次為風速，最小為氣溫(表1)。通過逐步回歸分析發現，在年尺度上風速、相對濕度和日照時數的組合可以預測ET0的年際變化(表2)。

圖5　參考作物蒸散量與氣象要素的相關關系

注：*通過0.01的顯著性檢驗，下同。

表2　ET0與氣象要素的逐步回歸模型

4　結論

(1) 研究區ET0平均值約為1 223.7 mm，通過M-K檢驗得知研究區平均ET0突變時間約為1982年，從1982年以后出現下降的趨勢，從1990年以后出現顯著下降的趨勢。

(2) 滇中地區ET0在空間上分布特征總體上表現為：中部高，東西低，高值區出現在元謀地區。在季節上，春季最大，夏季高于秋季，冬季最小。

(3) ET0與風速、氣溫和日照時數呈現顯著的正相關關系，與相對濕度呈現極顯著的負相關關系。在年尺度上風速、相對濕度和日照時數的組合可以預測ET0的年際變化。

[1]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.

[2]卓瑪蘭草,劉普幸,張亞寧,等.甘肅黃土高原區潛在蒸散量時空變化與成因研究[J].水土保持研究,2012,19(1):70-75.

[3]張麗,吳金亮,楊國范.基于BP神經網絡的東港灌區作物需水量預測研究[J].水土保持研究,2012,19(6):207-210.

[4]Brutsaert W, Parlange M B. Hydrologic cycle explains the evaporation paradox[J]. Nature, 1998,396(6706):30-30.

[5]高歌,陳德亮,任國玉,等.1956—2000年中國潛在蒸散量變化趨勢[J].地理研究,2006,25(3):378-387.

[6]Thomas A. Spatial and temporal characteristics of potential evapotranspiration trends over China[J]. International Journal of Climatology, 2000,20(4):381-396.

[7]Allen R G, Pereira L S, Raes D, et al. Crop Evapotranspiration Guidelines for Computing Crop Water Requirements [M]. FAO Irrigation and Drainage，1998.

[8]汪言在,茍詩薇,張述林.三峽庫區(重慶段)降雨侵蝕力變化趨勢及突變分析[J].水土保持研究,2012,19(3):1-6.

[9]于東平,張鑫,何毅,等.青海東部高原農業參考作物蒸散量的時空變化[J].農業工程學報2012,28(2):66-71.

Spatiotemporal Variations of Reference Crop Evapotranspiration in Central Yunnan

WANG Jie, WU Hao, HUANG Ying, WANG Shupeng

(Yunnan Hydraulic Research Institute, Kunming 650228, China)

The ET0is the basis of irrigation and irrigation scheme. This paper, based on the meteorological measurement of meteorological stations in central Yunnan, the yearly reference evapotranspiration (ET0) was estimated. The spatiotemporal characteristics and the influence of meteorological elements were also analyzed. The results showed that the abrupt change of ET0happened in 1982 approximately, the downward trend happened from 1982 to 2002, the increased trend happened from 2003 to 2012, the regional average ET0was 1 223.7 mm. the main space characteristics of ET0is high in central and low in east and west, the maximum value appeared in Yuanmou area in the seasons, the maximum value of ET0occurred in spring, and the value of ET0in summer was higher than that in autumn, the minimum value of ET0occurred in spring. The ET0increased with wind velocity, temperature and sunshine duration significantly and decreased with the relative humidity significantly. Partial correlative analysis and multiple regression analysis indicated that wind velocity, relative humidity and sunshine duration were the most significant predictors of the annual ET0.

reference crop evapotranspiration; central Yunnan; Penman-Monteith equation

2014-08-06

2014-10-15

云南省社會發展科技計劃項目(社會事業發展專項)(2012CA021);云南省水利廳水資源費項目(201203)

王杰(1979—),男,甘肅靜寧人,高工,博士,主要從事水文水資源研究。E-amil:wangjie@lzb.ac.cn

S161.4

A

1005-3409(2015)04-0041-04