基于Mann-Kendall檢驗和小波分析的南昌降水變化特征研究

林思平 高震東 王秋霜 胡耀丹 楊舒惠

摘要：近年來，氣候變化已成為國際研究的熱點問題，其中全球變暖在自然生態與社會經濟方面都給可持續發展帶來了重要影響。為了解全球變暖以及土地利用方式轉變下南昌地區的氣候變化規律，在相關氣候要素中，選擇降水量作為衡量指標，利用Mann-Kendall趨勢檢驗和Morlet小波分析方法對其近64年降水量的年際變化趨勢以及其周期性變化特征進行研究。結果表明：（1）近64年來南昌地區降水變化不明顯，但總體年降水量呈略增加趨勢。（2）通過小波分析得出：降水量序列存在5個周期，即7年、11年、28年、50年和63年。主振蕩周期為28年，且從主振蕩周期尺度上來看，南昌在未來一段時間里將處于降水偏多階段。

關鍵詞：南昌;Mann-Kendall趨勢檢驗;小波分析;降水

20世紀以來，全球大部分地區都經歷了不同程度的升溫[1]，加速變暖的氣候已經引起了社會各界的高度關注。降水量是各類氣候指標中，具有定量化分析某地區氣候變化特征和趨勢的主要因素之一。近年來，相關學者和機構進行了不同類型的區域降水變化趨勢和規律研究[2，3]。這其中，不同尺度的區域降水分析和局域降水變化研究可以顯著提高地方應對中長期氣候變化和旱澇災害的能力。

南昌市位于我國江西省中部偏北，與我國最大的淡水湖泊——鄱陽湖毗鄰，是我國主要的糧食生產基地之一。近年來，在全球變暖的背景下，疊加南昌市土地利用方式的轉變、農田轉為其他用地等因素，區域內的能量平衡發生了很大程度的轉變。隨之而來的是區域氣候的變化，以及環境本底的改變。因此，本文利用Mann-Kendall趨勢檢驗和Morlet小波分析等技術，對南昌地區近64年來的降水變化趨勢進行規律分析。相關研究結果可以為南昌地區應對氣候變化，提高該地區中長期適應降水變化的能力。同時，該研究結果也可以為區域農業生產應對未來旱澇災害提供基礎數據。

1 研究區概況

南昌市位于江西省的中部偏北區域，贛江和撫河的下游沖積平原，處在鄱陽湖西南岸。東經115°27′至116°35′，北緯28°10′至29°11′之間。南昌地區以平原為主，占國土面積的35%以上，東南地區相對平坦，西北則地勢較高。總體上，南昌及周邊水網密布，湖泊面積大，總水域達2200km2以上，占30%左右。南昌屬于亞熱帶季風氣候，氣候濕潤、雨熱同期，農業生產自古以來較為發達。1980年代以來，南昌市的城市化進程逐漸加速，城區面積逐年擴大，人口增速較快，由此引起了一系列的城市環境氣候問題。

2 研究方法

（1）Mann-Kendall分析法。Mann-Kendall分析法可以對各類時間尺度下的氣候要素進行變化趨勢檢驗，該方法屬于一種非參數統計檢驗方法，因此具有較好的檢測分析效果，常被應用于氣溫和降水等序列的變化檢測中[4，5]。該方法的優勢主要在于計算過程較為簡便，檢測的時間尺度范圍廣，不易受到異常值的控制，并且在檢測之初無須選擇某種特定的分布函數。因此，本文采用該方法檢驗南昌市近64年的降水序列，定量反映降水變化趨勢的顯著性。

（2）小波分析法。Morlet小波分析是一種具有時頻多層次分辨功能的方法，它能夠有效的獲取水文氣象序列不同時間尺度的周期性特征，并能評估序列未來一段時間內的發展趨勢[6-10]。本文通過運用小波分析檢驗方法，能夠很好的獲取南昌市降水的周期性變化規律。

3 結果與分析

（1）降水量趨勢分析。通過結合國家基本氣象站數據信息，統計了南昌市1953—2016年的日降水數據，分析其降水量的年代變差系數。如圖1所示，20世紀60年代的最大降水量在研究期間最小，20世紀50年代的最大降水量在研究期間最大。由此可見，20世紀60年代可能屬于降水低谷區、20世紀50年代可能屬于一個降水豐沛區。從平均值這一屬性來看，20世紀90年代的降水平均值最大，屬于降水相對較豐富的時期。20世紀60年代的降水平均值最低，這一時間段則屬于一個降水偏少期[11]。

如圖2所示，近64年來該地區的降水變化的總體趨勢并不顯著。由5年滑動平均的計算結果來看，1960s、1980s和2000年之后的降水處于相對低值區。50年代初、70年代和90年代降水均較多。總體平均年降水量呈略增加趨勢。MK趨勢檢驗的Z統計量為1.755，在90%的可信度范圍內，說明從長期來看，該地區的降雨量序列呈顯著的上升變化趨勢。總體而言，1950—1960年代、1970—1980年、2000年前后、2010年前后降水量較大，超過了近60年的平均值。而1960—1970年，以及1980—1990年屬于降水的低值區間。通過統計分析，南昌地區降水具有一定的變化趨勢性。

（2）降水變化周期。通過運用小波分析方法，從而得到南昌降水的周期變化特征，如圖3所示，南昌地區近64年的降水具有一定的周期變化特征。小波變換可以將降水周期變化過程中的精細結構顯示出來，一般小波實部變換實部值的低值區和高值區是周期變化，這便可以進行降水變化的周期識別。由圖2的小波變換實部時頻分布圖可以看出，在7—8年、11—12年、24—28年等時間尺度上可以較明顯的識別出來，其小波變換實部值具有高低變化的過程。同時，大部分周期變化可以通過95%的置信度檢驗。由此說明，南昌地區的降水具有一定的變化規律性。總體而言，在長時間尺度上，南昌地區降水的Morlet小波變換實部時頻分布值較高，但是大部分時間段超過檢驗的置信度區間。

小波方差圖可以反映年南昌地區降水變化趨勢的能量譜特征[5]，因此可以進行主周期和不同時間尺度的分析。小波方差分析的結果如圖4所示。

結合圖3、圖4分析得出降水量序列存在5個主周期，即7年、11年、28年、50年和63年。前三個振幅比較小，后面振幅逐漸增強，主振蕩周期為28a。從主振蕩周期上看，南昌在未來一段時間里將處于降水偏多階段。

4 結論

本文通過采用Mann-Kendall趨勢檢驗、Morlet小波分析方法分析了南昌市的降水變化特征。結果表明，64年間南昌市的降水量呈略增加趨勢，且在周期性的分析中表現出未來將繼續增加的趨勢。總體而言，南昌地區降水變化還處在一定的波動范圍之內。未來農業和城市發展可能會導致新的用水需求，因此從社會可持續發展的角度考慮應該進一步進行工農業和生活用水的綜合調度與規劃。通過統籌不同部門的用水規劃，以此促進南昌地區的用水合理化與社會經濟協調發展。

參考文獻：

[1]萬智巍，賈玉連，蔣梅鑫.華南北熱帶11.5—2.5ka B.P.溫度集成重建與特征分析[J].熱帶地理，2018，38（5）：641-650.

[2]袁喆，許繼軍，王永強，洪曉峰，周彥辰，陳進.基于水熱平衡的撫河流域地表徑流長期預估[J].工程科學與技術，2019，51（1）：60-67.

[3]萬智巍，賈玉連，章鳴，洪祎君，蔣梅鑫.1470-2014年鄱陽湖流域干濕變化與特征分析[J].干旱區資源與環境，2018，32（6）：114-119.

[4]萬智巍，蔣梅鑫，洪祎君，賈玉連.1470-2014年鄱陽湖流域降水量重建與旱澇災害診斷[J].災害學，2018，33（4）：93-98.

[5]簡虹，駱云中，謝德體.基于Mann-Kendall法和小波分析的降水變化特征研究——以重慶市沙坪壩區為例[J].西南師范大學學報（自然科學版），2011，36（4）：217-222.

[6]袁喆，楊志勇，鄭曉東，袁勇.近50年來淮河流域降水時空變化特征分析[J].南水北調與水利科技，2012，10（2）：98-103.

[7]蔣佑承，劉蛟，張越關，商瀅.長江源區降水變化的時空分布規律分析[J].水資源研究，2020，9（3）：335-343.

[8]楊潔，謝保鵬，張德罡.基于InVEST模型的黃河流域產水量時空變化及其對降水和土地利用變化的響應[J].應用生態學報，2020，31（8）：2731-2739.

[9]李娜，王寶榮，安韶山，焦峰，黃倩.黃土高原草地土壤細菌群落結構對于降水變化的響應[J].環境科學，2020，41（9）：4284-4293.

[10]方鶴楠，何真，陳茜，夏天.1985—2018年查干湖時空變化及其對降水變化響應分析[J].中國農業信息，2020，32（1）：64-73.

[11]馮敏玉，魏麗，胡逢喜，萬正，肖金香.南昌地區近55年的氣溫與降水變化分析[J].江西農業大學學報，2009，31（04）：781-784.