侯偉 馬明明 胡德強 施思 官滿元

摘 要 利用1951～2010年海口市國家基準氣候站逐月平均氣溫、降水量、日照時數(shù)等數(shù)據(jù)，采用線性傾向率、Mann-Kendall非參數(shù)檢驗、滑動T檢驗、累積距平法和Morlet小波分析等數(shù)理統(tǒng)計方法，分析了?？谑?951～2010年的氣候變化特征。結果表明：?？谑薪?0 a氣候變化表現(xiàn)為氣溫顯著上升、日照時數(shù)顯著減少，但降水量變化不明顯。20世紀90年代是氣溫升高、降水量增多、日照時數(shù)減少的顯著變化穩(wěn)定期，2000年以后氣候變化波動最強烈。旱季氣溫增幅明顯大于雨季，但降水量變化較小。此外，旱季氣候要素波動大于雨季。氣候突變和小波分析得出，1986年是氣溫由低到高的突變點，1967～1978、1998～2010年分別是冷期和暖期；1976年是降水量發(fā)生由多到少的突變點，1993年是日照時數(shù)發(fā)生由多到少的突變點。?？谑袣鉁睾徒邓孔兓嬖?～10 a左右的周期，日照時數(shù)存在5～8 a、25～28 a左右的周期。

關鍵詞 氣候變化；Mann-Kandall非參數(shù)檢驗；滑動T檢驗；累積距平；Morlet小波；?？谑?/p>

中圖分類號 P467 文獻標識碼 A

Abstract Haikou is one of the most important ecological tourist city in China， even throughout the World. The ecological environment and climate is the most precious natural resources. Carrying out the research of climate change characteristics in Haikou could be used to enhance its ability to adapt to the long-term effects of climate change， and delay or avoid the threats caused by climate change to regional ecological security. This paper used the data of monthly mean temperature， precipitation and sunshine duration during the period from 1951-2010 observed by national benchmark climate station and combined with some mathematical statistics methods such as linear regression method， Mann-Kendall nonparametric test， moving T-test， cumulative anomaly method and Morlet wavelet to analyze climate change characteristics in Haikou City. The results indicated that the climate change presented increasing annual mean temperature， decreasing annual sunshine duration， and climate tendency rate was 0.18 ℃/（10 a）and -83.1 h/（10 a）at the significance level α=0.01 respectively， the increasing rate of annual mean temperature was distinctly lower than that in whole China， and decreasing rate of annual sunshine duration significantly higher than in whole China. Moreover， no significant long-term change was observed in annual precipitation in the past 60 a. The stable period of climate changed significantly appeared in 1990s， mainly showed that increasing annual mean temperature and annual precipitation， decreasing annual sunshine duration， and fluctuations in climate change became the strongest since 2000 years. This indicated that the frequency and intensity of extreme weather events such as heat damage， floods and droughts will increase in the near future. In addition， the increasing rate of temperature in dry season was significantly higher than that in rainy season， but the change rate of precipitation was more obvious in rainy season than that in dry season. In addition， the fluctuation of climatic elements was more obvious in dry season than that in rainy season. The study of abrupt climate change and Morlet wavelet revealed that 1986 was an abrupt change point of annual mean temperature increasing， 1976 was an abrupt change point of annual precipitation reduction， 1993 was an abrupt change point of annual sunshine duration reduction， 1967～1978 and 1998～2010 correspond to cold period and warm period in Haikou respectively. Annual mean temperature and annual precipitation of Haikou existed evident characteristic time scale with 5-10 a， the periodic oscillation of annual sunshine duration were in 5-8 and 25-28 a.

Key words Climate change； Mann-Kendall test； Moving T-test； Cumulative anomaly method； Morlet wavelet； Haikou

doi 10.3969/j.issn.1000-2561.2016.10.028

氣候變化是21世紀全球面臨的一個重大問題[1]，IPCC在第五次評估報告（AR5）中指出，全球氣候系統(tǒng)變暖的事實毋庸置疑，1880～2012年全球平均溫度已升高0.85 ℃（0.65～1.06 ℃），在過去30 a中，每10 a溫度增幅均高于1850年以來任何時期，在北半球，1983～2012年可能是最近1 400 a氣溫最高的30 a。越來越多地區(qū)強降水頻率、強度和降水量在增加[2]。與全球氣候變暖相應，中國近100 a來平均溫度上升了0.79 ℃，尤其近50 a（1956～2005年）的增幅最明顯，氣候變暖對中國南澇北旱的分布格局產(chǎn)生了顯著影響[3-4]。

中國氣候類型多樣，時間和空間上的較大差異，氣候變化背景下的區(qū)域氣候響應差異性很大[3]。深入研究區(qū)域性長時間序列的氣候變化背景、特征、原因及影響不僅對目標區(qū)域自然生態(tài)保育和社會經(jīng)濟發(fā)展提供參考，更是對應對未來氣候變化的適應和應急能力，最大限度避免或減輕極端天氣氣候事件有重大意義。相較國外，中國在氣候變化研究方面起步較晚，觀測設備性能、觀測資料序列、研究結果可靠性等方面普遍不理想[3]，研究實力與國外發(fā)達國家仍有距離。目前，中國在區(qū)域性氣候變化方面的研究已較多，如葉正偉等[5]利用Mann-Kendall方法對廬山旅游區(qū)的氣溫和降水進行了變化趨勢分析和突變檢驗；張立偉等[6]應用Mann-Kendall突變檢驗和Kinging插值法分析了南嶺南北地區(qū)氣溫的時空變化特征；任永建等[7]考慮到氣溫序列的非均一性，通過對資料進行插補后對武漢區(qū)域近百年地表氣溫變化趨勢進行了研究；李淼等[8]對北京地區(qū)近300 a降水變化進行了趨勢和突變檢驗，并采用Morlet小波函數(shù)開展了降水變化的多時間尺度周期性變化分析；何萍等[9]利用常規(guī)距平法、三階多項式擬合趨勢法、5 a滑動平均法分析了云南楚雄市氣候年代際變化特征。目前針對海口市氣候變化研究報道較少，研究內(nèi)容多采用累積距平和線性傾向估計方法對氣溫、降水和日照等要素進行簡單變化趨勢分析，而對多要素的氣候突變多途徑驗證和多時間尺度周期性變化研究還屬空白。本文利用?？谑袊覛庀蠡鶞收咀?951～2010年逐月氣溫、降水量、日照時數(shù)等氣象資料，采用線性傾向率、Mann-Kendall非參數(shù)檢驗、滑動T檢驗、累積距平和Morlet小波分析等統(tǒng)計方法對?？谑薪?0 a來氣候變化特征進行初步研究，了解其氣候變化趨勢和周期規(guī)律，以期對區(qū)域自然生態(tài)保育、社會經(jīng)濟發(fā)展和氣象防災減災有指導作用。

1 材料與方法

1.1 資料來源

?？跉庀笥^測站是國家氣象基準觀測站，具有完整的氣象觀測資料序列，能代表海口市的氣候狀態(tài)和變化趨勢。本文所用1951～2010年逐月平均氣溫、降水量、日照時數(shù)等資料來自海南省氣象局，所有數(shù)據(jù)經(jīng)過均一性檢驗和質量控制。海口市地處低緯度熱帶北緣，屬熱帶海洋性季風氣候，一年中降水變化顯著，干濕分明，雨季降水量占全年的83.2%左右，故本文又分雨季（5～10 月）和旱季（11～次年4月）進行統(tǒng)計。

1.2 分析方法

利用海口氣象觀測站近60 a的年平均氣溫、年降水量、年日照時數(shù)等資料，采用滑動平均、氣候傾向率分析其氣候變化趨勢，對結果的顯著性使用F檢驗和Mann-Kendall（簡稱M-K）非參數(shù)趨勢檢驗。氣候突變采用M-K非參數(shù)突變檢驗[10]、滑動T檢驗[10]和累積距平法對時間序列突變點和突變區(qū)域進行分析，并檢驗突變點真?zhèn)危岣呓Y論的可信度，在進行氣候要素年代際變化趨勢分析中用到的變異系數(shù)（CV）是采用統(tǒng)計年代內(nèi)要素的標準差與平均數(shù)比值的百分率來計算，變異系數(shù)越大，表明統(tǒng)計年代的氣候要素波動越大。此外，本文采用Morlet小波分析[11]對?？谑袣夂蜃兓闹芷谛砸?guī)律進行了分析。

1.3 數(shù)據(jù)處理與分析

采用Microsoft Excel 2007對氣象數(shù)據(jù)資料進行統(tǒng)計分析。使用Compaq Visual Fortran 6.6進行M-K非參數(shù)突變檢驗和小波分析。Excel 2007和Adobe Illustrator CS5用于圖形制作。

2 結果與分析

2.1 氣候變化趨勢分析

?？谑心昶骄鶜鉁亍⒛杲邓亢湍耆照諘r數(shù)的變化趨勢見圖1。由圖1可看出，?？谑?951～2010年的年平均氣溫呈逐漸上升趨勢，其中最高出現(xiàn)在1998年（25.4 ℃），最低出現(xiàn)在1967、1968、1971、1974年（均為23.2 ℃），通過線性擬合的氣候傾向率為0.18 ℃/（10 a）。經(jīng)過F顯著性檢驗和M-K檢驗得知，近60 a來?？谑心昃鶜鉁卦鰷刳厔萃ㄟ^了α=0.01的顯著水平。

與氣溫變化不同，?？谑心杲邓孔兓厔莶幻黠@，波動有逐漸增大的趨勢，特別是2000年以后波動更為劇烈。降水量最多在2009年（2 628.2 mm），最少在1977年（871.5 mm），線性擬合的年降水量氣候傾向率為10.5 mm/（10 a）。年降水量增加趨勢未通過F檢驗，而在M-K趨勢檢驗中得知年降水量是減少的趨勢，同樣未通過α=0.05的水平。這說明?？谑心杲邓孔兓厔莶幻黠@，且波動較大的特征。

通過對日照時數(shù)的變化趨勢分析得出，海口市近60 a年來日照時數(shù)明顯下降，最高值在1963年（2 553.8 h），最低值在1953年（1 316.9 h），通過線性擬合的氣候傾向率為-83.1 h/（10 a）。由F檢驗和M-K檢驗得知，年日照時數(shù)增長趨勢通過了α=0.01的顯著水平。

從年代際變化來看（表1），20世紀70年代開始氣溫呈現(xiàn)上升趨勢，而變異系數(shù)則表現(xiàn)為下降（2000～2010年例外），90年代氣溫上升幅度最大，相比80年代上升0.57 ℃，但其變異系數(shù)較小，2000年以后變異系數(shù)最大，表明90年代是?？谑袣夂蝻@著增暖的平穩(wěn)期，而2000年以后氣溫開始出現(xiàn)較大波動。雨季、旱季氣溫變化與年代際氣溫變化趨勢類似，但旱季增溫幅度大于雨季。此外，旱季氣溫變異系數(shù)相比雨季大，但總體有減小趨勢。這表明旱季氣溫變化逐漸穩(wěn)定。

由表2看出，從20世紀90年代開始降水量顯著增加，但其變異系數(shù)較小，而2000年以后的10 a降水量增幅最大，相比90年代降水量增加124.33 mm，且變異系數(shù)最大，說明90年代降水量增加穩(wěn)定，2000年以后降水量變化異常。雨季和旱季降水量變化趨勢相反，即雨季降水量增加的年代正好是旱季降水量減少的年代。此外，雨季降水量變化幅度大于旱季，但波動有增大趨勢。

由表3看出，從20世紀70年代開始日照時數(shù)呈減少趨勢，變異系數(shù)略微增大，90年代減少幅度最大，相比80年代減少了193.91 h，其變異系數(shù)較小，說明90年代是日照時數(shù)減少穩(wěn)定期。雨季、旱季日照時數(shù)均減少，旱季日照時數(shù)波動明顯大于雨季，且波動有減緩趨勢。

2.2 氣候突變特征分析

?？谑心昶骄鶜鉁赝蛔兊腗-K和滑動T檢驗曲線見圖2。由圖2可知，經(jīng)M-K突變檢驗顯示，UF曲線分別在1967～1978年和1998～2010年超過了α=0.05的臨界線（2000～2010年超過了0.001的顯著水平），且UF與UB在1992～1993年有一個交點，交點位于α=0.05的臨界線范圍內(nèi)，說明其可能是突變點，1963～1978年和1998～2010年是可能的突變時間。為檢驗結果的可靠性，本文先后選取5、10、15 a滑動步長進行滑動T檢驗。5 a滑動檢驗結果為1978、1986和1997年均出現(xiàn)由低到高的突變，1981年出現(xiàn)由高到低的突變；10 a滑動檢驗得出，1971、1978、1986、1997年各出現(xiàn)由低到高的突變；15 a滑動檢驗得出，1986、1990年出現(xiàn)由低到高的突變。根據(jù)累積距平分析得出（圖略），60 a?？谑心昃鶜鉁乩鄯e距平呈現(xiàn)先降后升的“V”字形趨勢，以1986年為轉折點，其中1967～1978年、1998～2010年分別出現(xiàn)了明顯的下降和上升趨勢。綜合分析來看，1986年氣溫發(fā)生了由低向高的突變，1967～1978年、1998～2010年均通過α=0.05的檢驗，證實其為相應冷期和暖期。

?？谑心杲邓客蛔兊腗-K和滑動T檢驗曲線見圖3。由圖3看出，M-K突變檢驗得出?？诮?0 a降水量變化呈下降趨勢，但未達到α=0.05的顯著水平，UF與UB在1953年存在交點。經(jīng)過滑動T檢驗后得知，5 a滑動下的降水量在1976和2002年出現(xiàn)由多到少的突變；10 a滑動在1976年發(fā)生由多到少突變，15 a滑動沒有突變點。這也表明了降水量數(shù)值波動較大。經(jīng)累積距平分析得出，1951～1976年降水量累積距平大部分為正值，而1977～2010年大部分為負值，由此斷定1976年是降水量變化的突變點。綜合分析可知，海口市年降水量在1976年發(fā)生了由多到少的突變。

?？谑心耆照諘r數(shù)突變的M-K和滑動T檢驗曲線見圖4。由圖4可知，M-K突變檢驗中UF與UB在1990年有1個交點，但落在α=0.05的臨界線之外，故1990年是否為突變點需進一步檢驗，1982～1986年、1988～2010年是可能的突變時間。經(jīng)滑動T檢驗可知，5 a滑動在1965、1993年均發(fā)生由多到少的突變；10 a滑動結果與5 a滑動一致；15 a滑動下1969、1993年出現(xiàn)了由多到少突變。根據(jù)累積距平分析得出，1951～1993年日照時數(shù)呈現(xiàn)增加趨勢，1994～2010年轉為減少，1993年為轉折點。綜合分析，海口市年日照時數(shù)在1993年發(fā)生了由多到少的突變，僅1988～2010年突變時間順利通過0.05水平檢驗。

2.3 小波分析

基于?？谑心昶骄鶜鉁亍⒔邓?、日照時數(shù)時間序列，本文進行了連續(xù)Morlet小波變換，分析得到氣溫、降水量和日照時數(shù)的多時間尺度結構（圖5），圖中信號的強弱用小波系數(shù)的大小表示，正值區(qū)表示偏高或增多，負值區(qū)表示偏低或偏少。從各氣候要素的Morlet小波變換系數(shù)實部時頻變化分析來看，年平均氣溫小波系數(shù)等值線在5～10 a時間尺度上較為密集；降水量的小波系數(shù)等值線同樣在5～10 a時間尺度上密集，特別是在2000年以后最為顯著；日照時數(shù)的小波系數(shù)等值線在5～8、25～28 a時間尺度密集。

3 討論與結論

近60 a來?？谑心昶骄鶜鉁仫@著上升，線性擬合氣候傾向率為0.18 ℃/（10 a），氣溫上升速率明顯低于全國的0.25 ℃/（10 a）[12]和其他地區(qū)[3，13-18]，這可能與海口處于低緯度熱帶區(qū)域有關，史培軍等[19]對全國氣候開展區(qū)劃研究認為，近50 a中國氣溫上升速率與緯度分布有關，高緯度地區(qū)增溫大，低緯度地區(qū)增溫小。從年代際來看，20世紀90年代是?？跉夂蛟雠铒@著的穩(wěn)定期，且1998年是最暖的一年，這與任玉國等[12]研究所得出的中國20世紀后半葉最暖的十年及記錄中最暖的一年結論一致，這表明?？跉鉁厣仙厔菖c全國同步。雨季、旱季氣溫變化與年代際氣溫變化趨勢類似，但旱季氣溫增暖明顯大于雨季，這表明?？诙杭竟?jié)氣溫升高的潛力較大。經(jīng)突變分析得出，1986年是氣溫突變點，1967～1978、1998～2010年對應?？诘睦淦诤团?，氣候變化經(jīng)歷了一次先變冷再變暖的過程[20]。?？谑心昶骄鶜鉁夭▌又饾u增大，這造成未來氣溫變化異常風險性增大，極端天氣氣候事件發(fā)生頻率和強度也會越來越高。此外，海口市氣溫近60 a存在5～10 a左右的周期震蕩。

海口市年降水量變化趨勢不明顯，但其波動逐漸增大。任國玉等[12]研究1951～2000年全國降水量幾乎沒有變化；史培軍等[19]對中國1961～2010年氣候變化區(qū)劃中將?？趧潪槿A東-華中濕暖趨勢帶，認為?？谑墙邓龆嗟貐^(qū)，本文若以1961～2010年統(tǒng)計降水量可得到與他一致的結論。年代際變化分析顯示，最近20 a降水量明顯增加，這與中國在90年代以來降水量增多的結論一致[12]，這表明海口降水量增加趨勢與全國趨勢的同步性較好。此外，2000年以后降水量波動最不穩(wěn)定，這是未來降水量變動異常的表現(xiàn)，也可能是未來發(fā)生異常澇災和旱災的表征。雨季降水量變化幅度大于旱季，但波動明顯較小。經(jīng)突變和小波分析得出，?？谑薪邓吭?976年發(fā)生由多到少的突變，存在5～10 a的周期震蕩。

海口市年日照時數(shù)顯著減少，線性擬合氣候傾向率為-83.1 h/10 a，這與虞海燕等[21]研究的59 a來?？趨^(qū)域日照時數(shù)變化速率結論類似。但其減少速率要大于全國平均水平[12]。日照時數(shù)年代際波動表現(xiàn)輕微增大，這同吳名杰等[22]研究?？诮?0 a日照時數(shù)變化特征結論一致，這表明?？谀耆照諘r數(shù)變化趨勢相對平穩(wěn)。此外，旱季日照時數(shù)波動大于雨季。突變和小波分析得出，?？谑心耆照諘r數(shù)在1993年發(fā)生由多到少的突變，且存在5～8 a、25～28 a的周期震蕩。

綜上分析，海口市近60 a來氣候變化特征表現(xiàn)為年平均氣溫升高，年日照時數(shù)減少，年降水量變化趨勢不明顯。20世紀90年代是氣溫和日照時數(shù)顯著變化的穩(wěn)定期，這可能是全國氣候變化大背景和海南建省辦經(jīng)濟特區(qū)（1988年）的共同結果。?？谑凶鳛橐月糜伟l(fā)展為支柱產(chǎn)業(yè)的生態(tài)城市，應樹立起綠色、低碳的發(fā)展理念，在工業(yè)化和城鎮(zhèn)化快速發(fā)展過程中更進一步大力開展節(jié)能減碳和生態(tài)建設，以加強應對和適應氣候變化的能力。

參考文獻

[1] 秦大河， 丁一匯， 蘇紀蘭， 等. 中國氣候與環(huán)境演變評估（I）： 中國氣候與環(huán)境變化及未來趨勢[J]. 氣候變化研究進展， 2005， 1（1）： 4-9.

[2] 沈永平， 王國亞. IPCC第一工作組第五次評估報告對全球氣候變化認知的最新科學要點[J]. 冰川凍土， 2013， 35（5）： 1 068-1 076.

[3] 何云玲， 魯枝海. 近60 a昆明市氣候變化特征分析[J]. 地理科學， 2012， 32（9）： 1 119-1 124.

[4] 翟盤茂， 鄒旭愷. 1951-2003年中國氣溫和降水變化及其對干旱的影響[J]. 氣候變化研究進展， 2005， 1（1）： 16-18.

[5] 葉正偉， 吳 威. 廬山旅游區(qū)氣候變化特征及其影響因素分析[J]. 地理科學， 2011， 31（10）： 1 221-1 227.

[6] 張立偉， 宋春英， 延軍平. 秦嶺南北年極端氣溫的時空變化趨勢研究[J]. 地理科學， 2011， 31（8）： 1 007-1 011.

[7] 任永建， 陳正洪， 肖 鶯， 等. 武漢區(qū)域百年地表氣溫變化趨勢研究[J]. 地理科學， 2010， 30（2）： 278-282.

[8] 李 淼， 夏 軍， 陳社明， 等. 北京地區(qū)近300 a降水變化的小波分析[J]. 自然資源學報， 2011， 26（6）： 1 001-1 011.

[9] 何 萍， 李宏波， 黃 惠. 1960-2009年云南高原楚雄市氣候年代際變化特征及城市氣候分析[J]. 地理科學進展， 2011， 30（1）： 65-72.

[10] 符涂斌， 王 強. 氣候突變的定義和檢測方法[J]. 大氣科學， 1992， 16（4）： 482-493.

[11] 魏鳳英. 現(xiàn)代氣候統(tǒng)計診斷與預測技術（2版）[M]. 北京： 氣象出版社， 2007， 106-113.

[12] 任國玉， 郭 軍， 徐銘志， 等. 近50 a中國地面氣候變化基本特征[J]. 氣象學報， 2005， 63（6）： 942-956.

[13] 任朝霞， 楊達遠. 近50 a西北干旱區(qū)氣候變化趨勢研究[J]. 第四紀研究， 2006， 26（2）： 299-300.

[14] 袁晴雪， 魏文壽. 中國天山山區(qū)近40 a來的年氣候變化[J]. 干旱區(qū)研究， 2006， 23（1）： 115-118.

[15] 滿蘇爾·沙比提， 艾克拜爾·艾爾肯， 迪里夏提·司馬義. 近50年來新疆莎車縣氣候變化特征及其成因[J]. 干旱區(qū)研究， 2009， 26（5）： 621-627.

[16] 夏依木拉提·艾依達爾艾力. 近55 a來新疆伊寧市氣候變化過程分析[J]. 干旱區(qū)研究， 2010， 27（6）： 878-884.

[17] 滿蘇爾·沙比提， 米娜瓦爾·斯拉吉丁， 吐尼沙古麗·牙生. 近40 a來新疆策勒綠洲氣候變化特征分析[J]. 水文， 2009， 29（1）： 89-92.

[18] 李 帥， 李祥余， 何 清， 等. 阿勒泰地區(qū)近40 a的氣候變化研究[J]. 干旱區(qū)研究， 2006， 23（4）： 637-643.

[19] 史培軍， 孫 劭， 汪 明， 等. 中國氣候變化區(qū)劃（1961-2010年）[J]. 中國科學： 地球科學， 2014（10）： 2 294-2 306.

[20] Huang J， Yu H， Guan X， et al. Accelerated dryland expansion under climate change[J]. Nature Climate Change， 2015， 6： 166-171.

[21] 虞海燕， 劉樹華， 趙 娜， 等. 中國近59 a日照時數(shù)變化特征及其與溫度、風速、降水的關系[J]. 氣候與環(huán)境研究， 2011， 16（3）： 389-398.

[22] 吳名杰， 陳小敏. 近50 a?？谑腥照諘r數(shù)氣候變化特征分析[J]. 熱帶農(nóng)業(yè)科學， 2014， 34（1）： 94-97.