近40年黑龍江省降水量變化及趨勢特征分析

顏庭琦

摘要：為探究水文時間序列變化時空變化特征及趨勢，本文以黑龍江省各地區為例，基于線性回歸方法及Mann-Kendall突變檢驗法揭示1975—2021年

（564個月）降水變化特征及趨勢，結果表明：黑龍江省近40年內總體降水呈現較大波動，且趨勢性變化不明顯；大興安嶺地區擬合優度最優，統計檢驗值最接近-1.96，表明大興安嶺降水變化趨勢相較于其他地區更顯著；黑龍江省降水序列整體呈現下降趨勢，但降水曲線擬合效果差，整體降水變化無顯著規律，趨勢性不明顯，降水趨勢波動較大，降水趨勢預測較困難。黑龍江省降水變化特征及趨勢的研究為水資源管理及合理利用提供參考，并為農業經濟發展、城鎮供需水平衡、洪峰預報等提供現實依據。

關鍵詞：時間序列；降水；線性回歸方法；Mann-Kendall檢驗法

隨著全球氣候變暖趨勢的提升，降水作為社會-生態系統中重要的氣象要素，常表現為非線性、非確定性、突變性等特征[1]，并在長期演化過程中，受人為、自然等因素影響增多，造成降水預測越發困難[2]，并直接影響著農業管理、城鎮供需水矛盾等[3]實際性問題的解決。現階段，降水時空分布及變化趨勢特征的研究作為基礎且重要的工作，成為重點關注的科學問題。葉柏生等[4]采用線性回歸方法，分析1951—1998年全國范圍內的678個氣象站的降水變化趨勢；李淼等[5]采用Morlet小波函數分析北京地區近300年降水量時間序列，并揭示了降水變化的周期性變化規律；王若男等[6]采用累積距平法分析丹東地區降水量變化特征及演變趨勢。其中，Mann-Kendall檢驗法作為降水趨勢檢測的重要方法，樣本選擇可無序分布，并不受少部分異常信息點干擾[7]，可有效揭示降水序列時間變化趨勢。胡萌等[8]采用Mann-Kendall檢驗法分析青島市1956～2016年降水量變化趨勢，結果表明，降水量整體呈現下降趨勢，并在1986—2022年呈現顯著下降態勢；蘇萬磊[9]基于Mann-Kendall檢驗法對秦皇島周邊3個氣象站點近59年的年降水變化趨勢進行研究并確定突變檢驗，結果表明，秦皇島降水四季變化呈現下降趨勢，其中降水突變發生在1980年，并在1990—2010年呈現相對枯水狀態。綜上，以黑龍江省為例，采用線性回歸方法及Mann-Kendall檢驗法對降水時空分布及變化趨勢進行研究，為區域水資源管理及合理利用提供參考，并為農業經濟發展、城鎮供需水平衡、洪峰預報等提供現實依據。

1 材料與方法

1.1 研究區域

黑龍江省位于我國最東北部，以烏蘇里江為界毗鄰俄羅斯。全省共有13個市級行政區，包含哈爾濱、大興安嶺、黑河、齊齊哈爾、大慶、雙鴨山、鶴崗、雞西、牡丹江、七臺河、綏化、佳木斯、伊春。黑龍江省屬溫帶季風性氣候，四季分明，長期保持冬季寒冷干燥、夏季炎熱潮濕、春秋氣候適宜的特征。省內具有黑龍江、松花江、烏蘇里江三大水系，但人均水資源占有量低于全國水平。作為重要能源及原材料基地，黑龍江省具有豐富的能源資源及林業資源，森林面積及木材產量位居全國前列，并擁有全國最大油田——大慶油田；同時作為農業大省，黑龍江省已連續十三年位列全國糧食產量首位，為維護國家安全、糧食安全、能源安全等做出巨大貢獻。

1.2 數據來源

從“中國氣象數據網”收集到1975—2021年（47年）

黑龍江省內共計32個站點的日降水數據，并根據13個行政區劃分，整理分析數據得到各地區月平均降水

數據。

1.3 研究方法

1.3.1 線性回歸分析法

線性回歸分析法作為數理統計方法，通過近似確定一組數據的線性關系，定量反映兩組及多組變量之間的相互關聯的依賴度[10]。線性回歸模型方程為：

式中，y為因變量；a、b為方程參數，其中，a<0表示整體隨時間變化呈下降趨勢，a>0表示整體隨時間變化呈上升趨勢。

1.3.2 Mann-Kendall檢驗法

Mann-Kendall趨勢檢驗法（簡稱M-K法）作為氣象研究推薦算法，具有計算簡便、實用性強的特點，并常用于水文序列趨勢分析領域[11]。Mann-Kendall檢驗法計算原理如下：

假設時間序列樣本（x1，x2，...，xn）是n個獨立、隨機分布的數據，構造趨勢檢驗統計量：

式中，n為樣本個數；xj、xk分別為第j年、第k年序列值（j>k），sgn（xj-xk）為指示函數，依據xj-xk的正負號分別取值，即：

隨機變量序列近似服從正態分布，其均值為0，方差表示為：

（3）根據式（2-4）及式（2-5）構建統計檢驗值：

通過統計檢驗值的正負判定時間序列變化趨勢，若檢驗值ZS為正值，表明時間序列呈現整體上升趨勢，若ZS為負值，則呈現下降趨勢；ɑ為檢驗顯著水平，當ɑ取0.05時，查正態分布表得=1.96。

2 結果與分析

基于1975—2021年（564個月）黑龍江各地區月平均降水數據，采用線性回歸方法分析長期降水波動情況，并確定擬合參數及擬合優度，如表1所示。結果表明，黑龍江省各地區擬合曲線線性參數均幾乎趨近于0，表明黑龍江省降水整體變化趨勢不顯著；黑龍江省各地區降水擬合優度（R2）最大值所在地區為大興安嶺（0.0036），表明大興安嶺地區降水變化趨勢最明顯，但各地區擬合優度均遠小于1，當R2越趨近于0，表明擬合效果越差，說明黑龍江省整體降水變化無明顯規律，整體趨勢性不明顯。

根據Mann-Kendall檢驗法對黑龍江省13個地區進行趨勢性檢驗并繪制表格，如表2所示。當統計檢驗值|Z|<1.96時，表明降水趨勢不顯著；當Z為正值，表明降水序列呈現上升趨勢，反之呈現下降趨勢。如表2所示，黑龍江省各地區降水序列統計檢驗值均為負值，表明黑龍江省各地區降水序列均呈現下降趨勢；同時，各地區降水序列統計檢驗值均滿足|Z|<1.96，其中，大興安嶺地區統計檢驗值的絕對值最靠近1.96，表明大興安嶺地區降水變化趨勢相較于其他地區更為顯著，但黑龍江省整體降水變化趨勢不明顯。綜上，Mann-Kendall檢驗結論與線性回歸方法所得結論相同，所得結果具有合理性。

3 討論

根據線性回歸分析及Mann-Kendall檢驗法對黑龍江省各地區降水變化特征及趨勢進行探究，結果表明各方法均得到黑龍江省整體降水變化雖呈現下降趨勢，但整體變化態勢并不顯著特點，何學敏等[12]在“農業氣象災害和氣溫降水對東北三省糧食產量影響評估”一文所得研究結論相同。同時，基于線性回歸方法及Mann-Kendall檢驗法得到大興安嶺地區降水變化趨勢相較于其他地區更為顯著。作為黑龍江省森林覆蓋面積最廣的地區，大興安嶺地區城建面積較少，森林資源豐富，城市中人類生產生活導致的“熱島效應”降低，生態系統可維持自然狀態，使得降水變化趨勢更趨近于自然狀態。所得結論與孫繼松等[13]所得結論相同。綜上，本文所得結論具有合理性，可為黑龍江省水文序列變化趨勢探究提供數據支撐。

4 結論

本文基于線性回歸方法及Mann-Kendall檢驗法對黑龍江省降水時空變化特征及趨勢進行研究，并獲得結論：

基于1975—2021年（47年）黑龍江各地區月平均降水數據，采用線性回歸方法構成擬合曲線，分析長期降水波動情況，結果如下：黑龍江省總體降水呈現較大波動，并在47年內長期呈現波動幅度較大現象，并具有冬天干燥、夏天潮濕的氣候特征；其中降水擬合優度（R2）最大值所在地區為大興安嶺（0.0036），表明大興安嶺地區降水變化趨勢最明顯，但各地區擬合優度均遠小于1，當R2越趨近于0，表明擬合效果越差，說明黑龍江省整體降水變化無明顯規律，整體趨勢性不明顯。

根據Mann-Kendall檢驗法對黑龍江省13個地區進行趨勢性檢驗，結果表明：黑龍江省各地區降水序列均呈現下降趨勢，但整體變化趨勢不顯著，且Mann-Kendall檢驗結論與線性回歸方法所得結論相同，所得結果具有合理性。

參考文獻

[1] 柏林，臧淑英，張麗娟，等.扎龍濕地降水變化非線性特征研究[J].地理與地理信息科學，2014，30（3）：5.

[2] 周鑫，郭正堂，秦利.近百年來自然和人為因素對亞洲季風降水影響的時間序列分析研究[J].中國科學：地球科學，2010，40（12）：1718-1724.

[3] 孫宏勇，胡春勝，張喜英，等.三河市40年來溫度和降水變化及對農業生產的影響[J].中國生態農業學報，2006（2）：173-176.

[4] 葉柏生，李翀，楊大慶，等.我國過去50a來降水變化趨勢及其對水資源的影響（Ⅰ）：年系列[J].冰川凍土，2004（5）：587-594.

[5] 李淼，夏軍，陳社明，等.北京地區近300年降水變化的小波分析[J].自然資源學報，2011，26（6）：1001-1011.

[6] 王若男，馬林，郭宗凱，等.近60年丹東地區降水年際變化分析[J].現代農業科技，2015（6）：247+254.

[7] 王玉雪，李波，王槿妍，等.基于Mann-Kendall檢驗法的北運河流域降水和徑流變化趨勢分析[J].北京水務，2022（1）：24-28.

[8] 胡萌，盛英武.青島年降水量和水資源量變化特征研究[J].水文，2022，42（1）：103-108+28.

[9] 蘇萬磊.基于Mann-Kendall的秦皇島市降水多時間尺度分析[J].水利科技與經濟，2020，26（11）：23-26.

[10] 宋兵.基于Mann-Kendall檢驗的王瑤水庫降水、徑流變化趨勢及突變分析[J].陜西水利，2018（3）：77-78+81.

[11] 叢子成.東北地區村鎮住宅綠色模式建構研究[D].西安：西安建筑科技大學，2014.

[12] 何學敏.農業氣象災害和氣溫降水對東北三省糧食產量影響評估[D].沈陽：沈陽農業大學，2018.

[13] 孫繼松，舒文軍.北京城市熱島效應對冬夏季降水的影響研究[J].大氣科學，2007（2）：311-320.