京津冀地區(qū)近55年氣候演變特征分析

苗正偉 徐利崗 韓會玲

摘要：基于京津冀地區(qū)22個氣象站 1961-2015年的逐月氣候資料，采用氣候傾向率、MannKendall突變檢測、Morlet小波分析、反距離權重插值（IDW）等方法，研究了近55 a來京津冀地區(qū)氣象要素的時空演變規(guī)律。結果表明：（1）京津冀地區(qū)整體呈現(xiàn)出降水減少、氣溫上升、日照減少、相對濕度降低、風速下降的趨勢。（2）各氣候要素都在1979-1989年間發(fā)生了突變，年降水在1996年發(fā)生了第二次突變。（3）各氣候要素均表現(xiàn)出多時間尺度耦合的復雜特性，其第一主周期尺度介于8～25 a之間。（4）年均溫由南向北遞減，中、西部為升溫速率的高值區(qū)；降水及其變率均自東向西遞減；日照時數(shù)由南向北遞增，而其變率則與之相反；相對濕度由東南向西北遞減，其變率則由西南向東北遞減；平均風速及其變率的高值中心均為張北和塘沽。

關鍵詞：氣候變化；氣溫；降水；日照時數(shù)；相對濕度；風速；京津冀

中圖分類號：P467文獻標志碼：A文章編號：16721683（2018）03012510

Climate evolution characteristics of BeijingTianjinHebei Region for the last 55 years

MIAO Zhengwei1，XU Ligang2，HAN Huilingi1

（1.Hebei University of Water Resources and Electric Engineering，Cangzhou 061001，China；

2.Ningxia Institute of Water Resources Research，Yinchuan 750021，China）

Abstract：Based on the climatic data of 22 meteorological stations in BeijingTianjinHebei region from 1961 to 2015，by means of the climate tendency rate，MannKendall mutation test，Morlet wavelet analysis，and inverse distance weighting （IDW），we systematically studied the temporalspatial variations of climatic variables over the BeijingTianjinHebei region in the recent 55 years.The results showed that：1） In the past 55 years，the BeijingTianjinHebei region underwent an increase of temperature and a decrease of sunshine，precipitation，wind speed，and average relative humidity.2） The climatic variables all went through a sudden change in 19791989，and the annual precipitation went through a second sudden change in 1996.3） The climatic variables all showed a complex characteristic of coupling at multiple time scales，and their first main period scales ranged from 8 to 25 years.4） In BeijingTianjinHebei region，the average air temperature showed a decreasing trend from south to north，with the change rate peaking in the central and western area.Both precipitation and its change rate showed a decreasing trend from east to west.Sunshine hours showed an increasing trend from south to north，while its change rate showed the opposite trend.Relative humidity showed a decreasing trend from southeast to northwest，while its change rate showed a decreasing trend from southwest to northeast.The peaks of wind speed and its change rate both appeared in Zhangbei and Tanggu area.

Key words：climate change；air temperature；precipitation；sunshine hours；relative humidity；wind speed；BeijingTianjinHebei

在全球升溫的背景下，氣候變化具有顯著的區(qū)域差異，認識區(qū)域氣候變化特征，對于區(qū)域農業(yè)生產(chǎn)、防災減災等具有重要意義。京津冀地區(qū)歷來氣象災害多發(fā)，對人們的生產(chǎn)生活造成重大影響，甚至制約了社會、經(jīng)濟的可持續(xù)發(fā)展。在京津冀協(xié)同發(fā)展的國家戰(zhàn)略之下，更有必要對其氣候變化特征進行探討。近幾年，對于京津冀地區(qū)的氣候變化已經(jīng)有了一些研究。許穎等[1]基于懷來縣1954-2012年的月氣象數(shù)據(jù)研究了氣溫、降水的變化特征及其與水庫來水、糧食單產(chǎn)的關系。李鵬飛等[2]研究了京津冀地區(qū)1960-2010年氣溫、降水及潛在蒸散量的時空變化特征。劉芳圓等[3]分析了京津冀地區(qū)1956-2007年氣溫、降水的變化趨勢及其對農業(yè)的影響。張國華等[4]利用5月-8月的日最高氣溫數(shù)據(jù)研究了京津冀地區(qū)城市高溫的氣候特征及城市化效應。段麗瑤等[5]分析了天津1921-2010年降水和氣溫的多尺度特征。胡保昆等[6]研究了北京城區(qū)的氣溫變化特征。向亮等[7]研究了1961-2011河北省降水的演變情況。總的來看，這些研究都取得了一定的成果，但也存在一些局限：或者沒有探討氣候要素的周期性或突變性，或者在所分析氣候要素上各有側重，主要集中于對氣溫、降水的研究，或者資料不是最新，或者沒有對空間特征進行分析，等等，因此，本文對京津冀地區(qū)的降水、氣溫、日照時數(shù)、濕度、風速的趨勢性、突變性、周期性及空間分布特征進行研究，旨在進一步探究京津冀地區(qū)氣候變化特征，以期為合理利用氣候資源、區(qū)域產(chǎn)業(yè)結構調整、生態(tài)環(huán)境保護等提供理論參考。

第16卷 總第96期·南水北調與水利科技·2018年6月苗正偉等·京津冀地區(qū)近55年氣候演變特征分析1材料與方法

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

選取北京、天津、河北三地作為研究對象，從中國氣象科學數(shù)據(jù)共享服務網(wǎng)獲取京津冀地區(qū)連續(xù)性較好的22個站的1961-2015年的逐月觀測資料，包括氣溫、降水、日照時數(shù)、相對濕度、平均風速5個氣候要素，站點分布見圖1，個別缺測數(shù)據(jù)利用鄰近站點線性回歸方法進行插補。

1.2研究方法

將逐月觀測資料整理成各站年均值，基于ArcGIS102中的泰森多邊形法計算氣象要素的區(qū)域均值，用氣候傾向率法[8]估算各氣候要素的長期變化趨勢；MannKendall突變檢測與滑動t檢驗、累積距平相結合進行突變分析，原則是至少通過兩種方法檢驗的才可能被認定為突變點，三種方法的原理及具體步驟參見文獻[9]。

MannKendall突變檢測中，UF曲線為時間正序統(tǒng)計量曲線，其值小于0表示序列呈減?。ㄉ伲┶厔荩笥?表示序列呈增大（多）趨勢，當超過005的顯著性水平（臨界線為±196，下稱信度線）時，表示變化趨勢顯著，UB曲線為時間逆序統(tǒng)計量曲線，UF、UB曲線在信度線內的交點即為可能的突變點。

滑動t檢驗是通過比較兩組相鄰子樣本均值差異是否顯著來確定突變點，兩子樣本長度的選擇具有很大主觀性，但一般二者取相同長度，本文取兩子序列等長為10，且認為P≤005時即為通過了假設檢驗。累積距平曲線既可以用來描述序列的變化趨勢，還可以根據(jù)其升降轉折判斷大概的突變點。

小波分析是進行水文氣象要素周期性研究時的常用方法，且復數(shù)小波比實數(shù)小波效果更優(yōu)[10]。復數(shù)小波系數(shù)的實部大小可表征氣候要素的高低多少，正（負）表偏多（少）或偏高（低），為零則表示突變，絕對值越大，表示變幅越大。因此由實部等值線圖可判讀不同時間尺度下信號在時域中的周期性變化及其位相。小波方差圖顯示了波動能量隨時間尺度的變化情況，可用來判讀氣候序列的主周期尺度[11]。本文利用MATLAB的小波分析工具箱采用復Morlet小波進行周期性分析；借助ArcGIS102的反距離權重插值（IDW）獲取氣候要素及其氣候傾向率的空間分布特征。

2結果與分析

2.1氣候要素的年際變化

1961-2015年，京津冀地區(qū)年均溫、年降水、年日照時數(shù)、年相對濕度、年平均風速的均值分別為：1043 ℃、50881 mm、2 63568 h、5848%、229 m/s。其中，年均氣溫最低為875 ℃，發(fā)生于1969年，最高達1173 ℃，發(fā)生于1998年；1999年降水量最少，僅34626 mm，1964年最多，為78935 mm；1965年日照時數(shù)最多，高達3 0687 h，2003年最少，僅2 2627 h；相對濕度在1964年取得最大值為6850%，在2005年取得最小值為5408%； 1969年平均風速取得最大值，為294 m/s，2003年取得最小值，為180 m/s?？梢娊邓?、日照、濕度、風速等4要素的最大值都出現(xiàn)在20世紀60年代，最小值都出現(xiàn)在2000年前后，而年均氣溫正與之相反。各要素變差系數(shù)由大到小排序為：降水（0182）>風速（0116）>日照時數(shù)（0067）>年均溫（0066）>相對濕度（0043），可見，降水的年際變化最劇烈，相對濕度的最小。

圖2表明，近55 a來，京津冀地區(qū)年均氣溫（圖2（a））呈明顯上升趨勢，氣候傾向率為029 ℃/10a，通過了P<001的顯著性檢驗；日照時數(shù)（圖2（c））、相對濕度（圖2（d））、平均風速（圖2（e））均呈明顯下降趨勢，氣候傾向率分別為：-906 h/（10a）、-058%/（10a）、-012[（m/s）/（10a）]，且均通過了001的顯著性檢驗；年降水（圖2（b））減少趨勢不明顯，氣候傾向率為-1176 mm/（10a），未通過01的顯著性檢驗。圖2還顯示：年降水和風速的波動性變化很強烈，且平均風速在2003年后呈現(xiàn)明顯的上升趨勢，日照時數(shù)和氣溫的波動變化相對較弱，這與變差系數(shù)的結果基本一致。圖2近55年京津冀地區(qū)各氣候要素變化趨勢

Fig.2Trends of climatic variables over the BeijingTianjinHebei region in recent 55 years

2.2突變性分析

首先以降水為例，闡明本文確定突變點的過程，其余各要素突變點的確定與之類似。圖3（b1）顯示年降水的UF曲線全程都在信度線之內，因此其增減趨勢都不顯著，這與降水氣候傾向率的檢驗結果一致。UF、UB曲線在信度線內存在包括1979在內的多個交點，說明1979等可能是突變點。由圖3（b2）知，降水累積距平曲線總體上表現(xiàn)出了先升后降的趨勢，但曲線波動頻繁、形狀不規(guī)則，表明降水存在劇烈的年際變化，整體的下降趨勢并不明顯。曲線在1979年、1996年存在兩個明顯的峰值，1961-1979年以上升趨勢為主，1979-1996年波動中略有上升，1996-2015年呈急劇下降趨勢，因此1979年、1996年是兩個可能的突變點?；瑒觮檢驗結果顯示只有1996年通過了005的顯著性檢驗。根據(jù)三種方法的結果可判斷：京津冀地區(qū)年降水存在1979年、1996年兩個突變點，而且都是由增加突變?yōu)闇p少。

圖3顯示：年均溫UF、UB曲線（圖3（a1））在信度線內交于1988年，其累積距平曲線（圖3（a2））在1988年取得最小值，且滑動t檢驗表明1988年通過了001的顯著性檢驗，因此，年均溫在1988年發(fā)生了由升溫到顯著升溫的突變。日照時數(shù)UF、UB曲線（c1）在信度線內沒有交點，其累積距平曲線（圖3（c2））在1989年取得最大值，且1989年通過了顯著性水平為001的滑動t檢驗，因此，日照時數(shù)在1989年發(fā)生了由減少到顯著減少的突變。相對濕度UF、UB曲線（圖3（d1））在信度線內交于1994年，其累積距平曲線（圖3（d2））存在1979年和1991年兩個明顯峰值，僅1979年通過了005顯著性水平的滑動t檢驗，因此，相對濕度在1979發(fā)生了由增加到減少的突變。平均風速UF、UB曲線（e1）在信度線內交于1978-1979年，其累積距平曲線（圖3（e2））最大值發(fā)生于1981年，滑動t檢驗顯示，1978-1981年4個年份均通過了001的顯著性檢驗，進一步考慮到：滑動t檢驗時無論子序列取5、10還是12，只有1981年都通過了001的顯著性檢驗，因此判斷，年均風速在1981年發(fā)生了由減小到顯著減小的突變。圖3近55 a京津冀地區(qū)各氣候要素MannKendall檢測及累積距平曲線

Fig.3MannKendall test and cumulative anomaly of climatic variables over the Beijing-TianjinHebei region in recent 55 years

2.3周期性分析

由圖4（a2）知：年均溫的第1、2、3、4主周期尺度分別是21 a、14 a、10 a、6 a。圖4（a1）顯示：在6 a左右的尺度上，信號振蕩相對較強，影響范圍主要是1961年到20世紀70年代中后期。從1961-1973年氣溫存在相對比較規(guī)律的5次降升的準周期變化，平均周期為3 a左右。在10 a左右的尺度上，主要是在20世紀70年代中后期至21世紀初存在顯著的信號波動，1975-2003年間共存在7次比較規(guī)律的降溫升溫的準周期變化，平均周期為4 a左右。在14 a左右的尺度上，信號的強烈振蕩主要存在于1961年到21世紀初，該信號比較穩(wěn)定，存續(xù)時間也較長，在1961-1994年間經(jīng)歷了6次降溫升溫的準周期性變化，平均周期為6 a左右。在21 a左右的時間尺度上，周期性波動自1961年開始逐漸增強，20世紀90年代后，尤其進入21世紀后，振蕩變得極為顯著，整個研究時域內共經(jīng)歷了6～7次氣溫升降的準周期性變化，平均周期為9 a左右。圖4近55 a京津冀地區(qū)各氣候要素小波分析結果

Fig.4Morlet wavelet analysis of climatic variables over the BeijingTianjinHebei region in recent 55 years

以上分析表明：京津冀地區(qū)年均溫的年際變化中存在的典型周期包括9 a、6 a、4 a、3 a。

由圖4（b2）知：年降水的第1、2、3、4主周期尺度分別是8 a、11 a、21 a、26 a。圖4（b2）表明：在8 a左右的尺度上，振蕩最劇烈，但衰減很快，主要影響1961年到20世紀70年代中期，該時間尺度下，年降水量在1961-1977年經(jīng)歷了5次減增的準周期性變化，平均周期為3 a左右。在11 a左右的尺度上，振蕩強烈、穩(wěn)定、且基本貫穿全研究時域，1972-2006年存在8次降水增減的準周期性變化，平均周期為4 a左右。在21 a左右的尺度上，顯著的信號振蕩始于20世紀80年代，其后漸強，直到2015年，1982-2015年間信號的周期性變化很明顯，共有4次增減的準周期變化，平均周期為8 a左右，從該尺度上看，2015年時，負值等值線未閉合，表示在該尺度上將來一段時間京津冀地區(qū)年降水處于偏少階段。在26 a左右的尺度上，振蕩強烈，但衰減趨勢明顯，尤其是進入21世紀后，信號變得很微弱，1961-1991年共有3次降水增減的準周期變化，平均周期為10 a左右。以上分析表明：年降水的年際變化中存在的典型周期包括3 a、4 a、8 a、10 a。

圖4（c2）顯示：年日照時數(shù)的第1、2、3、4主周期尺度分別是8 a、16 a、20 a、4 a。由圖4（c1）可見：在8 a左右的時間尺度上，波動最強烈，主要影響范圍大致是1961年到80年代末期， 1961-1988年日照時數(shù)經(jīng)歷了9次增多減少的準周期性變化，平均周期為3 a左右。10 a時間尺度盡管不是主周期尺度，但該尺度上從20世紀90年代中后期開始出現(xiàn)一個強震蕩區(qū)，直到2015年，1995-2014年存在5次日照時數(shù)增多減小的準周期變化，平均周期為4 a左右。在16～20 a的尺度上，信號的強烈振蕩主要存在于20世紀70年代初到21世紀初，1973-2005年日照時數(shù)經(jīng)歷了5次比較強烈的減少增多的準周期性變化，平均周期為7 a左右。以上分析表明：日照時數(shù)的年際變化中存在的典型周期包括3 a、4 a、7 a。

圖4（d2）表明：年均相對濕度的第1、2、3、4主周期尺度分別是17 a、8 a、10 a、14 a。圖4（d1）顯示：在8 a左右的尺度上，振蕩能量非常強，但持續(xù)時間短，影響范圍主要是在1970年以前， 1961-1970年，存在3次減小增加的準周期性變化，平均周期為3 a左右。在10 a左右的尺度上，波動能量也比較強，主要影響范圍是20世紀80年代末到2005年前后，該時域內有4次減小增加的準周期性變化，平均周期為4 a左右。在10 a左右的尺度上，振蕩偏弱，但存在于全時域，1961-2015年都具有比較明顯的周期波動，平均周期為6 a左右。在17 a左右的尺度上，振蕩強烈且比較持久，在1961-1990年間經(jīng)歷了4次增加減小的準周期性變化，平均周期為7 a。以上分析表明：相對濕度的年際變化中存在的典型周期包括7 a、6 a、4 a、3 a。

圖4（e2）可見：年均風速的第1、2、3、4主周期尺度分別是25 a、9 a、13 a、4 a。由圖4（e1）可知：在9 a左右的尺度上，1973年之前信號波動很強烈，但衰減迅速，1987年之后變得極其微弱，進入21世紀后才稍有回升，1961-1975年存在4次風速增大減小的準周期性變化，平均周期為4 a左右。在13 a左右的尺度上，信號在1961-1971年和2001-2011兩個時域內振蕩較強，其中，1961-1971年間存在2次增大減小的周期變化，在2001-2011年間存在2次減小增大的周期變化，平均周期為6 a左右。在25 a左右的尺度上，波動非常強烈，但衰減快，1961-1991年間風速經(jīng)歷了3次減小增加的準周期性變化，平均周期為10 a左右，從該尺度上看，2015年時負值等值線尚未閉合，這表示在25 a時間尺度上京津冀地區(qū)的年均風速在未來一段時間仍會偏小。以上分析表明：京津冀地區(qū)年均風速的年際變化中存在的主要周期包括10 a、6 a、4 a。

2.4氣候要素的空間分布特征

普通克里金插值法與反距離權重法（IDW）是分析水文氣象要素空間分布特征時常用的方法，對兩者的適用性，學者們意見不一[1215]。本文經(jīng)對比發(fā)現(xiàn)，IDW更適于描述京津冀地區(qū)氣象要素的空間分布特征，因此本文借助ArcGIS軟件，通過IDW插值法分析京津冀地區(qū)氣候要素的空間特征。

圖5（a1）顯示多年平均氣溫基本呈現(xiàn)由南向北、由東向西遞減的趨勢，高溫區(qū)位于京津冀地區(qū)南部，低溫區(qū)位于西北部。其中，多年平均氣溫最高的是邢臺，達1393 ℃，次之為石家莊，達1357 ℃；多年平均氣溫最低的是張北，僅有334 ℃，次之是圍場，為526 ℃。由圖5（a2）可見京津冀地區(qū)存在一橫一縱兩個明顯的高幅增溫帶，其中橫貫京津冀地區(qū)的蔚縣北京遵化樂亭一帶增溫最明顯，它們占據(jù)了全區(qū)域增溫幅度的前四位，其中，增幅最大的是蔚縣，達到了050 ℃/（10a），其次是北京，達045 ℃/（10a）；另外一個增幅明顯的地帶為張北蔚縣石家莊邢臺一線，其中邢臺的氣溫傾向率達到了040 ℃/10a，石家莊為037 ℃/（10a），多年平均氣溫最低的張北地區(qū)其傾向率也達到了035 ℃/（10a）。京津冀地區(qū)北部增溫趨勢緩慢，其中，承德表現(xiàn)出了微弱的降溫趨勢。圖5近55 a京津冀地區(qū)各氣候要素及其氣候傾向率的空間分布

Fig.5The spatial distribution of climatic variables and their climate tendency rate over the BeijingTianjinHebei region in recent 55 years

圖5（b1）表明，多年平均降水量基本呈現(xiàn)由東向西遞減的趨勢，相對來說，京津冀地區(qū)東北一隅降水比較豐沛，其中多年平均降水量最大的是遵化，達到了71392 mm，其次是青龍，為69144 mm，整個京津冀地區(qū)只有遵化多年平均降水量在700 mm以上，在600 mm以上的除青龍外，還包括秦皇島、北京、唐山、樂亭四地。京津冀地區(qū)西部降水普遍相對較少，其中最少的是懷來，僅有38032 mm，次之是張北，也只有38242 mm；在廊坊以南，除南宮多年平均降水量為47650 mm以外，其他地區(qū)都在500～600 mm之間。圖5（b2）顯示除圍場（圖中藍色部分）以723 mm/（10a）的增幅呈現(xiàn)一定的降水增加趨勢外，其他地區(qū)均為減少趨勢，東部地區(qū)減幅最大，南部次之，西北部減幅最小。全區(qū)域減幅最大的是遵化，達到3028 mm/（10a），次之是黃驊，為2797 mm/（10a），青龍地區(qū)減幅為-2579 mm/（10a）；蔚縣、懷來、張家口、張北一帶減幅極其微弱。

圖5（c1）顯示，多年平均日照時數(shù)基本呈現(xiàn)以懷來-蔚縣兩地為峰向周圍遞減的趨勢，其中，懷來地區(qū)的多年平均日照時數(shù)最多，達2 970 h/a，次之是蔚縣，為2 888 h/a，最少的是邢臺地區(qū)，只有2 355 h/a，其次為石家莊，僅2 424 h/a。漸變過程中也存在波動，由北向南遞減的過程中，饒陽地區(qū)的日照時數(shù)是2 631 h/a，比其周圍保定、南宮、廊坊等地的都要高；由西向東遞減的過程中，塘沽的日照時數(shù)卻達到一個小高峰，為2 644 h/a。由圖5（c2）知全區(qū)域日照時數(shù)均呈減少趨勢，減幅以邢臺石家莊保定廊坊塘沽一線為嶺，向兩側地區(qū)逐漸變緩，其中減幅最大的是石家莊，達176 h/（10a），其次是塘沽，達149 h/（10a），再其次是邢臺，為139 h/（10a）。而在北京以北，除承德、遵化兩地外，其他地區(qū)減幅都在100 h/（10a）以下，全區(qū)域減幅最小的是豐寧，為3284 h/（10a），其次是蔚縣，為3440 h/（10a）。

由圖5（d1）可見京津冀地區(qū)多年平均相對濕度由南向北、由東向西逐漸降低。全區(qū)域最高的是樂亭地區(qū)，高達6588%，南宮稍低，為6404%，最低的是張家口，僅4745%，次低是懷來，為5049%。圖5（d2）中藍色表示相對濕度呈增加趨勢，可見，圍場承德秦皇島一帶相對濕度基本呈遞增趨勢，其中增幅最大的是秦皇島，達118%/（10a），其次為承德，增幅為105%/（10a），圍場增幅最小，為023%/（10a）。除該三地之外，其他均呈減少趨勢，由北向南減少速率越來越快，至北京達到第一個峰值，減幅為161%/（10a），到塘沽達到全區(qū)域減幅最大值，為204%/（10a），由北京塘沽一線向西南部基本呈減幅先降后增的趨勢，到石家莊、邢臺，年平均相對濕度減幅又形成兩個峰值，其中，石家莊減幅123%/（10a），邢臺減幅為184%/（10a）。

圖5（e1）顯示多年平均風速在京津冀地區(qū)的空間分布呈馬鞍型，有兩個峰值，一個是壩上張北，多年平均風速達405 m/s，一個是沿海的塘沽，多年平均風速達401 m/s，除此二者，其他地區(qū)多年平均風速均在31 m/s以下。張北到北京、塘沽到北京均呈遞減趨勢，北京多年平均風速為245 m/s，由張北北京塘沽一線分向西南、東北方向遞減，全區(qū)域風速最小的是承德，為128 m/s，次之是青龍，為127 m/s。北京以南，多年平均風速最低的是石家莊，為169 m/s。由圖5（e2）可見除圍場、豐寧、承德三地年平均風速呈微弱遞增趨勢之外，其他地區(qū)年均呈遞減趨勢，其中減幅最大的是壩上張北，達046[（m/s）/（10a）]，僅次于張北的是沿海的塘沽，為044[（m/s）/（10a）]，年平均風速減幅由此二地區(qū)均向北京方向遞減，北京年均風速傾向率為-009[（m/s）/（10a）]。年平均風速減幅由張北北京塘沽一線向西南、東北方向基本均呈遞減趨勢，全區(qū)域減幅最小的是蔚縣，僅001[（m/s）/（10a）]，其次是青龍，減幅為002[（m/s）/（10a）]，圍場、豐寧、承德三地年平均風速的增幅均低于009[（m/s）/（10a）]。

3討論

（1）研究發(fā)現(xiàn)，近55 a來京津冀地區(qū)日照時數(shù)的變化具有一定的特殊性，首先是其減少速率高達9060 h/（10a）（P<<001），明顯高于全國1961-2007年397 h/（10a）的減少速率[16]；再者，有研究表明中國地區(qū)日照時數(shù)減少趨勢在20世紀90年代后基本停止[1617]，而圖3（c1）顯示京津冀地區(qū)20世紀90年代后日照時數(shù)減少趨勢依然非常明顯，直到2015年也沒有變緩的傾向。本文認為京津冀地區(qū)嚴重的大氣污染和風速的顯著下降是導致該現(xiàn)象的重要原因。不少學者的研究成果[16，1819]表明：除自然因素外，人類活動導致的大氣污染是我國尤其是我國東部地區(qū)日照時數(shù)減少的主因，而京津冀地區(qū)大氣污染是極端嚴重的[2021]；另外，本文結果表明，近55 a來京津冀地區(qū)風速以011[（m/s）/（10a）]（P<001）的速率減小，而風速的減小不利于大氣污染物的消散，這進一步加劇了日照時數(shù)的減少趨勢。

（2）京津冀地區(qū)升溫顯著（P<<001），變率為029 ℃/（10a），這與全國1964-2007年028 ℃/（10a）的增溫速率[22]基本一致，略高于華北地區(qū)1960-2013年023 ℃/（10a）的升溫速率[23]。本文認為該區(qū)域在日照時數(shù)顯著減少的情況下氣溫卻顯著上升的原因之一仍然是大氣污染，大氣中煙、霾等微粒物質增多，導致氣溶膠含量升高，大氣逆輻射增強，進而致使溫度上升[23]。

（3）京津冀地區(qū)年降水存在1979年、1996年兩個突變點，而且都是由增加突變?yōu)闇p少，這與張一馳等[24]華北地區(qū)1951-2009年降水無突變的結論不一致。對此，時空差異是可能的原因之一，另外，本文用MannKendall突變檢測、滑動t檢驗、累積距平3方法相結合進行突變分析，而張一馳等[24]利用MannWhitney檢驗法進行突變識別，不同的突變檢測方法也可能導致南轅北轍的結論[9]。

（4）本文研究表明，承德地區(qū)氣候變化與眾不同，該地區(qū)涵蓋承德、圍場、豐寧3站，其中，承德與圍場2站的風速、濕度的變化趨勢與其他大多數(shù)站都背道而馳，此外，承德站的氣溫、圍場站的降水、豐寧站的風速變化趨勢也與其他多數(shù)站截然不同，可見承德地區(qū)的氣候變化與京津冀其他地區(qū)的變化趨勢迥異，但這其中的原因還需搜集更多的資料做進一步的分析。

4結論

（1）1961-2015年，京津冀地區(qū)年均氣溫以029 ℃/（10a）的速率顯著上升（P<<001）；年降水以1176 mm/（10a）的速率減少，但趨勢不顯著；日照時數(shù)、相對濕度、平均風速均呈顯著下降趨勢，氣候傾向率分別為-9060 h/（10a）、-058%/（10a）、-011[（m/s）/（10a）]。

（2）氣候要素的突變主要發(fā)生在80年代前后。最早發(fā)生突變的是降水和相對濕度，二者都在1979年由增加突變?yōu)闇p少，年平均風速在1981年由減小突變?yōu)轱@著減小，年均氣溫在1988年由上升突變?yōu)轱@著上升，年日照時數(shù)在1989年由減少突變?yōu)轱@著減少，年降水在1996年再次由增加突變?yōu)闇p少。

（3）研究時段內不同的子時域具有不同的時間尺度結構，而同一時間尺度上的信號振蕩在時域中也是不斷變化的，這就決定了氣候要素周期特征的復雜性，年均氣溫、年降水量、年日照時數(shù)、年平均相對濕度、年平均風速的第一主周期尺度分別是21 a、8 a、8 a、17 a、25 a，對應的平均周期分別是9 a、3 a、3 a、7 a、10 a。

（4）空間上，年均溫呈現(xiàn)出明顯的由南向北遞減的趨勢，而升溫速率則存在蔚縣北京樂亭和張北蔚縣石家莊邢臺一橫一縱兩個高值區(qū)；年降水及其變化率均表現(xiàn)為自東向西遞減的總體趨勢；日照時數(shù)由南向北遞增，而其減少速率則由南向北遞減；相對濕度由東南向西北遞減，其減小速率則由西南向東北遞減；平均風速及其變率均以張北和塘沽為兩個典型高值中心向其他地區(qū)遞減。

參考文獻（References）：

[1]許穎，唐海萍.河北懷來盆地近60a氣候變化特征及其影響[J].北京師范大學學報（自然科學版），2015，51（3）：293298.（XU Y，TANG H P.Characteristics and impacts of climate change in Huailai Basin during the past 60 years[J].Journal of Beijing Normal University （Natural Science），2015，51（3）：293298.（in Chinese）） Doi：10.16360/j.cnki.jbnuns.2015.03.013.

[2]李鵬飛，劉文軍，趙昕奕.京津冀地區(qū)近50年氣溫、降水與潛在蒸散量變化分析[J].干旱區(qū)資源與環(huán)境，2015，29（3）：137143.（LI P F，LIU W J，ZHAO X Y.The changes of atmospheric temperature，precipitation and potential evapotranspiration in BeijingTianjinHebei region in recent 50 years[J].Journal of Arid Land Resources and Environment，2015，29（3）：137143.（in Chinese）） Doi：10.13448 /j.cnki.jalre.2015.094.

[3]劉芳圓，肖嗣榮，劉寒，等.河北地區(qū)氣候變化及其對農業(yè)的影響[J].地理與地理信息科學，2014，30（4）：122126.（LIU F Y，XIAO S R，LIU H，et al.Research of impacts of climate change on agriculture in Hebei Region[J].Geography and GeoInformation Science，2014，30（4）：122126.（in Chinese）） Doi：10.3969 /j.issn.16720504.2014.04.025.

[4]張國華，張江濤，金曉青，等.京津冀城市高溫的氣候特征及城市化效應[J].生態(tài)環(huán)境學報，2012，21（3）：455463.（ZHANG G H，ZHANG J T，JIN X Q，et al.Climate characteristics and effects of urbanization of the urban high temperature of Beijing，Tianjin and Hebei，China[J].Ecology and Environmental Sciences，2012，21（3）：455463.（in Chinese）） Doi：10.16258/j.cnki.16745906.2012.03.010.

[5]段麗瑤，丁一匯，任雨.19212010年天津氣溫和降水量序列的多尺度分析[J].氣候與環(huán)境研究，2014，19（4）：515522.（DUAN L Y，DING Y H，REN Y.Multiscale examination on the temperature and precipitation series in Tianjin during 19212010[J].Climatic and Environmental Research，2014，19（4）：515522.（in Chinese）） Doi：10.3878/j.issn.10069585.2014.13020.

[6]胡保昆，竇以文，儲偉.北京城區(qū)氣溫變化特征[J].氣象科技，2014，42（5）：852855.（HU B K，DOU Y W，CHU W.Characteristics of temperature variation in Beijing urban area[J].Meteorological Science and Technology，2014，42（5）：852855.（in Chinese））[7]向亮，郝立生，安月改.等.51a河北省降水時空分布及變化特征[J].干旱區(qū)地理，2014，37（1）：5665.（XIANG L，HAO L S，AN Y G，et al.TimeSpatial distribution and variational characteristics of rainfall in Hebei province in 51 years[J].Arid Land Geography，2014，37（1）：5665.（in Chinese）） Doi：10.13826/j.cnki.cn651103/x.2014.01.008.

[8]蘆佳玉，延軍平，王鵬濤，等.全球變化背景下陜甘寧地區(qū)風速時空變化特征[J].中國沙漠，2017，37（3）：554561.（LU J Y，YAN J P，WANG P T，et，al.Wind speed change characteristics of ShanxiGansuNingxia Area in 19602014[J].Journal of Desert Research，2017，37（3）：554561.（in Chinese）） DOI：10.7522/j.issn.1000694X.2016.00114.

[9]魏鳳英.現(xiàn)代氣候統(tǒng)計診斷與預測技術（第2版）[M].北京：氣象出版社，2007.（WEI F Y.Modern climatic statistical diagnosis and prediction technology（2nd）[M].Beijing：China Meteorological Press，2007.（in Chinese））

[10]彭俊翔，伍永年，胡維平，等.長蕩湖近61a降水量演化特征[J].長江流域資源與環(huán)境.2016，25（2）：292299.（PENG J X，WU Y N，HU X P，et al.The evolution characteristics of precipitation in the past 61 years at the lake Changdanghu[J].Resources and Environment in the Yangtze Basin，2016，25（2）：292299.（in Chinese）） DOI：10.11870/cjlyzyyhj201602.015.

[11]胡毅鴻，李景保.19512015年洞庭湖區(qū)旱澇演變及典型年份旱澇急轉特征分析[J].農業(yè)工程學報.2017，33（7）：107115.（HU Y H，LI J B.Analysis on evolution of droughtflood and its abrupt alternation in typical year from 1951 to 2015 in Dongting Lake area[J].Transactions of the Chinese Society of Agricultural Engineering，2017，33（7）：107115.（in Chinese）） DOI：10.11975/j.issn.10026819.2017.07.014.

[12]趙雪雁，張欽，王亞茹，等.近55a來青藏高原東部氣候演變特征[J].干旱區(qū)研究，2015，32（6）：10881096.（ZHAO X Y，ZHANG Q，WANG Y R，et al.Climate evolution characteristics of eastern Tibetan Plateau for the last 55 years[J].Arid Zone Research，2015，32（6）：10881096.（in Chinese）） DOI：10.13866 /j.azr.2015.06.07.

[13]曹永強，劉佳佳，高璐.遼寧省大雨以上降水日數(shù)分布與趨勢分析[J].地理科學進展，2015，34（8）：10521060.（CAO Y J，LIU J J，GAO L.Distribution and trend analysis of heavy precipitation days in Liaoning Province[J].Progress in Geography，2015，34（8）：10521060.（in Chinese）） Doi：10.18306/dlkxjz.2015.08.012.

[14]史佳良，王秀茹，李淑芳，等.近50年來河南省氣溫和降水時空變化特征分析[J].水土保持研究，2017，24（3）：151156.（SHI J L，WANG X R，LI S F，et al.Variation characteristics of air temperature and precipitation in Henan province in recent 50 years[J].Research of Soil and Water Conservation，2017，24（3）：151156.（in Chinese）） Doi：10.13869/j.cnki.rswc.2017.03.028.

[15]王理萍，王樹仿，張劉東，等.19602013年云南省降水時空變化特征分析[J].水電能源科學，2016，34（12）：2024.（WANG L P，WANG S F，ZHANG L D，et al.Spatialtemporal characteristics of precipitation in Yunnan province from 1960 to 2013[J].Water Resources and Power，2016，34（12）：2024.（in Chinese））

[16]趙東，羅勇，高歌，等.1961年至2007年中國日照的演變及其關鍵氣候特征[J].資源科學，2010，32（4）：701711.（ZHAO D，LUO Y，GAO G，et al.LongTerm changes and essential climatic characteristics of sunshine duration over China during 19612007[J].Resources Science，2010，32（4）：701711.（in Chinese））

[17]李慧群，付遵濤，聞新宇.等.中國地區(qū)日照時數(shù)近50年來的變化特征[J].氣候與環(huán)境研究，2013，18（2）：203209.（LI H Q，F(xiàn)U Z T，WEN X Y，et al.Characteristic analysis of sunshine duration change in China during the last 50 years[J].Climatic and Environmental Research，2013，18（2）：203209.（in Chinese）） Doi：10.3878/j.issn.10069585.2012.11052.

[18]鄭小波，羅宇翔，段長春，等.云貴高原近45年來日照及能見度變化及其成因初步分析[J].高原氣象，2010，29（4）：992998.（ZHENG X B，LUO Y X，DUAN C C，et al.Change trends and causes in sunshine duration and visibility over YunnanGuizhou Plateau in recent 45 years[J].Plateau Meteorology，2010，29（4）：992998.（in Chinese））

[19]劉玉英，韋小麗，李宇凡.19612012年吉林省日照時數(shù)的變化特征及影響因素[J].自然資源學報，2015（8）：13671377.（LIU Y Y，WEI X L，LI Y F.Variation of sunshine duration and related driving forces in Jilin Province during 19612012[J].Journal of Natural Resources，2015（8）：13671377.（in Chinese）） Doi：10.10.11849/zrzyxb.2015.08.011.

[20]王躍思，張軍科，王莉莉，等.京津冀區(qū)域大氣霾污染研究意義、現(xiàn)狀及展望[J].地球科學進展，2014，29（3）：388396.（WANG Y S，ZHANG J K，WANG L L，et al.Researching significance，status and expectation of haze in BeijingTianjin Hebei region[J].Advance in Earth Science，2014，29（3）：388396.（in Chinese）） Doi：10.11867 /j.issn.10018166.2014.03.0388.

[21]趙普生，徐曉峰，孟偉，等.京津冀區(qū)域霾天氣特征[J].中國環(huán)境科學，2012，32（1）：3136.（ZHAO P S，XU X F，MENG W，et al.Characteristics of hazy days in the region of Beijing，Tianjin，and Hebei[J].China Environmental Science，2012，32（1）：3136.（in Chinese））

[22]范澤孟，岳天祥，陳傳法，等.中國氣溫與降水的時空變化趨勢分析[J].地球信息科學學報，2011，13（4）：526533.（FAN Z M，YUE T X，CHEN C F，et al.Spatian change trends of temperature and precipitation in China[J].Journal of GeoInformation Science，2011，13（4）：526533.（in Chinese）） Doi：10.3724/SP.J 1047.2011.00526.

[23]阿多，熊凱，趙文吉，等.1960～2013年華北平原氣候變化時空特征及其對太陽活動和大氣環(huán)境變化的響應[J].地理科學，2016，36（10）：15551564.（A D，XIONG K，ZHAO W J，et al.Temporal trend of climate change and mutation analysis of North China Plain during 1960 to 2013[J].Scientia Geographica Sinica，2016，36（10）：15551564.（in Chinese）） Doi：10.13249/j.cnki.sgs.2016.10.013.

[24]張一馳，吳凱，于靜潔，等.華北地區(qū)19512009年氣溫、降水變化特征[J].自然資源學報，2011，26（11）：19301941.（ZHANG Y C，WU K，YU J J，et al.Characteristics of precipitation and air temperature variation during 19512009 in North China[J].Journal of Natural Resources，2011，26（11）：19301941.（in Chinese））第16卷第3期