李忠賢，陳海山，曾剛，倪東鴻

(南京信息工程大學氣象災(zāi)害省部共建教育部重點實驗室，江蘇南京210044)

0 引言

大氣是一個多時空尺度的系統(tǒng)，這種特性決定了大氣的可預(yù)報性具有較強的時空依賴性。大氣運動表現(xiàn)為一個混沌系統(tǒng)，由于其內(nèi)部動力不穩(wěn)定和非線性相互作用抑制了對大氣環(huán)流的確定性預(yù)報，隨之出現(xiàn)了可預(yù)報性極限，導致天氣預(yù)報時效只有兩周左右(Lorenz，1982)。李建平和丑紀范(2003)指出大氣的可預(yù)報性遵循單調(diào)性原則，即在相同的初始特征和外源強迫特征條件下，時空尺度較大的系統(tǒng)具有較大的可預(yù)報性。季節(jié)—年際時間尺度的大氣環(huán)流變化主要有兩個方面的原因:一是由大氣內(nèi)部的動力不穩(wěn)定及其非線性相互作用所產(chǎn)生的影響，它在一定意義上也可歸結(jié)為初始場的影響;二是由氣候系統(tǒng)的其他子系統(tǒng)所產(chǎn)生的影響，如海溫、海冰、積雪、土壤濕度等異常下邊界條件的強迫作用。因此，氣候或大氣的可預(yù)報性是由大氣內(nèi)部動力過程和緩變的邊界強迫共同決定的。緩變的邊界強迫異常將產(chǎn)生潛在的可預(yù)報性，這構(gòu)成了短期氣候預(yù)測的科學基礎(chǔ)(李崇銀，2000;Collins and Allen，2002)。

中國地處東亞季風區(qū)，受東亞季風年際變化影響，氣候災(zāi)害發(fā)生頻率較高，旱、澇災(zāi)害極大地影響中國的經(jīng)濟建設(shè)和社會發(fā)展，對月、季時間尺度旱、澇氣候，尤其是汛期旱、澇趨勢的預(yù)測是中國氣象學家的一項重要課題(丑紀范和徐明，2001)。氣候的可預(yù)報性一直是氣象學家研究的核心問題。隨著人們對短期氣候預(yù)測需求的增強，關(guān)于短期氣候預(yù)測可預(yù)報性問題的研究也在不斷深入，可預(yù)報性研究不僅可以提高對氣候和氣候變化特征的認識，還有助于完善氣候模式，改進氣候預(yù)測(王會軍等，2007)。

短期氣候預(yù)測的可預(yù)報性表示較長時間尺度里大氣可預(yù)報的程度，即可預(yù)報的氣候信號超出不可預(yù)報的氣候噪音的程度。許多學者(Madden，1976;Shukla and Gutzler，1983;Trenberth，1984;樂群等，1999;柳艷菊等，2000)利用觀測資料(或者再分析資料)分析了實際大氣月、季平均的氣候可預(yù)報性問題。但是，利用大氣實際資料所得到的較長時間尺度可預(yù)報性實際上包含了依賴于初值條件的動力學可預(yù)報性和外界強迫影響下的可預(yù)報性兩個方面(李建平和丁瑞強，2008)，將它們嚴格地分離開是不可能的，因此，迄今還沒有一種完全嚴格的方法可以用來估計實際大氣的氣候噪音和氣候信號(趙曉川等，2008)。

隨著科學的發(fā)展和社會需求的增長，短期氣候數(shù)值預(yù)測越來越引起人們的關(guān)注(丑紀范和徐明，2001)，其中應(yīng)用氣候模式分析可預(yù)報性，已成為氣候研究非?；钴S的領(lǐng)域(趙彥等，2000;范曉青等，2003;Zhou et al.，2009;Zhou and Zou，2010)。國內(nèi)外學者往往利用數(shù)值模擬的辦法進行集合試驗，分析模式大氣對給定的外強迫源的響應(yīng)(Zwiers，1996)。由于海溫通常被認為是最主要的大氣下邊界強迫，因此許多學者(王會軍等，1997;Kusunoki et al.，2001;Quan et al.，2004;Lang and Wang，2005)利用觀測海溫資料強迫大氣環(huán)流模式(AGCM)來研究月、季平均的氣候可預(yù)報性，結(jié)果表明，短期氣候的可預(yù)報性隨地理位置和季節(jié)不同而不同，可預(yù)報性在熱帶地區(qū)高于熱帶外地區(qū)，冬半球的可預(yù)報性高于夏半球。施洪波等(2008)利用參加第2次季節(jié)預(yù)測模式比較計劃的5個大氣環(huán)流模式的輸出結(jié)果，討論了大氣環(huán)流模式在季節(jié)尺度上對亞洲夏季風的可預(yù)報性，研究表明，多數(shù)模式可以大致再現(xiàn)東亞、南亞和西太平洋夏季降水異常EOF1模態(tài)及其對應(yīng)的時間系數(shù)。趙彥等(2000)利用中國科學院大氣物理研究所兩層大氣環(huán)流模式IAP L2 AGCM 1-1進行17 a(1980—1996年)、每年9個單個積分的集合后報試驗，結(jié)果表明，在中國區(qū)域，春季和夏季的降水場、海平面氣壓場以及表面氣溫場的可預(yù)報性從中國南海向西北遞減。此外，研究指出模式大氣的短期氣候可預(yù)報性與ENSO事件存在密切聯(lián)系(Zhou et al.，2009)，與正常年份相比，模式在強El Ni～no年和強La Ni～na年表現(xiàn)出較高的可預(yù)報性(Brankovic and Palmer，1997;Yang et al.，1998;Hassan et al.，2004)。

綜上所述，國內(nèi)外許多學者研究了短期氣候預(yù)測的可預(yù)報性問題，但針對夏季東亞地區(qū)大氣環(huán)流潛在可預(yù)報性的研究還較為缺乏。本文主要利用NCAR CAM3模式分析海溫強迫作用下東亞地區(qū)夏季大氣環(huán)流的潛在可預(yù)報性，并探討其與夏季海溫異常的關(guān)系。

1 試驗設(shè)計和分析方法

1.1 試驗設(shè)計

利用大氣環(huán)流模式NCAR CAM3進行多初值的集合試驗，采用方差分析方法，研究觀測海溫(具有年際變化)強迫下的夏季東亞地區(qū)模式大氣的潛在可預(yù)報性。在1979—2000年的每年夏季(6—8月)的前1個月(5月)用8個初始場(即5月1日00時、2日00時、……、8日00時)進行積分(即積分時分別取5月1—8日00時NCEP/NCAR再分析的大氣資料作為初值，8個初值積分的結(jié)果作為一個集合)，積分時間為4個月，海溫資料為CAM3模式攜帶的1979—2000年逐月的觀測海溫資料，該資料是HadISST OI數(shù)據(jù)和Smith/Reynolds數(shù)據(jù)的結(jié)合。

1.2 方差分析方法

在模式大氣潛在可預(yù)報性分析中，廣泛采用了方差分析方法(Yang et al.，1998)。具體描述如下。

在某一格點上，假定φi，j是某一氣象要素第i個積分、第j年的月(或季)平均值，分別計算樣本總體平均值，以及第j年的集合平均值，其中＜φ＞j隨年份不同而變化。定義

為外部方差，即由海表溫度異常等外強迫變化引起的年際變率;定義

為內(nèi)部方差，即由大氣內(nèi)部動力原因引起的自然變率，它反映了單個積分成員之間的離散度。二者之和被稱作總方差，

其中:M為初始場的個數(shù)(M=8);N為積分的年數(shù)(N=22)。外部方差與總方差之比被記作R，作為模式大氣潛在可預(yù)報性的指標。它反映了外強迫引起的外部方差在變量總方差中所占的比例。R越大，說明外強迫對氣候變量的影響越大，當R＞0.5時，說明外強迫的作用超過了大氣內(nèi)部動力學原因的作用，外強迫對氣候變量具有較大的潛在可預(yù)報能力。

為了分析海溫強迫產(chǎn)生的模式大氣可預(yù)報性的年際變化，對特定年份(第j年)外部方差與內(nèi)部方差，可由下式(Shukla and Gutzler，1983;Quan et al.，2004)計算:

第j年的模式潛在可預(yù)報性可表示為:

2 夏季東亞大氣環(huán)流潛在可預(yù)報性分析

圖1給出了夏季東亞地區(qū)海平面氣壓場的潛在可預(yù)報性(R)的空間分布。如圖1所示，東亞地區(qū)的潛在可預(yù)報性總體偏低，基本小于0.6。在中國區(qū)域，潛在可預(yù)報性呈現(xiàn)出東南高、西北低的分布特征，中國華南地區(qū)R大于0.5，表現(xiàn)出一定的可預(yù)報性，而西北地區(qū)的可預(yù)報性低，R小于0.3。

圖1 東亞地區(qū)夏季海平面氣壓場的潛在可預(yù)報性分布Fig.1 Spatial pattern of potential predictability for sea level pressure over East Asia in summer

圖2為夏季東亞850 hPa緯向風場的潛在可預(yù)報性(R)的空間分布。如圖2所示，東亞地區(qū)850 hPa緯向風場的潛在可預(yù)報性在低緯度地區(qū)明顯高于中高緯度地區(qū)，R大于0.5的區(qū)域主要分布在20°N以南，中國區(qū)域R基本為0.2～0.5，可預(yù)報性較低。

圖2 東亞地區(qū)夏季850 hPa緯向風場的潛在可預(yù)報性分布Fig.2 Spatial pattern of potential predictability for summer zonal wind at 850 hPa over East Asia

用500～200 hPa平均的溫度場來代表對流層平均溫度。Zhou and Zou(2010)研究表明，夏季東亞地區(qū)對流層平均溫度變化具有較高的可預(yù)報性。如圖3所示，對流層平均溫度場的可預(yù)報性明顯高于海平面氣壓場和850 hPa緯向風場。20～30°N對流層平均溫度場的可預(yù)報性基本呈緯向分布，30°N以南地區(qū)R基本為0.5～0.9，中國華南部分地區(qū)R大于0.7，表現(xiàn)出了較高的可預(yù)報性。

圖3 東亞地區(qū)夏季500～200 hPa平均溫度場的潛在可預(yù)報性分布Fig.3 Distribution of potential predictability for the temperature averaged from 500 to 200 hPa over East Asia in summer

圖4給出了夏季東亞500 hPa位勢高度場的潛在可預(yù)報性(R)的空間分布。如圖4所示，R在東亞地區(qū)基本上呈緯向分布特征，從中高緯度地區(qū)向低緯度地區(qū)不斷增大，35°N以南地區(qū)R基本為0.5～0.8，中國華南地區(qū)R大于0.7，表現(xiàn)出了較高的可預(yù)報性?？傮w而言，夏季東亞500 hPa位勢高度場的潛在可預(yù)報性較高。

圖4 東亞地區(qū)夏季500 hPa位勢高度場的潛在可預(yù)報性分布Fig.4 Spatial pattern of potential predictability for summer geopotential height at 500 hPa over East Asia

3 夏季東亞500 hPa位勢高度場潛在可預(yù)報性的年際變化及其與海溫異常的聯(lián)系

圖5是夏季東亞地區(qū)(100～125°E，20～45°N)500 hPa位勢高度場的潛在可預(yù)報性的年際變化?？梢钥闯?，500 hPa位勢高度場的潛在可預(yù)報性具有較強的年際差異，如1998年R為0.91，而1996年R僅為0.10。

圖5 1979—2000年夏季東亞地區(qū)500 hPa位勢高度場潛在可預(yù)報性的年際變化Fig.5 Interannual variation of potential predictability for summer geopotential height at 500 hPa over East Asia from 1979 to 2000

為了分析影響東亞大氣環(huán)流潛在可預(yù)報性年際變化的夏季海溫關(guān)鍵區(qū)，計算了1979—2000年夏季東亞地區(qū)500 hPa位勢高度場的潛在可預(yù)報性與同期全球各點海溫距平的絕對值的相關(guān)系數(shù)(圖6)。圖6顯示，南海、西太平洋和北印度洋的相關(guān)系數(shù)較高，尤其是南海區(qū)域，相關(guān)系數(shù)值大于0.6，通過了0.01信度的顯著性檢驗。而在熱帶中東太平洋地區(qū)的相關(guān)系數(shù)值較小，為0.2左右。1979—2000年夏季東亞地區(qū)500 hPa位勢高度場的潛在可預(yù)報性與Nino3.4區(qū)海溫距平絕對值的相關(guān)系數(shù)僅為0.27，未通過0.1信度的顯著性檢驗。可見，夏季東亞地區(qū)500 hPa位勢高度場的潛在可預(yù)報性的年際變化與夏季南海海溫異常存在密切的聯(lián)系，而與夏季熱帶中東太平洋海溫的關(guān)系并不明顯。夏季東亞地區(qū)500 hPa位勢高度場的潛在可預(yù)報性與前期冬季和春季熱帶中東太平洋海溫存在較好的相關(guān)性(圖略)。

圖7給出了夏季南海區(qū)域(100～120°E，0°～20°N)海溫距平的年際變化。對比圖7和圖5可以看到，在南海區(qū)域海溫距平絕對值較大的年份，東亞地區(qū)500 hPa位勢高度場的潛在可預(yù)報性也較大，而南海區(qū)域海溫距平絕對值較小的年份，東亞地區(qū)500 hPa位勢高度場的潛在可預(yù)報性也較小。

圖6 夏季東亞地區(qū)500 hPa位勢高度場潛在可預(yù)報性與夏季海溫距平絕對值的相關(guān)系數(shù)(陰影區(qū)表示通過0.01信度的顯著性檢驗)Fig.6 Correlation coefficients of potential predictability of the summer 500 hPa geopotential height over East Asia and absolute value of summer SST anomalies(the shaded regions denote the CC exceeding the 0.01 significance level)

圖7 1979—2000年夏季南海區(qū)域海溫距平的年際變化(單位:℃)Fig.7 Interannual variation of summer SCS SST anomaly from 1979 to 2000(units:℃)

圖8為1979—2000年夏季南海區(qū)域海溫距平絕對值與東亞500 hPa位勢高度場的潛在可預(yù)報性的相關(guān)系數(shù)。如圖8所示，東亞地區(qū)的相關(guān)系數(shù)較高，基本大于0.4，30°N以南地區(qū)相關(guān)系數(shù)大于0.6，通過了0.01信度的顯著性檢驗?？梢?，夏季南海區(qū)域海溫異常的強度對東亞500 hPa位勢高度場的潛在可預(yù)報性具有重要影響。

為了進一步分析夏季南海區(qū)域海溫異常對同期東亞500 hPa位勢高度場的潛在可預(yù)報性的影響，分別選取了4個夏季南海海溫偏暖年(1987、1988、1995、1998年)、4個夏季南海海溫偏冷年(1982、1984、1985、1994年)和4個正常年份(1980、1981、1997、2000年)，以比較在不同的南海海溫狀況下，夏季東亞500 hPa位勢高度場的潛在可預(yù)報性的差異。

圖8 夏季南海海溫距平絕對值與東亞500 hPa位勢高度場潛在可預(yù)報性的相關(guān)系數(shù)(帶點號區(qū)域表示通過0.01信度的顯著性檢驗)Fig.8 Correlation coefficients of the absolute value of summer SCS SST anomaly and potential predictability of the summer 500 hPa geopotential height over East Asia(areas with dots denote the CC exceeding the 0.01 significance level)

圖9a和9b表明，夏季南海海溫偏暖年和偏冷年，東亞地區(qū)的潛在可預(yù)報性在30°N以南地區(qū)基本為0.6，表現(xiàn)出了較高的潛在可預(yù)報性。而在夏季南海海溫正常年份，東亞地區(qū)潛在可預(yù)報性基本小于0.3，明顯偏低?？梢?，夏季南海海溫異常對東亞500 hPa位勢高度場的可預(yù)報性具有非常重要的影響。

圖9 夏季南海海溫偏暖(a)、偏冷(b)與正常(c)年份的東亞500 hPa位勢高度場的潛在可預(yù)報性Fig.9 Spatial patterns of potential predictability for summer geopotential height at 500 hPa over East Asia for(a)warmer，(b)colder and(c)normal SCS SST years

4 結(jié)論和討論

利用大氣環(huán)流模式NCAR CAM3，進行了22 a(1979—2000年)、每年8個初值的集合試驗，并采用方差分析方法，研究了觀測海溫強迫下東亞夏季大氣環(huán)流的潛在可預(yù)報性，結(jié)論如下:

1)夏季東亞地區(qū)海平面氣壓場的潛在可預(yù)報性偏低，在中國區(qū)域基本呈現(xiàn)出東南高、西北低的分布特征，華南地區(qū)的可預(yù)報性明顯高于西北地區(qū)。850 hPa緯向風場、對流層500～200 hPa平均溫度場和500 hPa位勢高度場的潛在可預(yù)報性在低緯度地區(qū)高于中高緯度地區(qū)。500 hPa位勢高度場的潛在可預(yù)報性較高，東亞大部分地區(qū)大于0.5，尤其中國華南地區(qū)大于0.7。

2)夏季東亞500 hPa位勢高度場的潛在可預(yù)報性具有明顯的年際變化特征，并與夏季南海區(qū)域的海溫異常具有密切的聯(lián)系。南海區(qū)域海溫為正常年份時，夏季東亞地區(qū)500 hPa位勢高度場的潛在可預(yù)報性較低，基本小于0.3;南海區(qū)域海溫為偏暖或偏冷年份時，夏季東亞地區(qū)500 hPa位勢高度場的潛在可預(yù)報性較高，尤其是30°N以南地區(qū)的潛在可預(yù)報性大于0.6。

基于NCAR CAM3模式多初值集合試驗的模擬結(jié)果，分析了海溫對夏季東亞大氣環(huán)流可預(yù)報性的影響。分析發(fā)現(xiàn)，夏季東亞地區(qū)500 hPa位勢高度場的潛在可預(yù)報性的年際變化與夏季熱帶中東太平洋海溫的相關(guān)性并不高。這可能是由于ENSO的影響多在前冬至春季最顯著，隨著與ENSO相關(guān)的海溫異常的衰減，夏季熱帶中東太平洋海溫與東亞大氣環(huán)流異常的相關(guān)性減弱。此外，在集合試驗中，每年只積分了5—8月的情況，如果從前期冬季開始積分，則有可能會提高ENSO與東亞夏季大氣環(huán)流可預(yù)報性的相關(guān)關(guān)系。另外，受計算條件限制，在集合試驗中每年只選取了8個集合成員(集合成員的預(yù)報起始時間分別為5月1日00時、2日00時、……、8日00時)，由于不同預(yù)報起始時間對潛在可預(yù)報性可能存在影響，因此，如果設(shè)置的預(yù)報起始時間間隔更短，則有可能提高東亞夏季大氣環(huán)流的潛在可預(yù)報性。對此有必要設(shè)計更多的具有不同預(yù)報起始時間的多初值集合試驗來進一步驗證本文的結(jié)果。

丑紀范，徐明.2001.短期氣候數(shù)值預(yù)測的進展和前景［J］.科學通報，46(11):890-895.

范曉青，李維京，張培群.2003.目前短期氣候預(yù)測可預(yù)報性的研究概況［J］.氣象，29(2):3-6.

李崇銀.2000.氣候動力學引論［M］.2版.北京:氣象出版社:449-472.

李建平，丑紀范.2003.非線性大氣動力學的進展［J］.大氣科學，27(4):653-673.

李建平，丁瑞強.2008.短期氣候可預(yù)報期限的時空分布［J］.大氣科學，32(4):975-986.

柳艷菊，馬開玉，林振山.2000.我國月降水量氣候噪聲的估計［J］.應(yīng)用氣象學報，11(2):165-172.

施洪波，周天軍，萬慧，等.2008.SMIP2試驗對亞洲夏季風的模擬能力及其可預(yù)報性的分析［J］.大氣科學，23(1):36-52.

王會軍，薛峰，畢訓強.1997.氣候模式的年際變率和可預(yù)測性［J］.應(yīng)用氣象學報，8(增刊):217-222.

王會軍，陳麗娟，李維京，等.2007.中國區(qū)域月平均溫度和降水的模式可預(yù)報性分析［J］.氣象學報，65(5):725-732.

樂群，曹俊武，林振山，等.1999.中國月平均溫度的氣候噪聲和潛在可預(yù)報性［J］.氣象學報，57(5):604-612.

趙曉川，吳洪寶，何浪.2008.中國季降水量的氣候噪音和潛在可預(yù)報性估計［J］.南京氣象學院學報，31(6):819-827.

趙彥，郭裕福，袁重光，等.2000.短期氣候數(shù)值預(yù)報可預(yù)報性問題［J］.應(yīng)用氣象學報，11(增刊):64-71.

Brankovic C，Palmer T N.1997.Atmospheric seasonal predictability and estimates of ensemble size［J］.Mon Wea Rev，125(5):859-874.

Collins M，Allen M R.2002.Assessing the relatives roles of initial and boundary conditions in interannual to decadal climate predictability［J］.J Climate，15(21):3104-3109.

Hassan A S，Yang X-Q，Zao S-S.2004.Reproducibility of seasonal ensemble integrations with ECMWF GCM and its association with ENSO［J］.Meteor Atmos Phys，86(3/4):159-172.

Kusunoki S，Sugi M，Kitoh A，et al.2001.Atmospheric seasonal predictability experiments by the JMA AGCM［J］.J Meteor Soc Japan，79(6):1183-1206.

Lang Xianmei，Wang Huijun.2005.Seasonal differences of model predictability and the impact of SST in the Pacific［J］.Adv Atmos Sci，22(1):103-113.

Lorenz E N.1982.Atmospheric predictability experiments with a large numerical model［J］.Tellus，34:505-513.

Madden R A.1976.Estimates of the natural variability of time-averaged sea-level pressure［J］.Mon Wea Rev，104(7):942-952.

Quan X W，Webster P J，Moore A M，et al.2004.Seasonality in SST-forced atmospheric short-term climate predictability［J］.J Climate，17(6):3090-3108.

Shukla J，Gutzler D S.1983.Interannual variability and predictability of 500 mb geopotential heights over the Northern Hemisphere［J］.Mon Wea Rev，111(6):1273-1279.

Trenberth K E.1984.Potential predictability of geopotential heights over the Southern Hemisphere［J］.Mon Wea Rev，113(1):54-64.

Yang X-Q，Jeffrey L A，William F S.1998.Reproducible forced modes in AGCM ensemble integrations and potential predictability of atmospheric seasonal variations in the extratropics［J］.J Climate，11(11):2942-2959.

Zhou Tianjun，Wu Bo，Scaife A A，et al.2009.The CLIVAR C20C Project:Which components of the Asian-Australian monsoon circulation variations are forced and reproducible?［J］.Climate Dynamics，33:1051-1068.doi:10.1007/s00382-008-0501-8.

Zhou Tianjun，Zou Liwei.2010.Understanding the predictability of East Asian summer monsoon from the reproduction of land-sea thermal contrast change in AMIP-type simulation［J］.J Climate，23(22):6009-6026.doi:10.1175/2010JCLI3546.1.

Zwiers F W.1996.Interannual variability and predictability in an ensemble of AMIP simulations conducted with the CCC GCM2［J］.Climate Dynamics，12(12):825-847.