張小南
(長江水利委員會漢江水文水資源勘測局,湖北 襄陽 441000)
受人類活動及自然環境影響,氣候變化導致極端天氣事件發生幾率大幅增加[1]。并且氣候具有明顯的區域性特征,各區域間存在明顯的差異[2]。國內外中眾多學者基于不同空間尺度與時間尺度做出大量研究,Tabari[3]等對伊朗1966 年~2005 年的年降水量變化趨勢展開研究,并發現局部地區的年降水量呈增加趨勢。Gocic M[4]運用Mann-Kendall 趨勢檢驗對塞爾維亞局部地區多組氣象數據資料的年際和季節趨勢進行空間特征分析。朱鵬[5]利用Mann-Kendall 檢驗驗證了北疆地區降雨量存在明顯上升趨勢,并依賴小波分析對不同時間尺度降雨資料的分析發現其周期習性變化規律。張洪亮[6]采用了滑動平均的方法分析了年降雨的變化趨勢,張文對年降雨量時間變化趨勢研究中則利用一元線性回歸進行分析描述。
國內外學者在對局部降雨變化趨勢分析的研究中不僅改進了方法,并且對年降雨時序研究內容的擴充。類如,王鵬云[7]利用模比差積曲線和滑動平均曲線將降雨的時間序列劃分為枯水期和豐水期,并以此確定局部降雨周期。
漢江流域是南水北調中線工程的水源地,流域降水強度大、驟發性強、時空分布不均。降雨情況對于水資源有重大影響,充分了解降雨極值的變化能夠幫助我們進行流域水資源再分配,同時對于防汛抗旱有著重大意義。因此,就流域尺度而言,研究區域內降雨量時間變化趨勢,對預防或減輕其災害損失,具有十分重要的現實意義。
漢江是長江最大的支流,流經陜西、湖北、河南等六個省市,其中干流流經湖北、陜西兩省,全長1577 km。
中上游多以丘陵坡地為主,屬于亞熱帶季風氣候,氣候溫和,植被良好,雨量充沛。降水量年內分布不均,主要分布于5 月~10 月,年際變化較大,多年平均降雨量為800 mm~1200 mm。漢江中上游洪水期徑流多由降水補給,平水期徑流多由地下水補給。歷史上漢江流域水患頻發,洪水多由局部暴雨和持續連陰雨造成,典型洪水如1935 年、1954 年、1964 年、1983 年、1998 年和2005 年等。
本文選取漢江中上游流域,研究采用漢江中上游向家坪、長沙壩、白河、黃龍灘四個水文站的降水系列進行分析,這四個站能獲得的降水系列較長,能較好地代表漢江中上游降水情況,漢江流域各測站示意圖見圖1。

圖1 研究站點分布圖
Mann-Kendall 檢驗法受異常值的干擾較小,且不需要序列服從一定的分布,因此廣泛用于分析水文及氣象要素序列的變化趨勢。降水量是具有時間序列的自然樣本,Mann-Kendall檢驗定義一個時間序列Yi(i=1,2,3,…,n),Ci是樣本Yi大于Yj的累積數:

原序列為無趨勢自然序列,則此序列的數學期望為:

方差為:

Xk的標準化為:

Xk逆序排列后,再次標準化,得UB。本次研究中,由于降水時空尺度變異較大,其變化趨勢的顯著性相對較弱,故選取較低的趨勢顯著性檢驗標準。本文將顯著水平設為0.05,對應該顯著水平的Zα/2為1.96,對應的置信水平為95%。即若UF、UB 曲線在置信區間內,則此序列未發生顯著變化;若UF、UB在置信區間內有交點,則此交點可能為突變點;UF>0,則時間序列是上升趨勢,若大于1.96,此序列顯著上升;UF<0,則時間序列是下降趨勢,若小于-1.96,此序列顯著下降。
對漢江流域內1956 年~2018 年的汛期降雨量數據進行線性回歸以及相關統計分析,結果見圖2。


圖2 1956年~2018 年漢江中上游流域年降水量變化趨勢
圖中結果表明,漢江流域中上游地區除向家坪地區降雨量整體呈現上升趨勢,且上升趨勢較為平緩。長沙壩站以下流域降雨量整體呈現下降趨勢,且下降趨勢較為平緩。其中越靠近下游下降趨勢越為顯著。以向家坪流域為例,區域內多年平均降水量為795.0 mm,最大降水量為1982 年1154.3 mm,最小降水量為1966 年470.4 mm,豐枯極值比為2.45,降水量年際變化較大,表1 詳列了各測點的降雨量枯豐極值比。

表1 漢江流域中上游各水文站降雨量年際變化
對比圖1,不僅可以發現漢江流域降雨量年際變化較大,局部表現為規律性的漲落態勢,在長沙壩與黃龍灘地區,年際降雨量變化明顯高于上游向家坪。各測站均出現了部分年份降雨量連續增加或者連續下降的趨勢,并且越靠近下游,趨勢越為顯著。年際降雨趨勢中不僅存在連續的降雨變化趨勢,更存在“陡增”和“陡落”現象。表明漢江流域降雨年際變化具有顯著的隨機性。圖3 為利用Mann-Kendal 檢驗法對漢江流域上游各測站降雨量數據的檢驗結果。

圖3 漢江中上游流域年降水量MK 曲線
如圖3 所示,UF、UB 兩條趨勢線相互交織均有不止一個交點,降雨量變化趨勢發生過多次轉折。對比各測站,可見向家坪站與黃龍灘站交點相對稀少且明顯,其中向家坪站于1990 年(35 年)、2003 年(48 年)、2007 年(52 年)、2011 年(56 年)、2016 年(61 年)相交;黃龍灘站于1962 年、1963 年、1970 年、1971 年、1973 年、1985 年、1986 年、2011 年、2018 年出現明顯相交。而位于兩測站中間的白河與長沙壩站,除部分年份存在明顯相交,其中長沙壩站UF、UB 曲線在1956 年~1965 年、1973 年~1984 年頻繁相交,降雨量在這兩個階段呈現波動式變化,頻繁突變上升或突變下降。白河站于1956 年~1980 年交點密集出現,相比之下,具有更為連續的變化序列,降雨量呈現震蕩變化。
從UF 曲線值看,1956 年~1989 年區間內向家坪UF 曲線值均大于0,區域年降水量持續增加,但在1984 年UF 值超過臨界值1.96,表明此點上升趨勢顯著。
而對圖3(b)分析可得,在1966 年~1973 年,UF 曲線值均小于0,表明長沙壩站降雨量經歷五六十年代震蕩之后,進入持續減少期,但是減少不顯著。1973 年~1984 年降水量減少趨勢變緩,1983 年~1991 年UF 值大于0,這階段長沙壩降雨量呈上升趨勢,趨勢不顯著。1992 年~2018 年,UF 值絕大多數都小于0,長沙壩降雨量呈下降趨勢,其中2002 年、2005 年、2006年均為突變年份。
圖3(c)中白河站1982 年~2018 年UF 值均小于0,并且在1999 年~2003 年UF 曲線值顯著低于臨界值-1.96,故白河區域降水量呈明顯下降趨勢。
與上游測站對比,黃龍灘站降水量在1970 年~2018 年處于不斷減少的趨勢。而1968 年前,UF 曲線值在0 點上下波動變化,區域內降雨量呈現出不同的變化規律。
漢江中上游流域及其分區的年降水量受空間分布影響顯著,其上游向家坪區域呈現上升趨勢,而靠近下游的各測站則表現出下降趨勢,其顯著性水平東部高于西部。而對于降雨突變時段,上下游基本呈現一直性,但明顯下游突變時段更為連續與集中,呈現明顯的流域性、空間群發性。