顏鵬程1, 2 封國林2, 3 侯威2 吳浩2, 3

?

500 hPa溫度場時間序列的年代際突變過程統(tǒng)計特征

顏鵬程封國林侯威吳浩

1蘭州大學(xué)大氣科學(xué)學(xué)院，蘭州730000;2中國氣象局國家氣候中心，北京100081;3揚(yáng)州大學(xué)物理科學(xué)與技術(shù)學(xué)院，揚(yáng)州225002

本文利用Logistic模型，推導(dǎo)出描述均值突變的分段函數(shù)，用該函數(shù)對可能存在突變的時間序列進(jìn)行擬合，并結(jié)合概率分布理論，確定突變幅度最大的時段即為發(fā)生突變的過程，由此確定并分析序列中突變的開始時刻、突變幅度、突變變率、突變持續(xù)時間、系統(tǒng)不穩(wěn)定特性等參數(shù)。對全球500 hPa溫度場單點(diǎn)時間序列的突變持續(xù)過程展開研究：（1）對突變過程的開始時刻進(jìn)行統(tǒng)計，發(fā)現(xiàn)1956～1959、1970～1979、1986～1994、1994～2004年開始的突變所占比重較大，幾次突變過程中，平均來看除了1986～1994年開始的突變表現(xiàn)為降溫，其余幾次突變均表現(xiàn)為增溫；并且每次突變時，增溫幅度比較大的突變，其變率也較大；（2）從突變先后的空間分布上看，1956～1959、1970～1979年檢測到發(fā)生突變的格點(diǎn)在海洋上空的突變偏早、歐亞大陸上空偏晚，而1986～1994和1994～2004年突變期間這一情況正好相反；（3）針對每次突變過程中的突變幅度，低緯度區(qū)域溫度發(fā)生突變的變化幅度較小、高緯度區(qū)域較大；（4）對突變持續(xù)時間的檢測結(jié)果表明，在全球增暖背景下，完成一次突變的持續(xù)時間正在逐漸變短；（5）當(dāng)系統(tǒng)突變正在進(jìn)行時，系統(tǒng)的不穩(wěn)定性較強(qiáng)。

氣候突變 500 hPa溫度場 突變過程 持續(xù)時間 不穩(wěn)定性

1 引言

IPCC第四次評估報告（IPCC，2007）指出，過去一百年（1906～2005年），全球溫度升高0.56～0.92°C；海平面自1961年以來以平均速度1.8 mm a（1993年以來是3.1 mm a）速度抬升；而北極海冰面積自1993年以來也正以2.7% (10a)的速度消融。與此同時，越來越多有關(guān)地球氣候系統(tǒng)的突變證據(jù)正逐漸被發(fā)現(xiàn)（Overpeck and Cole，2006）并引起了各國決策者的密切關(guān)注（Will and Patrick.，2007）。有證據(jù)表明全球氣候突變正在改變北半球冰蓋的格局（Clark et al.，1999；Post et al.，2009），氣候突變還對海洋生態(tài)系統(tǒng)產(chǎn)生重大影響（Hoegh- Guldberg and Bruno，2010），而全球性的氣候增暖還引起了區(qū)域性干旱程度的增加（符淙斌和馬柱國，2008）。因此開展全球增暖背景下的氣候突變的研究顯得急切而重要（封國林等，2001，2002）。季飛等（2011）對全球環(huán)流場關(guān)聯(lián)分析后，認(rèn)為存在1978～1982和1996～1998兩次明顯的躍變；肖棟和李建平（2007）在對海表溫度（sea surface temperature，SST）的研究中發(fā)現(xiàn)，1867～2005年存在七次較為明顯的突變，并且突變的發(fā)生具有明顯的空間分布特征；嚴(yán)中偉等（1990a，1990b），嚴(yán)中偉（1992）在分析了60年代的500 hPa高度 場、SLP（Sea-Level Pressure）等要素后認(rèn)為突變存在一定的空間分布，并且不同要素突變的發(fā)生有一定的超前滯后特性；江志紅等（2004a，2004b）分析了溫度場的趨勢變化同樣具有一定的時空區(qū)域分布特征；丁一匯和張莉（2008）發(fā)現(xiàn)我國青藏高原與其他地區(qū)氣候突變時間相比具有明顯滯后現(xiàn)象，此外全球溫度突變背景下還引發(fā)了極端事件的頻繁發(fā)生（龔志強(qiáng)等，2009）。

目前學(xué)術(shù)界對于氣候突變尚無統(tǒng)一的定義（Thom，1972；符淙斌和王強(qiáng)，1992；李建平等，1996），加之突變的原因錯綜復(fù)雜，對于時間序列突變點(diǎn)位置的確定也有不同的方法（封國林等，2006，2008，2011），傳統(tǒng)的檢測方法有：低通濾波法、Mann-Kendall（M-K）、滑動?檢驗(yàn)（MTT）、Cramer、Yamamoto等（魏鳳英，1999）檢測方法；隨后一些基于啟發(fā)式分割（BG）算法（封國林等，2005）、復(fù)雜度（侯威等，2005）、去趨勢波動分析（楊萍等，2008；侯威等，2011a，2011b）、滑動移除重標(biāo)極差分析（MC-R/S）（何文平等，2010）、滑動移除近似熵（MC-ApEn）（金紅梅等，2012a，2012b）等發(fā)展的突變檢測新方法被有效地應(yīng)用于不同類型的突變檢測。近幾年，Scheffer et al.（2009），吳浩等（2012，2013）等進(jìn)一步研究了突變的早期預(yù)警信號。然而，對突變的認(rèn)識，絕大部分研究還僅僅認(rèn)為氣候突變是時間序列上的一個“點(diǎn)”，沒有考慮突變事件的發(fā)生、發(fā)展乃至消亡的過程，忽略了突變事件應(yīng)有的過程性。

一個關(guān)于氣候突變的完備性定義（李建平等，1996）描述為：對于某氣候統(tǒng)計量，一個性質(zhì)的氣候狀態(tài)變化到另一個性質(zhì)的氣候狀態(tài)，期間的過渡期遠(yuǎn)小于它們各自狀態(tài)的維持時間，則該統(tǒng)計量發(fā)生突變。注意到氣候系統(tǒng)在不同狀態(tài)間躍變時的“過渡期”是存在的，通過縮小觀察窗口可以捕捉到這一過渡期（突變持續(xù)過程）。對氣候突變過程的研究，尤其對突變幅度、突變變率（單位時間系統(tǒng)變化幅度）以及突變持續(xù)時間在全球區(qū)域分布的研究，有助于認(rèn)識氣候突變的發(fā)生、發(fā)展直至消亡的過程，也將大大有益于氣候突變事件的前兆信號捕捉研究和確定重點(diǎn)監(jiān)測的區(qū)域。本文提出一種對氣候突變過程進(jìn)行分析的方法，該方法從Logistic模型推導(dǎo)出滿足突變的分段函數(shù)，擬合出均值突變過程，并結(jié)合統(tǒng)計概率，對氣候突變的開始時刻、突變幅度、突變變率、突變持續(xù)時間、系統(tǒng)不穩(wěn)定性等參數(shù)進(jìn)行估計，以此為基礎(chǔ)，系統(tǒng)考察了1948～2012年500 hPa溫度場上突變時間序列的時空分布特征及突變過程中的規(guī)律性。

2 資料與方法

2.1 資料

本文使用的是美國國家海洋與大氣管理局（NOAA，http://www.esrl.noaa.gov/）公布的500 hPa月平均溫度場再分析資料，長度為1948年1月至2012年9月，分辨率為2.5°×2.5°，空間格點(diǎn)分布為經(jīng)向144個，緯向73個（其中±90°的2個極點(diǎn)數(shù)值在實(shí)際計算中不參與計算）。

圖1 Logistic模型系統(tǒng)狀態(tài)的演化曲線。虛線將曲線分為三段，其中斜線代表突變過程，分別是過程開始和結(jié)束狀態(tài)

對500 hPa溫度場各個格點(diǎn)上的資料進(jìn)行去除季節(jié)趨勢預(yù)處理（距平化），處理后的數(shù)據(jù)不包含月尺度的季節(jié)震蕩信號。處理方法為序列上的數(shù)據(jù)減去對應(yīng)月份的平均，可用方程描述為：，其中和是預(yù)處理前后的數(shù)據(jù)，角標(biāo)和分別表示年和月，顯然有是序列中第月62年的平均。

2.2 方法介紹

2.2.1 Logistic模型推導(dǎo)下的分段擬合函數(shù)

Logistic模型被用來描述蟲口數(shù)量由一個狀態(tài)擴(kuò)張（收縮）至另一個狀態(tài)的非線性過程（May，1976；劉秉正和彭建華，2004），該模型的類似形式出現(xiàn)在諸多氣候、生物學(xué)模型中（Carpenter and Brock，2006，Guttal and Jayaprakash，2008），可用方程表述如下：

稱之為突變變率，其物理意義在于，系統(tǒng)在兩個狀態(tài)（）之間轉(zhuǎn)換時系統(tǒng)狀態(tài)變量的變化幅度與持續(xù)時間的比，即單位時間系統(tǒng)的變化幅度，突變變率越大表明系統(tǒng)變化幅度越大，持續(xù)時間越短。方程（1）的解析解為

，（3）

將式（3）帶入定義式（2），得到：

據(jù)式（4），則可以將時間序列劃分三段，用函數(shù)表示為

。 （7）

2.2.2 觀察窗口的選取和突變識別標(biāo)準(zhǔn)

由于氣候突變在不同時間尺度上有不同的表現(xiàn)形式（李建平等，1996），需要選取合適的時間長度作為觀察窗口以考察該氣候尺度下突變的發(fā)生情況。設(shè)定原序列長度為，取觀察窗口，長度為，對該觀察窗口下的序列進(jìn)行方法2.2.1的參數(shù)提取，標(biāo)記擬合得到的突變幅度為，并對之進(jìn)行統(tǒng)計。圖2是任意選取格點(diǎn)（10，110）的溫度序列，觀察窗口長度取，統(tǒng)計突變幅度概率分布圖，黑色曲線是高斯擬合結(jié)果，可以看出對于溫度序列的突變幅度近似滿足正態(tài)分布。

圖2 突變幅度概率的統(tǒng)計分布，橫坐標(biāo)是突變幅度大小，縱坐標(biāo)是其統(tǒng)計概率

由于突變事件屬于小概率事件，則可以對圖2所對應(yīng)的概率分布函數(shù)進(jìn)行積分，認(rèn)為落點(diǎn)在積分區(qū)間之外（如圖中虛線標(biāo)記的位置外側(cè)）的屬于突變事件。查正態(tài)分布函數(shù)表（黃嘉佑，1990）發(fā)現(xiàn)當(dāng)＞98%時，即有

成立，此時通過式（8）可以給出判斷突變是否發(fā)生的標(biāo)準(zhǔn)，式（8）中，分別表示突變幅度序列的數(shù)學(xué)期望和標(biāo)準(zhǔn)差。

3 500 hPa溫度場突變分布特征研究

對全球500 hPa溫度場1948年1月到2012年9月的溫度再分析資料，利用第2.2節(jié)方法進(jìn)行氣候突變過程的識別，并估計得到突變開始時刻、突變幅度、突變變率、突變持續(xù)時間等參量，考察500 hPa溫度場全球范圍內(nèi)的溫度突變事件在時間和空間上的分布情況。

3.1 全球格點(diǎn)突變開始時刻統(tǒng)計規(guī)律及空間分布

考察突變持續(xù)時間小于10年的突變，本文取觀察窗口為10年，對每個格點(diǎn)的溫度序列進(jìn)行突變識別，統(tǒng)計所有格點(diǎn)上檢測到的突變開始時刻，其空間分布如圖3所示。在1948年1月到2012年9月，存在1956～1959、1970～1979、1986～1994、1994～2004年四次（圖3中加粗部分）較為明顯的氣候突變（為分析方便，1970～1979和1986～1994年兩次突變沒有進(jìn)一步細(xì)化），這與針對SST和SLP資料（肖棟和李建平，2007；肖棟，2008；Xiao and Li，2007；Xiao et al., 2012）的氣候突變檢測結(jié)果是相一致的。

圖3 全球格點(diǎn)突變開始時刻統(tǒng)計分布，x坐標(biāo)是突變開始時刻，y坐標(biāo)是依據(jù)極端閾值法（公式6）識別出來的突變的統(tǒng)計概率

對500 hPa溫度場的四次突變，圖4分別給出檢測到發(fā)生突變格點(diǎn)的溫度平均序列，并進(jìn)行突變過程擬合。從溫度平均來看，發(fā)現(xiàn)開始于1956～1959、1970～1979、1994～2004的突變前后平均溫度升高，上升幅度分別為0.37、0.16和0.23°C；而開始于1986～1994年的突變，全球溫度降低了0.37°C。從突變變率來看，1986～1994和1994～2004年的兩次突變變率分別為0.220和0.759°C a，明顯要大于1956～1959（0.060°C a）和1970～1979（0.034°C a）的兩次突變。就序列平均溫度而言，1986～1994和1994～2004年的兩次突變，序列平均溫度分別為0.10°C a和0.18°C a，要高于1956～1959（－0.20°C a）和1970～1979（－0.23°C a），表明在增暖背景下的突變變率是增加的。值得注意的是1970～1979和1986～1994年突變過程中，在整體性均值突變的基礎(chǔ)上還存在明顯的轉(zhuǎn)折型突變。

圖4 對四次突變過程平均溫度序列的擬合，橫坐標(biāo)是年份，縱坐標(biāo)是溫度

圖5是前述四次突變開始時刻的空間分布圖，整體上來看中高緯區(qū)域突變開始時刻偏早，低緯度區(qū)域突變開始時刻偏晚。

1956～1959年開始的突變，海洋上發(fā)生偏早的區(qū)域主要在北太平洋中北部、南太平洋中東部區(qū)域、印度洋西南部、大西洋西南部以及北冰洋西半球（為便于描述，本文涉及下墊面區(qū)域，均指對應(yīng)的500 hPa溫度層上的情況），陸地上主要以在南美洲中東部、北美洲北部部分區(qū)域以及亞洲中東部和歐洲西部偏早（1956～1957年）。發(fā)生突變偏晚（1958～1960年）的區(qū)域主要在海洋上，北太平洋東部、太平洋中東部、大西洋中部以及北冰洋東半球。

1970～1979年開始的突變，海洋上北太平洋東北部和中部、南太平洋南部、印度洋南部北冰洋西半球突變發(fā)生均偏早（1970～1974年），突變發(fā)生偏晚區(qū)域集中在太平洋中東部、印度洋北部和陸地上的非洲中部、南美洲北部、歐亞大陸部分區(qū)域突變發(fā)生偏晚（1976～1979年）。

1986～1994年開始的突變，突變偏早（1986～1989年）區(qū)域集中在北太平洋中部、印度洋東南部、北冰洋、北大西洋中部區(qū)域，陸地是主要在北美洲西北部、亞洲東南部以及非洲東南部區(qū)域。突變發(fā)生偏晚（1990～1994年）區(qū)域主要以南太平洋、印度洋南部、南大西洋和歐洲中部區(qū)域。

1994～2004年開始的突變，北太平洋中部區(qū)域、印度洋中南部，亞洲大陸、南美洲南部區(qū)域突變發(fā)生偏早（1994～1998年），而太平洋中東部、南太平洋，大西洋中部區(qū)域突變發(fā)生偏晚（2000～2004年）。

總體上來看，1979年之前的突變，高緯度區(qū)域、海洋上空的突變發(fā)生偏早，低緯度區(qū)域、亞歐大陸上空偏晚；而1979年之后的突變，東半球偏早一些，西半球太平洋上空突變偏晚。

圖5 不同時段的突變開始時刻空間分布：（a）1956～1959年；（b）1970～1979年；（c）1986～1994年；（d）1994～2004年。深藍(lán)色區(qū)域表示突變開始年份偏早，橙色區(qū)域表示突變開始年份稍晚，青綠區(qū)域突變開始年份則偏晚

圖6 不同突變時期溫度突變幅度的空間分布：（a）1956～1959年；（b）1970～1979年；（c）1986～1994年；（d）1994～2004年。品紅色表示增溫幅度超過3.5°C的區(qū)域，青綠色表示增溫幅度在3.5°C以內(nèi)的區(qū)域，橙色表示降溫幅度在3.5°C以內(nèi)的區(qū)域，深藍(lán)色表示降溫幅度超過3.5℃的區(qū)域

3.2 不同突變時期溫度突變幅度的分布規(guī)律

1956～1959年突變，全球降溫幅度較大區(qū)域集中北太平洋北部、北美大陸南太平洋南部、北極、南太平洋南部等區(qū)域；降溫幅度較小的區(qū)域集中在太平洋東海岸和西海岸，大西洋以及印度洋南部區(qū)域。增溫幅度較大區(qū)域集中在歐洲西部和南美洲南部區(qū)域；增溫幅度較小的區(qū)域在北美洲南部、北太平洋中部、南太平洋東南部、澳大利亞南部等區(qū)域。

1970～1979年突變?nèi)蚪禍胤容^大區(qū)域較少，主要集中在歐洲西北部、北太平洋、北美洲東北部、大西洋北部以及南極區(qū)域；降溫幅度較小的區(qū)域集中在非洲北部、印度洋中部、太平洋南部以及大西洋南部等區(qū)域。而增溫幅度較大的區(qū)域集中在歐亞、北美大陸、北太平洋以及南太平洋南部區(qū)域和南極。增溫幅度較小的區(qū)域在太平洋、大西洋和印度洋北部。

1986～1994年突變幅度較大區(qū)域集中在歐洲、太平洋北部區(qū)域；降溫幅度較小的區(qū)域集中在太平洋西海岸赤道區(qū)域、印度洋中部以及南太平洋南部等區(qū)域。增溫幅度較大區(qū)域主要集中在北美洲、北冰洋區(qū)域以及南極部分區(qū)域；增溫幅度較小的區(qū)域在太平洋中東部，大西洋以及印度洋部分區(qū)域。

1994～2004年突變幅度較大區(qū)域集中在亞洲中北部、美洲東北部以及南太平洋西南部等部分區(qū)域；降溫幅度較小的區(qū)域集中在太平洋中部區(qū)域以及南太平洋東南部。增溫幅度較大區(qū)域在主要集中在南極和北半球高緯度區(qū)域；增溫幅度較小的區(qū)域在太平洋東西海岸、印度洋和大西洋，以及南美洲北部和非洲東北部等部分區(qū)域。

整體來看，高緯度區(qū)域增降溫幅度明顯要高于低緯度區(qū)域，1956～1959年突變期間大部分區(qū)域表現(xiàn)為降溫，且南北緯40°以外區(qū)域的增降溫幅度超過3.7°C；1970～1979年（北緯30°以北及南緯40°以南區(qū)域的增降溫幅度超過4°C）、1986～1994年（北緯30°以北及南緯50°以南的增降幅度超過4°C）、1994～2004年（北緯30°以北及南緯50°以南的增降幅度超過3.5°C）突變增降溫幅度較大的格點(diǎn)向低緯度靠近。

3.3 突變持續(xù)時間的統(tǒng)計規(guī)律

式（5）給出了均值突變序列突變過程中的持續(xù)時間，一次完整的突變過程，應(yīng)該包含突變前后系統(tǒng)狀態(tài)的保持時間和完整的突變持續(xù)時間，對于持續(xù)過程在10年內(nèi)的突變，全球檢測得到的突變持續(xù)時間統(tǒng)計結(jié)果如圖7所示。

圖7 不同突變時期各格點(diǎn)溫度發(fā)生突變持續(xù)時間統(tǒng)計分布，橫坐標(biāo)是突變持續(xù)時間，縱坐標(biāo)是統(tǒng)計概率：（a）1956～1959年；（b）1970～1979年；（c）1986～1994年；（d）1994～2004年

對四次突變的持續(xù)時間進(jìn)行統(tǒng)計，其均在60個月前后出現(xiàn)極小值，而1970～1979、1986～1994和1994～2004年突變，在30個月處也出現(xiàn)極小值，據(jù)此可以將突變分為三種類型，一是突變持續(xù)時間不大于30個月的突變，定義為A類突變；二是大于30個月但不大于60個月的突變，定義為B類突 變；對于大于60個月的定義為C類突變。統(tǒng)計三類突變的持續(xù)時間，如表1，屬于C類突變的格點(diǎn)數(shù)目除了1970～1979年期間的突變53.04%外，其余的均約為60%，認(rèn)為這一類型的突變沒有變化；屬于A類突變的格點(diǎn)數(shù)目，則持續(xù)增加，在1970～1979年期間，這一類型的格點(diǎn)突增至29.42%，表明在全球增暖背景下，越來越多的格點(diǎn)上完成一次突變的持續(xù)時間在縮短，持續(xù)時間較長的格點(diǎn)數(shù)基本保持不變。

表1 突變過程的持續(xù)時間分類統(tǒng)計

3.4 突變檢測過程中不穩(wěn)定性探討

表2 不穩(wěn)定性參數(shù)隨突變幅度參數(shù)的統(tǒng)計規(guī)律

Table 2 Statistics distribution of the unstable parameter with the abrupt change amplitude parameter changing

表2 不穩(wěn)定性參數(shù)隨突變幅度參數(shù)的統(tǒng)計規(guī)律

突變時段＞0增溫情況＜0降溫情況 ＞0情況突變占總突變百分比＜0情況突變占總突變百分比＞0情況突變占總突變百分比＜0情況突變占總突變百分比 1956～1959年17.4%29.9%20.2%32.4% 1970～1979年30.6%38.0%19.3%12.1% 1986～1994年15.9%42.7%24.0%17.3% 1994～2004年30.1%38.0%20.3%11.6%

圖8 不同突變過程下的參數(shù)的統(tǒng)計分布，橫坐標(biāo)是參數(shù)取值，縱坐標(biāo)是其概率

圖9 參數(shù)與始—末狀態(tài)，x軸突變過程開始狀態(tài)，y軸是結(jié)束狀態(tài)

4 結(jié)論

利用Logistic模型推導(dǎo)出分段參數(shù)方程，對500 hPa溫度場各格點(diǎn)溫度序列進(jìn)行突變參數(shù)提取，從而實(shí)現(xiàn)對突變的過程性分析。分析結(jié)果發(fā)現(xiàn)自1948～2012年，500 hPa溫度場出現(xiàn)四次明顯的突變過程：1956～1959、1970～1979、1986～1994、1994～2004年。對四次突變過程分析，得到如下結(jié)論：

1）從平均溫度序列來看，1956～1959年全球平均溫度距平－0.20°C，1970～1979年－0.23°C，而1986～1994年和1994～2004年則分別為0.10°C和0.18°C。發(fā)生突變期間，1986～1994年突變期間，全球溫度降低0.33°C；1956～1959、1970～1979 和1994～2004年突變期間，全球增溫分別為：0.37、0.16和0.23°C。

2）從突變早晚來看1956～1959、1970～1979年突變發(fā)生偏早區(qū)域主要在北半球的北太平洋中北部、歐洲西北部和美洲東北部上空；南半球的南太平洋、印度洋南部以及南美洲南部上空。而1986～1994和1994～2004年突變，北半球歐亞大陸上空突變發(fā)生時間偏早、北太平洋相對偏晚；南半球南太平洋上空同樣突變偏晚。

3）從全球溫度突變幅度來看，在大約以南北緯40°附近分界，出現(xiàn)低緯度區(qū)域溫度變化幅度較小、高緯度區(qū)域較大的明顯分界線。并且四次突變，這一分界線有逐漸向低緯度靠近的趨勢。

4）從突變持續(xù)時間來看，約20%的格點(diǎn)完成一次突變持續(xù)時間小于30個月，近50%格點(diǎn)突變持續(xù)時間不超過60個月。屬于A類突變格點(diǎn)（突變持續(xù)時間少于30個月）數(shù)目的比重也由1956～1959年的19.71%，增加到1986～1994年的22.63%和1994～2004年的25.75%（其中1970～1979年為29.42%）。

5）對突變過程中的系統(tǒng)不穩(wěn)定性分析結(jié)果表明，檢測過程中不少突變正處于平衡態(tài)向增/減少或者正在增加/減少至平衡的狀態(tài)，在此期間系統(tǒng)的不穩(wěn)定性明顯較強(qiáng)。

可以看出，對于較多文獻(xiàn)中分析認(rèn)為的七十年代突變前后，500 hPa溫度場的分布型和序列的突變演化都有了較大的結(jié)構(gòu)性調(diào)整，具體表現(xiàn)為七十年代末之前的兩次突變和之后的兩次突變，在突變早晚區(qū)域的分布型由太平洋偏早、歐亞偏晚，轉(zhuǎn)變?yōu)闅W亞偏早、太平洋偏晚。同時在突變持續(xù)時間上，在全球增暖背景下，越來越多的格點(diǎn)完成一次突變的持續(xù)時間縮短。在進(jìn)一步的研究工作中將對更多溫度場的數(shù)據(jù)進(jìn)行突變過程性分析，以期找到溫度序列的突變事件在垂直方向上的演化結(jié)構(gòu)。

Carpenter S R, Brock W A. 2006. Rising variance: A leading indicator ecological transition [J]. Ecology Letters, 9: 311–318.

Clark P U, Alley R B, Pollard D. 1999. Northern Hemisphere Ice-Sheet influences on global climate change [J]. Science, 286: 1104–1111.

丁一匯, 張莉. 2008. 青藏高原與中國其他地區(qū)氣候突變時間的比較 [J]. 大氣科學(xué), 32 (4): 794–805. Ding Yihui, Zhang Li. 2008. Intercomparison of the time for climate abrupt change between the Tibetan Plateau and other regions in China [J]. Chinese Journal of Atmospheric Sciences (in Chinese), 32 (4): 794–805.

董文杰, 封國林, 龔志強(qiáng), 等. 2005. 基于啟發(fā)式分割算法的氣候突變檢測研究 [J]. 物理學(xué)報, 54 (11): 5494–5499. Dong Wenjie, Feng Guolin, Gong Zhiqiang, et al. 2005. Abrupt climate change detection based on heuristic segmentation algorithm [J]. Acta Phys. Sin. (in Chinese), 54 (11): 5494–5499.

封國林, 戴新剛, 王愛慧, 等. 2001. 混沌系統(tǒng)中可預(yù)報性的研究 [J]. 物理學(xué)報, 50 (4): 606–611. Feng Guolin, Dai Xin’gang, Wang Aihui, et al. 2001. On numerical predictability in the chaos system [J]. Acta Phys. Sin. (in Chinese), 50 (4): 606–611.

封國林, 董文杰, 賈曉靜, 等. 2002. ?！?dú)庹袷幾又械臉O限環(huán)解 [J]. 物理學(xué)報, 51 (6): 1 181–1185. Feng Guolin, Dong Wenjie, Jia Xiaojing, et al. 2002. On the dynamics behaviour and instability evolution of air–sea Oscillator [J]. Acta Phys. Sin. (in Chinese), 51 (6): 1181–1185.

封國林, 董文杰, 龔志強(qiáng), 等. 2006. 觀測數(shù)據(jù)非線性時空分布理論和方法 [M]. 北京: 氣象出版社, 1–227. Feng Guolin, Dong Wenjie, Gong Zhiqiang, et al. 2006. Nonlinear theories and methods on spatial-temporal distribution of the observational data (in Chinese) [M]. Beijing: China Meteorological Press 1–227.

封國林, 龔志強(qiáng), 支蓉. 2008. 氣候變化檢測與診斷技術(shù)的若干新進(jìn)展 [J]. 氣象學(xué)報, 66 (6): 892–905. Feng Guolin, Gong Zhiqiang, Zhi Rong. 2008. Latest advances of climate change detecting technologies [J]. Acta Meteorologica Sinica (in Chinese), 66 (6): 892–905.

封國林, 侯威, 支蓉, 等. 2011. 極端氣候事件的檢測、診斷和可預(yù)測性研究 [M]. 北京: 科學(xué)出版社, 350pp. Feng Guolin, Hou Wei, Zhi Rong, et al. 2011. Research on Detecting, Diagnosing and Predictability of Extreme Climate Events (in Chinese) [M]. Beijing: Science Press, 350pp.

符淙斌, 王強(qiáng). 1992. 氣候突變的定義和檢測方法 [J]. 大氣科學(xué), 16 (4): 482–493. Fu Congbin, Wang Qiang. 1992. The definition and detection of the abrupt climatic change [J]. Scientia Atmospherica Sinica 16(4): 482–493.

符淙斌, 馬柱國. 2008. 全球變化與區(qū)域干旱化 [J]. 大氣科學(xué), 32 (4): 752–760. Fu Congbin, Ma Zhuguo. 2008. Global change and regional aridification [J]. Chinese Journal of Atmospheric Sciences (in Chinese), 32 (4): 752–760.

龔志強(qiáng), 王曉娟, 支蓉, 等. 2009. 中國近58年溫度極端事件的區(qū)域特征及其與氣候突變的聯(lián)系 [J]. 物理學(xué)報, 58 (6): 4342–4352. Gong Zhiqiang, Wang Xiaojuan, Zhi Rong, et al. 2009. Regional characteristics of temperature changes in China during the past 58 years and its probable correlation with abrupt temperature change [J]. Acta Physica Sinica (in Chinese), 58 (6): 4342–4352.

Guttal V, Jayaprakash C. 2008. Changing skewness: An early warning signal of regime shifts in ecosystems [J]. Ecology Letters, 11: 450–460.

何文平, 吳瓊, 鄧北勝, 等. 2010. 一種基于重標(biāo)極差方法的動力學(xué)結(jié)構(gòu)突變檢測新方法[J]. 物理學(xué)報, 59 (11): 8264–8271. He Wenping, Wu Qiong, Deng Beisheng, et al. 2010. A new method of detecting abrupt dynamic change based rescaled range analysis [J]. Acta Phys. Sin. (in Chinese), 59 (11): 8264–8271.

侯威, 封國林, 董文杰. 2005. 基于復(fù)雜度分析logistic映射和Lorenz模型的研究 [J]. 物理學(xué)報, 54 (8): 3940–3946. Hou Wei, Feng Guolin, Dong Wenjie. 2005. Investigation about the Lorenz model and Logistic equation based on the complexity [J]. Acta Phys. Sin. (in Chinese), 54 (8): 3940–3946.

侯威, 章大全, 周云, 等. 2011a. 一種確定極端事件閾值的新方法: 隨機(jī)重排去趨勢波動分析方法 [J]. 物理學(xué)報, 60 (10): 109202. Hou Wei, Zhang Daquan, Zhou Yun, et al. 2011a. Stochastic re-sorting detrended fluctuation analysis: A new method to define the threshold of extreme event [J]. Acta Physica Sinica (in Chinese), 60 (10): 109202.

侯威, 章大全, 錢忠華, 等. 2011b. 基于隨機(jī)重排去趨勢波動分析的極端低溫事件研究及其綜合指標(biāo)的建立 [J]. 物理學(xué)報, 60 (10): 109203. Hou Wei, Zhang Daquan, Qian Zhonghua, et al. 2011b. Extremely low temperature and its composite index based on stochastically re-sorting detrended fluctuation analysis [J]. Acta Physica Sinica (in Chinese), 60 (10): 109203.

Hoegh-Guldberg O, Bruno J F. 2010. The impact of climate change on the world’s marine ecosystems [J] Science, 328: 1523–1528.

黃嘉佑. 1990. 氣象統(tǒng)計分析與預(yù)報方法[M]. 北京: 氣象出版社, 28–30, 130–139. Huang Jiayou. 1990. Meteorological Statistical Analysis and Prediction (in Chinese) [M]. Beijing: China Meteorological Press, 28–30, 130–139.

IPCC, 2007. Climate Change. 2007. Synthesis Report [M]. Contribution of Working Groups I, II and III to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge, UK: Cambridge University Press, 2007.

季飛, 支蓉, 龔志強(qiáng), 等. 2011. 全球環(huán)流系統(tǒng)關(guān)聯(lián)性的時空演化特征研究 [J]. 氣象學(xué)報, 69 (6): 1038–1045. Ji Fei, Zhi Rong, Gong Zhiqiang, et al. 2011. Analysis of temporal and spatial characteristics of correlations between the global circulation systems [J]. Acta Meteorologica Sinica (in Chinese), 69 (6): 1038–1045.

金紅梅, 何文平, 侯威, 等. 2012a. 不同趨勢對滑動移除近似熵的影響 [J]. 物理學(xué)報, 61 (6): 069201. Jin Hongmei, He Wenping, Hou Wei, et al. 2012a. Effects of different trends on moving cut data-approximate entropy [J]. Acta Phys. Sin. (in Chinese), 61 (6): 069201 (in Chinese)

金紅梅, 何文平, 張文, 等. 2012b. 噪聲對滑動移除近似熵的影響 [J]. 物理學(xué)報, 61 (12): 129202. Jin Hongmei, He Wenping, Zhang Wen, et al. 2012b. Effect of noises on moving cut data—Approximate entropy [J]. Acta Phys. Sin. (in Chinese), 61 (12): 129202.

江志紅, 李建平, 王梅華, 等. 2004a. 20世紀(jì)全球溫度場趨勢變化的區(qū)域特征分析 [J]. 氣候與環(huán)境研究, 9 (3): 422–434. Jiang Zhihong, Li Jianping, Wang Meihua, et al. 2004. Regional characteristics of the trend change for global temperature field during the last century [J]. Climatic and Environmental Research (in Chinese), 9 (3): 422–434.

江志紅, 李建平, 屠其璞, 等. 2004b. 20世紀(jì)全球溫度年代和年代際變化的區(qū)域特征 [J]. 大氣科學(xué), 28 (4): 545–558. Jiang Zhihong, Li Jianping, Tu Qipu, et al. 2004. Regional characteristics of the decadal and interdecadal variations for global temperature field during the last century [J]. Chinese Journal of Atmospheric Sciences (in Chinese), 28 (4): 545–558.

李建平, 丑紀(jì)范, 史久恩. 1996. 氣候突變的完備定義和類型 [J]. 北京氣象學(xué)院, (1): 7?12. Li Jianping, Chou Jifan, Shi Jiuen. 1996. Complete definition and types of abrupt climate change [J]. Journal of Beijing Meteorological College (in Chinese), (1): 7?12.

劉秉正, 彭建華. 2004. 非線性動力學(xué) [M]. 北京: 高等教育出版社, 215–218. Liu Bingzheng, Peng Jianhua. 2004. Nonlinear Dynamics (in Chinese) [M]. Beijing: Higher Education Press, 215–218.

May R M. 1976. Simple mathematical models with very complicated dynamics [J]. Nature (in Chinese), 261: 459–467.

Overpeck J T, Cole J E. 2006. Abrupt change in earth’s climate system [J]. Annu. Rev. Environ. Resour., 31: 1–31.

Will M, Patrick M. 2007. Why is it so difficult? Trade liberalization under the Doha Agenda [J]. Oxford Review of Economic Policy, 23, 347–366

Post E, Forchhammer M C, Bret-Harte M S, et al. 2009. Ecological dynamics across the arctic associated with recent climate change [J]. Science, 325: 1355–1358.

Scheffer M, Bascompte J, Brock W A, et al. 2009. Early-warning signals for critical transitions [J]. Nature, 461: 53–59.

雷內(nèi)·托姆. 1972. 結(jié)構(gòu)穩(wěn)定性與形態(tài)發(fā)生學(xué) [M]. 成都: 四川教育出版社. Thom R. 1972. Stability Structural and Morphogenesis (in Chinese) [M]. Chengdu: Sichuan Education Press.

魏鳳英. 1999. 現(xiàn)代氣候統(tǒng)計診斷與預(yù)測技術(shù) [M]. 北京: 氣象出版社, 62–76. Wei Fengying. 1999. Modern Climatic Statistical Diagnosis and Forecasting Technology (in Chinese) [M]. Beijing: China Meteorological Press, 62–76.

吳浩, 侯威, 錢忠華, 等. 2012. 基于氣候變化綜合指數(shù)的中國近50年來氣候變化敏感性研究[J]. 物理學(xué)報, 61 (14): 149205. Wu Hao, Hou Wei, Qian Zhonghua, et al. 2012. The research on the sensitivity of climate change in China in recent 50 years based on composite index [J]. Acta Physica Sinica (in Chinese), 61 (14): 149205.

吳浩, 侯威, 顏鵬程. 2013. 試用臨界慢化原理探討氣候突變 [J]. 物理學(xué)報, 62 (3): 039206. Wu Hao, Hou Wei, Yan Pengcheng. 2013. Try to use the principle of critical slowing down to discuss the abrupt climate change [J]. Acta Physica Sinica (in Chinese), 62 (3): 039206.

肖棟. 2008. 全球海氣系統(tǒng)年代際突變的研究 [D]. 蘭州大學(xué)大氣科學(xué)學(xué)院博士學(xué)位論文. Xiao Dong. 2008. Study of the decadal abrupt changes of global ocean-atmosphere system [D]. Ph. D. dissertation (in Chinese), College of Atmospheric Science, Lanzhou University.

肖棟, 李建平. 2007. 全球海表溫度場中主要的年代際突變及其模態(tài) [J]. 大氣科學(xué), 31 (5): 839–854. Xiao Dong, Li Jianping. 2007. Main decadal abrupt changes and decadal modes in global sea surface temperature field [J]. Chinese Journal of Atmospheric Sciences (in Chinese), 31 (5): 839–854.

Xiao D, Li J P. 2007. Spatial and temporal characteristics of the decadal abrupt changes of global atmosphere-ocean system in the 1970s [J]. J. Geophys. Res. (in Chinese), 112, D24S22, doi: 10. 1029/2007JD008956.

Xiao D, Li J P, Zhao P, 2012. Four-dimensional structures and physical process of the decadal abrupt changes of the northern extratropical ocean-atmosphere system in the 1980s [J]. Int. J. Climatol. (in Chinese), 32: 983–994, doi:10.1002/joc.2326.

嚴(yán)中偉, 季勁鈞, 葉篤正. 1990a. 60年代北半球夏季氣候躍變——I. 降水和溫度變化 [J]. 中國科學(xué)(B輯), 20:97–103. Yan Zhongwei, Ji Jinjun, Ye Duzheng. 1990a. Northern Hemispheric summer climatic jump in the 1960s, Part I—Rainfall and temperature [J]. Science in China Series B (in Chinese), 20: 97–103.

嚴(yán)中偉, 季勁鈞, 葉篤正. 1990b. 60年代北半球夏季氣候躍變——II. 海平面氣壓和500 hPa高度變化 [J]. 中國科學(xué)(B輯), 20: 879–885. Yan Zhongwei, Ji Jinjun, Ye Duzheng. 1990b. Northern Hemispheric summer climatic jump in the 1960s, Part II—Sea level pressure and 500 hPa height[J]. Science in China Series B (in Chinese), 20: 879–885.

嚴(yán)中偉. 1992. 60年代北半球夏季氣候躍變過程的初步分析 [J]. 大氣科學(xué), 16 (1): 111–119. Yan Zhongwei. 1992. A primary analysis of the process of the 1960s northern hemispheric summer climatic jump [J]. Chinese Journal of Atmospheric Sciences (in Chinese), 16 (1): 111–119.

顏鵬程, 侯威, 胡經(jīng)國. 2012. 基于Logistic模型的均值突變時間序列臨界預(yù)警研究 [J]. 物理學(xué)報, 61 (18): 189202. Yan Pengcheng, Hou Wei, Hu Jingguo. 2012. The critical warning research of the mean time series mutations based on Logistic model [J]. Acta Phys. Sin. (in Chinese), 61 (18): 189202.

閆冠華, 顏鵬程, 侯威, 等. 2012. 一種基于Logistic模型的過程性突變分析方法及其應(yīng)用[J]. 物理學(xué)報, 62 (7): 079202. Yan Guanhua, Yan Pengcheng, Hou Wei, et al. 2012. A method of abrupt change process analysis based on Logistic model and its applications [J]. Acta Phys. Sin. (in Chinese), 62 (7): 079202.

楊萍, 侯威, 封國林. 2008. 基于去趨勢波動分析方法確定極端事件閾值[J]. 物理學(xué)報, 57 (8): 5333–5342. Yang Ping, Hou Wei, Feng Guolin. 2008. Determining the threshold of extreme events with detrended fluctuation analysis [J]. Acta Phys. Sin. (in Chinese), 57 (8): 5333–5342.

顏鵬程，封國林，侯威，等. 2014. 500 hPa溫度場時間序列的年代際突變過程統(tǒng)計特征[J]. 大氣科學(xué), 38 (5): 861?873, doi:10.3878/j.issn.1006-9895. 1403.13106. Yan Pengcheng, Feng Guolin, Hou Wei, et al. 2014. Statistical characteristics on decadal abrupt change process of time sequence in 500 hPa temperature field [J]. Chinese Journal of Atmospheric Sciences (in Chinese), 38 (5): 861?873.

Statistical Characteristics on Decadal Abrupt Change Process of Time Sequence in 500 hPa Temperature Field

YAN Pengcheng, FENG Guolin, HOU Wei, and WU Hao

1,7300002,,100081,,225002

In this paper, we fit the time series, which likely has abrupt change processes (ACP), with piecewise function deducing for describing the mean abrupt change from the logistic model. Thus combined with the theory of probability distribution formed the basis for considering the process in which abrupt change amplitude is maximum, as is the ACP. The parameters reflecting the beginning moment of the abrupt change, abrupt change amplitude, abrupt change rate, persistence time, and instability are determined and analyzed. By testing the time series in a 500 hPa temperature field during the period 1948–2012, we determine the following results: (1) Abrupt changes started in 1956–1959, 1970–1979, 1986–1994, and 1994–2004. We considered more probabilities by testing the start moment of the ACP. In these changes, the temperatures of all changes in the means show increases except for that in 1986–1994, and the rate increases if the abrupt change amplitude is large. (2) In spatial distribution of the abrupt change moment, the moment of the grid points in which the abrupt change occurred in 1956–1959 and 1970–1979 above the sea occurred earlier than that in the grid points above Eurasia, whereas the abrupt changes in 1986–1994 and in 1994–2004 show opposite behavior. (3) The abrupt change amplitude of temperature occurring in low latitudes is lower than that occurring in high latitudes. (4) Under the context of global warming, the statistics distribution of the persistence time shows that an increasing number of grid points require shore time to complete an abrupt change. (5) During the abrupt change process, the instability of the system was enhanced significantly.

Abrupt climate change, 500 hPa temperature field, Abrupt change process, Persist time, Instability

1006–9895(2014)05–0861–13

P467

A

10.3878/j.issn.1006-9895.1403.13106

2013?01–10，2014–03–05收修定稿

國家自然科學(xué)基金項目41175067，國家重點(diǎn)基礎(chǔ)研究發(fā)展計劃項目2012CB955901、2013CB430204

顏鵬程，男，1987年出生，博士研究生，主要從事氣候突變早期預(yù)警信號研究。E-mail: pch.yem@gmail.com

封國林，E-mail: fenggl@cma.gov.cn