呼和浩特市降水的時(shí)空變化及干旱特征分析

［關(guān)鍵詞］降水；干旱；時(shí)空分布；呼和浩特

全球氣溫持續(xù)升高，對(duì)生態(tài)環(huán)境、水資源等諸多方面均產(chǎn)生較大影響，間接造成干旱、暴雨等極端事件頻發(fā)的不良后果。中國(guó)降水量分布的時(shí)空差異較大，對(duì)農(nóng)牧業(yè)經(jīng)濟(jì)效益具有直接影響，全面揭示降水與干旱的時(shí)空分布特征，摸清其變化規(guī)律，是當(dāng)前急需開(kāi)展的基礎(chǔ)性研究之一。

已有諸多學(xué)者在全球[1-4]，歐洲、美洲、亞洲等多個(gè)地區(qū)，利用克里金插值[5、6]、SDSM降尺度[7、8]、極點(diǎn)對(duì)稱模態(tài)分解（ESMD）[9-11]等方法開(kāi)展了降水的時(shí)空變化研究，北美洲、歐洲等地降水量呈增加趨勢(shì)[12]，塞爾維亞近十年降水量呈明顯的下降趨勢(shì)[13]。中國(guó)降水量整體呈增加趨勢(shì)，尤其是在20世紀(jì)中期至21世紀(jì)初，但各季節(jié)降水總量差異較大，夏季最多，春、秋兩季次之，冬季最少，在空間上表現(xiàn)為青藏高原[14]、西北地區(qū)[15]和南方部分地區(qū)[16]呈增加趨勢(shì)，東北、華北地區(qū)呈減少趨勢(shì)[17-20]。亦有學(xué)者以全球[2]、歐洲[3]、亞洲[4]等地為研究區(qū)，采用Mann-Kendall 突變檢驗(yàn)[21-22]、線性傾向估計(jì)[23]、Morlet小波分析等方法[23-25]，在年[26]、季[27]等的時(shí)間尺度下，針對(duì)干旱的時(shí)空變化特征展開(kāi)研究，結(jié)果表明，20世紀(jì)下半葉以來(lái)，全球干旱化程度呈上升趨勢(shì)[2]，歐洲東部、南部干旱化趨勢(shì)顯著[3]。中國(guó)極端干旱事件日益頻發(fā)[28]，東北地區(qū)發(fā)生嚴(yán)重干旱事件的頻率最高[29]；西北地區(qū)除新疆南部干旱程度較重外，其他地區(qū)多以中旱為主[30]。

本文利用呼和浩特及其周邊地區(qū)22 個(gè)氣象站點(diǎn)1951-2018年降水量年（月）數(shù)據(jù)資料，使用氣候傾向率等方法，全面揭示1951-2018年呼和浩特降水、干旱的時(shí)空變化特征，不僅豐富了研究成果，也為氣候變化歸因分析、水資源問(wèn)題應(yīng)對(duì)等提供了參考。

1. 資料與方法

1.1 研究區(qū)概況

研究區(qū)位于華北北部，內(nèi)蒙古自治區(qū)中部，地勢(shì)整體由東北向西南傾斜，主要的地貌類型有丘陵、山地、平原等，主要河流為大黑河、小黑河，總流域面積1.767km²，如圖1所示，處于北溫帶，屬于溫帶大陸性季風(fēng)氣候，夏季短暫，雨熱同期，冬季漫長(zhǎng)，寒冷干燥。

1.2 數(shù)據(jù)來(lái)源

本次選取了呼和浩特及其周邊地區(qū)22個(gè)測(cè)站1951-2018年的逐年（月）降水量資料，主要來(lái)源于國(guó)家氣象科學(xué)數(shù)據(jù)中心（http：//data.cma.cn/）。

1.3 數(shù)據(jù)分析方法

（1）對(duì)于缺測(cè)數(shù)據(jù)，利用臨近的數(shù)據(jù)完整的站點(diǎn)數(shù)據(jù)采用相關(guān)分析和回歸分析方法[12]進(jìn)行插補(bǔ)延展。

（2）采用了氣候傾向率法[31]對(duì)于降水量變化趨勢(shì)進(jìn)行分析。

（3）采用變異系數(shù)表征降水變化的劇烈程度[32]，計(jì)算公式如下：

CV=|SD/MN|

式中：SD為標(biāo)準(zhǔn)差，MN為平均值。

（4）對(duì)于降水量及其特征的空間變化，采用ArcGIS克里金插值法[6]進(jìn)行插值。

（5）采用標(biāo)準(zhǔn)化降水指數(shù)（SPI）表征研究區(qū)的干旱狀況，并根據(jù)前人研究成果對(duì)干旱等級(jí)進(jìn)行劃分，見(jiàn)表1。

2. 結(jié)果

2.1 呼和浩特地區(qū)年際降水量的時(shí)間變化特征

本文對(duì)研究區(qū)及周邊地區(qū)22個(gè)氣象站點(diǎn)1951-2018年年際降水的變化趨勢(shì)進(jìn)行了分析，圖2為研究區(qū)多年平均降水量、變異系數(shù)、傾向率的空間分布。研究區(qū)多年平均降水量整體呈由南向北的遞減趨勢(shì)，并呈現(xiàn)階梯狀分布。研究區(qū)大部分站點(diǎn)年平均降水量在350-500 mm 之間，如圖2a所示；土默特左旗、呼和浩特等站點(diǎn)變異系數(shù)明顯大于其他地區(qū)，變異系數(shù)大多介于0.322～0.304，說(shuō)明這些代表站點(diǎn)降水量波動(dòng)程度高于其他地區(qū)，而武川和托克托等站點(diǎn)次之（0.249～0.267），如圖2b所示；研究區(qū)降水多年變化率由東南向西北遞增，整體呈東南減少、西北增加的分布規(guī)律，土默特左旗降水多年變化率最大（0.071mm/10a），其次是土默特左旗、呼和浩特，呼和浩特郊區(qū)站（-0.048～0.071 mm/10a）降水多年變化率整體較小，其余站點(diǎn)降水多年變化率大多在-0.405～-0.048 mm/10a之間，如圖2c所示。

2.2呼和浩特地區(qū)季節(jié)降水量的時(shí)間變化特征

圖3為研究區(qū)季節(jié)平均降水量、變異系數(shù)、傾向率的空間分布。由圖可知，研究區(qū)季節(jié)平均降水量整體呈由南向北的遞減趨勢(shì)，與年際降水量的空間分布規(guī)律一致。夏季降水大多在205～317mm之間，占全年降水量的61.41%左右，秋季降水量（81～106mm）多于春季（58～75mm），冬季最少（12mm左右）。研究區(qū)不同季節(jié)變異系數(shù)空間分布變化較大，但總體呈現(xiàn)由西向東遞減的趨勢(shì)，且冬季波動(dòng)程度最大（0.50～0.88），春季（0.54～0.71）、秋季（0.45～0.64）次之，夏季最?。?.35～0.47）。春季降水傾向率空間變化較小，大多介于0.07～0.29mm/10a，夏季傾向率自東西兩側(cè)向中部遞減。秋季在托克托站點(diǎn)最大，并以此為中心，向其他方向減少。冬季降水量?jī)A向率的分布與春季降水量?jī)A向率的分布類似。

2.3呼和浩特地區(qū)年際、季節(jié)干旱事件的變化特征

表2為研究區(qū)不同時(shí)間尺度下的干旱的發(fā)生次數(shù)及干旱頻率。研究區(qū)各時(shí)間尺度下干旱事件總數(shù)分別為21次、17 次、22 次、23 次和15 次，發(fā)生干旱的頻率分別為30.9%、25%、32.4%、33.8%和22.1%；年際尺度下發(fā)生輕旱的年數(shù)為11a，中旱、重旱的年數(shù)則均為5a；春季發(fā)生輕旱的年數(shù)為5a，中旱的年數(shù)為7a，重旱的年數(shù)為3a，且在1955年和1994年發(fā)生了2次特旱，SPI值分別為-2.18和-2.16；夏季發(fā)生輕旱的年數(shù)最多（11a），其次是重旱（5a）和中旱（5a），在1965年發(fā)生了1次特旱，SPI值為-2.05；秋季發(fā)生干旱事件的年數(shù)整體多于其他時(shí)間尺度，且發(fā)生干旱事件的年數(shù)隨干旱程度遞減，即輕旱（13a）gt;中旱（6a）gt;重旱（3a）gt;特旱（1a）；冬季發(fā)生干旱事件的年數(shù)最少，但發(fā)生年數(shù)整體隨干旱程度遞增，即發(fā)生輕旱的年數(shù)最少（4a），其次是中旱（5a），重旱年數(shù)最多（6a），并未發(fā)生特旱事件。

3. 結(jié)論

（1）研究區(qū)多年平均降水量整體自東南向西北遞減，以下降趨勢(shì)為主，托克托下降速度最快（0.3752mm/10a），和林格爾等站點(diǎn)次之，清水河最慢（-0.022mm/10a）；亦有8個(gè)站點(diǎn)降水量呈上升趨勢(shì)，土默特左旗上升最快（0.165mm/10a），武川最慢（0.02mm/10a）。

（2）研究區(qū)夏季多年平均降水量最大，秋季、春季次之，冬季最??；春季、秋季降水量以上升趨勢(shì)為主，夏、冬兩季反之。就變化速度而言，春季降水量上升速度（最大值出現(xiàn)在清水河，最小值出現(xiàn)在呼和浩特）整體大于秋季，夏季降水量減少速度（呼和浩特下降速度最快，武川最慢）整體大于冬季。

（3）研究區(qū)秋季發(fā)生干旱事件的頻率最高（33.8%），夏季（32.4%）、年際（30.9%）和春季（25%）次之，冬季最低（25%）；研究區(qū)發(fā)生輕旱事件的頻率整體最高，其次是中旱和重旱，特旱事件頻率最低。