王永財(cái) ，孫艷玲 ，張 靜 ，王中良 ，1b

（1.天津師范大學(xué)a.城市與環(huán)境科學(xué)學(xué)院，b.天津市水資源與水環(huán)境重點(diǎn)實(shí)驗(yàn)室，天津300387；2.中國(guó)科學(xué)院植物研究所，北京100093）

在全球變暖的背景下，氣候變化受到國(guó)際普遍關(guān)注.IPCC第4次會(huì)議指出過去50年全球變暖趨勢(shì)為平均每10a上升1.3℃（1.0～1.6℃/10a），幾乎是過去100a的2倍，全球大氣平均溫度和海洋平均溫度均在增加[1].由于氣候變化對(duì)農(nóng)業(yè)、水資源和自然生態(tài)系統(tǒng)等均有很大影響，因此許多學(xué)者對(duì)區(qū)域氣候變化進(jìn)行了研究和探討.如楊特群等[2]統(tǒng)計(jì)分析了黃河流域1951年—2007年氣溫降水變化情況，指出黃河流域氣溫在1950—1970年代前比較穩(wěn)定，1980年代后各區(qū)間均呈上升趨勢(shì)，黃河流域降水在各區(qū)間變化特點(diǎn)不盡相同.吉奇等[3]對(duì)東北1951年—2000年氣溫和降水研究表明東北地區(qū)年平均氣溫在以0.6℃/10a的速率上升，降水量以0.27mm/a的速率下降.任朝霞等[4]對(duì)西北干旱區(qū)1951年—2000年的氣溫和降水情況進(jìn)行了研究，結(jié)果表明西北干旱區(qū)氣溫上升趨勢(shì)為0.22℃/10a，上升趨勢(shì)南北高、中間低，降水量增加速率為3.2mm/10a，并呈現(xiàn)出由東南向西北遞增的格局.祁威等[5]對(duì)珠穆朗瑪峰地區(qū)尼泊爾境內(nèi)1971年—2009年的氣候變化情況進(jìn)行研究，指出1971年—2009年，珠穆朗瑪峰南坡線性升溫速率為0.25℃/10a，年平均降水量以4.27mm/a的速率增加，但增加并不顯著，南坡降水量高于北坡，北坡增溫趨勢(shì)較南坡明顯.

作為我國(guó)重要的工農(nóng)業(yè)生產(chǎn)地，海河流域氣候趨于暖干，持續(xù)干旱事件經(jīng)常發(fā)生，給當(dāng)?shù)亟?jīng)濟(jì)和社會(huì)帶來了十分嚴(yán)重的影響[6].海河流域氣候的變化特征和氣候變化的響應(yīng)已成為海河流域氣候變化研究的主要內(nèi)容.褚健婷等[7]利用K-S檢驗(yàn)等方法分析了海河流域1995年—2004年的降水時(shí)空變化特征，結(jié)果表明海河流域以太行山和燕山為界，分為山前多雨帶、山前平原區(qū)和山后少雨帶3部分.郝春灃等[8]以海河流域1956年—2005年的氣溫和降水作為研究對(duì)象，對(duì)海河流域氣溫和降水的變化特征和規(guī)律進(jìn)行了分析，指出海河流域近50年氣溫呈顯著上升趨勢(shì)，而降水呈下降趨勢(shì).據(jù)王利娜等[9]對(duì)海河流域1950年—2010年降水的時(shí)空變化特征分析表明：時(shí)間上，海河流域總體降水呈減少趨勢(shì)；空間上，整個(gè)流域年降水量由南向西減少.

綜上所述，目前對(duì)海河流域氣候變化的研究成果較多，但主要集中在該區(qū)域氣溫和降水單一方面的時(shí)間或空間變化特征，沒有從時(shí)間和空間2個(gè)方面研究氣溫和降水的變化趨勢(shì)特征、突變特征和周期性特征，不利于全面綜合地把握海河流域氣候的時(shí)空變化特征和變化規(guī)律.已有研究將由氣溫和降水計(jì)算得到的相對(duì)濕潤(rùn)指數(shù)作為干旱指標(biāo)，對(duì)海河流域1960年—2010年作物生長(zhǎng)季干旱的時(shí)空特征進(jìn)行分析，取得了很好的結(jié)果[10].因此，本研究以氣溫和降水量2個(gè)氣候表征參數(shù)作為研究對(duì)象，利用海河流域1960年—2010年的氣候資料，采用線性傾向估計(jì)、Mann-Kendall非參數(shù)檢驗(yàn)、滑動(dòng)T檢驗(yàn)和小波分析等方法，分析海河流域氣溫和降水變化的趨勢(shì)性、突變性和周期性等特征，以期更全面地揭示其氣候變化特征和事實(shí).

1 研究材料與方法

本研究所用資料來源于中國(guó)氣象科學(xué)數(shù)據(jù)共享服務(wù)網(wǎng)，選取分布較均勻、遍布整個(gè)海河流域的30個(gè)具有完整性、可靠性和代表性的站點(diǎn)，取得每個(gè)站點(diǎn)1960年—2010年的月平均氣溫和月降水資料，統(tǒng)計(jì)得到年平均氣溫和年降水量.為了獲得可以代表全流域的均一氣溫和降水資料，利用反距離加權(quán)平均法和泰森多邊形法將氣溫和降水資料展布到全流域，得到近51年的氣溫和降水量序列.

本研究中對(duì)海河流域氣候的線性傾向估計(jì)以海河流域年平均氣溫和降水量的距平作為統(tǒng)計(jì)量，采用氣候傾向率和5年滑動(dòng)平均分析其氣候變化趨勢(shì)特征，計(jì)算結(jié)果用Mann-Kendall非參數(shù)檢驗(yàn).其中，Mann-Kendall非參數(shù)檢驗(yàn)可以很好地揭示時(shí)間序列的趨勢(shì)變化，是目前比較常用的分析方法[11]，在長(zhǎng)時(shí)間序列的數(shù)據(jù)趨勢(shì)檢驗(yàn)中應(yīng)用廣泛，方法如下：

式（1）和式（2）中：Z 為標(biāo)準(zhǔn)化檢驗(yàn)統(tǒng)計(jì)量；xj和 xi為序列數(shù)據(jù)；n是序列數(shù)據(jù)個(gè)數(shù)；S為檢驗(yàn)統(tǒng)計(jì)量，當(dāng)n≥8時(shí)，對(duì)S進(jìn)行近似正態(tài)分布，其均值和方差為

海河的氣候變化周期性采用小波分析，小波分析是在Fourier分析的基礎(chǔ)上發(fā)展起來的一種新的時(shí)頻局部化分析方法，在時(shí)間序列研究中，小波分析主要用于時(shí)間序列的消噪和濾波、信息量系數(shù)和分形維數(shù)的計(jì)算、突變點(diǎn)的監(jiān)測(cè)、周期成分的識(shí)別以及多時(shí)間尺度的分析等[15].本研究采用墨西哥小帽小波（Mexican hat function）[16]對(duì)海河流域1960年—2010年的年平均氣溫和年降水量時(shí)間序列進(jìn)行分析，得到海河流域在不同時(shí)間尺度上的氣溫和降水的變化特征，利用小波方差確定主要的振蕩周期.

2 結(jié)果與分析

2.1 氣溫變化特征分析

2.1.1 氣溫變化趨勢(shì)分析

海河流域1960年—2010年年平均氣溫變化趨勢(shì)如圖1所示.

圖1 1960年—2010年海河流域年平均氣溫距平變化Fig.1 Annual variation of air temperature anomalies in Haihe River Basin from 1960to 2010

由圖1可以看出，海河流域年平均氣溫呈逐年升高的趨勢(shì)，近51a來，年平均最低氣溫出現(xiàn)在1969年，為8.7℃；最高年平均氣溫出現(xiàn)在1998年，為11.3℃.1960年—2010年年平均氣溫的增溫速率為0.34℃/10a，與郝椿灃等研究的海河流域1956年—2005年0.30℃/10a的年平均氣溫增溫速率結(jié)果相似[8]，但明顯高于文獻(xiàn)[17]中全國(guó)1951年—2002年0.22℃/10a的年平均氣溫的增溫速率[17].通過Mann-Kendall非參數(shù)檢驗(yàn)，1960年—2010年海河流域年平均氣溫增溫趨勢(shì)通過了0.01的顯著性檢驗(yàn)，說明海河流域近51年增溫趨勢(shì)非常明顯.從5年滑動(dòng)平均來看，海河流域近51年氣溫變化波動(dòng)上升，60年代年平均氣溫呈下降趨勢(shì)，下降速率達(dá)到1.1℃/10a；70年代開始?xì)鉁爻什▌?dòng)上升；80年代初期氣溫變化比較平穩(wěn)，后期開始上升；90年代氣溫顯著升高，升高速率高達(dá)1.4℃/10a；2000年以來，海河流域增溫趨勢(shì)減緩.

海河流域氣溫的空間變化情況如圖2所示.由圖2可以看出，海河流域以五臺(tái)山為中心的西部地區(qū)增溫趨勢(shì)最為明顯，中部地區(qū)增溫趨勢(shì)較為明顯，而南部和東北部增溫趨勢(shì)較弱.經(jīng)Mann-Kendall非參數(shù)檢驗(yàn)可知，海河流域30個(gè)站點(diǎn)中，僅有榆社和承德2個(gè)站點(diǎn)未通過α=0.05的顯著性檢驗(yàn)，其他站點(diǎn)均通過α=0.05的顯著性檢驗(yàn)，占所有站點(diǎn)的93%，說明海河流域氣溫上升趨勢(shì)的顯著性.

圖2 海河流域氣溫空間變化Fig.2 Spatial variation of temperature in Haihe River Basin

2.1.2 氣溫突變特征分析

根據(jù)Mann-Kendall突變檢驗(yàn)結(jié)果分析（圖3），海河流域1960年—2010年年平均氣溫正序列UF在1970年前下降，之后開始波動(dòng)上升，90年代超過了0.05的顯著性水平，說明上升趨勢(shì)顯著.正序列UF和逆序列UB在0.05的顯著性水平臨界范圍內(nèi)有1個(gè)交點(diǎn)，因此可能的突變點(diǎn)為1988年.為了獲得正確的突變點(diǎn)和突變發(fā)生的時(shí)間區(qū)域，對(duì)上述Mann-Kendall突變檢驗(yàn)所得的突變點(diǎn)進(jìn)行滑動(dòng)T檢驗(yàn)，結(jié)果表明1988年的突變超過了0.05的顯著性水平.因此，1988年是海河流域年平均氣溫的突變點(diǎn)，是氣溫迅速升高的開始.

圖3 海河流域年平均氣溫突變的Mann-Kendall檢驗(yàn)Fig.3 Abrupt change of annual average temperature in Haihe River Basin using Mann-Kendall test

對(duì)海河流域30個(gè)站點(diǎn)進(jìn)行Mann-Kendall突變檢驗(yàn)，結(jié)果如圖4所示.

圖4 海河流域氣溫突變檢驗(yàn)Fig.4 Abrupt change test of temperature in Haihe River Basin

由圖4可知，1960年—2010年海河流域30個(gè)站點(diǎn)有18個(gè)年平均氣溫發(fā)生了突變，突變發(fā)生的年份多在1972年—1974年、1986年—1989年和1990年—1994年這幾個(gè)時(shí)間段內(nèi)；東部地區(qū)大部分站點(diǎn)氣溫發(fā)生突變，且突變點(diǎn)集中在1986年—1989年時(shí)間區(qū)域內(nèi).

2.1.3 氣溫變化周期性分析

從海河流域年平均氣溫小波圖（圖5）分析，海河流域近51年的年平均氣溫存在不同時(shí)間尺度的震蕩，期間經(jīng)過了低溫、高溫、相對(duì)低溫和相對(duì)高溫共4個(gè)階段，存在7a左右和10～15a左右2個(gè)周期，其中7a左右是主周期，因此可以預(yù)測(cè)2010年后，海河流域?qū)⒃?010年以后以7a左右或10～15a為尺度進(jìn)入偏冷期.

圖5 海河流域年平均氣溫小波分析Fig.5 Wavelets analysis of annual mean temperature in Haihe River Basin

2.2 降水變化特征分析

2.2.1 降水變化趨勢(shì)分析

海河流域年降水量變化情況如圖6所示.

圖6 海河流域年降水量距平變化Fig.6 Annual variation of air precipitation anomalies in Haihe River Basin from 1960to 2010

由圖6可以看出，近51年海河流域年降水量呈下降趨勢(shì)，年降水量最高的年份是1964年，為859.6mm；年降水量最小的是1997年，年降水量為370.9mm.從海河流域多年降水變化趨勢(shì)來看，年降水量緩慢波動(dòng)下降，1960年—2010年的下降速率為18.9mm/10a，這與任國(guó)玉等[16]研究的中國(guó)1956年—2002年全國(guó)年平均降水量呈上升趨勢(shì)相反，也低于郝椿灃[8]和王曉霞[18]等分析海河流域年降水量的下降速率22.0mm/10a和27.2mm/10a.經(jīng)Mann-Kendall非參數(shù)檢驗(yàn)，海河流域年降水量未通過0.05的顯著性檢驗(yàn)，進(jìn)一步說明海河流域降水量呈下降趨勢(shì)，但下降趨勢(shì)并不顯著.海河流域降水量波動(dòng)較大，60年代降水量波動(dòng)下降，但下降較緩，下降速率為7.9mm/10a，70年代和90年代降水量下降較快，而80年代和2000年以來降水量變化較緩.

海河流域降水量的空間變化情況如圖7所示.由圖7可以看出，海河流域內(nèi)天津、唐山等東部地區(qū)和西部山西地區(qū)降水量下降趨勢(shì)明顯，而北部和南部地區(qū)降水量下降趨勢(shì)較弱.經(jīng)Mann-Kendall非參數(shù)檢驗(yàn)，五臺(tái)山、榆社和遵化3個(gè)站年降水量通過了0.05的顯著性檢驗(yàn)，其他27個(gè)站點(diǎn)未通過0.05的顯著性檢驗(yàn)，說明海河流域降水量呈下降趨勢(shì)，但下降趨勢(shì)相對(duì)不明顯.

圖7 海河流域降水量的空間變化Fig.7 Spatial variation of precipitation in Haihe River Basin

2.2.2 降水突變特征分析

Mann-Kendall檢驗(yàn)結(jié)果如圖8所示.

圖8 海河流域年降水Mann-Kendall檢驗(yàn)Fig.8 Abrupt change of annual average precipitation in Haihe River Basin using Mann-Kendall test

由圖8可知，1960年—2010年海河流域年降水量的正序列UF在1965年之前上升，后波動(dòng)下降，1997年后顯著下降.正序列UF和逆序列UB在0.05的臨界范圍內(nèi)有一個(gè)交點(diǎn)，所以1997年可能為突變點(diǎn)，經(jīng)滑動(dòng)T檢驗(yàn)，結(jié)果超過了0.05的顯著性水平，因此1997年是海河流域年降水量的突變點(diǎn)，后降水量進(jìn)入相對(duì)較少期.

對(duì)海河流域30個(gè)站點(diǎn)的降水量進(jìn)行Mann-Kendall突變檢驗(yàn)，結(jié)果如圖9所示.

圖9 海河流域降水的突變檢驗(yàn)Fig.9 Abrupt change test in Haihe River Basin

由圖9可知，1960年—2010年海河流域30個(gè)站點(diǎn)中有7個(gè)站點(diǎn)年降水量發(fā)生了突變，發(fā)生突變時(shí)間主要在1980年—1982年和2002年—2005年2個(gè)時(shí)間區(qū)域內(nèi)，海河流域東北部的青龍、承德和圍場(chǎng)年均降水量均發(fā)生了突變，且發(fā)生突變的時(shí)間為1980年—1982年.

2.2.3 降水變化周期性分析

海河流域年降水量小波分析結(jié)果如圖10所示.由圖10可以看出，近50年海河流域降水量變化經(jīng)過了偏少、偏多、相對(duì)偏少和相對(duì)偏多共4個(gè)震蕩階段，存在5～7a和10～15a 2種周期，其中5～7a為主周期，即海河流域在2010年以后可能以5～7a或10～15a為尺度進(jìn)入降水偏豐期.

圖10 海河流域年降水小波分析Fig.10 Wavelets analysis of annual mean temperature in Haihe River Basin

3 結(jié)論

以氣候的表征參數(shù)氣溫和降水作為研究對(duì)象，從趨勢(shì)性、突變性和周期性等角度對(duì)海河流域1960年—2010年的氣溫和降水特征進(jìn)行分析，結(jié)果表明：

（1）從變化趨勢(shì)來看：近51年，海河流域氣溫顯著升高，降水呈下降趨勢(shì)，年平均氣溫增溫速率為0.34℃/10a，降水量的下降速率為18.9mm/10a；海河流域西部氣溫上升趨勢(shì)強(qiáng)于南部和北部，東西部降水量下降趨勢(shì)較強(qiáng)，而南北部較弱.

（2）從突變性來看：海河流域年平均氣溫突變點(diǎn)在1988年，降水量突變點(diǎn)在1997年，海河流域30個(gè)站點(diǎn)中有18個(gè)發(fā)生氣溫突變，且突變點(diǎn)集中在1986年—1989年和1990年—1994年2個(gè)時(shí)間區(qū)域內(nèi)；有7個(gè)站點(diǎn)的降水量發(fā)生了突變，突變點(diǎn)集中在1980年—1982年這個(gè)時(shí)間區(qū)域內(nèi).

（3）從周期性來看：海河流域年平均氣溫變化情況包含7a左右和10～15a左右2個(gè)周期，其中5年左右是主周期；降水量變化情況包含5～7a和10～15a 2個(gè)周期，其中5～7a是主周期.

氣候變化極具復(fù)雜性，本研究以氣溫和降水作為研究對(duì)象，從氣候變化的趨勢(shì)性、突變性和周期性等角度分析氣候變化的特征，只能從中短時(shí)期內(nèi)氣候變化的基本特征和規(guī)律，對(duì)短時(shí)期的氣候變化作出預(yù)測(cè)，無法揭示氣候變化的動(dòng)因并對(duì)長(zhǎng)期氣候變化作出預(yù)測(cè)，因此，氣候變化需要更多探索性的研究.

[1]IPCC.IPCCFourth Assessment Report（AR4）[M].Cambridge：Cambridge University Press，2007.

[2]楊特群，饒素秋，陳冬伶.1951年以來黃河流域氣溫和降水變化特點(diǎn)分析[J].人民黃河，2009，31（10）：76-77.

[3] 吉奇，宋冀鳳，劉輝.近50年東北地區(qū)溫度降水變化特征分析[J].氣象與環(huán)境學(xué)報(bào)，2006，22（5）：1-5.

[4] 任朝霞，楊達(dá)源.西北干旱區(qū)近50年氣候變化特征與趨勢(shì)[J].地球科學(xué)與環(huán)境學(xué)報(bào)，2007，29（1）：99-102.

[5] 祁威，張鐿鋰，高俊剛，等.1971-2009年珠穆朗瑪峰地區(qū)尼泊爾境內(nèi)氣候變化[J].地理學(xué)報(bào)，2013，68（1）：82-94.

[6] 高歌，李維京，張強(qiáng)，等.華北地區(qū)氣候變化對(duì)水資源的影響及2003年水資源預(yù)評(píng)估[J].氣象，2003，29（8）：26-30.

[7] 褚健婷，夏軍，許崇育，等.海河流域氣象和水文降水資料對(duì)比分析及時(shí)空變異[J].地理學(xué)報(bào)，2009，64（9）：1083-1089.

[8] 郝椿灃，賈仰文，龔家國(guó)，等.海河流域近50年氣候變化特征及規(guī)律分析[J].中國(guó)水利水電科學(xué)研究院學(xué)報(bào)，2010，8（1）：39-43.

[9] 王利娜，朱厚華，魯帆，等.海河流域近50年降水量時(shí)空變化特征分析[J].干旱地區(qū)農(nóng)業(yè)研究，2012，30（2）：242-246.

[10]孫艷玲，王中良，王永財(cái).近51年海河流域作物生長(zhǎng)季干旱時(shí)空分布特征研究[J].天津師范大學(xué)學(xué)報(bào)：自然科學(xué)版，2013，33（1）：59-63.

[11]于延勝，陳興偉.R/S和Mann-Kendall法綜合分析水文時(shí)間序列未來的趨勢(shì)特征[J].水資源與水工程學(xué)報(bào)，2008，19（3）：41-44.

[12]黃森旺，李曉松，吳炳方，等.近25年三北防護(hù)林工程區(qū)土地退化及驅(qū)動(dòng)力分析[J].地理學(xué)報(bào)，2012，67（5）：589-598.

[13]魏鳳英.現(xiàn)代氣候統(tǒng)計(jì)診斷與預(yù)測(cè)技術(shù)[M].北京：氣象出版社，1999：63-72.

[14]何云玲，魯枝海.近60年昆明市氣候變化特征分析[J].地理科學(xué)，2012，32（9）：1119-1124.

[15]趙紅怡，武夢(mèng)龍，曹淑琴.小波分析在突變信號(hào)檢測(cè)中應(yīng)用[J].北方工業(yè)大學(xué)學(xué)報(bào)，2004，16（1）：21-24.

[16]徐建華.現(xiàn)代地理學(xué)中的數(shù)學(xué)方法[M].北京：高等教育出版社，2002：417-425.

[17]任國(guó)玉，郭軍，徐銘志，等.近50年中國(guó)地面氣候變化基本特征[J].氣象學(xué)報(bào)，2005，63（6）：942-956.

[18]王曉霞，徐宗學(xué)，紀(jì)一鳴，等.海河流域降水量長(zhǎng)期變化趨勢(shì)的時(shí)空分布特征[J].水利規(guī)劃與設(shè)計(jì)，2010，1：35-38.