趙國永, 韓 艷, 劉明華, 侯俊麗, 時海萍, 劉維真, 郭伊利, 喬 強
(1.信陽師范學院 地理科學學院, 河南 信陽 464000; 2.信陽師范學院 氣候與環境演變重點實驗室,河南 信陽 464000; 3.信陽師范學院 南北過渡帶研究中心, 河南 信陽 464000)
隨著全球變暖,極端氣候事件發生的頻率和強度急劇增加[1]。相較于氣候平均狀況,極端氣候事件帶來不利影響更加嚴重[2],引起各國政府和國際機構高度重視,已成為當前學術界關注的一個熱點問題[3]。極端降水事件是極端氣候事件重要組成部分,國內外學者進行了大量研究。與極端氣溫事件相比,極端降水事件不具備全球一致性[4-5],如一些地區總降水量減少或不變,但是極端降水的頻率和強度具有增加趨勢[6]。我國不同地區極端降水事件變化趨勢不一致,長江中下游地區、西北地區北部和西南地區西部年極端降水頻數具有明顯增加趨勢,而華北地區呈現減少趨勢[2]。同時,學者對我國[7]及不同區域極端氣候事件進行研究,如:西藏[8]、我國南方[9]、西南地區[10]、我國半干旱區[11]、黃土高原[12]、淮河流域[13]、瀏陽河流域[14]、山東省[15]等。
河南省位于黃河中下游,地勢低平,是我國人口大省和農業大省,是我國糧食生產的核心產區。已有學者對河南省極端降水事件進行研究,1957—2005年,河南省暴雨日數和極端降水事件具有不顯著增加趨勢[16]。1961—2006年,河南省汛期極端降水事件空間變化特征為:降水量南多北少,極端降水事件發生頻次呈現北多南少;汛期極端降水事件發生頻次存在南北差異[17]。相較于其他區域極端降水事件研究,河南省研究程度較低,且已有河南省極端降水事件研究指數較少、研究時間尺度較短。
本文擬以1961—2013年河南省17個站點逐日降水量數據為研究對象,運用Sen′s傾向估計、Mann-Kendall顯著性檢驗和空間插值等方法,分析河南省11個極端降水指數時間和空間變化特征及其影響因素。通過本文研究,有助于理解河南省極端降水事件變化特征,為決策部門提供參考,盡可能降低極端降水事件帶來的影響和損失。
本文研究數據由中國氣象科學數據共享服務網提供,包括1961—2013年河南省17個站點逐日降水量,并經過極值檢驗、時間一致性檢驗和均一化檢驗等嚴格質量控制。運用R編輯器、RClimDex軟件對數據進行進一步質量控制,包括數據記錄日期是否與現實一致、錯誤值及異常值的篩選等;且要求觀測站連續觀測時間段為1961年1月1日—2013年12月31日。
本研究對極端降水指數定義和計算采用世界氣象組織氣候委員會(WMO-CCL)和氣候變化與可預測性計劃(CLIVAR)氣候變化檢測、監測和指數專家組(ETCCDMI)確定的指數體系[1],從中選取11個極端降水指數(表1)。采用基于R語言RClimDex 1.0軟件(來源于http:∥cccma.seos.uvi.ca/ETCCDMI)進行極端降水指數計算。
運用Sen′s傾向估計分析11個極端降水指數變化趨勢[18]。運用Mann-Kendall方法進行變化趨勢顯著性檢驗[19]。在ArcGIS 10.0環境下,通過反距離加權(Inverse Distance Weighting)插值繪制河南省11個極端降水指數及變化趨勢空間分布圖。
1961—2013年,河南省5個極端降水指數(PRCPTOT,RX5day,R10,R20和CWD)具有減少趨勢,且均沒有通過α=0.05顯著性檢驗;6個指數(RX1day,R95p,R99p,SDⅡ,R25和CDD)具有升高趨勢,只有SDⅡ通過α=0.05顯著性檢驗。同時,11個指數變化幅度不明顯(表2),共同說明河南省極端降水指數變化趨勢不顯著。其PRCPTOT變化趨勢與韓艷等[20]基于月觀測數據研究結果一致。
2.2.1 空間分布 自河南省東南部向西北部,10個極端降水指數(PRCPTOT,RX1day,RX5day,R95p,R99p,SDⅡ,R10,R20,R25,CWD)逐漸減少,而CDD逐漸增加,說明河南東南部降水量、降水日數、降水強度較大;西北部降水量、降水日數、降水強度較小。其原因是河南省東南部受東亞季風影響強度大于西北部。
2.2.2 變化趨勢空間分布 近53 a,河南省17個站點中PRCPTOT趨于減少站點占76.4%,主要分布在河南南部和西北部,最大減幅為21 mm/10 a;其余站點趨于增加,主要分布在商丘和西華等地,最大增幅為21 mm/10 a(圖1A)。RX1day趨于減少站點占47.1%,主要分布在三門峽、孟津、新鄉、駐馬店、固始等地區;其余站點趨于增加(圖1B)。RX5day趨于減少站點占58.8%,主要分布在安陽、新鄉、孟津、寶豐、許昌、固始等地,最大減幅為8 mm/10 a;其余站點趨于增加,主要分布在西華、商丘等地區,最大增幅為12 mm/10 a。其中,西華變化趨勢通過α=0.05顯著性檢驗(圖1C),說明西華具有顯著增加趨勢。
R95p趨于減少站點占41.2%,主要分布在三門峽、孟津、鄭州等地,最大減幅為12 mm/10 a;其余站點趨于增加,主要分布在西華、南陽、許昌、商丘等地區(其中西華變化趨勢通過α=0.05顯著性檢驗),最大增幅為30 mm/10 a(圖1D)。R99p趨于減少站點占47.1%,主要分布在西峽、三門峽、孟津、安陽等地;其余站點趨于增加(圖1E)。SDⅡ趨于減少站點占29.4%,主要分布在河南省東南部和西北部,每10 a最大減幅為0.2 mm/d;其余站點趨于增加,主要分布在西華、商丘等地,每10 a最大增幅為0.6 mm/d。其中,西華和商丘變化趨勢通過α=0.05顯著性檢驗(圖1F),說明西華和商丘具有顯著增加趨勢。

表1 極端降水指數定義

表2 1961-2013年河南省極端降水指數變化趨勢及顯著性檢驗
注:*表示通過α=0.05顯著性檢驗。
R10趨于增加站點占17.6%,主要分布在盧氏、西華和商丘,最大增幅為0.2 d/10 a;其余站點均趨于減少,最大減幅為1 d/10 a(圖1G)。R20趨于減少站點占41.2%,主要分布在信陽、固始、新鄉等地區;其余站點趨于增加(圖1H)。R25趨于減少站點占41.2%,主要分布在信陽、新鄉、安陽等地;其余站點趨于增加(圖1I)。
CDD趨于減少站點占41.2%,主要分布在安陽、商丘等地;其余站點趨于增加(圖1J)。CWD趨于增加站點占11.8%,主要分布在鄭州和固始,最大增幅為0.1 d/10 a;其余站點均趨于減少,最大減幅為0.3 d/10 a(圖1K)。
過去53 a,11個極端降水指數中,商丘和西華9個指數(PRCPTOT,RX5day,RX1day,R95p,R99p,SDⅡ,R10,R20和R25)具有升高趨勢,且增加幅度為每一個指數最大值;安陽、新鄉、孟津和三門峽9個指數(PRCPTOT,RX5day,RX1day,R95p,R99p,R10,R25,CDD和CWD)均具有降低趨勢,且降低幅度幾乎為每一個指數最大值。在河南省區域內,商丘和西華降水量、降水日數和降水強度顯著增加;安陽、新鄉、孟津、三門峽降水量、降水日數和降水強度顯著減少。
與全球同期7個極端降水指數相比,河南省4個指數(PRCPTOT,RX1day,RX5day和R10)具有相反變化趨勢,3個指數(SDⅡ,CDD和CWD)具有相同變化趨勢(表3),說明河南省和全球極端降水指數變化趨勢存在差異。
與我國其他區域同期指數相比。與西藏已有9個極端降水指數[8]均具有相同變化趨勢,與我國南方已有5個指數中[9]4個指數具有相同變化趨勢,與我國半干旱區已有7個指數[11]中3個指數具有相同變化趨勢(表3),同時河南省降水強度變化幅度均大于我國南方、青藏高原和半干旱區,說明河南省極端降水指數變化特征與我國南方和青藏高原地區相似,而與我國半干旱區不同,體現了極端降水指數區域性差異。
河南省11個極端降水指數與緯度相關性,8個指數相關系數達到0.8以上,同時10個指數通過了顯著性檢驗(9個α=0.01,1個α=0.05顯著性檢驗)(表4),說明河南省極端降水指數與緯度密切相關。同時,10個極端降水指數與緯度呈負相關,1個指數(CDD)與緯度呈正相關(表4),說明隨著緯度增加(自南向北),河南省降水量、降水日數、降水強度逐漸減少,而連續干旱日數逐漸增加。

注:黑色方框表示變化趨勢通過α=0.05顯著性檢驗。
圖11961-2013年河南省極端降水指數變化趨勢空間變化分布
河南省11個極端降水指數與經度和海拔相關性,大部分指數相關系數均較高,說明河南省極端降水指數與經度和海拔具有較好相關性。極端降水指數與經度和海拔相關系數中,只有3個指數(RX1day,RX5day,SDⅡ)通過顯著性檢驗(表4),說明河南省極端降水指數與經度和海拔密切程度低于緯度。同時,9個極端降水指數與經度呈正相關,2個指數(CDD,CWD)呈負相關且相關系數接近于0;9個極端降水指數與海拔呈負相關,2個指數(R10,CWD)呈正相關,說明隨著經度遞減(自東向西)和海拔升高(自東向西),河南省極端降水指數均具有減少趨勢。河南省降水量自東部向西部逐漸減少,與月觀測數據記錄結果一致。其原因是河南省處于東亞季風區范圍內,東部更加靠近海洋,東部受季風影響強于西部。
從表5可以看出,河南省9個極端降水指數與年總降水量呈正相關,1個指數(CDD)呈負相關。正相關系數在0.59以上,最高達到0.95,且均通過α=0.01顯著性檢驗,說明本文選取11個極端降水指數對河南省降水具有較好的指示作用,極端降水指數的增多/減少反映了降水量的增多和減少[12]。同時,河南省每一個站點4個指數(R95p,R10,R20和R25)與年總降水量正相關系數均通過α=0.05顯著性檢驗,說明極端降水量(R95p)、年降水日數降水量(R10,R20和R25)對年總降水量貢獻最大。

表3 河南省極端降水指數變化趨勢與其他區域對比
注:*,**分別表示通過α=0.05和α=0.01顯著性檢驗,青藏高原數據、全球數據來自于參考文獻[8],我國南方數據來自于文獻[9],我國北方半干旱區數據來自于文獻[11]。

表4 1961-2013年河南省極端降水指數與緯度、經度和海拔相關系數
注:*,**分別表示通過α=0.05和α=0.01顯著性檢驗。

表5 1961-2013年河南省極端降水指數與年總降水量相關系數
注:**表示通過α=0.01顯著性檢驗,SNo表示通過α=0.05顯著性檢驗站點數占總站點數比例。
(1) 時間尺度上,河南省極端降水指數變化趨勢不顯著。(2) 空間尺度上,河南東南部降水量、降水日數、降水強度強于其西北部。過去53 a,商丘和西華降水量、降水日數、降水強度顯著增加,而安陽、新鄉、孟津、三門峽顯著減少。(3) 緯度因素對河南極端降水指數變化影響強于經度和海拔因素,極端降水量(R95p)、年降水日數降水量(R10,R20和R25)對年總降水量貢獻最大。