基于CFSv2對夏季西太平洋副高逐日預測能力的檢驗評估

付炳秀， 王同美， 韋翠， 鄧志明

(1．仁化縣氣象局，廣東仁化　512300；2．中山大學大氣科學系，廣東廣州　510275；3.清遠市氣象局，廣東清遠　511510)

?

基于CFSv2對夏季西太平洋副高逐日預測能力的檢驗評估

付炳秀1， 王同美2， 韋翠3， 鄧志明1

(1．仁化縣氣象局，廣東仁化512300；2．中山大學大氣科學系，廣東廣州510275；3.清遠市氣象局，廣東清遠511510)

摘要：利用NCEP的氣候預報系統(tǒng)第2版(CFSv2)提供的2000—2010年夏季逐日500 hPa位勢高度歷史回報試驗結果，通過均方根誤差分析等方法，檢驗評估了該系統(tǒng)對夏季500 hPa西太平洋副高的預測能力，發(fā)現CFSv2對夏季西太平洋500 hPa位勢高度場的整體相關性較高。CFSv2對副高西脊點和脊線位置的1～7 d的預報準確率在61.8%以上，對副高面積指數、強度指數以及北界位置6～42 h預報時效的預報準確率在60%以上。此外，對副高面積指數、強度值和北界位置的預測結果作功率譜分析，得出CFSv2對500 hPa副高強度指數1～7 d預測能力的變化在6—8月存在著顯著的2.2～3.3 d的變化周期。

關鍵詞：氣象預報； 氣候預報系統(tǒng)第2版； 西太平洋副高； 預測檢驗評估

美國國家環(huán)境預報中心(NCEP)的氣候預測系統(tǒng)第2版(CFSv2)于2011年3月開始進行業(yè)務實時預測，并提供了1982—2010年的回報試驗結果[1]。目前基于CFSv2的研究工作主要集中于氣溫、降雨以及季節(jié)振蕩方面的預測能力研究[2-6]。在大尺度環(huán)流的回報檢驗方面，陳官軍等[7]曾利用CFS提供的歷史回報試驗結果檢驗并評估了其對夏季東亞地區(qū)大氣環(huán)流的預報技巧和系統(tǒng)誤差，而對一些在天氣預測方面具有指示意義的具體環(huán)流系統(tǒng)的預測評估，如西太平洋副熱帶高壓這一對我國夏季天氣及短期氣候預測都有著非常重要影響的系統(tǒng)，暫無相關工作報告發(fā)表。

本研究將采用CFSv2的預報產品，對夏季500 hPa副高的中短期歷史回報試驗結果進行檢驗，找出CFSv2預報模式對500 hPa副高的預測上存在的偏差，從而為預報員分析數值預報產品提供客觀依據，向模式研發(fā)人員提供更多的反饋信息。

1資料與方法

1.1資料

CFSv2所提供的2000—2010年的夏季(6—8月)逐日歷史回報試驗結果資料，預報時效從6 ～162 h，間隔6 h，本研究使用此歷史回報試驗結果資料作為副高預測場。NCEP提供的2000—2010年的夏季(6—8月)FNL 500 hPa位勢高度場再分析資料，作為CFSv2對500 hPa副高預測檢驗的再分析場。資料的水平網格分辨率均為1°×1°。

1.2方法

1)異常相關系數(abnormal correlation coefficient，ACC)：異常相關系數越趨近于1.0，表示2個變量之間相關關系密切程度越高，2個變量的變化趨勢也越相近。

2)均方根誤差(root．mean—square error，RMSE)：均方根誤差用于衡量觀測值同真值之間的偏差。

3)功率譜分析：用于了解對副高的預測能力在6—8月的低頻信號分布。

2預測能力檢驗

2.1夏季500 hPa位勢高度預測場和再分析場的相關分析

在天氣預報分析過程中，588 dagpm特征線是劃分副高脊線、西脊點的重要參考線。因此，有必要了解預測檢驗工作是否有位勢高度場的高相關性進行匹配。本研究對夏季500 hPa位勢高度場的預測場和再分析場進行相關分析，由于對500 hPa副高預測的檢驗主要區(qū)域在60°—180°E，10°N以北，故選取該區(qū)域進行相關分析。

圖1是從各預報時效對應的相關系數分布圖中選取的4個預報時效所對應的相關系數分布圖，在中緯度10°—50°N，80°—180°E區(qū)域的相關性都較高，在靠近赤道區(qū)域相關性都很差，而在30°N以北地區(qū)，相關性隨著預報時效的增長有明顯的下降。CFSv2歷史回報試驗結果對500 hPa位勢高度場的預報能力在空間上隨預報時效的增長而下降。在10°—50°N，80°—180°E的范圍內，預測場與再分析場的相關系數也基本達到0.65，通過α=0.05的可信度檢驗，CFSv2對西太平洋500 hPa副高的預測值得檢驗。

圖1　預測場與分析場的相關系數分布

2.2對500 hPa西太平洋副高的預測檢驗

1)副高指數定義。

目前我國業(yè)務預報的氣候預測中使用副高的5個指數(北界位置、脊線位置、面積指數、強度指數、西伸脊點)來表征副高。本研究參考中央氣象臺長期預報組在長期天氣預報技術經驗總結工作中副高指數的定義[8]，陳官軍等人在基于低頻振蕩特征的夏季江淮持續(xù)性降水延伸期預報方法工作中副高強度指數、脊線位置以及西脊點的定義[9]，并結合CFSv2提供的歷史回報試驗結果數據以及編程計算的特點，對引用的各種副高指數定義修改后定義如下：

(1)面積指數(Subtropical High Area index，SubHA)：在網格數據為1°×1°的500 hPa位勢高度圖上，110°—180°E，10°N以北范圍內，位勢高度大于588 dagpm網格點數稱為副高的面積指數。

(2)強度指數(Subtropical High Intensity index，SubHI)：10°N以北，110°—150°E范圍內，大于588 dagpm的編碼之和為副高的強度指數。(對≥588 并<589 dagpm的數值編碼為1，≥589 并<590 dagpm的數值編碼為2，其余類推)

(3)西脊點(Western end of the Subtropical High ridge，SubHW)：取90°—180°E范圍內，位勢高度值大于等于588 dagpm并且位置最西的格點所在經度定義為副高的西脊點。

(4)脊線位置(The position of Subtropical High Ridge，SubHR)：在110°—150°E范圍內，以赤道500 hPa位勢高度值為基準，求北半球各個維度與赤道500 hPa位勢高度值偏差最大的格點的緯度平均值，將此平均值定義為脊線位置。

(5)北界位置(The northern boundary of Subtropical High，SubHN)：在110°—150°E范圍內，分別求各個經度所對應的大于等于588 dagpm并且位于該經度最北位置的格點的緯度，取這些格點的緯度平均值定義為北界位置。

2)對副高指數的預測檢驗。

(1)副高指數均方根誤差檢驗。

均方根誤差用來衡量觀測值同真值之間的偏差，誤差越小，表明預報越準確。基于預測場與再分析場的相關系數基本達到0.65，CFSv2對西太平洋500 hPa副高的預測值得檢驗的結論，本研究對6—8月份各預報時段的預測場和再分析場的5個副高指數分別求RMSE，得到一個隨6—8月以及預報時效變化的二維均方根誤差分布圖(圖2)。

圖2　各副高指數RMSE隨時間分布

在圖2各副高指數RMSE隨時間分布圖中，各副高指數的RMSE值均隨著預報時效的增長而增大，表明預測場對5個副高指數的預測能力均隨著預報時效的增長而下降。副高面積指數、強度指數和北界位置的RMSE值分布在6—8月則呈現先增后減的變化，說明CFSv2對副高面積指數和強度指數以及北界位置的預測能力在6—8月先減后增。預測場對副高西脊點的預測能力在6—8月期間遞增，對脊線位置的預測能力在6—8月期間遞減。從各副高指數的RMSE值分布圖對比，預測場對副高西脊點和脊線位置的預測效果比對副高面積指數和強度指數以及北界位置的預測效果更好。

此外，CFSv2對副高面積指數、強度指數和北界位置的預測能力在6—8月的變化中有比較明顯的、有規(guī)律的波動，這種波動在72 ～144 h預報時效的預報中更為明顯。

(2)預測場對各副高指數的預測評分。

根據各副高指數的均方根誤差時間分布情況，選定面積指數的誤差范圍≤100，強度指數的誤差范圍為≤200，西脊點的誤差范圍≤5，脊線位置的誤差范圍≤3，北界位置的誤差范圍≤6，再分別對各副高指數在預測場中的值和再分析場中的值做差值后取其絕對值，若所得的絕對值在選定的相應誤差范圍內則認為預報成功。為更直觀了解預測評分變化，本研究分別將統(tǒng)計的各副高指數預報成功的次數與總預報次數求比值，并將所得比值作為預測評分，得到預測評分圖(圖3)。圖3是CFSv2對夏季500 hPa 5個副高指數6～162 h預報時效的預測評分，預報時效間隔為6 h，預測評分滿分為1分。此外，圖3中評分數值乘以100%后即為預報準確率。

圖3　預測場對各副高指數6～162 h預報時效的預測評分

從圖3中的數據可以看出，預測場對夏季500 hPa西太平洋5個副高指數的預測評分在整體上均隨著預報時效的增長呈現著下降的趨勢。副高西脊點和脊線位置的預測評分最低分為0.618和0.665，其預報成功的誤差范圍分別是5個經度和3個緯度。CFSv2對夏季副高西脊點和脊線位置在7 d內的預報中整體的預報評分較高，誤差范圍較小。CFSv2對西脊點和脊線位置的1～7 d的預報準確率在61.8%以上，對面積指數、強度指數以及北界位置6～42 h預報時效的預報準確率基本在60%以上。

(3)預測能力的頻譜分析。

基于在均方根誤差時間分布的分析中預測場對副高面積指數、強度指數以及北界位置的預測能力在6—8月份存在較明顯且有規(guī)律的波動特點，本研究對副高面積指數、強度指數以及北界位置6—8月的均方根誤差值進行功率譜分析，并將相應周期的譜密度除以α=0.05的紅噪音標準譜，繪制成6—8月CFSv2對500 hPa副高面積指數、強度指數和北界位置預測能力的譜比值分布(圖4)，當譜比值>1.0時代表該周期超過0.05顯著性檢驗水平，比值越大表明周期越顯著。

副高強度指數預測能力的譜比值分布圖(圖4)顯示，在6～162 h預報中，CFSv2對6—8月500 hPa副高強度指數的預測能力變化存在著通過檢驗的2.2～3.3 d的變化周期，其中最顯著的變化周期為2.5和2.9 d。在副高面積、北界位置指數預測能力的譜比值分布圖(圖略)中，CFSv2對副高面積指數和北界位置的預測能力在6—8月期間通過檢驗的變化周期零散分布于個別預報時效中，變化周期不明顯。CFSv2對副高面積指數的預測能力變化在42、66、78、90、96和138 h預報中有通過檢驗的2.5～3.3 d的變化周期。CFSv2對6—8月500 hPa北界位置的預測能力變化在42和66 h預報中都存在著通過檢驗的3.3 d左右的變化周期，在60和78 h預報中存在通過檢驗的2.5 d左右的周期。通過對副高面積指數、強度指數和北界位置預測能力的譜比值分布圖的對比分析，發(fā)現在副高面積指數、強度指數、北界位置的周期顯著性中，CFSv2對強度指數的預測能力變化周期最為明顯。

圖4　6—8月CFSv2對500 hPa副高強度預測能力的譜比值分布

3結論

1)在10°—50°N，80°—180°E的范圍內，預測場與再分析場的相關系數基本達到0.65，通過了α=0.05的可信度檢驗，CFSv2對西太平洋500 hPa副高的預測值得檢驗。CFSv2歷史回報試驗結果對500 hPa位勢高度場的預報能力在空間上隨預報時效的增長而下降。

2)CFSv2歷史回報試驗結果對5個副高指數的預測能力都隨著預報時效的增長而下降。CFSv2對副高面積指數、強度指數和北界位置的預測能力在6—8月期間先減后增，對副高西脊點的預測能力在6—8月期間遞增，對副高脊線位置的預測能力在6—8月期間遞減。

3)預測場對夏季西太平洋500 hPa 5個副高指數的預測評分在整體上均隨著預報時效的增長呈現著下降的趨勢。CFSv2對副高西脊點和脊線位置的1～7 d的預報準確率在61.8%以上，對副高面積指數、強度指數以及北界位置6～42 h預報時效的預報準確率在60%以上。

4)在6～162 h各個預報時段中，CFSv2對副高強度指數的預測能力在6—8月期間有2.2～3.3 d的變化周期，其中最顯著的變化周期為2.5和2.9 d。CFSv2對副高面積指數和副高北界位置的預測能力在6—8月期間通過檢驗的變化周期零散分布于個別預報時效中，變化周期不明顯。

參考文獻:

[1]Suranjana Saha，Shrinivas Moorthi，Xingren Wu，et al.The NCEP Climate Forecast System version 2[J].Journal of Climate，2014，27(6)：2185-2208.

[2] Wang Wanqiu， Hung Mengpai，Scott J Weaver，et al.MJO prediction in the NCEP Climate Forecast System version2[J].Climate Dynamics，2014，42(9)：2509-2520.

[3] Zhang Qin，Chul-Su Shin，Huug van den Dool，et al.CFSv2 prediction skill of stratospheric temperature anomalies[J].Climate Dynamics，2013，41(7)：2231-2249.

[4] Luo Lifeng，Tang Wei，Lin Zhaohui，et al.Evaluation of summer temperature and precipitation predictions from NCEP CFSv2 retrospective forecast over China[J].Climate Dynamics，2013，41(8)：2213-2230.

[5]Abhilash S，Sahai A K ，Borah N，et al.Prediction and monitoring of monsoon intraseasonal oscillations over Indian monsoon region in an ensemble prediction system using CFSv2[J].Climate Dynamics，2014，42(9)：2801-2815.

[6]陳官軍，魏鳳英，鞏遠發(fā).NCEP/CFS模式對東亞夏季延伸預報的檢驗評估[J].應用氣象學報，2010，21(6)：659-670.

[7]中央氣象臺長期預報組.長期天氣預報技術經驗總結(附錄)[Z].北京：中央氣象臺，1976.

[8]陳官軍，魏鳳英.基于低頻振蕩特征的夏季江淮持續(xù)性降水延伸期預報方法[J].大氣科學，2012，36(3)：633-644.

[9]李江南，蒙偉光，王安宇，等.西太平洋副熱帶高壓強度和位置的氣候特征[J].熱帶地理，2003，23(1)：35-39.

收稿日期：2015-08-22

作者簡介：付炳秀(1990年生)，女，本科，主要從事短期天氣預報工作和地面氣象觀測。E-mail：fubingx@126.com

中圖分類號：P456

文獻標識碼：A

doi:10.3969/j.issn.1007-6190.2016.02.004

Testing and Assessment of Capabilities of Day-to-Day Predicting ofSummertime West Pacific Subtropical High Based on CFS v2

FUBing-xiu1,WANGTong-mei2,WEICui3,DENGZhi-ming1

(1. Meteorological Bureau of Renhua County, Renhua 512300; Department of Atmospheric Sciences,Sun Yet-Sen University, Guangzhou 510275; Meteorological Bureau of Qingyuan City, Qingyuan 511510)

Abstract：Using historical hindcasts of day-to-day 500-hPa geopotential heights from 2000 to 2010 given by the second version of Climate Forecasting System of NCEP (CFSv2) and with the methods such as the root mean square analysis, we tested and assessed its capabilities of predicting summertime day-to-day 500-hPa geopotential heights, with the finding that CFSv2 has high overall correlation with the geopotential field for the same level. Its accuracy is more than 61.8% in predicting the westernmost point and position of the subtropical high ridge one to seven days ahead and more than 60% in predicting the area and intensity indexes as well as the position of the northern boundary of the subtropical high for the coming 6-42 h. Besides, CFSv2 varies with a significant period of 2.2 to 3.3 d in June to August when it is used to predict the intensity index at 500 hPa.

Key words：weather forecast; second version of Climate Forecasting System; western Pacific subtropical high; testing and assessment of prediction

付炳秀， 王同美， 韋翠， 等.基于CFSv2對夏季西太平洋副高逐日預測能力的檢驗評估[J].廣東氣象，2016,38(2)：15-19.