陳燕麗，謝 映，張 會(huì)，韋春霞

（1.廣西壯族自治區(qū)氣象科學(xué)研究所，廣西南寧 530022；2.保定市氣象局，河北保定 071030）

極端氣候事件具有破壞性強(qiáng)、可預(yù)測性差的特點(diǎn)，對生態(tài)環(huán)境、生產(chǎn)、生活及經(jīng)濟(jì)發(fā)展均有較大影響，已引起社會(huì)各界的廣泛關(guān)注[1-2]。學(xué)者們采用趨勢分析和空間分析等方法研究不同地區(qū)極端氣候事件變化規(guī)律，發(fā)現(xiàn)極端氣溫和降水事件存在次數(shù)增加、強(qiáng)度增大的趨勢[3-7]。喀斯特地貌區(qū)是廣西典型的生態(tài)脆弱區(qū)，石漠化現(xiàn)象突出，特殊的生境使得該地區(qū)植被對氣候變化響應(yīng)非常敏感且承災(zāi)能力弱，摸清該地區(qū)極端氣候事件變化規(guī)律可為其石漠化治理、植被保護(hù)和生態(tài)修復(fù)提供科學(xué)參考。

水熱分配不均導(dǎo)致的旱澇災(zāi)害是影響喀斯特地區(qū)植被和石漠化最重要的因素?？λ固氐貐^(qū)的植被比非喀斯特地區(qū)更易受干旱影響[8]，且喀斯特地區(qū)受旱程度越重，植物對干旱的敏感性越強(qiáng)[9]。暴雨對喀斯特植被的影響存在閾值效應(yīng)；在貴州省喀斯特地區(qū)，降水要素與石漠化呈顯著正相關(guān)，暴雨量級的增加會(huì)加速石漠化的發(fā)生和發(fā)展[10]。采用常規(guī)指數(shù)分析喀斯特地區(qū)氣象災(zāi)害時(shí)空規(guī)律已有研究[11-12]。近年來，西南喀斯特地區(qū)的極端災(zāi)害事件發(fā)生次數(shù)明顯增多，嚴(yán)重威脅喀斯特地區(qū)植被健康生長和石漠化治理。目前，有關(guān)喀斯特地區(qū)極端氣候事件時(shí)空演變規(guī)律的報(bào)道較少?？卒h[3]詳細(xì)報(bào)道了全國極端冷暖事件變化格局，針對廣西地區(qū)的研究僅分析極端氣溫的變化且選用的氣象站點(diǎn)較少，研究結(jié)論在廣西喀斯特地區(qū)應(yīng)用時(shí)有局限性。

為了更全面研究廣西喀斯特地區(qū)極端氣候事件時(shí)空變化規(guī)律，本研究基于長時(shí)間序列歷史氣象觀測資料，選擇世界氣象組織推薦的6個(gè)極端氣溫指數(shù)和2個(gè)極端降水指數(shù)，采用Mann-Kendall法和線性趨勢分析法詳細(xì)分析廣西喀斯特地區(qū)極端氣候事件時(shí)空變化特征，擬為該地區(qū)石漠化治理、植被保護(hù)修復(fù)和氣象防災(zāi)減災(zāi)等提供參考。

1 材料與方法

1.1 研究區(qū)概況

廣西壯族自治區(qū)（104°26'～112°04'E，20°54'～26°24'N）地處中國南部，形似盆地，山地多、平原少；喀斯特地貌區(qū)面積占廣西總面積的37.8%，土地貧瘠，石漠化現(xiàn)象突出，自然災(zāi)害嚴(yán)重，嚴(yán)重制約了當(dāng)?shù)亟?jīng)濟(jì)發(fā)展，屬于全國貧困地區(qū)。廣西喀斯特地區(qū)分布如圖1所示。

圖1 廣西喀斯特地區(qū)分布Fig.1 Distribution of Guangxi karst areas

1.2 氣象數(shù)據(jù)

氣象數(shù)據(jù)為廣西氣象信息中心提供的69個(gè)氣象站點(diǎn)1961—2019年逐日的最高、最低氣溫和降水?dāng)?shù)據(jù)。氣溫和降水是影響喀斯特地區(qū)植被生長最重要的氣象因子，本研究選取世界氣象組織（WMO）推薦的27個(gè)極端氣候指數(shù)中的6個(gè)極端氣溫指數(shù)和2個(gè)極端降水指數(shù)（表1）。

1.3 研究方法

采用RClimDex軟件計(jì)算極端氣候指數(shù)值，采用Mann-Kendall法對廣西喀斯特地區(qū)極端氣候指數(shù)進(jìn)行突變檢驗(yàn)，采用氣候傾向率分析極端氣候事件的年際變化趨勢。

Mann-Kendall法是一種非參數(shù)統(tǒng)計(jì)檢驗(yàn)方法，用于明確時(shí)間變化序列突變開始時(shí)間[13-14]，計(jì)算方法為：

假定時(shí)間序列為{xi}（i=1，2，3，……，n），定義統(tǒng)計(jì)量（dk）為：

式中，mi為第i個(gè)樣本xi＞xj[1≤j≤i]的累計(jì)數(shù)。

在時(shí)間序列獨(dú)立的假定下，定義統(tǒng)計(jì)量為：

式中，E(dk)和var(dk)分別為統(tǒng)計(jì)量dk的均值和方差。UFk為標(biāo)準(zhǔn)正態(tài)分布序列；當(dāng)UF＞0時(shí)，表示該時(shí)間序列呈上升趨勢；當(dāng)UF＜0時(shí)，為下降趨勢；UF和UB曲線在臨界線內(nèi)存在交點(diǎn)，該交點(diǎn)對應(yīng)的時(shí)間為突變發(fā)生的時(shí)間。當(dāng)UF超過95%置信水平（α=1.96）時(shí)，說明指數(shù)上升或下降趨勢顯著。

利用氣候傾向率分析極端氣候指數(shù)（y）的年際變化趨勢和大小，一般采用一次線性方程表示，計(jì)算公式為[15]：

式中，t為年份，t=1，2，…，n；b×10為氣候傾向率；a為常數(shù)。氣候傾向率為正時(shí)，表示極端氣候指數(shù)隨時(shí)間呈上升趨勢，反之呈下降趨勢。

2 結(jié)果與分析

2.1 極端氣溫指數(shù)時(shí)空變化趨勢

2.1.1 時(shí)間變化趨勢

1961—2019 年，廣西喀斯特地區(qū)TX10p呈減少趨勢，變化不顯著，氣候傾向率約為-0.6 d/10a（圖2a）；TN10p呈顯著減少趨勢（P＜0.05），氣候傾向率約為-2.9 d/10a（圖2b）；TX90p和TN90p均呈顯著增加趨勢（P＜0.05），兩者氣候傾向率分別為2.9和4.7 d/10a（圖2c～d）；TX10p減少幅度小于TN10p，TX90p增加幅度小于TN90p，說明廣西喀斯特地區(qū)夜間氣溫的增加幅度大于白天。CSDI呈顯著減少趨勢（P＜0.05），氣候傾向率為-1.8 d/10a（圖2e）；WSDI呈增加趨勢，變化不顯著，氣候傾向率為1.0 d/10a（圖2f）。說明廣西喀斯特地區(qū)熱天氣持續(xù)日數(shù)在增加。6個(gè)極端氣溫指數(shù)中，TN10p、TX90p、TN90p和CSDI的年際變化趨勢均能較好地反映喀斯特地區(qū)近60年的時(shí)間變化趨勢。廣西喀斯特地區(qū)對全球變暖為正響應(yīng)，其氣溫變化趨勢以變暖為主。

圖2 廣西喀斯特地區(qū)極端氣溫指數(shù)變化趨勢Fig.2 Trends of extreme temperature indices in Guangxi karst areas

MK檢驗(yàn)結(jié)果表明，近60年廣西喀斯特地區(qū)TX10p未發(fā) 生 顯著 突變，UF和UB曲線 在1962—1964年和1996—1997年相交，但未對整體變化趨勢產(chǎn)生較大影響（圖3a）。TN10p的UF和UB曲線在1988—1990年相交，TN10p發(fā)生突變，其變化前后分別為41.3和31.1天；UF曲線在1997年后超過了0.05顯著性水平線，說明TN10p呈顯著下降趨勢（圖3b）。TX90p的UF和UB曲線在2004—2005年相交，TX90p發(fā)生突變，變化前后分別為31.6和46.5天；UF曲線在1970—1980年超過了0.05顯著性水平線，說明TX90p呈顯著下降趨勢；2009年后超過0.05顯著性水平線，說明TX90p呈顯著上升趨勢（圖3c）。TN90p的UF和UB曲線在1996—1998年相交，TN90p發(fā)生突變，變化前后分別為29.0和44.3天；UF曲線在1982年后超過0.05顯著性水平線，說明TN90p呈顯著上升趨勢（圖3d）。CSDI的UF和UB曲線在1990—1994年相交，CSDI發(fā)生突變，突變前后分別為21.3和14.6天；UF曲線在2017年超過0.05顯著性水平線，說明CSDI呈顯著下降趨勢（圖3e）。WSDI的UF和UB曲線在2004—2006年、2010—2013年相交，WSDI在該期間波動(dòng)劇烈，除1970—1977年外，其他時(shí)段UF曲線均在0.05顯著性水平線內(nèi)，總體變化趨勢不明顯；在1961—1977年，WSDI呈下降趨勢，其UF曲線在1970—1977年超過0.05顯著性水平線，說明WSDI呈顯著下降趨勢（圖3f）。

圖3 廣西喀斯特地區(qū)極端氣溫指數(shù)突變特征Fig.3 Mutation characteristics of extreme temperature indices in Guangxi karst areas

2.1.2 空間變化特征

1961—2019 年，TX10p在空間上多呈下降趨勢，氣候傾向率為-20.6～2.8 d/10a，下降站點(diǎn)數(shù)占比為95.7%，約2.9%站點(diǎn)通過0.05顯著性水平（圖4a）。河池市中部和崇左市中偏東部地區(qū)，TX10p增加趨勢較明顯；桂林市中部、東部和百色市東北部、西南部，TX10p減少趨勢較明顯。

TN10p在空間上全部呈下降趨勢，氣候傾向率為-45.1～-13.7 d/10a，下降站點(diǎn)數(shù)占比達(dá)100%，約94.2%站點(diǎn)通過0.05顯著性水平（圖4b）。桂林市中偏南部、北部和百色市西北部，TN10p減少趨勢較明顯；河池市南部、西部及崇左市西南部，TN10p減少趨勢較弱。

TX90p在空間上多呈上升趨勢，氣候傾向率為-3.8～88.9 d/10a，上升站點(diǎn)數(shù)占比達(dá)97.1%，約68.1%站點(diǎn)通過0.05顯著性水平（圖4c）。TX90p由西南向東北呈減少趨勢；百色市東北部、西南部，TX90p增加趨勢較明顯；桂林市東北部、河池市中部、柳州市西南部、來賓市西北部和賀州市中部，TX90p增加趨勢較弱。

TN90p在空間上多呈上升趨勢，氣候傾向率為-10.5～100.0 d/10a，上升站點(diǎn)數(shù)占比98.6%，約89.9%站點(diǎn)通過0.05顯著性水平（圖4d）。TN90p變化趨勢空間分布規(guī)律與TX90p相似。

CSDI在空間上全部呈下降趨勢，氣候傾向率為-32.5～-4.8 d/10a，下降站點(diǎn)數(shù)占比100%，約42.0%站點(diǎn)通過0.05顯著性水平（圖4e）。西北部和東北部CSDI減少趨勢高于中部和東南部；桂林市南部和百色市西部、北部，CSDI減少趨勢較明顯。

WSDI在空間上多呈上升趨勢，氣候傾向率為-7.7～37.6 d/10a，上升站點(diǎn)數(shù)占比84.1%，約31.9%站點(diǎn)通過0.05顯著性水平（圖4f）。河池東南部與來賓市交界處、桂林市東北部、百色市東北部及崇左市中南部部分地區(qū)，WSDI減少趨勢較明顯；百色市西南部、東北部和桂林市中部，WSDI增加趨勢較明顯。

圖4 廣西喀斯特地區(qū)極端氣溫指數(shù)變化趨勢空間分布Fig.4 Spatial distributions of change trends of extreme temperature indices in Guangxi karst areas

2.2 極端降水指數(shù)時(shí)空變化趨勢

2.2.1 時(shí)間變化趨勢

1961—2019 年，廣西喀斯特地區(qū)R95p呈增加趨勢，變化不顯著，氣候傾向率約為19.9 mm/10a；R99p呈顯著增加趨勢（P＜0.05），氣候傾向率約為14.6 mm/10a（圖5a～b）。廣西喀斯特地區(qū)降水變化以增加趨勢為主，極強(qiáng)降水量增加趨勢比強(qiáng)降水量更明顯。MK檢驗(yàn)結(jié)果表明，59年間廣西喀斯特地區(qū)R95p和R99p均并未發(fā)生顯著變化，R95p、R99p的UF和UB曲線在多個(gè)時(shí)段相交，但UF曲線仍在0.05顯著性水平線內(nèi)，說明盡管部分時(shí)期強(qiáng)降水量和極強(qiáng)降水量發(fā)生較大波動(dòng)，但并未對整體變化趨勢產(chǎn)生較大影響（圖6a～b）。

圖5 廣西喀斯特地區(qū)極端降水指數(shù)變化趨勢Fig.5 Change trends of extreme precipitation indices in Guangxi karst areas

圖6 廣西喀斯特地區(qū)極端降水指數(shù)突變特征Fig.6 Mutation characteristics of extreme precipitation indices in Guangxi karst areas

2.2.2 空間變化特征

R95p在空間上多呈上升趨勢，氣候傾向率為-51.8～707.0 mm/10a，上升站點(diǎn)數(shù)占比89.9%，但僅11.6%站點(diǎn)通過0.05顯著性水平（圖7a）。R95p由西南向東北呈增加趨勢，桂林市大多數(shù)地區(qū)增加趨勢明顯高于其他地區(qū)，崇左市中部、北部及其與南寧市和百色市的交界地區(qū)減少趨勢較明顯。R99p在空間上多呈上升趨勢，氣候傾向率為-216.0～651.0 mm/10a，上升站點(diǎn)數(shù)占比87.0%，但僅18.8%站點(diǎn)通過0.05顯著性水平（圖7b）。R99p的空間分布與R95p相似，由西南向東北呈增加趨勢，桂林市中部增加趨勢明顯高于其他地區(qū)，崇左市南部減少趨勢較明顯。

圖7 廣西喀斯特地區(qū)極端降水指數(shù)變化趨勢空間分布Fig.7 Spatial distributions of change trends of extreme precipitation indices in Guangxi karst areas

3 討論

研究發(fā)現(xiàn)廣西喀斯特地區(qū)極端氣溫和極端降水總體表現(xiàn)為增加趨勢，極端氣溫和極端降水增加、減少較明顯的地區(qū)與平均氣溫和平均降水的高低值中心較符合[16]，說明本身較暖的地方趨于更暖，降水多的地方更容易發(fā)生強(qiáng)降水事件。本研究時(shí)間尺度為年，從季節(jié)尺度上研究極端氣候指數(shù)演變特征可進(jìn)一步研究冬季暖冬、夏季洪澇的變化趨勢。已有研究利用月平均降水和氣溫計(jì)算標(biāo)準(zhǔn)化降水蒸散指數(shù)研究廣西喀斯特地區(qū)干旱特征，發(fā)現(xiàn)研究區(qū)冬旱和夏旱呈波動(dòng)減弱趨勢，春旱和秋旱呈增強(qiáng)趨勢[11]，如能結(jié)合極端氣溫和極端降水指數(shù)時(shí)空演變特征，可為該地區(qū)重旱、特旱的發(fā)展趨勢預(yù)測提供更科學(xué)的參考依據(jù)。

近60年廣西喀斯特地區(qū)極端氣溫和降水的變化呈現(xiàn)較大的隨機(jī)性和不連續(xù)性。6個(gè)極端氣溫指數(shù)中，冷晝?nèi)諗?shù)、冷持續(xù)日數(shù)和暖持續(xù)日數(shù)突變并不明顯。對于極端降水指數(shù)而言，盡管強(qiáng)降水量和極強(qiáng)降水量均呈增加趨勢，但兩者的突變趨勢均不明顯，說明廣西喀斯特地區(qū)極端氣溫和降水事件變化的不規(guī)律性給該地區(qū)氣象災(zāi)害防御帶來極大挑戰(zhàn)。研究中采用的極端氣候指數(shù)為國際標(biāo)準(zhǔn)[17]；其中，極端氣溫指數(shù)為10%和90%分位數(shù)，極端降水指數(shù)為5%和1%分位數(shù)。研究區(qū)某些指標(biāo)的變化趨勢并不明顯，如對于極端降水指數(shù)，R95p（5%分位數(shù)）為不顯著增加，R99p（1%分位數(shù)）為顯著增加，進(jìn)一步研究中，針對極端氣候指數(shù)可補(bǔ)充計(jì)算更多指標(biāo)體系的極端氣候指數(shù)，更全面分析極端氣候事件變化情況，為定量評估極端氣候事件對喀斯特地區(qū)植被影響提供更科學(xué)的參考依據(jù)。

4 結(jié)論

本研究分析廣西喀斯特地區(qū)極端氣溫和極端降水時(shí)間的時(shí)空演變規(guī)律。廣西喀斯特地區(qū)對全球變暖為正響應(yīng)，極端氣溫變化以變暖趨勢為主，白天和夜間均趨于更暖，且夜間氣溫的增加幅度大于白天，表明廣西喀斯特地區(qū)近60年來白天和夜間的高溫天氣事件明顯增多。極端降水指數(shù)總體呈增加趨勢，極強(qiáng)降水量表現(xiàn)尤其明顯，表明研究區(qū)極強(qiáng)降水事件增加，洪澇災(zāi)害發(fā)生概率增大。不同極端氣溫指數(shù)空間分布上存在較大異質(zhì)性，總體而言研究區(qū)中部增溫趨勢較小，西北部（百色市）增溫趨勢較明顯。極端降水指數(shù)空間異質(zhì)性比極端氣溫小，極端降水增加趨勢最明顯地區(qū)主要分布在東北部（桂林市），減少趨勢最明顯的地區(qū)主要分布在西南部（崇左市）。