關穎紅
(遼寧省東港市水利技術推廣中心 遼寧 東港 118300)
蒸發(fā)蒸騰是植物生長的重要生理活動之一,同時也是地球水文循環(huán)的重要組成部分,直接關系著地表能量與水量平衡[1],落在地面的所有降水中,約有62%被蒸騰消耗掉了[2]。因此,它的精確計算對水資源的有效利用及水文研究有著重要意義,尤其對作物需水量的計算、灌溉制度的確定、干旱管理等有著顯著影響。近幾十年來,通過參考作物騰發(fā)量(ET0)計算實際作物騰發(fā)量進而確定作物需水量的方法在國內外受到廣泛關注,不少學者對其進行了深入研究并形成了諸多理論。其中,penman-montieth(P-M)公式被FAO推薦為計算 ET0的唯一標準公式[2~3]。
丹東市位于遼寧省東南部,地理坐標東經124°23′,北緯40°07′。屬暖溫帶亞濕潤季風型氣候,年平均雨量多在800mm~1200mm,占遼寧省降水量的70%,是我國北方雨量最多的地區(qū)。年平均氣溫南部在8℃~9℃,北部6℃~7℃,受季風影響,季節(jié)變化明顯,四季分明,是東北地區(qū)最溫暖最濕潤的地方,素有“北方江南”、“東北蘇杭”之稱。
本文利用丹東市4個氣象站1960年~2006年的逐日氣象資料,采用P-M公式計算各站逐日ET0,并對各站計算結果求平均得到丹東地區(qū)逐日ET0,再通過分類匯總得到整個地區(qū)逐月及逐年ET0,在此基礎上分析研究區(qū)不同時間尺度下ET0的演變特征,并對引起該變化的主要氣象影響因素進行了探討,以期為研究丹東地區(qū)未來ET0的變化和用水管理提供依據(jù)。逐日ET0的計算公式如下[4]:
(1)ET0式中:ET0是參考作物騰發(fā)量,Rn是到達作物表面的凈輻射[MJ/(m2·d)],G是土壤熱通量 [MJ/(m2·d)],T是作物冠層2m高處的空氣溫度(℃),U2是作物冠層2m高處的風速 (m/s),ea是飽和水汽壓(kPa),ed是實際水汽壓(kPa),Δ是水汽壓對溫度的斜率(kPa/℃),γ 是干濕球常數(shù)。
采用M-K秩次相關法對丹東市ET0及各氣象因素的變化趨勢進行顯著性檢驗,其有關理論參見文獻[5];對于ET0變化的影響因素分析,采用基于偏導數(shù)的敏感性分析法。其計算公式如下[3,6]:
(2)式中:SXi為參數(shù) Xi的敏感系數(shù)(無量綱)。SXi取正值表明ET0隨著Xi的增加而增加,反之則表明ET0隨著Xi的增加而減?。粅SXi|越大表明Xi變化對ET0的影響越大,反之亦然。
氣象要素及ET0的變化趨勢用一次線性方程表示,即[7]
(3)式中:xi為氣象要素的擬合值,10a1為氣候傾向率,表示氣象要素每10年的變化率。氣候傾向率描述的是序列的升降程度,其小于0說明研究對象在統(tǒng)計時段內呈下降趨勢,反之呈增長趨勢[8]。
圖1 丹東地區(qū)年均日ET0年際變化及其累計距平曲線
圖2 各站點及全區(qū)平均ET0年際變化曲線
圖3 丹東地區(qū)ET0月際變化情況
表1 各站點及全區(qū)逐月日ET0變化的氣候傾向率 /mm·(10y)
表2 各站點及全區(qū)平均各氣象因子變化的氣候傾向率(1960年~2006年) /mm·(10y)
表3 各站點及全區(qū)年均日ET0與各氣象因子的相關系數(shù)及敏感系數(shù)
丹東地區(qū)年均日ET0年際變化及其累計距平曲線如圖1。全區(qū)ET0總體上緩慢減少,但趨勢不明顯(P=0.783),最大值出現(xiàn)在 1965 年,其值為 2.414mm/(y·d),最小值出現(xiàn)在2005年,其值為 1.921mm/(y·d)。結合累計距平圖(圖1b)可知,研究時段內丹東地區(qū)ET0變化大致可分為四個階段:第一階段和第三階段分別對應于1960年~1968年、1996年~2003年,該階段累計距平值持續(xù)上升,表明這兩個階段的ET0比平均ET0大,屬于騰發(fā)較強的階段;第二階段和第四階段分別對應于1968年~1996年、2003年~2006年,該階段累計距平值持續(xù)下降,說明這個階段內ET0比平均ET0小,屬于騰發(fā)較弱的階段;其中第二階段內,ET0的累計距平曲線出現(xiàn)一定幅度的波動,說明在這個時段內,丹東地區(qū)ET0不穩(wěn)定。
由圖2可知,丹東地區(qū)各氣象站ET0與全區(qū)平均ET0年際變化趨勢相似,其中丹東站與全區(qū)平均值最接近,可近似代表丹東地區(qū)近幾十年來ET0年際變化情況,而與全區(qū)平均值差異最大的是寬甸站,最大差值為0.612mm/(y·d)。
丹東地區(qū)ET0月際變化比較大,各站點ET0年內分布情況相似,最高值一般發(fā)生在5月份~6月份,最低值出現(xiàn)在12月份~1月份。由圖3可知,全區(qū)平均連續(xù)最大6個月ET0值出現(xiàn)在4月份~9月份,約占全年比例的73.2%,其中5月份~6月份約占27.5%;而12月份~1月份僅占全年比例的5.0%。由圖也可以看出,一年四季中,丹東地區(qū)冬季ET0最低,夏季最高,這可能與冬季溫度低、日照短而夏季溫度高、日照長有關。
分析各站點及全區(qū)ET0的氣候傾向率(如表1)可知,丹東地區(qū)ET0年內各月變化趨勢不盡相同,其中12月份~3月份及8月份~9月份ET0有上升趨勢,4月份~7月份及10月份~11月份ET0有下降趨勢,除5月份外其余各月變化趨勢均未達到顯著水平。各站點及全區(qū)2月份ET0均呈上升趨勢,5月份ET0均呈下降趨勢,其余各月規(guī)律沒有一致性。其中,氣候傾向率值最大的是東港站5月份的ET0,其值為-0.114mm/(10y),最小的是丹東站8月份的 ET0,其值為 0.00002mm/(10y)。
選擇平均溫度、相對濕度、風速和太陽輻射4個氣象因素進行分析。由表2、表3可知,絕大部分站點及全區(qū)平均太陽輻射和風速均有顯著的下降趨勢,除寬甸站外,相對濕度也呈下降趨勢,但大多數(shù)趨勢欠顯著,各站點及全區(qū)平均溫度均呈顯著上升趨勢,這與全球氣候變暖的客觀事實相符。丹東地區(qū)年均ET0與各氣象因素的相關性由強到弱依次為:太陽輻射>相對濕度>平均溫度>風速,其中,太陽輻射、平均溫度及風速與ET0呈正相關,相對濕度與ET0呈負相關,這與常理相符。另外,通過敏感性分析得,丹東地區(qū)年ET0對太陽輻射的變化最敏感,其次是相對濕度的變化和平均溫度的變化,最后是風速的變化,這與相關分析結果一致,由此可以說明,丹東地區(qū)ET0變化受太陽輻射的影響最大,但還不能確定其余3個氣象因子對ET0變化影響強弱的先后順序,需要從各氣象因子擾動對ET0變化的貢獻值角度做進一步研究。
用同樣的分析方法分析影響月ET0值變化的主要氣象因素。結果表明,不同月份對丹東地區(qū)ET0變化影響最大的氣象因子不盡相同,12月份~3月份相對濕度對ET0影響較大,4月份~7月份、10月份~11月份太陽輻射對ET0影響較大??傮w來看,溫度和風速對各月ET0影響不顯著。
研究表明,丹東地區(qū)多年參考作物騰發(fā)量呈不明顯下降趨勢,太陽輻射是影響丹東地區(qū)多年參考作物騰發(fā)量變化的最主要因素,這一結論印證了太陽輻射是蒸發(fā)蒸騰最基本的能量來源。
月際分析表明,研究區(qū)參考作物騰發(fā)量最高值發(fā)生在5月份~6月份,最低值出現(xiàn)在12月份~1月份,各站點變化情況具有相似性。從季節(jié)上看,除冬季ET0有上升趨勢外,春、夏、秋3季ET0均呈緩慢減少趨勢,但這些趨勢均未達到顯著水平。相對濕度對冬季ET0影響較大,其余各季太陽輻射對ET0影響最顯著。陜西水利
[1]劉小莽,鄭紅星,劉昌明,等.海河流域潛在蒸散發(fā)的氣候敏感性分析 [J].資源科學,2009,31(9):1470-1476.
[2]V.Th.AMBAS,E.BALAS.Sensitivity Analysis of Different Evapotranspiration Methods Using A New Sensitivity Coefficient[J].Global Nest Journal,2012,14(3):335-343.
[3]曾麗紅,宋開山,張柏,等.東北地區(qū)參考作物蒸散量對主要氣象要素的敏感性分析[J].中國農業(yè)氣象,2010,31(1):11-18.
[4]張淑杰,張玉書,隋東,等.東北地區(qū)參考蒸散量的變化特征及其成因分析 [J].自然資源學報,2010,25(10):1750-1760.
[5]K.H.HAMED.Exact distribution of the Mann Kendall trend test statistic for persistent data[J].JournalofHydrology,2009,365:86-94.[6]GONGLB,CHONGYX,CHEN DL,etal.Sensitivity of the Penman-Monteith reference evapotranspirationtokeyclimatevariablesinthe Changjiang(YangtzeRiver)basin[J].Journalof Hydrology,2006,32(9):620-629.
[7]孫蘭東,劉德祥.西北地區(qū)熱量資源對氣候變化的響應特征[J].干旱氣象,2008,26(1):8-12.
[8]李春友,張勁松,孟平,等.黑龍港流域參考作物蒸散量的時序變化和分形特征[J].中國農業(yè)氣象,2011,32(1):81-88.