程胡華 王益柏 武 帥 趙 亮 智茂林

1.63729 部隊，太原，030027

2.61741 部隊，北京，100094

3.32021 部隊，北京，100094

4.中國科學院大氣物理研究所大氣科學和地球流體力學數(shù)值模擬國家重點實驗室（LASG），北京，100029

1 引 言

重力波對全球大氣的動力、熱力結(jié)構(gòu)變化及調(diào)整具有重大影響（Holton，1983；Fritts，et al，2003），與多種不同尺度天氣現(xiàn)象均有密切聯(lián)系。中外研究均表明，大氣重力波影響著臺風（洪軍等，2015；陳丹等，2011；Chen，et al，2013）、雨雪（鄧少格等，2012；程胡華，2016；孫艷輝等，2015）、冰雹（李啟泰等，1993；覃衛(wèi)堅等，2010）等其他天氣（Thomas，et al，1999；林永輝等，2003；陳丹等，2014；Plougonven，et al，2014；Kramer，et al，2015；李馳欽等，2018）的產(chǎn)生、發(fā)展和消亡；因此，為提高大氣模式對不同尺度天氣現(xiàn)象的預報能力，需要在大氣模式中使用重力波參數(shù)化的方法加入重力波影響（徐國強等，2010；王宇虹等，2017），利用高垂直分辨率探空資料獲取全球范圍內(nèi)大氣重力波參數(shù)氣候特征及其分布，已成為研制大氣模式中重力波參數(shù)化方案的重要前提條件之一（Yamamori，et al，2006；Zhang，et al，2010）。

20世紀80年代，科學家們使用無線電探空資料中的溫度廓線研究重力波能量及譜斜率（van Zandt，1982），隨后，Kitamura等（1989）、Hamilton（1991）、Allen等（1995）利用風場、溫度場廓線計算重力波參數(shù)，經(jīng)過不斷的完善（Eckermann，1996；Eckermann，et al，1999），形成了依據(jù)風場和溫度場廓線計算重力波參數(shù)的一套完整方法，該方法已被認可并廣泛應用（Guest，et al，2000；卞建春等，2004；Zhang，et al，2005，2007；Ki，et al，2010；白志宣 等， 2016；Kafando， et al， 2016； Marlton，et al，2016；Jackson，et al，2018；Chen，et al，2019；吳泓錕等，2019）。在利用此方法獲取大氣重力波參數(shù)時，首先需獲取溫度場、緯向風和經(jīng)向風的垂直擾動場，目前，不同階曲線擬合法是獲取擾動場的常用方法（Guest，et al，2000；卞建春等，2004；Zhang，et al，2005，2007；Ki，et al，2010；白志宣等，2016；Kafando，et al，2016；Marlton，et al，2016；Jackson，et al，2018；Chen，et al，2019；吳泓錕等，2019），例如，白志宣等（2016）利用2階曲線擬合法得到中國南北向8個站兩年多高垂直分辨率無線電探空資料下平流層擾動場，基于該擾動場獲取下平流層慣性重力波，并對該重力波活動的季節(jié)變化及其隨緯度分布特征進行研究；Ki等（2010）基于韓國10個站高垂直分辨率探空資料，通過3階曲線擬合擾動場獲取下平流層大氣重力波參數(shù)，并對重力波源及其特征進行研究；Guest等（2000）利用4階曲線擬合擾動場獲取麥夸里島下平流層重力波參數(shù)，并對重力波參數(shù)的季節(jié)變化及波源特征進行了研究。

目前，不同階曲線擬合擾動場獲取的重力波參數(shù)氣候特征是否與實況一致，尚未得到大量觀測資料或?qū)嶒灧椒ǖ挠行炞C，且針對不同階曲線擬合擾動場之間的差異特征并分析依此計算得到大氣重力波參數(shù)間氣候差異特征的研究較少，程胡華等（2013，2016）利用Morlet小波對不同階曲線擬合擾動場進行分析發(fā)現(xiàn)，大多數(shù)情況下擾動場中包含背景場信號，其含有的背景場信號強弱對大氣重力波參數(shù)有較大影響；程胡華（2017）以2013年8月6—8日晉西北地區(qū)出現(xiàn)的一次雷陣雨天氣過程為例，利用不同方法獲取高垂直分辨率探空資料擾動場，發(fā)現(xiàn)所得擾動場垂直廓線之間存在一定差異，該差異對獲取大氣重力波參數(shù)的大小及變化趨勢存在影響，但對不同重力波參數(shù)的影響程度不同；文中在已有研究（程胡華等，2013，2016；程胡華，2017）的基礎上，利用2014—2017年太原地區(qū)高垂直分辨率探空資料，進一步針對不同階曲線擬合擾動場對下平流層（17—24 km高度）大氣重力波參數(shù)氣候特征的影響進行統(tǒng)計研究。

2 資料及計算方法

所用研究資料為山西太原氣象臺（37.78°N，112.55°E）提供的2014—2017年高垂直分辨率探空資料，每日08時（北京時）的溫度、濕度和氣壓的垂直分辨率約為10 m，風速的垂直分辨率約為300 m，由于原始資料的垂直間隔不等，在此利用3次樣條插值方法得到垂直分辨率為50 m的氣象要素值。

對于高垂直分辨率探空資料中的緯向風速（u）、經(jīng)向風速（v）、溫度（T），由于對流層和平流層大氣間的浮力頻率差異較大，并且為了避開對流層頂?shù)臉O低溫度和對流層上部緯向風急流的影響，選擇下平流層高度范圍為海拔17 —24 km（馬蘭夢等，2012）；分別采用2—4階曲線擬合法得到下平流層背景場(uˉ、vˉ 、Tˉ )，用原始數(shù)據(jù)(u、v、T)減去背景場（uˉ、vˉ 、Tˉ ），即得到擾動場（u′、v′、T′）。獲取重力波參數(shù)所用的程序包為澳大利亞阿德萊德大學R.Vincent教授提供，線性重力波理論（卞建春等，2004；鄧少格等，2012；Sawyer，1961）和多種觀測資料均表明，由重力波產(chǎn)生的水平風速矢量(u′(z)，v′(z))端點的連線，大致表現(xiàn)為橢圓形狀，若重力波的水平傳播方向為自北順時針旋轉(zhuǎn)角度（α），則 α方向即為水平風擾動矢量偏振橢圓的長軸方向（卞建春等，2004；鄧少格等，2012；Hines，1989）。通過利用 Stokes參數(shù)（卞建春等，2004；鄧少格等，2012；Eckermann，1996）和主成分分析法（卞建春等，2004；鄧少格等，2012），借助Hilbert變換后的歸一化溫度擾動和速度擾動分量的關系（卞建春等，2004；鄧少格等，2012），可確定出重力波水平傳播方向（φ）、橢圓偏振度（d）和偏振橢圓的長短軸，由偏振橢圓長短軸之比（卞建春等，2004；鄧少格等，2012），即可求出波固有頻率（）和固有周期通過歸一化溫度擾動（）的垂直波數(shù)功率譜可以估算出波動的垂直波數(shù)（m）（Hoffmann，et al， 2006； Moffat-Griffin， et al，2011），依據(jù)波數(shù)與波長的關系，可求取垂直波長（λz， λz=2π/m），將上述得到的和m代 入由線性重力波理論得到的頻散關系式，可得到重力波的水平波長（λh）、固有相速（ci）、對地水平群速（cg）等一系列的波參數(shù)。

基于最新研究成果（白志宣等，2016；Jackson，et al，2018；Chen，et al，2019；吳泓錕等，2019）采用2階曲線擬合擾動場，程胡華（2017）研究發(fā)現(xiàn)，利用2階曲線擬合擾動場獲取大氣重力波參數(shù)的演變特征，與2013年8月6—8日的大氣動力演變特征基本吻合，故以2階曲線擬合擾動場及計算得到的重力波參數(shù)為基準，統(tǒng)計其與其他階曲線擬合擾動場的相關系數(shù)，并對相關系數(shù)在不同區(qū)間內(nèi)的占有率進行統(tǒng)計，隨后對計算得到的重力波參數(shù)間的相關、偏差特征及占有率進行統(tǒng)計分析，其中，占有率為物理量數(shù)值在某區(qū)間內(nèi)的個數(shù)占總數(shù)的百分比。

3 擾動場特征對比

記Rank2、Rank3、Rank4分別代表2階、3階、4階曲線擬合的擾動場，Upert、Vpert、Tpert分別為緯向風擾動場、經(jīng)向風擾動場、溫度擾動場。表1列出了2—4階曲線的擬合優(yōu)度（Goodness of Fit）參數(shù)值統(tǒng)計結(jié)果，其中，參數(shù) SSE （The sum of squares due to error） 、 RMSE （Root mean squared error）、R2（Coefficient of determination）、Radj（Degree of freedom adjusted coefficient of determination）分別表示誤差平方和、均方根誤差、決定系數(shù)、自由度調(diào)整決定系數(shù)，SSE和RMSE值越小、R2和Radj值越接近1，表明曲線擬合優(yōu)度越高，從表1可以看出，4階曲線的擬合優(yōu)度最高，其次為3階曲線擬合方法，而2階曲線擬合方法的擬合優(yōu)度較差。

表1 2—4 階曲線的擬合優(yōu)度參數(shù)值統(tǒng)計結(jié)果Table 1 Statistical results of goodness of fit parameter values for 2—4 degree curve fitting

從2014—2017年太原地區(qū)下平流層（17—24 km）2—4 階曲線擬合得到的Upert、Vpert、Tpert隨高度變化特征（圖1）可見，平均緯向風擾動場隨高度變化趨勢基本一致，但振幅存在一定差異（圖1a），在整個高度范圍內(nèi)，2—4階平均緯向風擾動場范圍分別為 ?0.19 （24.00 km） — 0.12 m/s （17.00 km） 、?0.08（20.95 km）—0.08 m/s （19.55 km）和?0.07（18.45 km）—0.08 m/s （19.60 km），2 階與 3 階平均經(jīng)向風擾動場隨高度變化基本重合，但在高度17—18、23—24 km均與4階擾動場的變化趨勢相反（圖1b），在高度 17—24 km，2—4 階平均經(jīng)向風擾動場范圍分別為?0.10 （22.45 km）—0.12 m/s（20.95 km）、?0.10 （22.45 km）—0.12 m/s （20.95 km）和?0.08 （17.00 km）—0.06 m/s （20.95 km），2—4 階平均溫度擾動場之間的振幅及隨高度變化特征均存在一定差異（圖1c），其中，2階平均溫度擾動場振幅明顯偏大，在高度17—24 km，2—4階平均溫度擾動場范圍分別為?0.48 （18.95 km）—0.79℃（17.00 km）、?0.17（22.45 km）—0.24℃ （24.00 km）和?0.20 （17.00 km）—0.11℃ （20.45 km）。

為反映不同階曲線擬合擾動場間的相關，通過計算每個時次的3階、4階曲線擬合與2階曲線擬合擾動場垂直廓線的相關系數(shù)，得到2014—2017年樣本數(shù)的數(shù)組；圖2給出3階、4階曲線擬合與2階曲線擬合的緯向風擾動場的相關系數(shù)、經(jīng)向風擾動場相關系數(shù)、溫度擾動場相關系數(shù)在不同區(qū)間內(nèi)占有率的分布特征，除相關系數(shù)在[0.0，0.1]內(nèi)占有率均為0外，在其他區(qū)間占有率均大于0，其中，3階曲線擬合擾動場相關系數(shù)值沿橫坐標均呈現(xiàn)出增大特征，而4階曲線擬合擾動場相關系數(shù)值沿橫坐標呈現(xiàn)出增大、基本維持不變的趨勢特征。3階、4階曲線擬合緯向風擾動場相關系數(shù)分別在（0.9，1.0 ]、（0.7，0.8 ]內(nèi)占有率最大（圖2a），分別為45.46%和16.54%，相關系數(shù)超過0.8的占有率分別為62.54 %和27.38%，而相關系數(shù)不超過0.5的占有率分別為5.31%和22.31%，類似緯向風擾動場相關系數(shù)分布特征，3階、4階曲線擬合經(jīng)向風擾動場相關系數(shù)分別在（0.9，1.0]、（0.8，0.9]內(nèi)占有率最大（圖2b），分別為45.46%和17.77%，相關系數(shù)超過0.8的占有率分別為65.38%和30.23%，相關系數(shù)不超過0.5的占有率分別為4.31%和18.69%，3階、4階曲線擬合溫度擾動場相關系數(shù)均在（0.9，1.0]內(nèi)的占有率最大（圖2c），分別為47.62%和20.46%，而相關系數(shù)超過0.8的占有率分別為67.54%和40.08%，相關系數(shù)不超過0.5的占有率分別為2.31%和9.46%。

從上述分析可知，與4階曲線擬合擾動場相比，2階與3階曲線擬合擾動場的相似度更高；因此，下面針對2階、3階曲線擬合擾動場獲取的下平流層大氣重力波參數(shù)結(jié)果進行對比研究。

4 重力波參數(shù)特征對比

4.1 群 速

2階、3階曲線擬合擾動場的群速值范圍、年內(nèi)變化趨勢存在一定差異（圖3a），2—5月，2階曲線擬合擾動場的群速表現(xiàn)為增大、減小、增大的變化特征，而 3階曲線擬合擾動場的群速表現(xiàn)為減少、增大的變化特征；7—9月，2階曲線擬合擾動場的群速表現(xiàn)為減小、增大的變化特征，而3階曲線擬合擾動場的群速為一致減??；1—12月，2階曲線擬合 擾 動 場 的 群 速 值 為 ?1.52（12月 ） — 8.97 m/s（7月），平均值為 3.51 m/s，1月、11—12月為負值；而3階曲線擬合擾動場的群速為?3.64（12月）—7.30 m/s（7 月），平均為 2.10 m/s，1 月、10—12 月為負值；2階與3階曲線擬合擾動場之間的群速相關系數(shù)年內(nèi)變化特征（圖3b）可見，除4月超過0.8外，其他月份的群速相關系數(shù)均小于0.8，1—12月，群速相關系數(shù)的值為 0.57（6月）—0.81（4月），平均為0.71。

2階、3階曲線擬合擾動場的群速在不同區(qū)間內(nèi)占有率分布特征（圖4）可見，不同區(qū)間占有率沿橫坐標均呈現(xiàn)出增大、減小的變化，其最大值位于[0，10）（圖4a），分別為 40.05%、43.35%，而在（?∞，?20]、[30，+∞）內(nèi)的占有率均較小，合計分別為1.65%、1.65%；3階與2階曲線擬合擾動場之間的群速偏差（以2階曲線擬合擾動場的群速為基準，下同）在不同區(qū)間內(nèi)占有率的分布特征（圖4b）表明，其占有率最大值出現(xiàn)在[?10，0），數(shù)值為56.88%，主要出現(xiàn)在 [?0，10），占有率為90.87%，而在其他區(qū)間內(nèi)的占有率不超過10%，合計為9.13%。

4.2 水平波長

2階、3階曲線擬合擾動場的水平波長及水平波長相關系數(shù)的年內(nèi)變化特征（圖5）表明，1—12月，2階曲線擬合擾動場的水平波長均明顯大于對應的3階曲線擬合擾動場（圖5a），2階、3階曲線擬合擾動場的水平波長最大值均出現(xiàn)在8月、最小值均出現(xiàn)在12月；1—3月，2階曲線擬合擾動場的水平波長呈現(xiàn)出增大、減小的變化特征，而3階曲線擬合擾動場的水平波長為一致減??；6—8月，2階曲線擬合擾動場的水平波長表現(xiàn)為減小、增大的變化特征，而3階曲線擬合擾動場的水平波長為一致增大；1—12月，2階曲線擬合擾動場的水平波長范圍為 279.37 （12 月）—444.55 km （8 月），平均為346.36 km，而3階曲線擬合擾動場的水平波長范圍為 217.07（12月）—306.93 km（8月），平均為261.83 km；2階與3階曲線擬合擾動場之間的水平波長相關系數(shù)年內(nèi)變化（圖5b）表明，1—12月，水平波長相關系數(shù)范圍為0.28（6月）—0.74（11月），平均值為0.55。

圖2 2014—2017年太原地區(qū)下平流層 3階、4階曲線擬合與2階曲線擬合的緯向風擾動場相關系數(shù) （a）、經(jīng)向風擾動場相關系數(shù) （b） 和溫度擾動場相關系數(shù) （c） 的分布Fig.2 Distribution characteristics of zonal wind disturbance field correlation coefficient （a），meridional wind disturbance field correlation coefficient （b） and temperature disturbance field correlation coefficient （c） between 3rd，4th order curve fitting methods and 2nd order curve fitting method in the lower stratosphere above Taiyuan during 2014—2017

圖3 2014—2017 年太原地區(qū)下平流層 3 階與 2 階曲線擬合擾動場的群速值 （a） 和群速相關系數(shù) （b） 的月變化Fig.3 Annual variations of group velocity from disturbance fields obtained by second and third order curve fitting （a） and group velocity correlation coefficient （b） in the lower stratosphere above Taiyuan during 2014—2017

圖4 2014—2017年太原地區(qū)下平流層3階與2階曲線擬合擾動場的群速 （a）和群速偏差（b）在不同區(qū)間占有率的分布Fig.4 Occupancy distribution characteristics of group velocity from disturbance fields obtained by second and third order curve fitting （a） and group velocity deviation （b） in different interval ranges in the lower stratosphere above Taiyuan during 2014—2017

2階、3階曲線擬合擾動場的水平波長在不同區(qū)間內(nèi)占有率沿橫坐標方向呈現(xiàn)出減小趨勢，但最大值均未超過 50 %（圖6a），且均在（0，200）占有率最高，分別為 39.50 % 和 49.72 %，其次為 [200，400），分別為 28.56%和 28.64%，而在 [600，800）、 [800，+∞）內(nèi)的占有率均小于10%。3階與2階曲線擬合擾動場間的水平波長偏差在不同區(qū)間內(nèi)的占有率分布特征（圖6b）表明，其在 [?200，0）內(nèi)占有率最高（42.64%），主要出現(xiàn)在[?200，200]，占有率為73.09%，而在 [200，+∞）內(nèi)占有率為 6.06%，在（?∞，?200）內(nèi)占有率為20.85%。

4.3 垂直波長

2階、3階曲線擬合擾動場的垂直波長及垂直波長相關系數(shù)的年內(nèi)變化特征（圖7）表明，1—12月，2階曲線擬合擾動場的垂直波長均大于對應的3階曲線擬合擾動場（圖7a），2—4月，2階曲線擬合擾動場的垂直波長為一致減小，而3階曲線擬合擾動場的垂直波長表現(xiàn)為減小、增大變化特征，5—7月，2階曲線擬合擾動場的垂直波長為一致增大，而3階曲線擬合擾動場的垂直波長表現(xiàn)為減小、增大變化特征；1—12月，2階曲線擬合擾動場的垂直波長范圍為 2.64（5月）—2.94 km（12月），平均為2.76 km，3階曲線擬合擾動場的垂直波長范圍為2.20（6月）—2.55 km（12月），平均為 2.38 km；2階與3階曲線擬合擾動場垂直波長的相關系數(shù)年內(nèi)變化特征（圖7b）表明，1—12月，垂直波長相關系數(shù)范圍為 0.56（6月）—0.86（1月），平均為 0.72。

圖5 2014—2017 年太原地區(qū)下平流層 3 階與 2 階曲線擬合擾動場的水平波長 （a）和水平波長相關系數(shù) （b） 的月變化Fig.5 Annual variations of horizontal wavelength from disturbance fields obtained by second and third order curve fitting （a） and horizontal wavelength correlation coefficient （b） in the lower stratosphere above Taiyuan during 2014—2017

圖6 2014—2017 年太原地區(qū)下平流層 3 階與 2 階曲線擬合擾動場的水平波長 （a） 和水平波長偏差 （b） 在不同區(qū)間占有率的分布Fig.6 Occupancy distribution characteristics of horizontal wavelength from disturbance fields obtained by second and third order curve fitting （a） and horizontal wavelength deviation （b） in different interval ranges in the lower stratosphere above Taiyuan during 2014—2017

2階與3階曲線擬合擾動場的垂直波長及其之間的垂直波長偏差，在不同區(qū)間內(nèi)占有率的分布特征（圖8）表明，其占有率沿橫坐標方向呈現(xiàn)出增大、減小的變化特征（圖8a），最大值均出現(xiàn)在[2，3），分別為50.67%和60.19%，其次分別出現(xiàn)在 [3，4）和[1，2），分別為 32.02% 和 26.51%，而在（0，1）內(nèi)均為0，且3階曲線擬合擾動場的垂直波長在 [4，+∞）占有率為0，而相應的2階曲線擬合擾動場垂直波長為4.09%；垂直波長偏差在不同區(qū)間內(nèi)的占有率分布特征如圖8b所示，其占有率沿橫坐標軸方向呈現(xiàn)出增大、減小的變化特征，且在 [?1，0）內(nèi)占有率最大（70.97%），在 [0，1）占有率為次大（18.57%），而在 [1，+∞）、（?∞，?3）占有率均為 0.00%。

4.4 周期

2階、3階曲線擬合擾動場的周期及周期相關系數(shù)的年內(nèi)變化（圖9）表明，1—12月，2階曲線擬合擾動場的周期值均大于對應的3階曲線擬合擾動場，且兩者年內(nèi)變化趨勢基本一致，均呈現(xiàn)出增大、減小的趨勢（圖9a），1—12月，2階曲線擬合擾動場的周期為6.20（12 月）—8.80 h（8 月），平均為7.4012 h，3階曲線擬合擾動場的周期值范圍為5.68（12月）—7.92 h（6 月），平均為 6.86 h；2 階與 3 階曲線擬合擾動場的周期相關系數(shù)1—12月的變化范圍為 0.38（6月）—0.75（1月），平均為 0.63（圖9b）。

圖7 2014—2017 年太原地區(qū)下平流層 3 階與 2 階曲線擬合擾動場的垂直波長 （a） 和垂直波長相關系數(shù) （b） 的月變化Fig.7 Annual variations of vertical wavelength from disturbance fields obtained by second and third order curve fitting （a） and vertical wavelength correlation coefficient （b） in the lower stratosphere above Taiyuan during 2014—2017

圖8 2014—2017 年太原地區(qū)下平流層 3 階與 2 階曲線擬合擾動場的垂直波長 （a） 和垂直波長偏差 （b） 在不同區(qū)間的占有率Fig.8 Occupancy distribution characteristics of vertical wavelength from disturbance fields obtained by second and third order curve fitting （a） and vertical wavelength deviation （b） in different interval ranges in the lower stratosphere above Taiyuan during 2014—2017

2階、3階曲線擬合擾動場的周期在不同區(qū)間內(nèi)的占有率沿橫坐標軸方向的變化較?。▓D10a），其最大值分別位于 [6，9）和 [3，6），分別為 22.19%和23.84%，次大值均位于（0，3），分別為21.09%和23.45%，而最小值分別位于 [9，12）和 [12，+∞），分別為16.84%和14.48%；3階與2階曲線擬合擾動場之間的周期偏差，在不同區(qū)間內(nèi)的占有率（圖10b）沿橫坐標軸方向呈現(xiàn)出增大、減小變化特征，最大值位于 [0，3）（34.30%），次大值占有率位于 [?3，0）（32.97%），而在（?∞，?12 ]和 [12，+∞）的占有率均較小，合計為1.02%。

4.5 固有相速

2階、3階曲線擬合擾動場的固有相速及固有相速相關系數(shù)的年內(nèi)變化（圖11）表明，1—12月，2階與3階曲線擬合擾動場的固有相速年內(nèi)變化趨勢存在一定差異，且2階曲線擬合擾動場的固有相速值均大于對應的3階曲線擬合擾動場（圖11a）。1—3月，2階曲線擬合擾動場的固有相速表現(xiàn)為增大、減小變化特征，而3階曲線擬合擾動場的固有相速為一致減??；9—12月，2階曲線擬合擾動場的固有相速表現(xiàn)為減小、增大變化特征，3階曲線擬合擾動場的固有相速為一致增大變化特征；1—12月，2階曲線擬合擾動場的固有相速為10.72（1 月）—12.45 m/s（8 月），平均為 11.67 m/s，3 階曲線擬合擾動場的固有相速為9.00（5月）—10.05 m/s（8月），平均為 9.65 m/s；2階與 3階曲線擬合擾動場之間的固有相速相關系數(shù)（圖11b）1—12月變化為 0.38（6月）—0.75（1月），平均為 0.58。

圖9 2014—2017 年太原地區(qū)下平流層 3 階與 2 階曲線擬合擾動場的周期 （a） 和周期相關系數(shù) （b） 的月變化Fig.9 Annual variations of period from disturbance fields obtained by second and third order curve fitting （a） and period correlation coefficient （b） in the lower stratosphere above Taiyuan during 2014—2017

圖10 2014—2017 年太原地區(qū)下平流層 3 階與 2 階曲線擬合擾動場的周期 （a） 和周期偏差 （b） 在不同區(qū)間占有率的分布Fig.10 Occupancy distribution characteristics of period from disturbance fields obtained by second and third order curve fitting （a） and period deviation （b） in different interval ranges in the lower stratosphere above Taiyuan during 2014—2017

2階與3階曲線擬合擾動場的固有相速在不同區(qū)間內(nèi)的占有率沿橫坐標軸方向的變化趨勢特征存在差異，2階曲線擬合擾動場的固有相速呈現(xiàn)一致增大趨勢，而3階曲線擬合擾動場的固有相速表現(xiàn)為增大、減小變化特征（圖12a）。2階、3階曲線擬合擾動場的固有相速占有率最大值分別位于[12，+∞）和 [8，10），分別為 38.71%和30.45%；次大值分別位于 [10，12）和 [6，8），分別為 23.76%和24.86%；2階、3階曲線擬合擾動場的固有相速在（0，2）的占有率均為0.00%，2階曲線擬合擾動場的固有相速在[2，4）的占有率為0，而相應的3階曲線擬合擾動場的固有相速占有率為0.08%。固有相速偏差在不同區(qū)間內(nèi)的占有率（圖12b）沿橫坐標軸方向呈現(xiàn)出增大、減小變化特征，最大值位于 [?2，0）（41.31%），次大值位于 [?4，?2）（20.61%），固有相速偏差在 [2，+∞）的占有率較小，合計為2.99%。

4.6 能量上傳率

2階與3階曲線擬合擾動場的能量上傳率（FWUPE，fraction of wave upward propagation energy）的年內(nèi)變化趨勢基本一致，且兩者的數(shù)值比較接近（圖13a），1—12月，2階曲線擬合擾動場的能量上傳率為 63.79（12月）—75.92%（8月），平均為70.01%；3階曲線擬合擾動場的能量上傳率為65.37（12月）—74.80%（8月），平均為 70.44%。2階與3階曲線擬合擾動場的能量上傳率相關系數(shù)（圖13b）1—12月為 0.78（6月）—0.89（5月），平均為0.85。

圖11 2014—2017 年太原地區(qū)下平流層 3 階與 2 階曲線擬合擾動場的固有相速 （a） 和固有相速相關系數(shù) （b） 的月變化Fig.11 Annual variations of intrinsic phase velocity from disturbance fields obtained by second and third order curve fitting （a） and intrinsic phase velocity correlation coefficient （b） in the lower stratosphere above Taiyuan during 2014—2017

圖12 2014—2017 年太原地區(qū)下平流層 3 階與 2 階曲線擬合擾動場的固有相速 （a） 和固有相速偏差 （b） 在不同區(qū)間占有率的分布Fig.12 Occupancy distribution characteristics of intrinsic phase velocity from disturbance fields obtained by second and third order curve fitting （a） and intrinsic phase velocity deviation （b） in different interval ranges in the lower stratosphereabove Taiyuan during 2014—2017

2階與3階曲線擬合擾動場的能量上傳率及其之間的能量上傳率偏差在不同區(qū)間內(nèi)的占有率分布特征（圖14）表明，2階與3階曲線擬合擾動場之間的能量上傳率的占有率在各區(qū)間內(nèi)基本相同，且均沿橫坐標軸方向呈現(xiàn)出增大特征（圖14a），最大值均位于 [80，100]，分別為39.02%和39.81%，最小值均位于 [0，20），分別為1.65%和1.57%；能量上傳率偏差在不同區(qū)間內(nèi)的占有率沿橫坐標軸方向呈現(xiàn)出增大、減小變化（圖14b），最大值位于 [0，10）（36.66%），次大值位于 [?10，0）（32.02%），能量上傳率偏差在（?∞，?50）的占有率為 0.00%，在 [?50，?40）和 [40，+∞）的占有率分別為0.08%和0.24%。

5 結(jié)論和討論

大氣重力波在大氣中普遍存在，其對全球大氣的動力、熱力等變化具有重要影響，且與多種不同尺度的天氣現(xiàn)象密切相關，為提高大氣模式對多種不同尺度天氣現(xiàn)象的預報能力，需通過重力波參數(shù)化方案引入重力波的影響；目前，確定重力波參數(shù)氣候特征是研制大氣模式中重力波參數(shù)化的一個重要前提條件，這些參數(shù)的選取需要根據(jù)觀測資料來確定。中外常采用2—4階曲線擬合方法擾動場來計算得到大氣重力波參數(shù)，并以此結(jié)果分析大氣重力波的氣候統(tǒng)計特征，但較少針對不同階曲線擬合方法擾動場對重力波參數(shù)氣候特征的影響研究；本文對2—4階曲線擬合方法擾動場進行對比，統(tǒng)計發(fā)現(xiàn)2階與3階曲線擬合方法擾動場的相似程度更高；因此，選取相似程度更高的2階、3階曲線擬合方法擾動場分別獲取大氣重力波參數(shù)，并對大氣重力波參數(shù)之間的氣候差異特征進行分析，得到如下初步結(jié)論：

圖13 2014—2017 年太原地區(qū)下平流層 3 階與 2 階曲線擬合擾動場的能量上傳率 （a） 和能量上傳率相關系數(shù) （b） 的月變化Fig.13 Annual variations of FWUPE from disturbance fields obtained by second and third order curve fitting （a） and FWUPE correlation coefficient （b） in the lower stratosphere above Taiyuan during 2014—2017

圖14 2014—2017 年太原地區(qū)下平流層 3 階與 2 階曲線擬合擾動場的能量上傳率 （a） 和能量上傳率偏差 （b） 在不同區(qū)間的占有率Fig.14 Occupancy distribution characteristics of FWUPE from disturbance fields obtained by second and third order curve fitting （a） and FWUPE deviation （b） in different interval ranges in the lower stratosphericabove Taiyuan during 2014—2017

（1）基于2—4階曲線擬合方法，得到高垂直分辨率探空資料緯向風擾動場、經(jīng)向風擾動場及溫度擾動場的值變化范圍及隨高度變化趨勢均存在較明顯差異。

（2）3階、4階曲線擬合與2階曲線擬合方法擾動場之間的相關性存在較明顯差異；其中緯向風擾動場相關系數(shù)值大于0.8的占有率分別為62.54%和27.38%，相應的經(jīng)向風擾動場分別為45.46%和17.77%，相應的溫度擾動場分別為67.54%和40.08%。

（3）不同階曲線擬合方法擾動場對大氣重力波參數(shù)大小、變化趨勢均具有重大影響；相對3階曲線擬合方法擾動場，除能量上傳率外，基于2階曲線擬合方法擾動場獲取的大氣重力波參數(shù)均偏大（表2）。

（4）3階與2階曲線擬合方法擾動場，獲取的不同重力波參數(shù)之間相關性存在較明顯差異；相關性最強為能量上傳率（0.85），最弱為水平波長（0.55）（表2）。

（5）3階與2階曲線擬合方法擾動場計算得到的重力波參數(shù)及其偏差，在不同區(qū)間內(nèi)的占有率存在較明顯差異；其占有率最大值及對應區(qū)間見表3。

目前，已有研究成果尚未證實哪階曲線擬合擾動場獲取的重力波參數(shù)氣候特征與實況更接近，并且針對不同階曲線擬合擾動場對重力波參數(shù)氣候特征影響的研究較少，Dutta等（2017）利用2012年4月30日—5月4日（時間分辨率6 h）共20次高垂直分辨率探空資料，同樣發(fā)現(xiàn)不同階曲線擬合得到的擾動場廓線存在差異，但沒給出相應的重力波參數(shù)平均值、年內(nèi)變化趨勢以及重力波參數(shù)在不同區(qū)間內(nèi)占有率的差異特征；與此同時，一方面，重力波參數(shù)氣候特征是確定大氣模式中重力波參數(shù)化方案的基礎數(shù)據(jù)，通過該基礎數(shù)據(jù)有利于研制更合適的重力波參數(shù)化方案，從而進一步提高大氣模式對多種不同尺度天氣現(xiàn)象的預報能力；另一方面，本研究發(fā)現(xiàn)，雖然2階與3階曲線擬合方法擾動場的相似程度較高，但基于它們擾動場得到大氣重力波參數(shù)的數(shù)值及變化趨勢均存在較明顯差異；考慮到目前中外常采用不同階曲線擬合方法擾動場獲取大氣重力波參數(shù)氣候特征，且較少關注不同階曲線擬合方法擾動場對其氣候特征的影響，雖然本文對不同階曲線擬合擾動場對重力波參數(shù)氣候特征影響進行了對比研究，并發(fā)現(xiàn)它們之間存在較明顯差異，如何確定幾階曲線擬合擾動場得到的重力波參數(shù)氣候特征更接近實況，需要考慮以各階曲線擬合方法擾動場得到的重力波參數(shù)氣候特征為試驗方案，研制大氣模式重力波參數(shù)化方案，隨后對多種不同尺度天氣系統(tǒng)進行數(shù)值模擬，最后，選取總體上模擬準確率最高的大氣模式重力波參數(shù)化方案為最佳方案。

表2 基于2階、3階曲線擬合擾動場計算得到的重力波參數(shù)平均值及其之間的相關系數(shù)Table 2 The average gravity wave parameters and their correlation coefficients based on the second-order and third-order curve fitting disturbance field

表3 在不同區(qū)間范圍內(nèi)，3階、2階曲線擬合擾動場計算得到重力波參數(shù)值和偏差值的占有率最大值及對應的區(qū)間Table 3 In different intervals，The maximum occupancy rate of gravity wave parameter values and deviation values，and the corresponding interval range obtained obtained by the second-order and third-order curves fitted disturbance fields