孫密娜，管兆勇，張蓬勃，3，曹舒婭，4

(1.南京信息工程大學(xué)氣象災(zāi)害省部共建教育部重點實驗室，江蘇南京210044;2.天津市氣象局，天津300074;3.南京市氣象局，江蘇南京210009;4.江蘇省氣象局，江蘇南京210008)

0 引言

人們很早就注意到海表溫度(sea surface temperature，SST)的異常變化對全球天氣氣候的影響，并對此進行了大量的研究。黃榮輝和孫鳳英(1994)研究了熱帶西太平洋暖池上空的對流活動對我國江淮流域夏季降水的影響。余貞壽等(2005)討論了太平洋SST對中國東部夏季降水的影響。吳勝安等(2005)研究了太平洋SST變化與中國夏季降水年代際變化的相互關(guān)系，發(fā)現(xiàn)20世紀50年代至70年代后期，熱帶西太平洋SST變化為正異常時，長江以南地區(qū)夏季降水偏多，而長江以北地區(qū)降水偏少;之后則反之。陳烈庭和吳仁廣(1998)的研究結(jié)果就表明前期太平洋主要洋流區(qū)的各種海溫距平不同配置的共同影響，是導(dǎo)致中國東部各種雨帶分布類型的一個重要因素。嚴華生等(2002)指出北太平洋12月SST與中國5月降水有較好的相關(guān)關(guān)系。

近年來，許多工作將研究擴展到南太平洋。Hsuing and Newell(1983)指出ENSO(El Ni～no-Southern Oscillation)模態(tài)是太平洋SST場的最主要模態(tài):赤道太平洋和美洲西岸SST偏高，其余地區(qū)偏低，這種分布型顯示了ENSO的年際特征。Garreaud and Battisti(1999)指出南太平洋SST場中存在著比ENSO頻率低的類ENSO結(jié)構(gòu)。Luo and Yamagata(2001)討論了南太平洋類ENSO變化的長期趨勢。Holbrook and Bindoff(1997)發(fā)現(xiàn)西南太平洋海溫存在與南極副極地環(huán)狀海流相聯(lián)系的模態(tài)。而Guan and Yamagata(2001)、盧楚翰等(2008)的研究表明南北半球間存在著南北濤動(inter-hemisphere oscillation，IHO)，這表明北半球包括中國夏季降水與氣候預(yù)測可能與南北半球間的海氣信號存在內(nèi)在聯(lián)系。

然而，整個南太平洋SST變化對中國降水的影響，尤其是南半球夏季的SST變化對中國降水影響的研究還不是很多。本文首先分析了南太平洋12—2月SST變化的主要模態(tài)及其表征的含義，然后分析了它們對中國夏半年降水的影響，并給出了回歸預(yù)報方程。

1 資料與方法

采用NCEP/NCAR 1951—2006年2°×2°分辨率的SST資料，以及NCEP/NCAR 2.5°×2.5°分辨率的高度場、風(fēng)場、溫度場資料。中國160站1951—2006年降水資料。采用經(jīng)驗正交函數(shù)(empirical orthogonal function，EOF)展開、小波分析、相關(guān)分析、多元回歸等方法。

南半球環(huán)狀模(Southern-Hemisphere annual mode，SAM)指數(shù)采用40°S與70°S緯向月平均海平面氣壓標準化差值(南極濤動指數(shù)采用的是40°S與65°S差值)。這是因為緯向平均的海平面氣壓(sea level pressure，SLP)在40°S與70°S的相關(guān)系數(shù)高于在40°S與65°S的相關(guān)系數(shù)(南素蘭和李建平，2005a，2005b)。

2 12—2月南太平洋SST的主要模態(tài)

為了研究南太平洋SST異常的變化規(guī)律，利用EOF分析提取SST異常變化的主要模態(tài)(前3個特征向量通過了North準則檢驗(North et al.，1982))。前兩個特征向量的累積方差貢獻超過50%(表1)，能基本反映相應(yīng)年的12—2月SST異常的空間和時間變化特征。

表1 12—2月南太平洋SST異常EOF分解各個特征向量的方差貢獻率和累積方差貢獻率Table 1Variance contribution ratios of the first three EOFs of the South Pacific SST anomalies from December to February%

2.1 EOF1

EOF1特征向量場(方差貢獻占總方差的41.7%)，反映了12—2月南太平洋SST異常變化的主要空間分布(圖1a)。整個南太平洋海域，SST異常呈現(xiàn)西低東高的分布特征。從南美洲沿岸向西有一暖舌結(jié)構(gòu)伸向赤道西太平洋，從澳大利亞沿岸向東南太平洋有一冷舌結(jié)構(gòu)。30°S的中部海區(qū)也存在一個異常中心。

圖1 12—2月南太平洋SST變化EOF1及其時間序列的小波功率譜分析a.第1特征向量;b.第1特征向量所對應(yīng)的時間系數(shù);c.時間序列的小波功率譜分析(陰影區(qū)域表示通過95%置信水平的檢驗)Fig.1 EOF1 of the DJF South Pacific SST anomalies and its power spectrum analysis of the Mor let wavelet transform from December to Februarya.EOF1;b.its time coefficient;c.power spectrum analysis(the shaded areas denote the significance at 95%confidence level)

當EOF1的時間系數(shù)(圖1b)處于正位相時，赤道海域呈現(xiàn)西冷東暖型分布，赤道外中西部海域海溫降低。反之，當時間系數(shù)處于負位相時，赤道海域呈現(xiàn)東冷西暖型分布，赤道外中西部海域海溫升高。

為了了解其周期變化特征，圖1c給出了EOF1時間系數(shù)的Morlet小波功率譜，結(jié)果表明EOF1所表示的特征場存在著3～5 a的周期，1968—1973年存在準3 a的周期，1980—1988年存在著準4 a的周期，到1995年周期變?yōu)闇? a。研究表明ENSO的主周期為3～7 a(Lau and Shen，1988)，說明EOF1的周期與ENSO的周期相接近。

2.2 EOF2

圖2 12—2月南太平洋SST變化EOF2及其時間序列的小波功率譜分析a.第2特征向量;b.第2特征向量所對應(yīng)的時間系數(shù);c.時間序列的小波功率譜分析(陰影區(qū)域表示通過95%置信水平的檢驗)Fig.2 EOF2 of the DJF South Pacific SST anomalies and its power spectrum analysis of the Morlet wavelet transform from December to Februarya.EOF2;b.its time coefficient;c.power spectrum analysis(the shaded areas denote the significance at 95%confidence level)

EOF2特征向量場(方差貢獻占總方差的12.2%)反映的SST異常的分布型(圖2a)與EOF1的結(jié)構(gòu)不同，赤道太平洋呈現(xiàn)三極型分布，這種分布是否與Ashok et al.(2007)提出的El Ni～no-Modoki有聯(lián)系，目前尚不清楚。同時在南太平洋中部存在兩個與赤道中部相反的中心。

小波功率譜分析表明EOF2所表示的特征場在20世紀60—70年代初以及1995年以后存在著準2 a的周期。

2.3 EOF1、EOF2與ENSO及SAM的聯(lián)系

2.3.1 與ENSO的聯(lián)系

圖3 EOF1對應(yīng)的時間系數(shù)(擴大10倍，空心圈)與標準化的CPC(Climate Prediction Center)冬季Ni～no3指數(shù)(實心圈)Fig.3 The time coefficients of EOF1(×10，hollow circles)and the DJF standardized(Climate Prediction Center)Ni～no3 index in winter(solid circles)

圖4 EOF1對應(yīng)的時間系數(shù)回歸表面(a)、200 hPa(b)u、v風(fēng)場的回歸系數(shù)分布及其與全球SST(a)、全球SLP(b)的相關(guān)分布(陰影區(qū)通過了95%置信水平的檢驗)Fig.4 The series of the global reanalyzed(a)surface and(b)200 hPa(u，v)wind fields regressed upon the time coefficient of EOF1 and correlations of the latter with(a)SST and(b)SLP anomalies(The shaded areas denote the significance at 95%confidence level)

EOF1對應(yīng)的時間系數(shù)和Ni～no3指數(shù)存在顯著的正相關(guān)關(guān)系，相關(guān)系數(shù)達0.89(圖3)。EOF1時間序列與全球海溫的相關(guān)系數(shù)空間分布(圖4a)顯示出顯著的偶極型特征，即赤道中東太平洋是正相關(guān)區(qū)，在赤道西太平洋是負相關(guān)區(qū)。北太平洋中部存在一個負相關(guān)區(qū)，而北美西海岸存在正相關(guān)區(qū)，這種分布與PDO(Pacific Decade Oscillation)的海溫分布類似，反映了南太平洋與北太平洋SST也存在著聯(lián)系(呂俊梅等，2005)。與全球SLP的相關(guān)系數(shù)空間分布(圖4b)，表現(xiàn)為澳大利亞、印度洋和西太平洋為正相關(guān)區(qū)，東太平洋為負相關(guān)區(qū)，與El Ni～no相對應(yīng)。

用EOF1時間系數(shù)回歸表面風(fēng)場，回歸系數(shù)分布(圖4a)顯示，赤道太平洋地區(qū)出現(xiàn)偏西風(fēng)異常，南太平洋上存在一個氣旋，意味著ENSO暖位相時，南太平洋副高減弱，引起太平洋中部風(fēng)暴軸氣候平均位置的改變(Sinclair et al.，1997)。在200 hPa上，赤道太平洋上空存在一對關(guān)于赤道對稱的反氣旋對，這與赤道東太平洋地區(qū)的海表加熱有關(guān)。Gill理論可以解釋這一現(xiàn)象，西風(fēng)帶與Kelvin波相聯(lián)系，東風(fēng)帶與Rossby波相聯(lián)系(Gill，1980)。

2.3.2 與南半球環(huán)狀模的聯(lián)系

EOF2對應(yīng)的時間系數(shù)與3—5月南半球環(huán)狀模指數(shù)的相關(guān)達到－0.59(圖5)，去掉線性趨勢后相關(guān)為－0.29。EOF2與南太平洋1 000 hPa和500 hPa高度場的相關(guān)(圖6)顯示:圍繞極地有一個條帶狀的相關(guān)區(qū)，并且與低緯地區(qū)的符號相反，反映了南半球環(huán)狀模的特征。

EOF2與全球SST的相關(guān)分布(圖7a)表現(xiàn)為:南太平洋中部是負相關(guān)區(qū)，在大陸附近負相關(guān)區(qū)范圍更大。此外，在印度洋也存在著負相關(guān)區(qū)。

EOF2與全球6—8月海平面氣壓的相關(guān)系數(shù)分布(圖7b)顯示出顯著的南北蹺蹺板特征:中緯地區(qū)是負相關(guān)區(qū)，高緯是正相關(guān)區(qū)，在60°S以南形成繞極的帶狀分布。在EOF2負相位時，繞極低壓加深，SAM加強。

用EOF2時間系數(shù)回歸表面風(fēng)場，回歸系數(shù)分布(圖7a)顯示:當EOF2處于正位相時，赤道東太平洋地區(qū)出現(xiàn)偏東風(fēng)異常，在北太平洋存在一個位于北美大陸西岸的氣旋，中心位于50°N。在南太平洋高緯地區(qū)存在著一個以70°S為中心的反氣旋。在北半球由于受陸地的影響風(fēng)向比較散亂。

圖5 EOF2對應(yīng)的時間系數(shù)(擴大10倍，空心圈)與南半球環(huán)狀模指數(shù)(實心圈)Fig.5 The time coefficients of EOF2(×10，hollow circles)and the DJF Southern hemisphere Annual Mode index(solid circles)

圖6 南太平洋12—2月EOF2與南半球1 000 hPa(a)和500 hPa(b)高度場的相關(guān)(淺、深陰影分別表示通過95%、99%置信水平的檢驗)Fig.6 Correlations of time coefficient of EOF2 with(a)1 000 hPa and(b)500 hPa geopotential height(The light and dark shaded areas denote the significance at 95%and 99%confidence levels respectively)

圖7 EOF2對應(yīng)的時間系數(shù)回歸表面(a)、200 hPa(b)u、v風(fēng)場的回歸系數(shù)分布及其與全球SST(a)、全球SLP(b)的相關(guān)分布(陰影區(qū)通過了95%置信水平的檢驗)Fig.7 The series of the global reanalyzed(a)surface and(b)200 hPa(u，v)wind fields regressed upon the time coefficient of EOF2 and correlations of the latter with(a)SST and(b)SLP anomalies(The shaded areas denote the significance at 95%confidence level)

200 hPa上(圖7b)，在南太平洋上存在著與PSA(Pacific South American)型相對應(yīng)的波列，北太平洋上的波列與PNA(Pacific Nouth American)型相對應(yīng)。與EOF1的回歸結(jié)果相似，但是波列的分布更清晰。在印度洋上分布著一對關(guān)于30°S對稱的氣旋對。

2.4 與我國夏季降水的聯(lián)系

圖8 南太平洋EOF1對應(yīng)的時間系數(shù)去除線性趨勢前(a)、后(b)與中國5月降水的相關(guān)和EOF2對應(yīng)的時間系數(shù)去除線性趨勢前(c)、后(d)與中國6—8月降水的相關(guān)(陰影區(qū)表示通過90%置信水平檢驗)Fig.8 Correlations of time coefficient of EOF1(a)before and(b)after removing the linear trend with the rainfall of China in May;correlations of time coefficient of EOF2(c)before and(d)after removing the linear trend with the JJA rainfall of China(The shaded areas denote the significance at 90%confidence level)

EOF1、EOF2時間系數(shù)與降水的聯(lián)系顯示出顯著的滯后性。圖8a顯示了EOF1時間系數(shù)與中國5月降水的相關(guān)空間分布，顯示出在華北地區(qū)有正相關(guān)區(qū)，最大相關(guān)達0.41，說明在EOF1時間系數(shù)為正值時，5月華北夏季的降水偏多。研究(楊芳林和袁重光，1995;張秉倫等，1998;勵申申和壽紹文，2000;孫淑清和馬淑杰，2003;高輝，2006)顯示:當前期ENSO為暖位相狀態(tài)時，夏季中國主要雨帶位置偏南，長江流域降水偏多;反之，當前期ENSO為冷位相狀態(tài)時，夏季中國主要雨帶位置偏北，長江流域降水偏少。圖8c顯示EOF2時間系數(shù)與同年我國夏季降水的相關(guān)空間分布，在江南地區(qū)存在負相關(guān)，在華北存在正相關(guān)區(qū)。去掉EOF1和EOF2的線性趨勢后(圖8b、8d)，這種相關(guān)仍然存在，只是范圍有所改變。

為了驗證南太平洋SST對中國降水的影響，圖9分別給出了12—2月南太平洋SST與中國降水的相關(guān)。南太平洋12—2月平均SST與華北5月降水的相關(guān)區(qū)域集中在熱帶東太平洋地區(qū)，與華北6—8月平均降水的相關(guān)區(qū)域集中在中緯度地區(qū)，與長江中下游夏季平均降水的相關(guān)區(qū)域集中在東南太平洋地區(qū)。

2.5 對中國5月降水的預(yù)報

圖9 華北5月(a)、華北夏季(b)以及長江中下游夏季(c)降水序列與南太平洋12—2月平均SST的相關(guān)Fig.9 Correlations of time coefficient of DJF South Pacific mean SST with(a)North China rainfall in May，(b)North China summer rainfall and(c)summer rainfall over the middle and lower reaches of the Yangtze River(The shaded areas denote the significance at 90%confidence level)

為了更好地對中國降水做出預(yù)報，比較南北太平洋對中國降水的影響，對北太平洋12—2月SST變化做EOF分析，前兩個模態(tài)如圖10所示，研究顯示:北太平洋冬季EOF1主要反映了ENSO信號，EOF2反映的主要是12—2月北太平洋海溫的異常模態(tài)，代表了北太平洋海溫中緯度與低緯度及高緯度的局地差異。即東亞沿海、西太平洋暖池區(qū)SST與北太平洋西風(fēng)漂流區(qū)SST呈反位相變化(解思梅等，1986;高慶九和屠其璞，2003)。對應(yīng)EOF1時間序列與中國5月降水的相關(guān)(圖11a)與南太平洋EOF1時間序列的相關(guān)范圍相近，都反映了ENSO對中國5月降水的影響。12—2月EOF2時間系數(shù)與同年中國夏季月降水在江南地區(qū)存在負相關(guān)(圖11b)。利用南、北太平洋的EOF1時間序列(南太平洋:x1;北太平洋:x2)與中國華北地區(qū)5月降水建立回歸預(yù)報方程(y=35.22+66.35x1－2.15x2)，其方差貢獻達到20.3%。圖11c反映的是預(yù)報序列與實際降水序列的相關(guān)，在華北地區(qū)的相關(guān)很高，這一結(jié)果顯示太平洋的SST能夠影響中國夏季降水并能做出預(yù)測。

圖10 北太平洋SST變化EOF1(a)、EOF2(b)及其所對應(yīng)的時間序列(c、d)Fig.1 0EOFs of the North Pacific SST anomalies in DJF(a，c.EOF1 and its time coefficient;b，d.EOF2 and its time coefficient)

圖11 北太平洋EOF1對應(yīng)的時間系數(shù)與中國5月(a)、EOF2對應(yīng)的時間序列與6—8月(b)降水的相關(guān)以及華北5月降水預(yù)報序列與實際降水量的相關(guān)(c)(陰影區(qū)表示通過90%置信水平檢驗)Fig.1 1Correlations of time coefficient of EOF1，EOF2 with the rainfall of China in(a)May，and(b)JJA，and(c)correlation between the prediction series and observed rainfall in May(The shaded areas denote the significance at 90%confidence level)

3 結(jié)論

1)EOF1主要反映了ENSO信息，EOF2與南半球環(huán)狀模聯(lián)系比較緊密。通過小波功率譜分析，發(fā)現(xiàn)EOF1的周期和ENSO的周期相近，主要是年際變化，EOF2的周期為準兩年的周期。EOF1對應(yīng)的時間序列與Ni～no3指數(shù)的相關(guān)達到0.89，反映了其與ENSO的聯(lián)系。EOF2時間序列與3—5月SAM指數(shù)的相關(guān)明顯，與500 hPa和200 hPa高度場的相關(guān)在高低緯呈反相的環(huán)狀分布，反映了其與南半球環(huán)狀模之間的聯(lián)系。EOF2與太平洋的行星波的遙相關(guān)有關(guān)，反映出EOF2與高緯(60°S以南)地區(qū)的大氣活動聯(lián)系密切。

2)通過分析EOF1與我國5月降水的相關(guān)場，得出當南太平洋SST出現(xiàn)西低東高偶極型分布的年份華北地區(qū)5月降水偏多，反之則降水偏少。EOF2與中國6—8月平均降水的相關(guān)場反映出:赤道中太平洋SST偏高而其他海域SST偏低的年份我國江南地區(qū)夏季降水偏少而華北降水偏多，反之則江南地區(qū)降水偏多，華北降水偏少。去掉EOF1和EOF2的線性趨勢后，這種相關(guān)仍然存在，只是范圍有所縮小。

3)南太平洋的EOF1與中國5月降水的相關(guān)與北太平洋EOF1的相關(guān)分布一致，這是由于它們都反映了ENSO信號。利用南、北太平洋的EOF1時間序列建立回歸方程，可為華北5月降水的預(yù)測提供線索。

致謝:NOAA-CIRES Climate Diagnostics Center和國家自然科學(xué)基金委員會地球科學(xué)部大氣資料服務(wù)中心提供資料服務(wù)，謹致謝忱!

陳烈庭，吳仁廣.1998.太平洋各區(qū)海溫異常對中國東部夏季雨帶類型的共同影響［J］.大氣科學(xué)，22(5):718-722.

高輝.2006.淮河夏季降水與赤道東太平洋海溫對應(yīng)關(guān)系的年代際變化［J］.應(yīng)用氣象學(xué)報，17(1):1-9.

高慶九，屠其璞.2003.20世紀70年代前后北太平洋海溫場氣候特征的比較［J］.南京氣象學(xué)院學(xué)報，26(2):223-251.

黃榮輝，孫鳳英.1994.熱帶西太平洋暖池的熱狀況及其上空的對流活動對東亞夏季氣候異常的影響［J］.大氣科學(xué)，18(2):141-151.

勵申申，壽紹文.2000.赤道東太平洋海溫與我國江淮流域夏季旱澇的成因分析［J］.應(yīng)用氣象學(xué)報，11(3):331-338.

盧楚翰，管兆勇，梅士龍，等.2008.大氣質(zhì)量南北濤動的季節(jié)循環(huán)［J］.科學(xué)通報，53(21):2638-2645.

呂俊梅，琚建華，張慶云，等.2005.太平洋年代際振蕩冷、暖背景下ENSO循環(huán)的特征［J］.氣候與環(huán)境研究，10(2):238-249.

南素蘭，李建平.2005a.春季南半球環(huán)狀模與長江流域夏季降水的關(guān)系:Ⅰ基本事實［J］.氣象學(xué)報，63(6):837-846.

南素蘭，李建平.2005b.春季南半球環(huán)狀模與長江流域夏季降水的關(guān)系:Ⅱ印度洋和南海海溫的“海洋橋”作用［J］.氣象學(xué)報，63(6):847-856.

孫淑清，馬淑杰.2003.海溫異常對東亞夏季風(fēng)及長江流域降水影響的分析及數(shù)值試驗［J］.大氣科學(xué)，27(1):36-51.

吳勝安，江志紅，劉志雄，等.2005.中國夏季降水與太平洋SSTA年代際變化關(guān)系的初步研究［J］.熱帶氣象學(xué)報，21(2):153-162.

解思梅，青木孝，吉野正敏.1986.北太平洋海表面溫度的年際變化特征及其與東亞氣候的相關(guān)關(guān)系［J］.海洋學(xué)報，8(6):547-557.

嚴華生，呂俊梅，琚建華，等.2002.冬季太平洋海溫變化對中國5月降水的影響［J］.氣象科學(xué)，22(4):410-415.

楊芳林，袁重光.1995.夏季赤道東太平洋海溫異常對全球及東亞短期氣候變化影響的數(shù)值試驗［J］.大氣科學(xué)，19(5):535-544.

余貞壽，孫照渤，曾剛.2005.太平洋SSTA對中國東部夏季降水的影響I:觀測分析［J］.熱帶氣象學(xué)報，21(2):467-477.

張秉倫，王成興，曹永忠.1998.厄爾尼諾與江淮流域旱澇災(zāi)害的關(guān)系［J］.自然雜志，20(5):289-293.

Ashok K，Behera S K，Rao S A，et al.2007.El Ni～no Modoki and its possible teleconnection［J］.J Geophys Res，112，C11007，doi:10.1029/2006JC003798.

Garreaud R D，Battisti D S.1999.Interannual(ENSO)and interdecadal(ENSO-like)variability in the Southern Hemisphere tropospheric circulation［J］.J Climate，12(7):2113-2122.

Gill A E.1980.Some simple solutions for heat induced tropical circulations［J］.Quart J Roy Meteor Soc，106(5):447-462.

Guan Z Y，Yamagata T.2001.Interhemispheric oscillations in the surface air pressure field［J］.Geophys Res Lett，2001，28(2):263-266.

Holbrook N J，Bindoff N L.1997.Interannual and decadal temperature variability in the Southwest Pacific Ocean between 1955 and 1988［J］.J Climate，10(5):1035-1049.

Hsuing J，Newell R E.1983.The principal nonseasonal modes of global sea surface temperature［J］.Pap Phys Oceanogr Meteor，13(10):1957-1967.

Luo J J，Yamagata T.2001.Long-term El Ni～no-Southern Oscillation(ENSO)-like variation with special emphasis on the South Pacific［J］.Geophys Res Lett，106(22):22211-22227.

Lau K M，Shen P.1988.Annual cycle，QBO，and southern oscillation in global precipitation［J］.J Geophys Res，93:10975-10988.

North G R，Bell T L，Cahalan R F.1982.Sampling errors in estimation of empirical orthogonal function［J］.Mon Wea Rev，110(7):699-706.

Sinclair M R，Renwick J A，Kidson J W.1997.Low-frequency variability of Southern Hemisphere sea level pressure and weather system activity［J］.Mon Wea Rev，125(10):2531-2542.