閆彩霞, 金榮花
(1.成都信息工程學院大氣科學學院,四川成都610225;2.國家氣象中心,北京100081)
在南、北半球副熱帶地區(qū)對流層上部各存在著一支強而窄的高速氣流帶,稱之為副熱帶西風急流,其位置的南北移動以及強度的異??梢砸鸫蠓秶沫h(huán)流異常,進而導致中國乃至東亞地區(qū)天氣氣候異常,與東亞大氣環(huán)流的季節(jié)轉(zhuǎn)換、亞洲夏季風的爆發(fā)、中國東部地區(qū)夏季降水等有著密切的聯(lián)系[1-8]。因此,東亞副熱帶西風急流監(jiān)測與分析在預報預測中至關重要。在中國的預測預報業(yè)務中,建立了東亞夏季風、西太平洋副熱帶高壓以及歐亞阻塞高壓等重要大型環(huán)流系統(tǒng)的監(jiān)測指標[9],但是,過去受高空資料質(zhì)量和數(shù)值模式性能的限制,對東亞副熱帶急流及其物理機制的認識還比較局限,尚未能建立東亞副熱帶急流監(jiān)測指標。近些年,隨著大氣科學研究向三維結構更深入更細致方向發(fā)展,東亞副熱帶急流已經(jīng)成為普遍關注的重要系統(tǒng),為實時獲取系統(tǒng)的、規(guī)范的東亞副熱帶急流的長序列資料,滿足業(yè)務監(jiān)測和科學研究的需要,建立東亞副熱帶急流的監(jiān)測指標是目前迫切需要解決的問題。
以客觀、定量和自動方式實現(xiàn)對東亞副熱帶急流的監(jiān)測,首先要研究東亞副熱帶急流各種特征參數(shù)的表征方法。以往研究工作中,提出了很多表征副熱帶急流的方法[10-20],但是,研究的目的和分析的角度不同,所構造的方法各不相同,應用于業(yè)務存在局限性。因此有必要對以往東亞副熱帶西風急流客觀定量表征方法進行系統(tǒng)分析和客觀評價,并通過分析指出以往各種客觀定量化表征西風急流的方法的優(yōu)點與適用范圍,以期為科研工作者和業(yè)務技術人員合理應用東亞副熱帶急流指數(shù)提供參考,論文通過總結歸納以往研究工作中的東亞副熱帶急流特征參數(shù)表征方法,基于各種方法的性能分析結果,結合東亞副熱帶急流特征及其對中國降水的影響,對比分析東亞副熱帶急流特征參數(shù)表征方法的優(yōu)缺點,為實現(xiàn)對東亞副熱帶急流的客觀實時監(jiān)測和診斷分析提供參考和依據(jù)。
分析所使用資料為NCEP/NCAR逐日再分析等壓面資料,要素為200hPa水平緯向風和水平經(jīng)向風,水平分辨率為2.5°×2.5°,序列長度為1960~2009年。國家氣象信息中心提供中國714站20時(北京時)日降水量資料,文中采用Cressman方法將降水量插值到1°×1°網(wǎng)格,序列長度為1960~2009年。
為了便于比較,在對各種東亞副熱帶急流客觀定量表征方法進行數(shù)值反演時,研究時間統(tǒng)一定為夏季6~8月,2005年降雨量接近常年同期或偏多,西風急流北跳時間顯著,對雨帶推移有較大影響[21]。2005年天氣實況分析表明中國東部地區(qū)三次降水帶北進較為顯著,因此選擇2005年作為逐日演變分析個例,季節(jié)尺度年際變化分析序列統(tǒng)一為1960~2009年。以往研究工作中提出了很多客觀定量表征東亞副熱帶急流的方法,按照其功能劃分主要有位置、強度兩類特征參數(shù)。下面對原作者的方法逐一介紹,并進行反演和性能分析。
位置參數(shù)是用以描述副熱帶西風急流位置移動的參數(shù),其演變主要反映在急流位置南北移動變化上,以往客觀定量表征方法歸納為3種。
況雪源等[10]采用200hPa等壓面上(70°E~120°E,30°N~50°N)區(qū)域內(nèi)最大西風所在緯度的平均值表示東亞副熱帶急流的位置,用參數(shù) A1表示。黃安麗等[11]采用動力學方法分析對流層高、低空急流耦合作用時,確定以105°E~145°E,45°N 以南作為統(tǒng)計區(qū),300hPa等壓面上風速值大于30m·s-1的地區(qū)作為高空急流區(qū),用高空急流區(qū)風速最大值位置表示東亞副熱帶急流的位置,用參數(shù) A2表示。如圖1所示采用急流軸線方法,計算了2005年6~8月急流位置變化(圖2)及指數(shù)與同期東部地區(qū)夏季主要降水中心位置(100°E~130°E)的關系,對比結果,A1方法采用了更符合經(jīng)典的東亞副熱帶急流定義的200hPa等壓面資料,能夠準確描述夏季急流位置的南北變化。指數(shù)與中心雨帶位置相關系數(shù)達0.52,能較好地反映急流對雨帶南北移動的影響,定義合理且有明確的物理意義,方法簡單易行、意義清晰、表現(xiàn)直觀,適用于逐日實時監(jiān)測和診斷分析。但是,由于采用統(tǒng)計區(qū)內(nèi)各經(jīng)度上最大西風的位置平均來表征急流軸位置,對于南北兩支西風急流并存并且北支強于南支情況,對于盛夏(春季)急流軸北抬(南落)至統(tǒng)計區(qū)以北(以南)的情形,如2005年7月19日,該參數(shù)無法捕捉到真實信息。A2指數(shù)變化趨勢跟雨帶位置移動具有較好的一致性,五點滑動平均后相關系數(shù)達0.56,呈現(xiàn)顯著的相關,方法簡單易行、意義清晰,對中國東部地區(qū)降水位置的南北移動具有較好的指示意義。但是此方法采用的副熱帶西風急流為300hPa層面上,與經(jīng)典高空急流定義不相符合。
由于海陸分布不均及青藏高原的影響,急流軸的季節(jié)變化具有東西向的不一致性[10]。兩種定義所選擇的急流軸關鍵區(qū)不同,A1方法確定的關鍵區(qū)主要位于中國中西部,A2方法選擇區(qū)域側(cè)重于中國東部地區(qū)(江淮地區(qū)),2個指數(shù)相同點在于對急流表示本質(zhì)上都采用急流最大風速代表的軸線表示。
圖1 軸線法表征的東亞副熱帶西風急流軸位置參數(shù)與同期中國東部降水帶位置演變(參數(shù)A1、A2,多年平均位置以及降水中心位置逐日演變,單位:緯度)
圖2 2005年7月19日200hPa緯向風等值線(單位:m·s-1)
劉飛等[12]對1958~2000年43年東亞地區(qū)夏季200hPa緯向風場進行了EOF分解,第一特征向量反映西風急流繞夏季平均急流軸的南北變動,因此定義第一模態(tài)對應的標準化時間系數(shù)為東亞夏季副熱帶急流位置參數(shù)。況雪源等[13]同樣采用此方法對北半球冬季副熱帶西風急流指數(shù)進行了定義及探討,結果表明此指數(shù)對副熱帶西風急流變動有較好的指示意義,并與長江中下游地區(qū)降水相關最為顯著。故以劉飛方法為代表,用參數(shù) A3表示。
廖清海等[14]根據(jù)東亞地區(qū)急流月平均位置的方差分布,選取各月南北兩個活動中心對應的正規(guī)化的200hPa緯向風的差作為當月東亞副熱帶西風急流指數(shù),并以此獲得東亞地區(qū)夏季(6~8月)對應的副熱帶急流位置參數(shù)。楊蓮梅等[15]和Lin等[16]均采用指定范圍內(nèi),西風急流氣候平均位置南北各5°或10°的兩個區(qū)域內(nèi)的200hPa緯向風的差反映東亞西風急流位置,正(負)值表明東亞西風急流偏南(北)。Liang等[17]對冬季200hPa緯向風場進行EOF分解,第一模態(tài)呈現(xiàn)出南北不同性質(zhì)的正負中心結構,在30°N以北是正值中心,中心范圍為(100°E ~ 130°E,30°N ~ 40°N),以南是負值中心,中心范圍為(90°E ~ 120°E,10°N ~ 25°N),因此定義南北活動中心的緯向風距平的標準化值相減,正(負)值表明東亞西風急流偏南(北)。以廖清海方法為代表,用參數(shù) A4表示。
采用上述2種急流位置參數(shù)表征方法,以1960~2009年6~8月平均為例進行了西風急流夏季平均位置逐年參數(shù)數(shù)值反演(圖3-圖6),對比兩種方法性能發(fā)現(xiàn):2種參數(shù)的優(yōu)點在于抓住了西風急流典型分布形式的時間變化特征,定性表示出西風急流位置較平均態(tài)偏南或偏北特征,對于氣候評估和診斷季節(jié)內(nèi)副熱帶急流的位置特征不失為好方法。但是,不能用以定量形式給出急流軸具體位置,此外,這兩種方法對統(tǒng)計區(qū)域比較敏感,區(qū)域變化將可能導致結果發(fā)生較大變化。從圖5,圖6兩種指數(shù)與降水的關系圖中可以看出 A3方法與長江中下游地區(qū)降水呈顯著的負相關,最高相關系數(shù)達-0.53,與華北東南沿海地區(qū)降水呈現(xiàn)顯著的正相關,通過0.01顯著性水平檢驗。A4指數(shù)指示意義與 A3正好相反,與長江中下游地區(qū)降水呈現(xiàn)顯著正相關,但相關性不如 A3方法顯著。因此在對降水的指示作用方面 A3方法更勝一籌。東亞副熱帶西風急流位置參數(shù)綜合對比分析結果見表1。
考慮兩種方法的優(yōu)點,技術規(guī)避缺點,建議使用此種方法時,構造指數(shù)如下:可采用EOF分解,得出夏季西風急流活動評價位置的南北中心,用兩中心緯向風的差作為當月東亞副熱帶西風急流指數(shù)。
圖3 參數(shù) A3表征的東亞副熱帶急流位置參數(shù)的演變情況
圖4 表征的東亞副熱帶急流位置參數(shù)的演變情況
圖5 參數(shù) A3表征的演變情況與中國夏季同期降水的相關
圖6 參數(shù) A4表征的演變情況與中國夏季同期降水的相關
表1給出了上述3種方法的主要特征、分析尺度及適用范圍,由表可見,表征位置指數(shù)的方法主要有軸線法、EOF分解法以及區(qū)域風速差法。實況反演及降水相關分析結果表明,急流軸線法能較好地反映急流的經(jīng)向移動,與中國東部雨帶進退有較為一致的對應關系,定義合理且有明確的物理意義,方法簡單易行、意義清晰、表現(xiàn)直觀,適用于逐日實時監(jiān)測和診斷分析。EOF分解法與區(qū)域風速差法抓住了西風急流典型分布形式的時間變化特征,定性表示出西風急流位置較平均態(tài)偏南或偏北特征,對于氣候評估和診斷季節(jié)內(nèi)副熱帶急流的位置特征不失為好方法;方法分析尺度不同,適用范圍不同,研究者應根據(jù)研究需要,恰當?shù)倪M行選擇。
表1 東亞副熱帶西風急流位置參數(shù)表征方法比較
強度參數(shù)是表征東亞副熱帶西風急流強度的特征量,急流強弱反映在雨量大小上,其客觀定量表征方法歸納為兩種。
史有瑜等[18]定義200hPa等壓面上給定區(qū)域內(nèi)(30°N~50°N,90°E~110°E)緯向風速的區(qū)域平均值為西風強度指數(shù)。楊蓮梅等[15]也采用同樣的方法,給定區(qū)域選為(30°N~45°N,100°E~140°E),并對區(qū)域平均的 200hPa緯向風之和進行標準化處理。Yang等[19]通過統(tǒng)計分析1968~2000年冬季東亞副熱帶西風急流極大值的位置、頻次和平均值,得出西風急流核通常位于(30°N~35°N,130°E~160°E)區(qū)域,所以定義冬季西風急流強度參數(shù)為該區(qū)域緯向風值的平均值。毛睿等[20]統(tǒng)計各年冬季緯向風極大值的出現(xiàn)范圍,范圍為(30°N~35°N,127.5°E~155°E),因此定義西風急流強度參數(shù)為這個區(qū)域冬季200hPa緯向風平均值的標準化值。
上述方法均采用區(qū)域平均風速來表示急流強度,以2005年6~8月為例,對急流強度參數(shù)B1進行逐日數(shù)值反演,為進一步驗證強度指標對于降水監(jiān)測的意義,計算了該指數(shù)與同期中國降水的相關,結果如圖7-圖9所示:6~8月區(qū)域急流強度圍繞多年平均值呈先減小后上升的趨勢,7月急流強度最弱,8月次之,6月最強,基本趨勢符合急流變化規(guī)律。指數(shù)與長江中下游地區(qū)降水及華北地區(qū)降水有顯著的相關,均通過0.01信度檢驗。對降水的指示意義較好。從參數(shù)B1的物理意義上看,區(qū)域平均西風風速的變化,與所選取的區(qū)域關系密切,與區(qū)域內(nèi)是否存在大于30m·s-1的急流中心范圍有密切聯(lián)系。例如2005年4月13日東亞副熱帶急流中心區(qū)有一半面積落在統(tǒng)計區(qū)域(30°N~50°N,90°E~110°E)之外,此種情況下用區(qū)域西風平均并不能真正反映東亞副熱帶急流的強度。不同范圍的選擇結果有很大的不同。(如何選擇范圍)關鍵區(qū)的選擇與所研究的區(qū)域有很大的關系,選擇不當,則不能準確捕捉急流的主要特征,因此在實際應用中,要針對研究的區(qū)域準確的選擇急流強度關鍵區(qū),方能取得滿意的結果。
圖7 2005年與相應多年平均6~8月急流強度參數(shù)逐日演變
圖8 演變與同期中國夏季降水相關
譚杰麗等[3]選取200 hPa上30m·s-1風速作為副熱帶急流特征風速值,大于30m·s-1特征風速值區(qū)域為急流中心范圍,針對區(qū)域(20°N ~ 60°N,90°E ~ 150°W)定義了 2 種急流中心強度,一種是急流中心絕對強度,即急流中心范圍內(nèi)所有格點的風速值與特征風速值的差值之和;另一種是急流相對強度,即急流中心范圍內(nèi)所有格點的風速值與特征風速值的差值之和與急流面積之比。急流面積為急流中心范圍內(nèi)大于特征風速值的格點總數(shù)。分析指出,東亞副熱帶急流的中心強度、絕對強度和相對強度的變化具有一致性,體現(xiàn)了急流強度風速變化的整體性和一致性。急流絕對強度參數(shù)用B2表示,急流相對強度參數(shù)用B3表示。
以2005年6~8月為例,對急流強度參數(shù)B2和B3進行逐日數(shù)值反演,并對結果進行標準化處理,結果如圖10~11所示。B2、B3指數(shù)呈先快速減小后緩慢增加的趨勢,與多年同期相比變化幅度不同,B2指數(shù)強度下降趨勢比B3要緩慢得多。這是2種方法的物理含義不同造成的,反映了2種方法的不同特性。從與降水的相關分布上來看,急流絕對強度參數(shù)B2和急流相對強度參數(shù)B3均能反映實際狀況,具有客觀定量、計算簡便、物理意義清晰、表現(xiàn)直觀等特點。絕對強度參數(shù)表征急流的整體絕對強度,與急流面積有關;相對強度參數(shù)表征急流的平均狀況,消除了面積大小的影響,兩者僅是物理含義不同,性能上不具有可比性。
圖9 2005年4月13日200hPa緯向風等值線(單位:m·s-1)
圖10 2005年與相應多年平均6~8月急流強度參數(shù)逐日演變
圖11 各參數(shù)與中國夏季降水相關
表2給出了上述3種方法的主要特征、分析尺度及適用范圍,由表可見,表征急流強度的3種方法均適應于逐日實時監(jiān)測與環(huán)流觀測,區(qū)域風速法適用于研究特定區(qū)域內(nèi)急流強度,關鍵在于選取急流活動影響的關鍵區(qū)。絕對強度參數(shù)表征急流的整體絕對強度,與急流面積有關。相對強度參數(shù)表征急流的平均狀況,消除了面積大小的影響,兩者僅是物理含義不同,性能上不具有可比性。參數(shù)物理意義不同,主要特征不同,研究者應根據(jù)研究需要,恰當?shù)倪M行選擇。
表2 東亞副熱帶西風急流強度參數(shù)表征方法比較
如何客觀定量表征東亞副熱帶急流是目前科研和業(yè)務領域迫切需要解決的問題,論文總結歸納了以往的東亞副熱帶急流客觀定量表征方法,基于對各種方法的性能分析結果,結合東亞副熱帶急流特征及其對中國夏季降水的影響,研究了東亞副熱帶西風急流軸位置、急流強度特征參數(shù)的客觀定量表征方法。利用NCEP/NCAR再分析資料,以2005年6~8月為例,實現(xiàn)各種方法的數(shù)值反演和性能分析,得出以下結論和有待進一步討論的問題:
(1)定義位置指數(shù)的方法主要有軸線法、EOF分解法及區(qū)域風速差法。軸線法定義合理且有明確的物理意義,方法簡單易行、意義清晰、表現(xiàn)直觀,適用于逐日實時監(jiān)測和診斷分析。EOF分解法與區(qū)域風速差法適用于氣候評估和季節(jié)診斷;表征急流強度的方法主要有區(qū)域風速法及中心強度法,均適應于逐日實時監(jiān)測與環(huán)流觀測,對華北地區(qū)及長江中下游雨量多寡有較好的指示意義,方法計算簡便、意義清晰、指示性好。每一種指數(shù)都有側(cè)重點,定義的方法有各自的優(yōu)點和局限性,在實際應用中,依據(jù)研究的目的和分析的角度,采用適當?shù)姆椒?才能取得比較滿意的效果。
(2)由于影響氣候短期變化的因子眾多,且期間又有復雜的相互作用,而目前尚未對其中任何一個因子的詳細物理過程及其影響有明晰的認識,因此僅以某種標準或公式的客觀方法實現(xiàn)對變化多樣的東亞副熱帶西風急流的表征是極其困難的,在實際的預報過程中,要綜合考慮各種氣象條件的因素,根據(jù)實際情況合理恰當?shù)倪x擇副熱帶西風急流指數(shù),以做出較為準確的預報。
下一步將利用在上述指數(shù)對東亞副熱帶急流與中國東部主要雨季的關系作統(tǒng)計分析,并研究這種對應關系通過何種機理影響中國東部夏季降水,希望研究成果能夠深化對東亞副熱帶急流的認識,豐富中國東部夏季降水的中期預報提供理論基礎。數(shù)值試驗基于2005年6~8月時段,由于不同年份間的副熱帶西風活動存在年際差異,1年時段試驗對比結果是否有普遍性尚待驗證。
[1] 丁一匯,趙深銘,傅秀琴.5~10月全球熱帶和副熱帶200hPa多年平均環(huán)流的研究(二)—行星風系[J].大氣科學,1988,12(3):242-249.
[2] Nakamura H,T Sampe.Trapping of synoptic scales disturbanve into the North-Paciffic subtropical jet core in midwinter[J].Geophys Res Lett,2002,29:1761-1776.
[3] 譚杰麗,江靜,袁俊鵬.副熱帶高空急流各中心強度時間變化分析[J].氣象科學,2009,29(4):482-489.
[4] 葉篤正,陶詩言,李麥村.在6月和10月大氣環(huán)流的突變現(xiàn)象[J].氣象學報,1958,29(4):250-263.
[5] 陶詩言,趙煜佳,陳曉敏.東亞的梅雨與亞洲上空大氣環(huán)流季節(jié)變化的關系[J].氣象學報,1958,29(2):119-134.
[6] 李崇銀,王作臺,林士哲,等.東亞夏季風活動與東亞高空西風急流位置北跳關系的研究[J].大氣科學,2004,28(5):641-658.
[7] 張耀存,況雪源.東亞副熱帶西風急流季節(jié)變化突變特征及其與東亞夏季風建立的聯(lián)系[J].第六次全國動力氣象學術會議,2005,18.
[8] 金榮花,矯梅燕,徐晶,等.2003年淮河多雨期西太平洋副高活動特征及其成因分析[J].熱帶氣象學報,2006,22(1):60-66.
[9] 金榮花,李艷,王式功.四種客觀定量表征阻塞高壓方法的對比分析[J].高原氣象,2009,28(5):1121-1128.
[10] 況雪源,張耀存.東亞副熱帶西風急流位置異常對長江中下游夏季降水的影響[J].高原氣象,2006,25(3):382-389.
[11] 黃安麗,高坤.對流層高低空急流耦合作用的動力學分析[J].杭州大學學報,1982,9(3):256-364.
[12] 劉飛,何金海,董國青.東亞夏季西風急流變化與中國夏季降水的相關分析[M].中國氣象學會2006年年會,138-152.
[13] 況雪源,張耀存,劉健.對流層上層副熱帶西風急流與東亞冬季風的關系[J].高原氣象,2008,27(4):701-712.
[14] 廖清海,高守亭,王會軍,等.北半球夏季副熱帶西風急流變異及其對東亞夏季風氣候異常的影響[J].地球物理學報,2004,47(1):10-18.
[15] 楊蓮梅,張慶云.夏季東亞西風急流Rossby波擾動異常與中國降水[J].大氣科學,2007,31(4):586-595.
[16] Lin Z D,Lu R Y.International meridional displacement of the East Asian upper-tropospheric jet stream in summer[J].ADV Atmos Sci,2005.22:199-211.
[17] Liang X Z,Wang W C.Asscociation between China monsoon rainfall and tropospheric jets[J].Q.J.R.Meteo.Soc.,1998,124(6):2597-2623.
[18] 史有瑜,張耀存.東亞副熱帶西風急流變化與華北夏季降水之間的關系[M].中國氣象學會2008年年會,625-634.
[19] Yang S,P J Webster.The effect of summer tropical heating on the location and intensity of the extratropical westerly jet streams[J].J Geophys Res,1990,95(D11):18705-18721.
[20] 毛睿,龔道溢,房巧敏.冬季東亞中緯度西風急流對我國氣候的影響[J].應用氣象學報,2007,18(2):137-146.
[21] 孫建華,衛(wèi)捷,趙思雄,等.2005年夏季的主要天氣及其環(huán)流分析[J].氣候與環(huán)境研究,2006,11(2):138-154.