歐亞季節(jié)環(huán)流氣候平均態(tài)和氣候變率：20CR再分析資料在歐亞及中國東部地區(qū)的評估

劉輝1, 2 范可1

?

歐亞季節(jié)環(huán)流氣候平均態(tài)和氣候變率：20CR再分析資料在歐亞及中國東部地區(qū)的評估

劉輝范可

1中國科學院大氣物理研究所竺可楨—南森國際研究中心，北京100029 2中國科學院大學，北京100049

本文評估了美國國家大氣海洋局（NOAA）新發(fā)布的20世紀再分析資料（20CR）對歐亞季節(jié)環(huán)流氣候平均態(tài)和氣候變率及中國東部氣溫降水的刻畫能力。結(jié)果表明：20CR再分析資料對歐亞地區(qū)四季環(huán)流氣候平均態(tài)刻畫能力與NCEP2資料的相比，均呈現(xiàn)北部中高緯度系統(tǒng)性偏高，南部中低緯度系統(tǒng)性偏低的特點，導致描繪的東亞冬季風偏弱，夏季風偏強。這可能與20CR資料在極地海岸地區(qū)海冰資料處理時產(chǎn)生的差錯有關(guān)。與中國東部的站點資料對比則顯示20CR對我國東部氣溫的刻畫偏低，而對降水的刻畫偏高，站點相關(guān)性氣溫好于降水，東南沿海地區(qū)優(yōu)于內(nèi)陸地區(qū)。平均場和空間相關(guān)場結(jié)合來看，秋季氣溫和降水20CR與站點觀測資料吻合最好。20CR資料較好地刻畫近百年北半球冬夏季的氣候指數(shù)（北極濤動、北大西洋濤動、北太平洋濤動、東亞冬季風、阿留申低壓等）的年際變率及年代際變化特征，很好地刻畫了阿留申低壓1970年代末的年代際增強，西伯利亞高壓1970年代末的下降和1990年后的上升趨勢及北太平洋濤動、北大西洋濤動和北極濤動指數(shù)1970年代末期由負位相到正位相的年代際轉(zhuǎn)變。

20CR再分析資料 歐亞氣候平均態(tài) 氣候變率 年代際變化 中國東部

1 引言

隨著資料同化技術(shù)和數(shù)值預報模式的不斷發(fā)展，再分析資料已經(jīng)廣泛應用于現(xiàn)代氣候研究的各個方面，成為模式驗證、氣候變率及長期變化趨勢診斷的有力工具。因此這種經(jīng)過同化加工的再分析資料能在多大程度上真實描述大氣的三維狀態(tài)就成了氣象工作者們十分關(guān)心的問題。再分析資料的質(zhì)量問題主要受兩個方面影響，首先受預報模式及資料同化分析系統(tǒng)的影響，如NCEP/NCAR 聯(lián)合發(fā)布的第一套再分析資料（NCEP1）中，由于其同化分析系統(tǒng)出現(xiàn)解碼錯誤，致使東亞地區(qū)的海平面氣壓在1967年之前出現(xiàn)了系統(tǒng)性偏低的現(xiàn)象（徐影和丁一匯，2001；Yang et al., 2002；Wu et al., 2005；趙天保和符淙斌，2009）。其次，所同化的觀測資料的不均一性也可能導致再分析資料中出現(xiàn)虛假變化趨勢（Marshall，2002；Marshall，2003；Hodges et al., 2003；Bromwich and Fogt，2004）。如歐洲中期數(shù)值預報中心（ECMWF）的近45年再分析資料（ERA40）計算的全球年平均總動能，在1970年代前后，隨著衛(wèi)星資料的引進有一個明顯的上升趨勢，而在考慮了衛(wèi)星資料的出現(xiàn)的敏感性試驗中，這種上升趨勢就會大大減弱（Bengtsson et al., 2004）。然而無論NCEP1資料還是ERA40資料，其時間自20世紀中后期開始，其覆蓋范圍都較短，無法滿足研究百年尺度年代際氣候變化的要求。

近年，美國國家大氣海洋局（NOAA）利用最先進的同化資料并結(jié)合觀測資料，提供了一套從地面到高空20世紀大氣再分析資料（20th Century Reanalysis Version 2，20CR）。這套資料不但時間跨度長（1871～2010年），且為避免觀測系統(tǒng)變更和資料內(nèi)部不協(xié)調(diào)所帶來的虛假趨勢，它僅同化地面和海平面的氣壓數(shù)據(jù)，用觀測的月平均海表面溫度和海冰作為下邊界驅(qū)動條件（Compo et al., 2011）。20CR資料在長期氣候變化趨勢的研究中已經(jīng)有了廣泛應用，但有些結(jié)果也受到質(zhì)疑。20CR資料計算的緯向平均降水量與蒸發(fā)量之差與20世紀氣候模式描繪的結(jié)果有顯著差異（Compo et al., 2011）；其氣溫和降水資料在美國中部也出現(xiàn)不連續(xù)性（Ferguson and Villarini，2012）；關(guān)于東北部大西洋風暴20CR資料得出的結(jié)果與觀測結(jié)果在20世紀前半葉呈相反趨勢（Krueger et al.，2012）。這些結(jié)果表明20CR資料還存在一定的缺陷，因此有必要在使用前對20CR資料進行評估。宋豐飛和周天軍（2012）對20CR資料在東亞季風區(qū)做了質(zhì)量評估，指出盡管20CR資料對東亞夏季氣候態(tài)和年際變率的刻畫與此前公布的其他再分析資料有較高的一致性，但其未能再現(xiàn)1970s末期東亞夏季風年代際減弱的信號，說明在年代際尺度上可靠性較低。Fan and Liu（2013）對20CR資料在南半球大氣環(huán)流的描述能力上也做了評估，指出在1979～2010年 間，20CR資料刻畫的南半球的氣候態(tài)、遙相關(guān)模態(tài)及越赤道氣流與對比資料有較高的一致性。但1920年代之前20CR資料與觀測資料差異明顯。

由于觀測站點的相對較少，再分析資料在歐亞地區(qū)的可信度要弱于北半球其他地區(qū)，而為了能在歐亞地區(qū)更合理的使用這套新的再分析資料，有必要對其描繪的歐亞季節(jié)環(huán)流氣候平均態(tài)和氣候變率做一個完整的評估。本文選取NCEP/NCAR聯(lián)合完成的NCEP2再分析資料作為對照資料，主要因為它是采用改進的預報模式和同化系統(tǒng)，修正了NCEP1中的人為誤差，被認為是一種校正的較好的全球再分析資料，也是應用較為廣泛的再分析資料之一。本文首先從四個季節(jié)氣候平均態(tài)的角度對比20CR與NCEP2再分析資料季節(jié)演變的異同。并在此基礎上分別對20CR的2 m氣溫和降水資料與中國站點的氣溫和降水資料進行了四季氣候平均態(tài)的對比，進一步驗證20CR再分析資料在中國東部地區(qū)的可信性。其次，通過與Hadley氣候中心的百年海平面氣壓場（SLP）資料做對比，考察了北半球主要氣候指數(shù)在近百年的年際及年代際變率的刻畫能力。

2 資料介紹及方法

20CR再分析資料是美國大氣海洋局（NOAA）最新發(fā)布的一套綜合性全球大氣環(huán)流數(shù)據(jù)，時間覆蓋范圍從1871年到2010年，水平分辨率為2°×2°。作為參照的NCEP2（Kanamitsu et al.，2002）再分析資料時間覆蓋范圍為1979年至今，分辨率為2.5°× 2.5°。取兩套資料1979～2010年相應變量的月平均數(shù)據(jù)，在做比較前，首先將NCEP2再分析資料對應變量水平方向上用雙線性插值到2°×2°。

在做20CR再分析資料月中國站點資料的對比時，選取20CR的2 m氣溫和降水資料，均為高斯格點分布，全球共192×94格點。將其降水率轉(zhuǎn)換為季節(jié)降水總量，同樣取1979年1月～2010年12月32年的資料。站臺資料是同期的中國160臺站月平均溫度和降水總量的觀測資料。對比之前將格點資料用雙線性插值方法插值到站點，得到中國160個站臺的20CR資料，由于對比區(qū)域選為105°E以東的中國東部地區(qū)，地形起伏不明顯，因此插值過程中沒有考慮地形因素。

此外，為了檢驗百年尺度海平面氣壓，本文還用到英國Hadley中心提供的HadSLPv2全球海平面氣壓月平均數(shù)據(jù)，該數(shù)據(jù)完全以海平面氣壓的觀測資料為基礎，水平分辨率為5°×5°，選用時間段為1871～2010年（Allan and Ansell，2006）。

分析20CR與NCEP2兩套再分析資料的差異場時，均為20CR減去NCEP2；分析20CR資料與中國東部站點資料的差異場時，均為20CR資料減去站點資料。

3 氣候平均態(tài)的季節(jié)對比

3.1 海平面氣壓場（SLP）的四季對比

圖1為20CR與NCEP2再分析資料海平面氣壓（SLP）的四季氣候平均態(tài)及其差異的空間對 比。如圖1a、b所示，冬季歐亞大陸上存在一強大的反氣旋環(huán)流，即西伯利亞高壓（SH）。海洋上 則對應為一較強的氣旋性低壓，即阿留申低壓（AL）。結(jié)合差異場（圖1c）可以看出，20CR的SH位置偏北，導致歐亞大陸45°N以北及中國東北地區(qū)20CR的SLP都顯著偏高。中國西南部、青藏高原、印度半島則顯著偏低。AL沒有顯著的位置和強度差異。結(jié)合850 hPa風場的差異（圖1c），歐亞大陸45°N以南是顯著的氣旋性差值環(huán)流，東亞沿岸盛行偏南風異常，說明20CR刻畫的東亞冬季風比NCEP2弱。選取（20°N～40°N，110°E～125°E）范圍加權(quán)平均的850 hPa經(jīng)向風為冬季風指數(shù)（Yang et al., 2002），定量計算可知1979～2010年時段20CR資料所計算的冬季風比NCEP2的計算結(jié)果高0.507 m/s。因冬季風符號為負值，說明20CR的冬季風更弱。春季時（圖1d、e）隨著歐亞大陸和西太平洋海陸熱力差異的減小，大陸上的反氣旋和海洋上的氣旋都相應變?nèi)酢H西撤到內(nèi)陸地區(qū)，亞洲低壓在印度半島附近開始形成。AL向西北轉(zhuǎn)移其范圍則延伸到大陸，夏威夷高壓漸漸增強。結(jié)合差異場（圖1f）可以看出，20CR所刻畫的西伯利亞高壓和亞洲低壓在春季更強一些，表現(xiàn)為大陸上45°N以北，我國東北華北及日本海地區(qū)大范圍正距平，青藏高原及南亞次大陸印度半島地區(qū)呈現(xiàn)顯著負距平。海洋上北太平洋地區(qū)兩套資料的差異不明顯，15°N以南的熱帶太平洋地區(qū)20CR的氣壓場偏低。從夏季SLP平均場（圖1g、h）可以看出，海洋上為一個強大的夏威夷高壓，歐亞大陸上呈現(xiàn)熱低壓，另一低壓中心在印度半島附近。我國東北有一氣旋性低壓。結(jié)合差值場（圖1i），20CR資料在45°N以南，120°E以西的歐亞大陸地區(qū)為顯著負距平，在45°N以北，120°E以東的歐亞大陸為顯著正距平，海洋上20°N以南的熱帶太平洋地區(qū)為顯著負距平，20°N以北的北太平洋地區(qū)差異不明顯。這說明20CR所刻畫的亞洲低壓偏強，東北低渦偏弱，夏威夷高壓兩套資料差異不大。結(jié)合850 hPa風場的差異（圖1i，箭頭），可看出我國東南華北地區(qū)都有顯著的南風異常，說明20CR描繪的夏季風要強于NCEP2。選取（20°N～40°N，110°E～125°E）范圍加權(quán)平均的850 hPa經(jīng)向風為夏季風指數(shù)（Wang, 2002），定量計算可知1979～2010時段20CR資料所計算的夏季風比NCEP2的計算結(jié)果高0.933 m/s。秋季時（圖1j、k、l）亞歐大 陸上西伯利亞高壓開始形成，海洋上阿留申低壓變強范圍擴大，夏威夷高壓變?nèi)醪⑾驏|南撤退，差異場上45°N以南的亞歐大陸地區(qū)負距平，以北正距平，850 hPa風場上顯示我國西南華北地區(qū)偏南風異常。

3.2 500 hPa高度場（H500）的四季對比

圖2為20CR與NCEP2再分析資料500 hPa高度場（H500）的四季氣候平均態(tài)及其差異的空間對比。冬季時（圖2a、b）鄂霍茨克海附近存在一明顯的東亞大槽，貝加爾湖地區(qū)和烏拉爾山附近分別有一弱脊，在20CR資料場中（圖2a）西太平洋副高不明顯，NCEP2資料場（圖2b）副高位于20°N以南，脊線大約在10°N左右，處于一年中最南的位置。差異場中（圖2c），60°N以北顯著正距平，40°N以南的大陸及海洋上為顯著的負距平。鄂霍茨克海處顯著正距平。可看出20CR的資料所描繪的東亞大槽偏弱，極渦也明顯偏弱，貝加爾湖高壓脊和烏拉爾山高壓脊都略強，副高明顯偏弱。春季時（圖2d、e、f）亞歐大陸上的槽脊系統(tǒng)都相應減弱且東移，副高的位置北移，脊線大概在15°N。差異場中我國東北、華北、及蒙古地區(qū)上空存在一反氣旋差值環(huán)流，20°N以南負異常顯著。夏季時（圖2g、h、i）東亞大槽已東移入海，貝加爾湖出現(xiàn)一淺槽，伊朗高壓和北太平洋副高都十分明顯，兩個高壓的中心值20CR都要略高于NCEP2，副高的脊線位置大約在30°N。差異場中（圖2i）20°N以北的歐亞大陸大范圍的顯著正距平，我國東北地區(qū)為正距平中心，20°N的印度洋，南亞次大陸及熱帶太平洋上空為顯著負距平。對比平均態(tài)，可知20CR所描繪的伊朗高壓要略強于NCEP2的結(jié)果，且位置偏北。副高5880和5860線的西伸端點較NCEP2的結(jié)果更偏西。 90°E附近的貝湖槽偏弱。秋季（圖2j、k、l）伊朗高壓和太平洋副高都有所減弱，脊線位置南撤到20°N左右，中心最高值NCEP2略高于20CR，亞歐大陸上鄂霍茨克海沿岸東亞大槽開始形成。結(jié)合差異場，可以看到20CR的資料30°N以南大范圍負距平，說明20CR所刻畫的秋季副高比NCEP2，明顯偏低。

3.3 200 hPa緯向風場（U200）的四季對比

圖3為200 hPa緯向風的四季平均態(tài)及差異對比圖。可看出高空西風急流的位置和強度隨季節(jié)變化有很大差異。冬季最強位置最南，急流中心在日本上空，中心風速最大達70 m/s以上，急流中心位置約在30°N，夏季時日本急流最弱位置最北，東移到北太平洋。另一西風急流中心出現(xiàn)我國新疆上空，中心風速僅30 m/s以上，急流中心位置約在42°N左右。從四季的差異場（圖3c、f、i、l）可以看出，以45°N為界，以北為正距平，以南為負距平，分別以貝加爾湖和我國南部地區(qū)為正負距平中心。其中春季差異值最大，夏秋季其次，冬季差異值及顯著性范圍都較小。同時注意到急流中心位置處均無顯著性差異，說明兩套資料對西風急流強度和位置的刻畫比較一致。同時，由于45°N以北的正異常和45°N以南的負異常在高度場上會體現(xiàn)為一個反氣旋高度異常，南亞高壓在夏季時最為明顯，其脊線位置大概位于20°N，45°N附近的反氣旋異常意味著20CR資料所刻畫的南亞高壓位置要更偏北，這點從200 hPa高度場的差值場也可以看出（圖略）。

總體來看，1979～2010年期間，歐亞地區(qū)20CR資料的SLP、H500、U200的四季氣候平均態(tài)與NCEP2資料的相比，呈現(xiàn)北部中高緯度系統(tǒng)性偏高，南部中低緯度偏系統(tǒng)性低的特點，這可能與20CR資料在極地海岸地區(qū)海冰資料處理時產(chǎn)生的差錯有關(guān)（Compo et al., 2011）。這種北正南負的異常會導致其在描述東亞氣候時相比NCEP2資料冬季風偏強，夏季風偏弱。

4 與中國站點資料的比對

一般來說，相對于高空再分析產(chǎn)品來說，再分析地表氣溫與降水資料的可信度相對較低，因此有必要對20CR的地表氣溫與降水資料與我國站點器測資料相對比進行評估。圖4、圖5分別為四季20CR插值結(jié)果與中國站點氣溫和降水資料的差值對比，菱形圈出的為差異通過0.05顯著性檢驗的站點，即存在顯著差異的站點。圖6、圖7分別為四季20CR插值結(jié)果與中國站點氣溫和降水的空間相關(guān)，紅色點為通過0.05顯著性相關(guān)檢驗的站點，藍色點為非顯著相關(guān)站點及負相關(guān)站點。時間范圍均為1979～2010年。由于105°E以西的地區(qū)地形變化較大，站點也較少，我們只選取105°E以東的區(qū)域進行比較。

從溫度的四季差異圖（圖4）可以看出，冬春季20CR的氣溫資料都偏低，夏季35°N以北的北方地區(qū)大部分偏高，南方地區(qū)則偏低。秋季中部偏高東北西南及東南溫度偏低。四個季節(jié)相比較，冬季無顯著差異的站點相對較多，多集中在華北、華南地區(qū)。冬季華南地區(qū)氣溫平均差值較小，大多在2°C以內(nèi)，東北地區(qū)差異較大。夏季時則東北地區(qū)差異較小，南方地區(qū)差異相對較大，但大多差異溫度在4°C以內(nèi)。

圖5為降水的四季差異結(jié)果，比較季節(jié)降水總量，可以看出20CR刻畫的降水在四季都普遍偏高。冬季的差異最為明顯，只有南部沿海少數(shù)幾個站點無顯著差異。夏秋季無顯著差異的站點較多，集中在山東半島及我國東南部地區(qū)。降水無顯著差異的站點數(shù)明顯多于氣溫的，這可能是由于降水的變率比較大，方差比較大，導致多年平均的差異易通過檢驗。

圖6是20CR再分析2 m氣溫資料與觀測氣溫資料的相關(guān)系數(shù)（）的空間分布特征，紅色的為相關(guān)性通過0.05信度檢驗的站點。由圖可見四季中再分析資料與大部分站點資料顯著相關(guān)，但相關(guān)系數(shù)達到0.9以上的站點很少，總體來看，相關(guān)性東部沿海及東北地區(qū)要好于內(nèi)陸地區(qū)，秋季的相關(guān)性為四季中最好的。

圖7為20CR再分析降水資料與器測降水資料的相關(guān)系數(shù)的空間分布特征。可以看出降水的相關(guān)性要明顯低于氣溫的相關(guān)性，不顯著相關(guān)的站點增多，甚至個別站點出現(xiàn)負相關(guān)。同樣東南沿海地區(qū)的相關(guān)性要好于內(nèi)陸地區(qū)，春秋季的相關(guān)性好于冬夏季。

圖4 中國105°E以東地區(qū)1979～2010年（a）冬、（b）春、（c）夏、（d）秋季節(jié)20CR氣溫插值結(jié)果與站點氣溫的差值。紅色為正，藍色為負；菱形圈出的為差異未通過了顯著性0.05檢驗的站點

5 氣候變率的比較

選取英國Hadley中心的海平面氣壓場和20CR的海平面氣壓場作對比，考察北半球主要氣候指數(shù)的年際和年代際變率（圖8、圖9）。圖8為冬季的氣候指數(shù)，分別為西伯利亞高壓指數(shù)（SHI=SLP）、阿留申低壓指數(shù)（ALI=SLP）、東亞冬季風指數(shù)EAWMI（郭其蘊，1994）、北大西洋濤動指數(shù)（NAOI=SLP－SLP）、北太平洋濤動指數(shù)（NPOI= SLP－SLP）、南方濤動指數(shù)（SOI=SLP－SLP）和北極濤動指數(shù)（AOI，20°N以北SLP場EOF第一模態(tài)的時間系數(shù)）。其中指數(shù)的計算公式中下標為區(qū)域的都表示所選取的SLP區(qū)域的加權(quán)平均。紅線為20CR資料的結(jié)果，藍線為HadSLP的結(jié)果，黑線為HadSLP結(jié)果的9年濾波曲線。右上角分別為1871～1920年和1921～2010年兩套資料原始序列指數(shù)的相關(guān)系數(shù)。圖9為夏季的氣候指數(shù)，分別為東亞夏季風指數(shù)（郭其蘊，1983）、夏季NAO（Greatbatch and Rong，2006）、NPO、SO指數(shù)，其中NPOI、SOI兩個指數(shù)與冬季時計算方法相同。所有指數(shù)均進行了標準化處理。東亞冬夏季風指數(shù)采用郭其蘊（1983，1994）的計算方法，具體為取10°～50°N范圍內(nèi)（每10度），110°～160°E之間的氣壓差小于（大于）或等于－5 hPa（5 hPa）的所有數(shù)值之和，再把各年的值與氣候平均態(tài)求比值，作為夏季風（冬季風）指數(shù)。

圖5 中國105°E以東地區(qū)1979～2010年冬（a）、春（b）、夏（c）、秋（d）季節(jié)20CR降水插值結(jié)果與站點降水的差值。紅色為正，藍色為負；菱形圈出的為差異未通過了顯著性0.05檢驗的站點

圖6 中國105°E以東地區(qū)1979～2010年（a）冬、（b）春、（c）夏、（d）秋季節(jié)20CR氣溫插值結(jié)果與站點氣溫的空間相關(guān)。紅色為相關(guān)性通過了顯著性0.05檢驗的站點

從各氣候指數(shù)的相關(guān)系統(tǒng)來看，20CR資料的結(jié)果與HadSLPv2的結(jié)果有很高的一致性，幾乎全部通過了0.01的顯著性檢驗。1920前的50年學 生檢驗99%和95%的置信水平分別為0.28和0.36，1920年后的90年為0.21和0.27。1920年前的夏季NPO指數(shù)（圖9c）相關(guān)性略低，但也通過了0.05的顯著性檢驗。同時對比1920年前后的相關(guān)系數(shù)，可以看到盡管1920年后的序列有90年，長于1920年前的50年，但1920年后兩套資料的相關(guān)系數(shù)反而更高，說明1920年后兩套資料的一致性要遠強于1920年前的，這從圖中也可以直觀看出，到20世紀后半葉，隨著觀測站點的增多，特別是到了1970年代之后衛(wèi)星資料的引入，指數(shù)的相關(guān)性更有明顯提高。

圖7 中國105°E以東地區(qū)1979～2010年（a）冬、（b）春、（c）夏、（d）秋季節(jié)20CR氣溫插值結(jié)果與站點氣溫的空間相關(guān)。紅色為相關(guān)性通過了顯著性0.05檢驗的站點

觀察冬夏季各指數(shù)的近百年序列，可以看到明顯的年際和年代際振蕩，地球氣候系統(tǒng)作為一個整體，大氣濤動和大氣活動中心的變化不僅會影響局地氣候，還可通過遙相關(guān)產(chǎn)生非局地的影響。因此關(guān)注各氣候指數(shù)的百年序列對研究近百年的氣候變化有很好的指示作用。

冬季各指數(shù)中西伯利亞高壓是影響我國北方大部分地區(qū)冬季氣溫的重要氣候因子，在近百年中，其變化趨勢呈現(xiàn)低（約1871～1920年）、高（約1921～1980年）、低（1981～2010年）的年代際震蕩（圖8a，黑線）。與中國區(qū)域近百年冬季氣溫的三個主要冷暖期基本一致（范可和劉輝，2013）。同時也很好的刻畫出西伯利亞高壓1970年代末，全球變暖背景下的下降趨勢（Panagiotopoulos et al., 2005）和1990年后由于北極海冰減少等因素導致的上升趨勢（Wu et al., 2011），西伯利亞高壓90年后與極地環(huán)流聯(lián)系更為密切（范可和劉輝，2013）。從圖8b可以看出，兩套資料所刻畫的AL指數(shù)在1920年后有非常高的一致性（相關(guān)系數(shù)達0.96），1980年以來AL指數(shù)有明顯的年代際下降，反映了阿留申低壓的加深，這與前人的研究結(jié)果一致（Trenberth and Hurrell，1994；Rodionov et al., 2007）。本文按照郭其蘊（1994）定義冬季風指數(shù)（圖8d），用海陸氣壓差作為指標，主要反映冬季風在大陸東岸向南擴展的程度，與長江上游、華南和東部沿海的U形地帶的氣溫密切相關(guān)。1980年代中國南方氣溫略有下降，可以用冬季風指數(shù)在1980年代后有所加強來解釋。同時注意到1980年代后20CR描繪的冬季風指數(shù)略低于HadSLPv2資料。需要指出的是，若用西伯利亞高壓或東亞大槽為指標所定義冬季風指數(shù)（琚建華 et al., 2004），則指數(shù)主要與中國北方及中部大范圍地區(qū)的氣溫密切相關(guān)，在1980年代后為年代際減弱的趨勢，最近研究表明東亞東季風指數(shù)在1980年代后與ENSO的關(guān)系減弱（Wang and He，2012）。AO、NAO、NPO（圖8c、e、f）這三個大氣濤動聯(lián)系緊密，往往有著一致的變化或滯后相關(guān)（Delworth et al., 2000；孫建奇和王會軍，2005；Yuan and Sun， 2009；趙南和王啟祎，2010），如圖中可以看出三個指數(shù)在1920年、1970年附近都出現(xiàn)了較低的負位相，且1980年代都出現(xiàn)上升趨勢，而在1990年代又開始下降。三個指數(shù)1990年代的下降是我國及歐洲出現(xiàn)寒冬暴雪事件的關(guān)鍵因子（Wang et al., 2007；Cohen et al., 2010；Seager et al., 2010；Liu et al., 2012）。南方濤動表明東太平洋和印度洋的氣壓變化有東西向蹺蹺板式震蕩特征，其年代際變化主要反映了熱帶大氣環(huán)流年代際變化。SOI指數(shù)（圖 8g）近百年來兩套資料比較一致，1980年代后長時間位于負距平，近年來指數(shù)又有明顯回升。

夏季各指數(shù)兩套資料的相關(guān)性都整體低于冬季的。夏季風指數(shù)（圖9a）有明顯的年代際變 化，值得注意的是20CR資料在1950～1980年間描繪的夏季風略低于HadSLPv2，而在1980年代后又略高，因此從20CR資料不能很好地描繪1970年代末東亞夏季風的減弱（Wang, 2001）。這是可能是由于20CR資料不能很好地刻畫出對流層中上層的年代際變冷引起的（宋豐飛和周天軍，2012）。夏季NAO指數(shù)近百年來經(jīng)歷了兩正兩負的四個年代際的變化（圖9b）。1871～1920年和1940～1960年處于負距平，而1920～1940年和1960～2000年處于正距平。根據(jù)Greatbatch and Rong（2006）的研究，1960年之后持續(xù)較高的正距平僅是夏季NAO內(nèi)部低頻振蕩特征引起的。但1990年代末期，SNAO指數(shù)有一個很明顯的下降，這可能是導致地中海地區(qū)近些年干旱的原因之一，其突變的原因很可能與溫室氣體排放有關(guān)（Sun et al., 2008；Sun and Wang，2012；Bladé et al., 2012）。夏季NPO和SOI指數(shù) （圖9c、d）在1970年代末都出現(xiàn)一個急劇的下降趨勢，此后長期處在負距平，近年來又有所回升。這可能與1970年代末北半球夏季大氣環(huán)流的突變有關(guān)（Chang et al., 2000；Mantua and Hare，2002）。同時注意到夏季SOI指數(shù)在1970年代之后，20CR資料的結(jié)果要高于HadSLP的結(jié)果，而夏季SOI指數(shù)又和南亞夏季風關(guān)系緊密，因此20CR資料描繪南亞夏季風的能力還有待考察。

6 結(jié)論與討論

本文首先對比分析了歐亞地區(qū)20CR與NCEP2再分析資料的冬、春、夏、秋四季的氣候平均態(tài)及其差異場，并將20CR的氣溫和降水資料插值到中國站點上與中國臺站資料進行了對比。最后在百年尺度上與HadSLPv2資料進行指數(shù)對比，得出以下結(jié)論：

（1）冬季時，相比NCEP2資料，20CR資料所刻畫的極渦偏弱，東亞大槽偏弱偏西，海平面氣壓場上西伯利亞高壓位置偏北，我國冬季風偏弱。夏季時，相比NCEP2資料，20CR資料描繪的南亞高壓偏北，副高西伸點偏東且副高偏弱，伊朗高壓略強且位置偏北。海平面氣壓場上亞洲低壓更強，東北低渦較弱，夏季風偏強。歐亞地區(qū)20CR資料的SLP、H500、U200的四季氣候平均態(tài)與NCEP2資料的相比，均呈現(xiàn)北部中高緯度系統(tǒng)性偏高，南部中低緯度偏系統(tǒng)性低的特點，這可能與20CR資料在極地海岸地區(qū)海冰資料處理時產(chǎn)生的差錯有關(guān)。

（2）與中國東部站點資料的比對顯示，20CR再分析資料與我國站點資料存在一些差異。總體來說，對氣溫的刻畫能力偏低，而對降水的刻畫能力偏高，與站點氣溫的相關(guān)性好于與站點降水的相關(guān)性，東南沿海地區(qū)優(yōu)于內(nèi)陸地區(qū)。平均場和空間相關(guān)場結(jié)合來看，秋季氣溫和降水20CR與站點觀測值吻合最好。注意到最近有作者發(fā)展了一些中國區(qū)域的高分辨率格點觀測資料（Xu et al., 2009；吳佳和高學杰，2013），這些資料能更好地檢驗模式資料在中國地區(qū)的適用性，在未來的工作中將予以使用。

（3）20CR資料較好地刻畫近百年北半球冬、夏季的氣候指數(shù)北極濤動、北大西洋濤動、北太平洋濤動、東亞冬季風、阿留申低壓等的變率特征，很好地刻畫了阿留申低壓1970年代末的年代際增強，西伯利亞高壓1970年代末的下降和1990年后的上升趨勢，及北太平洋濤動、北大西洋濤動和北極濤動指數(shù)1970年代末期由負位相到正位相的年代際轉(zhuǎn)變和1990年代的變化，但不能再現(xiàn)1970年代末東亞夏季風的減弱。

Allan R, Ansell T. 2006. A new globally complete monthly historical mean sea level pressure data set (HadSLP2): 1850–2004 [J]. J. Climate, 19: 5816–5842.

Bengtsson L, Hodges K, Hagemann S. 2004. Sensitivity of the ERA-40 reanalysis to the observing system: Determination of the global atmospheric circulation from reduced observations [J]. Tellus, 56A: 456–471.

Bladé I, Liebmann B, Fortuny D, et al. 2012. Observed and simulated impacts of the summer NAO in Europe: Implications for projected drying in the Mediterranean region [J]. Climate Dyn., 39: 709–727.

Bromwich D, Fogt R. 2004. Strong trends in skill of the ERA-40 and NCEP/NCAR reanalyses in the high and midlatitudes of the Southern Hemisphere, 1958–2001 [J]. J. Climate, 17: 4603–4619.

Chang C P, Zhang Y S, Li T. 2000. Interannual and interdecadal variations of the East Asian summer monsoon and tropical Pacific SSTs. Part I: Roles of the subtropical ridge [J]. J. Climate, 13 (24): 4310–4325.

Cohen J, Foster J, Barlow M, et al. 2010. Winter 2009–2010: A case study of an extreme Arctic Oscillation event [J]. Geophys. Res. Lett., 37 (17), doi: 10.1029/2010GL044256.

Compo G P, Whitaker J S, Sardeshmukh P D, et al. 2011. The twentieth century reanalysis project [J]. Quart. J. Roy. Meteor. Soc., 137 (654): 1–28.

Delworth T L, Dixon K W. 2000. Implications of the recent trend in the Arctic/North Atlantic Oscillation for the North Atlantic thermohaline circulation [J]. J. Climate, 13 (21): 3721–3727.

范可, 劉輝. 2013. 近百年東亞冬季氣溫及其大氣環(huán)流變化型態(tài)[J]. 大氣科學, 37 (2): 383–394. Fan Ke, Liu Hui. 2013. Winter temperature over East Asia and atmospheric circulation patterns in the last 100 years [J]. Chinese Journal of Atmospheric Sciences (in Chinese), 37 (2): 383–394.

Fan K, Liu H. 2013. Evaluation of atmospheric circulation in the Southern Hemisphere in 20CRv2 [J]. Atmos. Oceanic Sci. Lett., 6 (5): 337–342.

Ferguson C R, Villarini G. 2012. Detecting inhomogeneities in the twentieth century reanalysis over the central United States [J]. J. Geophys. Res., 117 (D5): doi:10.1029/2011JD016988.

Greatbatch R J, Rong P P. 2006. Discrepancies between different Northern Hemisphere summer atmospheric data products [J]. J. Climate, 19 (7): 1261–1273.

郭其蘊. 1983. 東亞夏季風強度指數(shù)及其變化的分析[J]. 地理學報, 38 (3): 207?217. Guo Qiyun. 1983. The summer monsoom intensity index in East Asia and its variation [J]. Acta Geographica Sinica (in Chinese), 38 (3): 207?217.

郭其蘊. 1994. 東亞冬季風的變化與中國氣溫異常的關(guān)系[J]. 應用氣象學報, 5 (2): 218?225. Guo Qiyun. 1994. Relationship between the variation of East Asian winter monsoon and temperature anomalies in China [J]. Journal of Applied Meteorological Science (in Chinese), 5 (2): 218?225.

Hodges K I, Hoskins B J, Boyle J, et al. 2003. A comparison of recent reanalysis datasets using objective feature tracking: Storm tracks and tropical easterly waves [J]. Mon. Wea. Rev., 131: 2012–2037.

琚建華, 任菊章, 呂俊梅. 2004. 北極濤動年代際變化對東亞北部冬季氣溫增暖的影響[J]. 高原氣象, 2 (4): 429–434. Ju Jianhua, Ren Juzhang, Lü Junmei. 2004. Effect of interdecadal variation of arctic oscillation on temperature increasing in north of East Asian winter [J]. Plateau Meteorology (in Chinese), 23 (4): 429–434.

Kanamitsu M, Ebisuzaki W, Woollen J, et al. 2002. NCEP-DOE AMIP-II Reanalysis (R?2) [J]. Bull. Amer. Meteor. Soc., 83: 1631–1643.

Krueger O, Schenk F, Feser F, et al. 2012. Inconsistencies between long-term trends in storminess derived from the 20CR reanalysis and observations [J]. J. Climate, 26: 868–874.

Liu J, Curry J A, Wang H, et al. 2012. Impact of declining Arctic sea ice on winter snowfall [J]. Proceedings of the National Academy of Sciences of the United States of America, 109 (11): 4074–4079.

Mantua N J, Hare S R. 2002. The Pacific decadal oscillation [J]. J. Oceanography, 58 (1): 35–44.

Marshall G J. 2002. Trends in Antarctic geopotential height and temperature: A comparison between radiosonde and NCEP/NCAR reanalysis data [J]. J. Climate, 15: 659–674.

Marshall G J. 2003. Trends in the Southern Annular Mode from observations and reanalysis [J]. J. Climate, 16: 4134–4143.

Panagiotopoulos F, Shahgedanova M, Hannachi A, et al. 2005. Observed trends and teleconnections of the Siberian high: A recently declining center of action [J]. J. Climate, 18 (9): 1411–1422.

Rodionov S N, Bond N A, Overland J E. 2007. The Aleutian low, storm tracks, and winter climate variability in the Bering Sea [J]. Deep Sea Res. Part II, 54 (23): 2560–2577.

Seager R, Kushnir Y, Nakamura J, et al. 2010. Northern Hemisphere winter snow anomalies: ENSO, NAO and the winter of 2009/10 [J]. Geophys. Res. Lett., 37 (14), doi:10.1029/2010GL043830.

宋豐飛, 周天軍. 2012. 20CR再分析資料在東亞夏季風區(qū)的質(zhì)量評估[J]. 大氣科學, 36 (6): 1207–1222. Song Fengfei, Zhou Tianjun. 2012. Reliability of the 20CR reanalysis data in measuring the East Asian summer monsoon variability [J]. Chinese Journal of Atmospheric Sciences (in Chinese), 36 (6): 1207–1222.

Sun J Q, Wang H J. 2012. Changes of the connection between the summer North Atlantic Oscillation and the East Asian summer rainfall [J]. J. Geophys. Res., 117 (D8), doi:10.1029/2012JD017482.

孫建奇, 王會軍. 2005. 北極濤動與太平洋年代際振蕩的關(guān)系[J]. 科學通報, 50 (15): 1648–1653. Sun Jianqi, Wang Huijun. 2005. Relationship between Arctic oscillation and Pacific decadal oscillation on decadal timescale [J]. Chinese Science Bulletin, 51 (1): 75?79.

Sun J Q, Wang H, Yuan W. 2008. Decadal variations of the relationship between the summer North Atlantic Oscillation and middle East Asian air temperature [J]. J. Geophys. Res., 113 (D15), doi: 10.1029/2007JD009626.

Trenberth K E. 1990. Recent observed interdecadal climate changes in the Northern Hemisphere [J]. Bull. Amer. Meteor. Soc., 71 (7): 988–993.

Trenberth K E, Hurrell J. 1994. Decadal atmosphere–ocean variations in the Pacific [J]. Climate Dyn., 9 (6): 303–319.

Wang H J. 2001. The weakening of the Asian monsoon circulation after the end of 1970’s [J]. Adv. Atmos. Sci., 18 (3): 376–386.

Wang H J. 2002. The instability of the East Asian summer monsoon–ENSO relations [J]. Adv. Atmos. Sci., 19 (1): 1–11.

Wang H J, He S P. 2012. Weakening relationship between East Asian winter monsoon and ENSO after mid-1970s [J]. Chin. Sci. Bull., 57 (27): 3535– 3540.

Wang L, Chen W, Huang R H. 2007. Changes in the variability of North Pacific Oscillation around 1975/1976 and its relationship with East Asian winter climate [J]. J. Geophys. Res., 112 (D11), doi:10.1029/2006JD008054.

Wu B Y, Su J Z, Zhang R H. 2011. Effects of autumn–winter Arctic sea ice on winter Siberian high [J]. Chin. Sci. Bull., 56 (30): 3220–3228.

吳佳, 高學杰. 2013. 一套格點化的中國區(qū)域逐日觀測資料及與其它資料的對比[J]. 地球物理學報, 56 (4): 1102–1111. Wu Jia, Gao Xuejie. 2013. A gridded daily observation dataset over China region and comparison with the other datasets [J]. Chinese Journal of Geophysics (in Chinese), 56 (4): 1102–1111.

Wu R G, Kinter J L, Kirtman B P. 2005. Discrepancy of interdecadal changes in the Asian region among the NCEP?NCAR reanalysis, objective analyses, and observations [J]. J. Climate, 18: 3048–3067.

徐影, 丁一匯. 2001. 美國NCEP/NCAR近50年全球再分析資料在我國氣候變化研究中可信度的初步分析[J]. 應用氣象學報, 12 (3): 337–347. Xu Ying, Ding Yihui. 2001. Confidence analysis of NCEP /NCAR 50-year global reanalyzed data in climate change research in China [J]. Journal of Applied Meteorological Science (in Chinese), 12 (3): 337–347.

Xu Y, Gao X J, Shen Y, et al. 2009. A daily temperature dataset over China and its application in validating a RCM simulation [J]. Adv. Atmos. Sci., 26 (4): 763–772.

Yang S, Lau K, Kim K. 2002. Variations of the East Asian jet stream and Asian–Pacific–American winter climate anomalies [J]. J. Climate, 15: 306–325.

Yuan W, Sun J Q. 2009. Enhancement of the summer North Atlantic Oscillation influence on Northern Hemisphere air temperature [J]. Adv. Atmos. Sci., 26(6): 1209–1214.

趙南, 王啟祎. 2010. 一個觀察北極濤動與北大西洋濤動關(guān)系的典型個例[J]. 氣象學報, 68: 847–854. Zhao Nan, Wang Qiyi. 2010. A typical case for observation of the relationship between the Arctic Oscillation and the North Atlantic Oscillation [J]. Acta Meteorologica Sinica (in Chinese), 68: 847–854.

趙天保, 符淙斌. 2009. 應用探空觀測資料評估幾類再分析資料在中國區(qū)域的適用性[J]. 大氣科學, 33 (3): 634–648. Zhao Tianbao, Fu Congbin. 2009. Applicability evaluation for several reanalysis datasets using the upper-air observations over China [J]. Chinese Journal of Atmospheric Sciences (in Chinese), 33 (3): 634–648.

Eurasia Seasonal Circulation Climatology and Variability: Evaluation of 20CR Reanalysis Data in Eurasia and East China

LIU Huiand FAN Ke

1,,100029 2,100049

This study evaluates the characterization capabilities of seasonal Eurasia climatology and climate variability of the newly released twentieth century reanalysis (20CR) data by the U.S. National Oceanic and Atmospheric Administration. The results shows that in the Eurasia area, the 20CR data exhibited higher values in the northern high latitudes and lower values in the southern low latitudes than those recorded by version II of the National Centers for Environmental Prediction reanalysis data, resulting in weaker winter monsoon and a stronger summer monsoon circulations. This may be partly attributed to the handling of specified sea-ice concentration in the polar coastal area in the 20CR data. Comparison with the station observation data of eastern China indicates that the 20CR data represents lower air temperature and a higher precipitation in the four seasons whereas the correlation of temperature between 20CR reanalysis data and the station observation is higher than that of precipitation. In addition, the variability of temperature and precipitation over the southeastern coastal area were captured better than those of the inland areas in 20CR reanalysis data. In particular, the highest consistency of the temporal–spatial feature of autumn temperature and precipitation in eastern China was between 20CR data and station observations. Major climate variability indices such as Arctic oscillation (AO), North Atlantic oscillation (NAO), North Pacific oscillation (NPO), Aleutian low, Siberian high, East Asian summer monsoon, East Asian winter monsoon, and Southern oscillation (SO) recorded during the last hundred years have also been evaluated in 20CR reanalysis data. The 20CR reanalysis data and the Hadley Centre sea level pressure dataset both showed the decadal change characteristics that the winter Aleutian low deepened during the late of 1970s, whereas the Siberian high weakened during the same period and intensified in the 1990s. In addition, negative phases of winter AO, NAO, and NPO transferred to positive phases during the late 1970s.

20CR reanalysis data, Eurasian climatology, Climate variability, Decadal change, East China

1006?9895(2014)03?0469?15

P468

A

10.3878/j.issn.1006-9895.2013.13182

2013?05?27，

2013?12?02收修定稿

中國科學院戰(zhàn)略性先導科技專項XDA05090426，全球變化研究國家重大科學研究計劃2010CB950304，中國科學院知識創(chuàng)新工程重要方向項目青年人才類KZCX2-YW-QN202

劉輝，女，1987年出生，博士研究生，主要從事氣候變化研究。E-mail: liuh@mail.iap.ac.cn

劉輝, 范可. 2014. 歐亞季節(jié)環(huán)流氣候平均態(tài)和氣候變率：20CR再分析資料在歐亞及中國東部地區(qū)的評估[J]. 大氣科學, 38 (3): 469–483, doi:10.3878/j. issn.1006-9895.2013.13182. Liu Hui, Fan Ke. 2014. Eurasia seasonal circulation climatology and variability: Evaluation of 20CR reanalysis data in Eurasia and East China [J]. Chinese Journal of Atmospheric Sciences (in Chinese), 38 (3): 469–483.