任 黎，蔡文美，梁士奎，朱佳晨，宋姝萍，劉雨桐，楊文嬌

(1.河海大學水文水資源學院，江蘇 南京 210098；2.華北水利水電大學水資源學院，河南 鄭州 450046)

1 概述

一直以來，降水作為水文循環(huán)的必要環(huán)節(jié)，在水資源中占有舉足輕重的地位，對其進行特征分析有非常重要的意義。憑借對經濟及社會產生的極大影響，成為全球科學家的研究熱點。在全球變暖的大環(huán)境下，極端氣候事件發(fā)生概率大幅度增加，短時間強降水會造成城市內澇，給社會及人身安全造成極大的威脅，極端氣候正在“常態(tài)化”，成為人類社會必須時時面對的生存危機之一，對水資源的高效管理及水災害的治理是現如今人類面臨的重要挑戰(zhàn)。

位于云貴高原的貴州省，由于處于高原斜坡上導致大江大河地勢陡峭，其經濟發(fā)達區(qū)大多處于分水嶺及小型流域，因此隨著經濟地快速發(fā)展，地區(qū)資源性缺水問題日益突出。福泉市位于貴州省中部，作為亞洲磷都以及西南出海的交通要塞，近年來隨著經濟社會的快速發(fā)展，社會需水總量升高，水資源利用率較低。因此研究福泉市的降水時空變化特征，為洪澇、干旱災害的防治提供依據，保障全市人民生命財產安全，促進研究區(qū)經濟發(fā)展具有重要意義。

2 研究區(qū)概況及方法

2.1 研究區(qū)概況

福泉市地處于貴州省中部，黔南布依族苗族自治州北部。東鄰黔東南的凱里市和黃平縣，南邊與黔東南的麻江縣接壤，西與貴定縣、龍里縣、開陽縣相鄰，北與甕安縣相連。福泉市地處云貴高原東部，屬亞熱帶季風氣候華中濕潤區(qū)，熱量豐富、雨量充沛，無霜期長，年均溫14℃左右，無霜期245～278d。年均降水量1033～1220mm，雨季一般在四月中旬開始，十月下旬結束。地勢西高南低，從西北到東南地形逐漸傾斜，其位于烏蒙山脈至苗嶺山脈的中段。全境河流隸屬于長江流域，以武夷山脈為分水嶺，劃分為沅江流域和烏江流域[7]。福泉市全境干、支流之數多達39條，其總河長為633.7km，河網密度為0.38km/km2，即每1平方公里就有0.38km的河流。

2.2 研究方法

M-K檢驗方法是一種非參數檢驗方法，也稱無分布檢驗，它無需遵循某一特定分布，且不會受到少量的異常影響的特點，對類型、次要變量具有更好的適應性，在降水、氣溫、徑流等非正態(tài)分布數據方面被廣泛應用。該方法的優(yōu)勢在于無需采樣遵循一定的分配，且不會受到少量的異常影響，而且適用于類型變量和連續(xù)變量，且易于進行運算。當Z>0時，則表示幅度增加，當Z<0時，則表示幅度減少。當|Z|≥1.28時，表示顯著性檢驗大于90%，當|Z|≥1.64時，表示顯著性檢驗大于95%，當|Z|≥2.32時，表示顯著性檢驗大于99%。

M-K突變檢驗方法以氣候、下墊面一致性為前提，要求序列隨機獨立且同分布，能夠判定氣候序列中的氣候突變，若有，則能判定其何時出現。該方法不僅計算簡單，還能確定變異起始時間，并能確定變異區(qū)，因此，是一種常用的突變檢測方法。在坐標軸上繪制UFk和UBk兩條曲線，當它們超過臨界置信水平直線時(檢驗置信水平α=0.05時，置信水平線為±1.96)，即表示上升或下降的趨勢顯著，突變的時間為超過臨界線的區(qū)間，若UFk和UBk兩條曲線出現交點且交點位于臨界線之間，則該點為便是突變開始的時間。

R/S分析方法又稱為重標極差分析法，它是英國水利學者赫斯特在尼羅河大壩項目中提出的一種分析方法，是可以對時間序列進行分析的方法。在描述自然界非線性系統(tǒng)的復雜性，揭示其系統(tǒng)內部本質規(guī)律及找出自身特點有著良好的效果，可以用來檢驗時間序列是否遵循獨立，Hurst指數也能很好地揭示出時間序列中的趨勢性成分，并且能夠代表趨勢性成分強度，目前已廣泛用于對不同的時間序列進行分析。Hurst指數分級見表1。

表1 Hurst指數分級表

RClimDex模型是一種用于計算極端氣候指數的軟件，基于R語言編輯器的開發(fā)，以Excel為基礎，旨在提供容易使用的軟件包去計算并監(jiān)測氣候變化的氣候極值。RClimDex模型會通過仔細篩查數據中是否存在無效數據，如果存在，就可以通過用戶以標準差的形式去定義閾值去除異常值。最后通過計算將會獲得27個核心極端氣候指數[10]。

3 降水特征分析

3.1 顯著性分析

將福泉市60年的降水數據利用Mann-Kendall檢驗法進行檢驗，結果見表2。

由表2可知，1959—2018年間，年降水量|Zc|=0.64，小于1.64，未通過顯著性檢驗，故福泉市年降水量呈現出不顯著的下降趨勢；春、夏兩季降水量|Zc|=0.59、|Zc|=0.03，均小于1.64，均未通過顯著性檢驗，故福泉市春、夏兩季降水量呈現出不顯著的下降趨勢；秋季降水量|Zc|=1.75，大于1.64，通過顯著性檢驗，故福泉市秋季降水量表現出明顯的下降趨勢；冬季降水量|Zc|=0.63，大于1.64，未通過顯著性檢驗，故福泉市冬季降水量表現出不明顯的上升趨勢。

表2 M-K指數結果

3.2 突變性分析

根據M-K突變檢驗法，將福泉市60年降水量數據繪制成M-K突變檢驗圖，如圖1所示。從圖1中可以看出，在1959—1985年間，除1961—1963年外，UFk曲線取值基本上位于0值以上，表明該期間，福泉市年降水量一直呈現出上升的趨勢。自1985年后，UFk曲線取值均位于0值以下，表明自20世紀80年代后，福泉市年降水量一直呈現出下降的趨勢。圖中UFk和UBk曲線存在較多交叉點，再結合變化趨勢分析的結果，認為福泉市年降水量在1959—2018年不存在顯著性突變時間點。

圖1 年降水量M-K突變檢驗

將福泉市季度降水量數據繪制成M-K突變檢驗圖，如圖2所示。圖2(a)中UFk曲線在1959—1975年間波動較大，自1975年后UFk曲線均低于0值，且在1996—1998、2001—2002年間超出信度線，表明自20世紀70年代后研究區(qū)春季降水量一直呈現出下降的趨勢，尤其在1996—1998、2001—2002年間下降的趨勢明顯。圖中UFk和UBk曲線存在較多交叉點，再結合變化趨勢分析結果，認為研究區(qū)年降水量在1959—2018年不存在顯著性突變時間點。圖2(b)中UFk曲線除1961—1962年外，均處于0值以上，在1964—1965年超出信度線，表明在1959年來研究區(qū)夏季降水量處于不明顯上升趨勢，在1964—1965年間為突變區(qū)間呈現出顯著性上升趨勢。圖2(c)中，在1959—1975年間UFk曲線在0值上下波動，變化趨勢不明顯。自1975年后，UFk曲線均處于0值以下且在1992年后UFk曲線在信度線上下波動，表明20世紀70年代后，研究區(qū)秋季降水量一直處于下降趨勢，特別從1992年后為突變區(qū)間呈現出顯著性下降趨勢。圖2(d)中，在1959—1970年間UFk曲線位于0值以下，表明該期間研究區(qū)冬季降水量處于不顯著的下降趨勢。自1970年后，除1983—1986、1988—1989年外UFk曲線均處于0值以上，且一直處于信度線之間，表明自20世紀70年代以來，研究區(qū)冬季降水量呈現出不顯著的下降趨勢。UFk和UBk曲線相交多點，結合趨勢結果分析可得出冬季降水量于1970年發(fā)生突變。

3.3 R/S分析方法

根據研究區(qū)逐日降水資料，利用Matlab軟件進行分析可得，Hurst指數H=0.1562，介于0～0.5之間，未來與過去的信息在時間尺度上表現出長程負相關性，時間序列具有較強反持續(xù)性，表示過去增加意味著未來減少，而過去減少意味著未來增加。年降水量呈現出下降趨勢意味著未來年降水量將會呈現出增加趨勢，并且其反持續(xù)性較強，所以未來的增加趨勢較明顯。

運用Matlab軟件使用R/S分析方法可得出各季度Hurst系數及相關指數，見表3。

表3 R/S分析相關系數表

從表3中可以看出四季度的Hurst系數均處于0～0.5之間，表示未來與過去的信息在時間尺度上表現出長程負相關性，時間序列具有反持續(xù)性。春季Hurst系數為0.3426，具有較弱反持續(xù)性，意味著未來春季降水趨勢與現在相反但關聯性較弱。夏季Hurst系數為0.2403，具有較強反持續(xù)性，意味著未來夏季降水趨勢與現在相反且關聯性較強。秋季Hurst系數為0.3861，具有弱反持續(xù)性，意味著未來秋季降水趨勢與現在相反但是關聯性弱。冬季Hurst系數為0.4202,具有弱反持續(xù)性且接近于0.5，意味著未來冬季降雨與現在相反,但關聯性弱。由于Hurst系數接近于0.5，即盡管過去60年冬季降水量增加，但未來增加趨勢減緩還有可能會降低，或者說它更接近于隨機變化，即冬季的降水趨勢與過去無關，是隨機的。

圖2 福泉市季度降水量M-K突變檢驗

3.4 極端降水變化特征分析

在1998—2001年間，世界氣象組織(WMO)在全球氣候變化監(jiān)測大會上，推出了一系列的極端氣候指標，并且成為氣候變化研究的統(tǒng)一標準。這些極端氣候指數是由日氣溫和日降水數據計算而得。本文選取降水強度(SDII)、強降水量(R95p)、日最大降水量(Rx1Day)、5日最大降水量(Rx5Day)、連續(xù)無雨天數(CDD)、連續(xù)有雨日數(CWD)這6個極端降水指數進行研究區(qū)極端降水特征分析，見表4。

表4 極端降水指數表

運用Rclimdex軟件計算福泉市6個極端降水指數，各項指數年際變化如圖3所示，從圖3中可以

圖3 福泉市極端氣候分析

整體看出，60年間福泉市6個極端降水指數(除連續(xù)有雨天數外)總體趨勢變動較為平緩。從圖3(a)可以看出，降水強度(SDII)反映平均降雨情況，總體呈現出下降趨勢，平均每10年下降約0.16mm，整體圍繞12mm上下波動。從圖3(b)可以看出，強降水量(R95p)反應極端降水量在總降水量中所占的比例，總體呈現出下降趨勢，平均每10年下降約15.41mm，整體波動幅度較大，表明福泉市極端降水發(fā)生頻率差異較大，極端現象發(fā)生概率不穩(wěn)定。從圖3(c)、(d)可以看出，日最大降水量(Rx1Day)、5日最大降水量(Rx5Day)總體呈現出下降趨勢，日最大降水量(Rx1Day)平均每10年下降約6.4mm，5日最大降水量(Rx5Day)平均每10年下降約11.8mm，在1965、1967、2012年出現了較為明顯的極大值，說明這3年日最大降水量(Rx1Day)出現了較為極端的情況。從圖3(e)可以看出，連續(xù)無雨天數(CDD)總體呈現出上升趨勢，平均每10年上升約0.99d，2009—2012年間波動幅度劇烈，表明這幾年間福泉市降水頻率較為極端，干旱現象頻繁。從圖3(f)可以看出，連續(xù)有雨天數(CWD)總體呈現出下降趨勢，平均每10年下降約0.04d，1959—2018年間波動幅度劇烈，極大值、極小值分布不均勻，1992—1996年連續(xù)有雨天數低于均值。

4 結論

(1)福泉市60年間年降水量沒有顯著的下降的趨勢，無明顯突變年份。春夏兩季降水量呈不顯著下降趨勢，秋季呈顯著性下降趨勢，冬季呈不顯著上升趨勢；秋季降水量突變時段較長，其余季度突變年不明顯。

(2)6個極端降水指數除連續(xù)無雨天數(CWD)呈上升趨勢外，其余均呈下降趨勢。強降水量(R95p)整體波動幅度較大，連續(xù)有雨天數(CDD)、連續(xù)無雨天數(CWD)波動劇烈；上述極端降水特征表示福泉市極端降水情況頻率較高，近年來，極端氣候事件增加，干旱現象頻發(fā)。