三江平原多年降水演變分析

王洪麗

(哈爾濱市雙城區水務局，黑龍江 哈爾濱 150100)

1 研究區域概況

三江平原位于黑龍江省東部，北起黑龍江、南臨興凱湖、西鄰小興安嶺、東至烏蘇里江，總面積約10.89萬km2，總人口862.5萬人，人口密度約為79人/km2。屬溫帶濕潤、半濕潤大陸性季風氣候，降水季主要集中于夏秋季，光照時數適中，雨熱同期，適合一年一熟的水稻、大豆，尤其是優質水稻和高油大豆等農作物的生長。三江平原有著非常豐富的自然資源儲備，是我國重要的商品糧生產基地，被中國國家地理雜志評選為中國十大“新天府”之一。

2 資料與數據來源

早在二十世紀八九十年代，已經有學者開始關注極端氣候變化所帶來的嚴重后果[1]，近年來關于極限降水演變特征的研究也越來越多。地球表面平均溫度在過去的百年間以0.74 ℃/10 a的速度升高[2]，表明全球氣溫增長持續加快。在未來全球持續變暖情況下，極端高溫事件將會增加，具體來說，熱浪持續時間更長，頻率更高；干旱的地區干旱加劇，濕潤的地區降水增加[3]，導致發生干旱和洪澇災害的頻率、強度總體上均呈上升趨勢[4]。

本文所采用的氣象數據資料來源于國家氣象科學數據中心中國地面氣候資料日值數據集(V3.0)。研究區域內選取了4個時間序列長、數據記錄完整、空間代表性強的氣象站(佳木斯站、富錦站、寶清站、和虎林站)1959—2018年的逐日降水資料。

所選各站基本資料見表1。

表1 氣象站概況

3 極限降水指標

表2 極限降水指標

4 結果與分析

4.1 一日最大降水量

基于該地區4個測站60年間的逐日降水資料可以計算出每年的一日最大降水量，畫出其隨時間變化的折線圖(如圖1所示)，進行線性擬合后可以看出變化趨勢：四個測站中佳木斯站的回歸系數為正值，說明佳木斯站60年間一日最大降水量呈上升趨勢；富錦站、寶清站、虎林站的回歸系數均為負值且絕對值均大于佳木斯站，說明佳木斯站的上升幅度遠遠小于其他3個站的下降幅度，富錦站的擬合直線相較來說更陡。因此，綜上4個測站結果來看，該地區一日最大降水量呈現波動下降趨勢，其中富錦站下降趨勢最為明顯，下降速度為1.60 mm/10 a。

圖1 不同站年一日最大降水量隨時間變化趨勢

計算出一日最大降水量之后進行了時間序列的趨勢分析和突變檢驗，運用的是M-K非參數檢驗方法，其結果見表3。

表3 不同站一日最大降水量傾向率和M-K顯著性檢驗 mm·10a-1

M-K趨勢檢驗表明，三江平原地區除佳木斯站之外的三個氣象站一日最大降水量均呈下降趨勢且佳木斯站的降水傾向率的絕對值并不算大，因此總體來看一日最大降水量呈波動下降趨勢，這與線性傾向估計法檢驗結果相同。突變檢驗以富錦站為主，見圖2。

圖2 富錦站一日最大降水量M-K突變檢驗

結果表明，UF和UB兩條曲線在1959年前后相交，在1987年前后出現明顯的突變點，表示富錦站一日最大降水量下降趨勢從1959年開始顯現，在1987年出現突變點。

4.2 強降水量

本文將一年中日降水量從小到大進行排序，取大于等于95%分位值的降水量之和作為這一年的強降水量。

以富錦站1959年為例，計算其強降水量，見表4。

按照以上方法分別計算出4個測站60年間的強降水量，畫出其隨時間變化的折線圖，并進行線性擬合可以看出變化趨勢，得到圖3。分析得出強降水量隨時間變化的趨勢：四個測站中只有佳木斯站的回歸系數為正值，說明該站強降水量逐年增加，但上升趨勢不明顯；富錦站、寶清站及虎林站的回歸系數均為負值且絕對值遠大于佳木斯站，寶清站的擬合直線最陡。因此，綜上4個測站結果來看，三江平原地區強降水量總體上呈下降趨勢，其中寶清站下降趨勢最為明顯，下降速度為8.55 mm/10 a。

表4 富錦站1959年強降水量

圖3 不同站強降水量隨時間變化趨勢

計算出強降水量之后進行了時間序列的趨勢分析和突變檢驗，運用的是M-K非參數檢驗方法，其結果見表5和圖4。M-K趨勢檢驗表明，三江平原地區4個氣象站的強降水量均呈現下降趨勢，佳木斯站的降水傾向率最小即下降趨勢最微弱。這與線性傾向估計法檢驗結果基本相同。結果表明，寶清站強降水量的下降趨勢從1962年開始顯現，在1987年出現突變點。

熊老師您好！我在一家專修歐系大眾車的修理廠工作，請您指導一下，如何利用五氣分析儀來確診故障？我特別弄不明白的是，混合汽稀時，為什么氮氧化合物反而增加？請您不吝賜教！

表5 不同站強降水量傾向率和M-K顯著性檢驗 mm·10a-1

4.3 最大一次降水量

計算出60 a的年最大一次降水量之后畫出其隨時間變化的折線圖，進行線性擬合后發現4個測站的回歸系數均為負值，可以明顯地看出變化趨勢，寶清站的擬合直線下降得最快。由圖5可以看出三江平原地區4個測站近60年來的最大一次降水量均呈波動下降趨勢。其中寶清站下降最為明顯，以6.83 mm/10 a的速度下降。

圖4 寶清站強降水量M-K突變檢驗

圖5 不同站年最大一次降水量隨時間變化趨勢

利用1959—2018年三江平原地區4個測站的降水資料，對單站年最大一次降水量指標進行時間序列的M-K趨勢分析和M-K突變檢驗，其結果見表6。

表6 不同站最大一次降水量傾向率和M-K顯著性檢驗 mm·10a-1

趨勢檢驗表明，4個站的降水傾向率均為負值，說明三江平原地區4個測站的年最大一次降水量均呈現下降趨勢，這與線性傾向估計法分析結果相同。突變檢驗結果以佳木斯站和寶清站為主，見圖6。

圖6 不同站年最大一次降水量M-K突變檢驗

佳木斯站最大一次降水量下降趨勢從1961年開始顯現,在1978年出現突變點，寶清站最大一次降水量的下降趨勢從1965年開始顯現,在1966年出現突變點。

4.4 連續干旱日數

連續干旱日數表征干旱程度，用來判斷一個區域是否容易受到旱災侵襲，本文將日降水量小于1 mm的最長連續天數定義為連續干旱日數。

得到連續干旱日數隨時間變化的折線圖后進行線性擬合，可以看出其隨時間的變化趨勢，由圖7可以看出：4個測站的回歸系數均為負值，說明該地區的連續干旱日數在4個站均呈現波動下降趨勢，其中虎林站的回歸系數絕對值最大即下降趨勢最明顯，擬合直線最陡峭，以2.24 d/10 a的速度下降。可見，三江平原地區連續干旱日數整體以緩慢的速度下降。

圖7 不同站連續干旱日數隨時間變化趨勢

分別對1959—2018年三江平原地區4個測站的單站連續干旱日數進行時間序列的M-K趨勢分析和突變檢驗，結果見表7。

趨勢檢驗表明，三江平原地區4個氣象站的降水傾向率均為負值，說明連續干旱日數均呈現下降趨勢，此結果與上述線性傾向趨勢分析結果是一致的。突變檢驗結果以寶清站和虎林站為主，見圖8。

圖8 不同站連續干旱日數M-K突變檢驗

寶清站連續干旱日數的減少趨勢從1996年開始顯現,在2006年出現突變點，虎林站連續干旱日數的減少趨勢從2003年開始顯現,在2015年出現突變點。

4.5 暴雨日數

三江平原地區暴雨多集中在7—8月，本文把日降水量大于50 mm的天數定義為暴雨日數。

畫出4個測站的暴雨日數隨時間變化的折線圖，并進行線性擬合，可看出該地區暴雨日數的變化趨勢，直線趨于平穩，見圖9，三江平原地區佳木斯站的暴雨日數逐年增長，以0.06 d/10 a的速度上升；寶清站下降趨勢較為緩慢，以0.12 d/10 a的速度下降。其余兩個測站的回歸系數非常接近0，說明暴雨日數無明顯變化趨勢。綜上4個測站的結果來看，該地區60年間暴雨日數并無明顯變化。

基于該地區的日降水資料可以得出4個測站每一年的暴雨日數，之后運用M-K非參數檢驗方法進行了時間序列的趨勢分析和突變檢驗，其結果見表8。趨勢檢驗表明，三江平原地區四個測站的暴雨日數均無明顯變化趨勢，說明三江平原地區近60年暴雨日數變化逐漸趨于平穩，這與線性傾向估計的結果大致相同。突變檢驗表明，60年中，三江平原地區暴雨日數尚未發生突變。

5 結 論

運用三江平原4個氣象站1959—2018年的逐日降水資料，計算一日最大降水量等五個極限降水指標，從而對該地區近60年極限降水事件的演變規律進行研究。得出一日最大降水量在除佳木斯站之外的3個測站均隨時間下降，表明該地區的一日最大降水量在西部有下降趨勢，在東部有上升趨勢；強降水量在4個測站均呈下降趨勢即該地區強降水量逐漸減少；最大一次降水量在4個測站均呈下降趨勢即該地區最大一次降水量逐漸減少；連續干旱日數在該地區呈緩慢下降趨勢；暴雨日數無明顯變化趨勢，整體來看變化趨于平穩，沒有發生突變。并且該地區極限降水具有一定的季節差異性。

圖9 不同站暴雨日數隨時間變化趨勢