張涵丹，衛(wèi) 偉，薛 萐

（1.中國科學(xué)院 生態(tài)環(huán)境研究中心城市與區(qū)域生態(tài)國家重點實驗室，北京100085；2.中國科學(xué)院大學(xué) 資源與環(huán)境學(xué)院，北京100049；3.中國科學(xué)院 水利部水土保持研究所，陜西楊凌712100）

近幾十年來，高溫、熱浪、冰川消融、海平面上升、永久凍土層融化等事件頻繁發(fā)生，氣候變化問題引起了人們的普遍關(guān)注。眾多觀測事實表明，全球氣候正在經(jīng)歷一次以增溫為主要特征的變化趨勢：20世紀(jì)的升溫速率是近1000年來最高的時期，而近30a又是20世紀(jì)升溫率最高的時期［1－2］。全球氣候變暖勢必會對自然生態(tài)系統(tǒng)和人類的生產(chǎn)生活產(chǎn)生直接或間接的影響，尤其是生態(tài)敏感性和生態(tài)脆弱性極高的干旱半干旱地區(qū)。研究表明，半干旱區(qū)在近100年來冬季升溫趨勢超過全球陸地平均氣溫升溫趨勢的35.4%，對全球陸地變暖的貢獻(xiàn)率達(dá)44.46%［3］。不同時期、不同地區(qū)的干旱半干旱區(qū)干濕變化趨勢呈現(xiàn)顯著差異性［4－5］：南非、南美洲、北美洲、中亞、中國西北等全球主要干旱區(qū)在1860—1995年這百年呈現(xiàn)暖干化趨勢［6］；澳洲和南美洲在1951—2002年溫度和降水均升高，而依然呈現(xiàn)干旱化趨勢［4，7］；中國西北地區(qū)氣溫持續(xù)增加［8－9］，降水量變化的波動性較大［10－12］，東北地區(qū)依然以暖干化趨勢為主［13－14］。

定西市位于黃土高原西南緣，處于半干旱半濕潤向干旱荒漠化轉(zhuǎn)變的過渡地帶，該地區(qū)水資源匱乏，降水量少而集中，水蝕與風(fēng)蝕嚴(yán)重，是黃土高原生態(tài)脆弱地區(qū)的典型代表。天然降水是主要的水分來源，農(nóng)林業(yè)的生產(chǎn)對氣候條件的依賴性極強(qiáng)［15］，氣候變暖及生態(tài)環(huán)境的干旱化程度對區(qū)域的經(jīng)濟(jì)影響重大［16－18］。作為西部大開發(fā)及扶貧攻堅的重點地區(qū)之一，為改善其生態(tài)環(huán)境現(xiàn)狀，植樹造林、退耕還林還草、輸引水源等諸多工程措施相繼得到開展［19－21］，在長期的歷史進(jìn)程中，區(qū)域氣候也在相應(yīng)的改變。毋庸置疑，基于長期氣象數(shù)據(jù)，對區(qū)域氣候變化特征進(jìn)行綜合分析，將為實現(xiàn)該地區(qū)的生態(tài)恢復(fù)與建設(shè)，深入了解西北半干旱地區(qū)的氣候變化特征提供重要科學(xué)依據(jù)。

1 材料與方法

1.1 數(shù)據(jù)來源

數(shù)據(jù)來源于中國氣象科學(xué)數(shù)據(jù)共享服務(wù)網(wǎng)（http：∥cdc.cma.gov.cn），采用1960—2004年安定區(qū)、華家?guī)X、通渭、隴西、渭源、臨洮、漳縣、岷縣8個站點的逐日氣象數(shù)據(jù)。除漳縣站1961—1966年的部分氣象數(shù)據(jù)缺少外，其余站點的觀測資料都比較齊全。此外，華家?guī)X、臨洮和岷縣為國家基本氣象站，其氣象數(shù)據(jù)持續(xù)至2013年。本文基于定西市45a的逐日氣象資料，分別計算該地區(qū)月平均氣溫、月極端最高氣溫、月極端最低氣溫、月降水量等氣象指標(biāo)，分別從年際尺度、季節(jié)差異以及極端天氣的變化等方面進(jìn)行系統(tǒng)分析，深入探討該地區(qū)多年氣候變化特征。季節(jié)劃分標(biāo)準(zhǔn)為：春季為3—5月，夏季為6—8月，秋季為9—11月，冬季為12—次年2月。

1.2 分析方法

采用線性回歸、5a滑動平均、氣候傾向率［22］、Mann－Kendall檢驗［23］和 R/S分析［24］等數(shù)學(xué)方法，對所搜集到的氣象資料進(jìn)行詳細(xì)分析。其中，線性回歸和氣候傾向率能夠從直觀角度反映定西地區(qū)長時間序列的氣候變化趨勢；5a滑動平均平滑處理了年際間的干擾，且用時間序列的平滑值來顯示氣候變化趨勢。

線性氣候傾向率計算公式為：

其趨勢變化率方程為：

式中系數(shù)a0，a1，可通過最小二乘法或經(jīng)驗正交多項式來確定，本文采用最小二乘法，即滿足式（3）：

Mann－Kendall檢驗法是世界氣象組織推薦并廣泛使用的非參數(shù)檢驗方法，該方法既可以檢測序列的變化趨勢，也可以進(jìn)行突變點檢驗，主要用于氣候要素在時間序列趨勢中的突變性檢測。該方法的優(yōu)點在于不需要樣本遵從一定的分布，不受少數(shù)異常值的干擾，且計算簡便［22］。

Mann－Kendall檢驗中，原假設(shè)H0為原序列（x1，x2，…，xn）無變化趨勢，備擇假設(shè)H1是原序列有增大或減小變化趨勢，為雙邊檢驗，構(gòu)造秩序列Sk。

定義統(tǒng)計變量：

UFk為標(biāo)準(zhǔn)正態(tài)分布，UF1＝0。對于給定顯著性水平α，若，則表明原序列存在明顯的趨勢變化。將上述方法逆向使用，時間序列逆序排列為xn，xn－1，…，x1，重復(fù)相同的操作，得到 UBk，滿足UBk＝－UFk，k＝1，2，…，n，UBn＝0。

通過UFk和UBk的曲線走勢可進(jìn)一步分析序列x的變化趨勢，同時也能夠明確突變的時間及顯著性變化的范圍。R/S分析方法又稱為重標(biāo)極差分析法（Rescaled Range Analysis），最早由被稱為尼羅河之父的H.E.Hurst提出。近年來，R/S方法在時間序列分析中得到廣泛應(yīng)用，其Hurst指數(shù)能很好的揭示時間序列中的趨勢成分［24－25］。

Hurst指數(shù)H（0＜H＜（1），存在以下幾種情況：

① 當(dāng)H＝0.5時，序列是一個隨機(jī)游動序列，反映在氣候指標(biāo)上則是各觀測結(jié)果之間完全獨立，氣候指標(biāo)隨機(jī)變化；

② 當(dāng)0≤H＜0.5時，系統(tǒng)是反持久性的或遍歷性的時間序列，具有比隨機(jī)序列更強(qiáng)的突變性或易變性。反應(yīng)在氣候指標(biāo)上則是氣候變化的總體趨勢與過去相反；

③ 當(dāng)0.5＜H≤1時，系統(tǒng)是一個持久性的或趨勢增強(qiáng)的序列，反應(yīng)今后的氣候變化總趨勢與過去相同。

2 結(jié)果與分析

2.1 氣溫變化

2.1.1 平均氣溫年際變化 1960—2004年定西市年平均氣溫呈遞增趨勢，線性升溫率為0.252℃/10a（圖1），略高于全國平均水平（0.25℃/10a）［26］，其中安定區(qū)、通渭縣、漳縣的年平均氣溫線性升溫率超過市平均升溫率，分別為0.338，0.255，0.272℃/10a，以安定區(qū)升溫最明顯（表1）。定西市年平均氣溫介于5.39～7.74℃之間。1967年是定西市60年代以來最冷的一年，年平均氣溫為5.39℃，比常年平均值低1.0℃；該市氣溫以增溫為主，氣溫距平值呈現(xiàn)遞增趨勢，1990年之后增溫明顯，以正向距平為主（圖1b）。1998年為定西市歷史上最暖年，年平均氣溫為7.74℃，高于常年平均值1.44℃。

圖1 1960－2004年定西市平均氣溫和氣溫距平年際變化

表1 定西市各站點平均氣溫對比

2.1.2 平均氣溫季節(jié)變化 各季平均氣溫的變化趨勢見圖2，如圖所示，定西市四季氣溫都呈現(xiàn)增暖趨勢，但不同季節(jié)的氣溫升高速率不同，其中以冬季增溫最明顯（0.377℃/10a），其次為秋季（0.271℃/10a），夏季（0.181℃/10a），春季的增溫幅度最?。?.156/10a）。20世紀(jì)90年代起各季節(jié)溫度均高于常年均值，20世紀(jì)初的4a變化尤為顯著（圖3），其中春季和冬季氣溫分別增長11.2%和14.6%。1990—1999年冬季平均氣溫為－4.83℃，比常年均值高0.5℃，21世紀(jì)初期進(jìn)一步升高為－4.55℃。2004年春季平均氣溫達(dá)8.74℃，比常年偏高1.4℃，是45a來春季最暖的一年；2002年夏季平均氣溫達(dá)18.34℃，比常年偏高1.46℃；秋季氣溫最高的年份為1998年8.39℃，高于常年均值1.81℃。

這些變化說明，定西市暖冬趨勢明顯，而且冬季增溫程度仍在加強(qiáng)，秋季的增溫幅度也有增大的趨勢，與全球氣候變暖明顯加速的結(jié)論相一致［23］。區(qū)域氣候變化一方面與大氣環(huán)流密切相關(guān)［27］，另一方面也受到該地區(qū)地形條件的影響。定西市地處中緯度地帶，屬典型大陸性氣候，同時深受印度洋、太平洋和歐亞大陸三大環(huán)流影響［28］，這些特點決定了其易受全球環(huán)境變化的影響。

2.1.3 極端氣溫變化特征 45a來該地區(qū)極端氣溫變化曲線如圖4所示，增暖趨勢明顯，且極端最高氣溫的氣候傾向率（0.338℃/10a）略大于極端最低氣溫的氣候傾向率（0.268℃/10a）。

定西市各站點極端氣溫的年際氣候傾向率和各季的氣候傾向率變化情況如表2所示。極端最高氣溫以秋季增溫最明顯（0.6℃/10a），為年際極端最高氣溫氣候傾向率的1.76倍；極端最低氣溫的氣候傾向率以冬季最高為0.5℃/10a，達(dá)年際極端最低氣溫氣候傾向率的1.85倍。極端氣溫的增溫幅度均以安定區(qū)的變化最明顯，且極端最低氣溫增溫幅度高于極端最高氣溫。極端最低氣溫的氣候傾向率除了秋季低于極端最高氣溫的氣候傾向率之外，春季、夏季、冬季均高于極端最高氣溫的氣候傾向率。進(jìn)一步表明定西市氣侯呈現(xiàn)明顯的變暖趨勢，且以秋冬季的增溫為主。

圖2 1960－2004年定西市平均氣溫季節(jié)性變化趨勢

圖3 1960－2004年定西市平均氣溫季節(jié)性變化

2.2 降水量變化特征

2.2.1 降水量年際變化特征 1960—2004年定西市年降水量呈減少趨勢（圖5），氣候傾向率為－22.62 mm/10a。年降水量變幅范圍為293.4～717mm，常年平均降水量為476.3mm，降水量最多的年份是最少年份的2.44倍，分別以1997年的降水量最少，1967年的降水量最高。從20世紀(jì)80年代以來，年降水量普遍減少，其年際間的波動相較于六七十年代更平緩。

圖4 1960－2004年定西市極端氣溫變化趨勢

表2 1960－2004年定西市極端氣溫氣候傾向率

圖5 定西市1960－2004年降水量變化趨勢

為了更清晰地探討降雨量的年際變化，繪制定西市降水量分布圖（圖6）。從圖6可知，降水量主要集中在7月，8月份，占年總降水量的38.6%。年降水量以定西市北部的安定區(qū)最低，東部華家?guī)X的年降水量最高，總體呈現(xiàn)南部高于北部，四周高于中心的變化趨勢。定西市的降水變化很大程度上受到地形條件與季風(fēng)變化的影響。

2.2.2 降水量季節(jié)變化特征 降水量季節(jié)變化趨勢如圖7所示：定西市各季節(jié)的降水量總體呈減少趨勢，其中以秋季的下降幅度最大，氣候傾向率為－15.30 mm/10a。各年代間的變化也不盡相同（表3），春季降水量以80年代最高，超過均值的7.6%；夏季降水量70年代最高，2000年以來明顯減少，比平均降水量低8.1%；秋季、冬季降水量分別以60年代和80年代最高，年平均降水量隨著時間梯度持續(xù)減少，20世紀(jì)初有所回升。

定西市處于東南季風(fēng)邊緣地區(qū)，是季風(fēng)變化的敏感地帶，季風(fēng)的強(qiáng)弱變化直接影響降水的多少。定西市降水量的減少主要發(fā)生在夏秋季節(jié)，近半個世紀(jì)以來東南季風(fēng)一直保持減弱趨勢，與區(qū)域降水量減少趨勢相一致［18，28］。此外，定西市地處印度洋、太平洋和歐亞大陸三大環(huán)流的匯聚位置，一方面受到秦嶺焚風(fēng)效應(yīng)和青藏高原下沉氣流影響，同時還受到全球尺度的大氣環(huán)流影響，這也是影響定西大范圍加水差異的可能性原因之一。

圖6 降水量分布狀況全市月平均以及各站點年降水量

圖7 定西市1960－2004年降水量季節(jié)性變化趨勢

2.3 氣候突變性檢驗

氣候突變泛指氣候從一種狀態(tài)到另一種狀態(tài)的較迅速轉(zhuǎn)變的現(xiàn)象。定西市氣候干旱少雨，氣溫總體呈現(xiàn)增溫趨勢，降水量呈遞減趨勢，其年平均氣溫和年平均降雨量的檢測結(jié)果見圖8。

根據(jù)Mann－Kendall趨勢檢驗原理，UF曲線值大于0，表明序列呈上升趨勢，小于0則表明呈下降趨勢。1987年之前UF曲線呈現(xiàn)波動，未達(dá)到顯著上升（或下降）趨勢，1987年之后呈現(xiàn)上升趨勢，1998—2004年升溫趨勢超過0.05顯著性水平臨界線（μ0.05＝±1.96），甚至超過0.01顯著性水平（μ0.01＝±2.56），表明定西氣溫呈現(xiàn)從冷到暖的顯著上升趨勢。UF曲線與UB曲線僅在1994年出現(xiàn)了一個交點且位于兩臨界線之間，表明定西市年平均氣溫的突變時間從1994年開始，此后，氣溫呈現(xiàn)明顯的遞增趨勢。相比而言，降水量的UF曲線和UB曲線有多處交點，沒有出現(xiàn)顯著的突變現(xiàn)象，降水量總體呈現(xiàn)減少趨勢，且在1997年低于0.05顯著性水平置信線，降水量進(jìn)入短期顯著減少階段。

表3 定西市各年代平均降水量季節(jié)性差異 mm

圖8 平均氣溫和降水量的Mann－Kendall檢測曲線

2.4 R/S分析

通過上述分析可知定西市氣溫具有上升趨勢而降水量呈下降趨勢。在此基礎(chǔ)上，對定西市年平均氣溫和年平均降水量分別進(jìn)行R/S分析（見圖9），得出定西市各站點氣候指標(biāo)的Hurst指數(shù)（見表4）。結(jié)果表明，定西市氣溫指標(biāo)中，各站點H值都大于0.5，存在明顯的Hurst現(xiàn)象，說明定西市氣溫的變化將延續(xù)過去的趨勢，持續(xù)上升。除了安定區(qū)、臨洮、漳縣降水指標(biāo)的H值略大于0.5，其余站點的H值均小于0.5，說明降水變化具有一定的波動性，有呈現(xiàn)反持續(xù)性趨勢，即降水量可能會有所增加。定西市氣候變化將可能由暖干型向暖濕型轉(zhuǎn)變。此外，華家?guī)X、臨洮、岷縣這三個站點2004—2013年的氣溫降水變化與R/S分析結(jié)果相一致。

圖9 定西市氣溫降水的R/S分析

當(dāng)前，對西北干旱半干旱地區(qū)氣候變化的相關(guān)結(jié)論并不一致，其中，Zhang等［4］指出暖干向暖濕型轉(zhuǎn)變是西北地區(qū)氣候變化的主要趨勢，而姚曉軍等［29］則表明，西北大部分地區(qū)干旱化程度增加。這主要與本研究的結(jié)果與前者的研究結(jié)果一致，定西市將呈現(xiàn)由暖干型向暖濕型轉(zhuǎn)變的趨勢。

表4 定西市各站點氣候指標(biāo)的Hurst指數(shù)

3 結(jié) 論

氣溫和降水是反映氣候變化的兩大要素，長時間序列分析能準(zhǔn)確呈現(xiàn)區(qū)域氣候的變化趨勢和規(guī)律［9，30］。定西市自1960—2004年的氣候變化具有以下特征：

（1）年平均氣溫整體呈上升趨勢，檢驗結(jié)果表明1994年是氣溫變化的突變年，之后進(jìn)入顯著增溫期。不同季節(jié)的平均氣溫均呈現(xiàn)增加趨勢，尤以冬季增溫最明顯（0.377℃/10a），春季的增溫幅度最?。?.156℃/10a）。

（2）極端最高氣溫和極端最低氣溫上升趨勢明顯，極端最高氣溫的氣候傾向率（0.338℃/10a）大于極端最低氣溫（0.268℃/10a）。從季節(jié)角度看，秋季和冬季的極端最高氣溫增溫趨勢大于年際間的增溫趨勢；極端最低氣溫以冬季增溫最顯著。

（3）年降水量呈下降趨勢，常年平均降水量為476.3mm，不存在明顯的突變現(xiàn)象。降水量總體呈現(xiàn)南部高于北部，四周高于中心的變化趨勢。降水量的減少主要發(fā)生在夏秋季，且秋季的下降幅度最大（－15.30mm/10a）。

（4）對定西市1960—2004年氣溫、降水的Mann－Kendall檢驗與R/S分析進(jìn)行綜合分析，結(jié)果表明，定西市平均氣溫呈現(xiàn)增長趨勢，且存在明顯的Hurst現(xiàn)象，即定西市氣候變化存在持續(xù)性，未來氣溫將持續(xù)上升；降水變化具有一定的波動性，有呈現(xiàn)反持續(xù)性趨勢，即未來降水量將有所增加。

［1］ 王紹武，謝志輝，蔡靜寧，等.近千年全球平均氣溫變化的研究［J］.自然科學(xué)進(jìn)展，2002，12（11）：27－31.

［2］ 錢維宏，陸波.千年全球氣溫中的周期性變化及其成因［J］.科學(xué)通報，2010，55（32）：3116－3121.

［3］ Huang J P，Guan X，Ji F.Enhanced cold－season warming in semi－arid regions［J］.Atmospheric Chemistry and Physics，2012，12（12）：5391－5398.

［4］ ZhangBaoqing，Wu Pute，Zhao Xining，et al.Drought variation trends in different subregions of the Chinese Loess Plateau over the past four decades［J］.Agricultural Water Management，2012，115：167－177.

［5］ 黃建平，季明霞，劉玉芝，等.干旱半干旱區(qū)氣候變化研究綜述［J］.氣候變化研究進(jìn)展，2013，09（1）：9－14.

［6］ 李新周，劉曉東，馬柱國.近百年來全球主要干旱區(qū)的干旱化特征分析［J］.干旱區(qū)研究，2004，21（2）：97－103.

［7］ Walther G R，Post E，Convey P，et al.Ecological responses to recent climate change［J］.Nature，2002，416（6879）：389－395.

［8］ 郭敏杰，張亭亭，張建軍，等.1982—2006年黃土高原地區(qū)植被覆蓋度對氣候變化的響應(yīng)［J］.水土保持研究，2014，21（5）：35－40.

［9］ 張曉曉，張鈺，徐浩杰.1961—2010年白龍江流域氣溫和降水量變化特征研究［J］.水土保持研究，2014，21（4）：238－245.

［10］ 呂妍，王讓會，蔡子穎.我國干旱半干旱地區(qū)氣候變化及其影響［J］.干旱區(qū)資源與環(huán)境，2009，23（11）：65－71.

［11］ 劉德坤，王軍邦，齊述華.基于濕潤指數(shù)的近35a青海省干濕狀況變化分析［J］.水土保持研究，2014，21（2）：246－250.

［12］ 卓瑪蘭草，劉普幸，張亞寧，等.甘肅黃土高原區(qū)潛在蒸散量時空變化與成因研究［J］.水土保持研究，2012，19（1）：70－75.

［13］ 賈建英，郭建平.東北地區(qū)近46a氣候變化特征分析［J］.干旱區(qū)資源與環(huán)境，2011，25（10）：109－115.

［14］ 高峰，隋波，孫鴻雁，等.1951—2008年東北地區(qū)冬季氣溫變化及環(huán)流場特征［J］.氣象與環(huán)境學(xué)報，2011，27（4）：12－16.

［15］ 吉奇，徐璐璐，宋冀鳳.1953—2005年本溪地區(qū)氣候變化及其對農(nóng)業(yè)生產(chǎn)的影響［J］.氣象與環(huán)境學(xué)報，2008，24（1）：36－39.

［16］ CookB I，Wolkovich E M，and Parmesan C.Divergent responses to spring and winter warming drive community level flowering trends［J］.Proceedings of the National Academy of Sciences，2012，109（23）：9000－9005.

［17］ Yu H，Luedeling E，and Xu J.Winter and spring warming result in delayed spring phenology on the Tibetan Plateau［J］.Proceedings of the National Academy of Sciences，2010，107（51）：22151－22156.

［18］ 張秀云，董覲華，姚玉璧，等.定西市氣候資源特點與開發(fā)利用［J］.成都信息工程學(xué)院學(xué)報，2004，19（3）：447－450.

［19］ Zhang Zhihua，Li Xiaoyan，Jiang Zhiyun，et al.Changes in some soil properties induced by re－conversion of cropland into grassland in the semiarid steppe zone of Inner Mongolia，China［J］.Plant and Soil，2013，373（1）：89－106.

［20］ Li Xiaoyan.Hydrology and biogeochemistry of semiarid and arid regions［J］.Forest Hydrology and Biogeochemistry，2011，216：285－299.

［21］ 周長進(jìn)，董鎖成，李岱.定西地區(qū)主要生態(tài)環(huán)境問題與對策［J］.干旱區(qū)資源與環(huán)境，2004，18（6）：32－37.

［22］ 林學(xué)椿，于淑秋.近40a我國氣候趨勢［J］.氣象，1990，16（10）：16－22.

［23］ 曲學(xué)斌，吳昊.呼倫貝爾市53a氣候變化特征分析［J］.水土保持研究，2014，21（1）：178－182.

［24］ 顧維慧，吳洪軍，喬曉波.R/S分析法在漣水縣降水趨勢預(yù)測中的應(yīng)用［J］.農(nóng)業(yè)與技術(shù)，2013，33（10）：204－206.

［25］ 何毅，楊太保，杜鵑，等.基于GIS和RS的中亞阿拉套山脈近22a來冰川變化［J］.水土保持研究，2013，20（6）：130－134.

［26］ 何云玲，魯枝海.近60a昆明市氣候變化特征分析［J］.地理科學(xué)，2012，32（9）：1119－1124.

［27］ 李雙雙，延軍平，萬佳.全球氣候變化下秦嶺南北氣溫變化特征［J］.地理科學(xué)，2012，32（7）：853－858.

［28］ 王文艷，張麗萍，吳東平，等.近19a干旱半干旱過渡區(qū)氣溫降雨變化特征分析［J］.環(huán)境科學(xué)與技術(shù)，2012，35（61）：57－61.

［29］ 姚曉軍，張曉，孫美平，等.1960—2010年中國西北地區(qū)水分盈虧量時空特征［J］.地理研究，2013，32（4）：607－616.

［30］ 周曉宇，趙春雨，張新宜，等.1961—2009年遼寧省氣溫、降水變化特征及突變分析［J］.干旱區(qū)資源與環(huán)境，2013，27（10）：87－93.