王曉婷 李森 史小康 黃江平 胡艷冰

(北京航空氣象研究所，北京 100085)

引言

平流層是介于對(duì)流層與中層之間的大氣層，層內(nèi)氣流平穩(wěn)，沒(méi)有強(qiáng)烈的對(duì)流運(yùn)動(dòng)[1]。在平流層冬夏季節(jié)轉(zhuǎn)換期間，緯向東西風(fēng)區(qū)之間會(huì)存在低風(fēng)速的緯向風(fēng)轉(zhuǎn)換層[2]，被稱為零風(fēng)層、弱風(fēng)層，它是部署平流層飛艇、高空氣球等臨近空間飛行器執(zhí)行對(duì)地觀測(cè)、網(wǎng)絡(luò)通信、環(huán)境監(jiān)測(cè)、防災(zāi)減災(zāi)、預(yù)警探測(cè)等任務(wù)的理想環(huán)境[3-5]。呂達(dá)仁等[6]針對(duì)臨近空間飛行器的特點(diǎn)，提出了準(zhǔn)零風(fēng)層(Quasi-Zero Wind Layer，QZWL)的概念，即平流層準(zhǔn)零風(fēng)層一般指平流層下層20 km高度附近的大氣層，上、下層緯向風(fēng)風(fēng)向相反，同時(shí)南北風(fēng)分量亦很小。平流層準(zhǔn)零風(fēng)層的低風(fēng)速特點(diǎn)不僅對(duì)于平流層飛艇、高空氣球等臨近空間浮空器有重要的應(yīng)用價(jià)值，同時(shí)也是影響臨近空間浮空器姿態(tài)及動(dòng)力控制的重要因素[7]。因此，加強(qiáng)平流層準(zhǔn)零風(fēng)層的規(guī)律及預(yù)報(bào)方法研究，對(duì)于加強(qiáng)臨近空間浮空器的開(kāi)發(fā)及應(yīng)用具有十分重要的意義。

在準(zhǔn)零風(fēng)層的概念提出后，國(guó)內(nèi)一些學(xué)者相繼開(kāi)展了有關(guān)準(zhǔn)零風(fēng)層的特征分析和影響機(jī)理研究。肖存英等[8]利用ERA-40再分析風(fēng)場(chǎng)資料首次分析了中國(guó)上空平流層準(zhǔn)零風(fēng)層的特點(diǎn)及其隨季節(jié)和地理位置的變化特征，指出低緯地區(qū)一般在冬季和初春存在準(zhǔn)零風(fēng)層；中高緯地區(qū)一般在春末和夏季存在準(zhǔn)零風(fēng)層；而中低緯過(guò)渡區(qū)域準(zhǔn)零風(fēng)層的存在與準(zhǔn)兩年震蕩(QBO)有關(guān)。陶夢(mèng)初等[9]根據(jù)不同的形成機(jī)制將準(zhǔn)零風(fēng)層劃分為2類，即由各高度處于不同準(zhǔn)兩年振蕩QBO位相期間形成的第1類準(zhǔn)零風(fēng)層和由平流層低層經(jīng)向溫度梯度逆轉(zhuǎn)形成的第2類準(zhǔn)零風(fēng)層，并分別研究了2類準(zhǔn)零風(fēng)層的變化規(guī)律。楊婷婷[10]分析了北半球平流層準(zhǔn)零風(fēng)層的分布特征，并分別探討了QBO、平流層經(jīng)向溫度梯度和行星波渦動(dòng)通量輸送對(duì)準(zhǔn)零風(fēng)層的影響。陳伯青[11]利用ERA-Interim再分析資料與庫(kù)爾勒站探空資料，分析了中國(guó)上空準(zhǔn)零風(fēng)層的時(shí)空分布特征及其影響因素。劉杰等[12]利用NCEP再分析風(fēng)場(chǎng)資料和探空資料研究了新疆地區(qū)平流層底層風(fēng)場(chǎng)的特征，結(jié)果表明從春季到秋季準(zhǔn)零風(fēng)層存在明顯的空間變化特征，其起始高度、穩(wěn)定高度呈“U”型分布。前人的研究結(jié)果表明，中國(guó)上空準(zhǔn)零風(fēng)層的分布特征隨季節(jié)和地理位置而變化。中緯度地區(qū)的準(zhǔn)零風(fēng)層存續(xù)時(shí)間最長(zhǎng)，其中新疆地區(qū)上空準(zhǔn)零風(fēng)層的累計(jì)天數(shù)可超過(guò)100 d[11]，且在夏季穩(wěn)定存在。

目前可用于風(fēng)場(chǎng)建模和預(yù)測(cè)的方法主要包括2類：數(shù)值天氣預(yù)報(bào)方法[13-15]和統(tǒng)計(jì)模型方法[16-18]。有些學(xué)者開(kāi)展了針對(duì)平流層風(fēng)場(chǎng)預(yù)測(cè)的研究。陳伯青[11]研究了基于Weibull函數(shù)的風(fēng)速分布擬合模型對(duì)于中低緯度地區(qū)準(zhǔn)零風(fēng)層分析的適用性，表明根據(jù)Weibull 分布計(jì)算特定的累積概率風(fēng)速值，可以作為選取適宜平流層飛艇運(yùn)行的低風(fēng)速風(fēng)場(chǎng)的判據(jù)。李魁等[19]提出一種基于本征正交分解(POD)的風(fēng)場(chǎng)數(shù)據(jù)降階方法，并提出一種可以對(duì)平流層風(fēng)場(chǎng)進(jìn)行預(yù)測(cè)的Fourier模型。龍遠(yuǎn)等[20]研究了基于PSO-BP神經(jīng)網(wǎng)絡(luò)的平流層區(qū)域風(fēng)場(chǎng)建模與快速預(yù)測(cè)方法。前人的工作大多側(cè)重于平流層風(fēng)場(chǎng)的統(tǒng)計(jì)建模，有關(guān)平流層風(fēng)場(chǎng)數(shù)值預(yù)報(bào)的研究尚不多見(jiàn)。鑒于此，本文以夏季新疆地區(qū)的準(zhǔn)零風(fēng)層為代表，開(kāi)展平流層低層風(fēng)場(chǎng)的特征分析與數(shù)值模擬研究，探索利用高分辨率數(shù)值模式開(kāi)展平流層風(fēng)場(chǎng)預(yù)報(bào)的可行性，為臨近空間浮空器的風(fēng)場(chǎng)預(yù)測(cè)提供手段。

1 資料

新疆地區(qū)夏季平流層低層風(fēng)場(chǎng)特征分析采用美國(guó)國(guó)家環(huán)境預(yù)報(bào)中心(NCEP)提供的全球分析資料(FNL)，水平分辨率為0.25°×0.25°，垂直方向31層(1000～1 hPa)，時(shí)間分辨率6 h。數(shù)值模擬采用NCEP發(fā)布的全球預(yù)報(bào)系統(tǒng)資料(GFS)，水平分辨率為0.25°×0.25°，垂直方向31層(1000～1 hPa)，時(shí)間分辨率3 h。采用新疆烏魯木齊站的探空資料用于模擬結(jié)果的檢驗(yàn)，該資料包括垂直各層的氣壓、高度、氣溫、露點(diǎn)、風(fēng)向、風(fēng)速等要素，時(shí)間分辨率12 h。

2 新疆地區(qū)夏季平流層低層風(fēng)場(chǎng)的特征分析

2.1 緯向風(fēng)隨緯度時(shí)間的變化特征

中緯度地區(qū)準(zhǔn)零風(fēng)層的形成主要跟平流層緯向風(fēng)的季節(jié)轉(zhuǎn)換有關(guān)，因此從緯向風(fēng)轉(zhuǎn)換的規(guī)律可以分析得到準(zhǔn)零風(fēng)層的變化規(guī)律。本文采用FNL資料分析新疆地區(qū)緯向風(fēng)場(chǎng)的轉(zhuǎn)換規(guī)律，為進(jìn)一步說(shuō)明FNL資料的可信度，選取烏魯木齊站為代表，對(duì)2017年6—8月探空資料與FNL資料平流層緯向風(fēng)隨時(shí)間的演變趨勢(shì)進(jìn)行對(duì)比分析(圖1)，二者的演變趨勢(shì)非常接近，70 hPa、50 hPa、30 hPa二者的相關(guān)系數(shù)分別為0.89，0.88和0.93，均通過(guò)0.05顯著性檢驗(yàn)，表明利用FNL再分析資料開(kāi)展平流層低層風(fēng)場(chǎng)的特征分析是可行的。在此基礎(chǔ)上，利用2017年6—8月的FNL再分析資料分析了沿85°E的30 hPa、50 hPa、70 hPa緯向風(fēng)隨緯度和時(shí)間的變化(圖2)。可以看到，在30 hPa高度上新疆大部分地區(qū)整個(gè)夏季盛行東風(fēng)。在50 hPa高度上，緯向風(fēng)隨時(shí)間的演變存在明顯的經(jīng)向差異。具體表現(xiàn)為，6月初至8月中旬前后，新疆地區(qū)自北向南西風(fēng)減弱，東風(fēng)增強(qiáng)；8月下旬開(kāi)始，南部的東風(fēng)減弱，北部的西風(fēng)增強(qiáng)。在70 hPa高度上，除40°N以南地區(qū)在6月下旬至8月初受東風(fēng)控制以外，其余地區(qū)和時(shí)段主要盛行西風(fēng)。從夏季50 hPa緯向風(fēng)隨時(shí)間的演變情況也可以看出，整個(gè)夏季緯向風(fēng)呈現(xiàn)先西風(fēng)轉(zhuǎn)東風(fēng)，再東風(fēng)轉(zhuǎn)西風(fēng)的特征。此外，緯向風(fēng)的轉(zhuǎn)換存在經(jīng)向差異，西風(fēng)轉(zhuǎn)東風(fēng)的起始時(shí)間新疆南部地區(qū)早于北部地區(qū)，而東風(fēng)轉(zhuǎn)西風(fēng)的起始時(shí)間北部地區(qū)早于南部地區(qū)。同時(shí)也說(shuō)明，新疆地區(qū)準(zhǔn)零風(fēng)層出現(xiàn)的時(shí)間南部地區(qū)早于北部地區(qū)，而準(zhǔn)零風(fēng)層結(jié)束的時(shí)間則是北部地區(qū)早于南部地區(qū)。

圖1 2017年6—8月烏魯木齊站探空資料(紅色)與FNL再分析資料(藍(lán)色)平流層低層緯向風(fēng)隨時(shí)間演變 (圖中實(shí)線：70 hPa；長(zhǎng)虛線：50 hPa；短虛線：30 hPa)

圖2 2017年6—8月沿85°E的緯向風(fēng)隨緯度-時(shí)間的分布:(a)30 hPa,(b)50 hPa,(c)70 hPa

2.2 緯向風(fēng)隨時(shí)間高度的變化特征

為進(jìn)一步揭示準(zhǔn)零風(fēng)層隨時(shí)間的變化特征，以烏魯木齊站為代表利用FNL資料繪制站點(diǎn)上空平流層緯向風(fēng)隨高度和時(shí)間的變化圖(圖3)。圖3中零線為緯向風(fēng)的零風(fēng)速線，可以表征準(zhǔn)零風(fēng)層高度?？梢钥吹剑瑴?zhǔn)零風(fēng)層在整個(gè)夏季均穩(wěn)定存在。準(zhǔn)零風(fēng)層下方為弱西風(fēng)層，準(zhǔn)零風(fēng)層上方為弱東風(fēng)層。整個(gè)夏季準(zhǔn)零風(fēng)層處于70～40 hPa之間，準(zhǔn)零風(fēng)層高度隨時(shí)間存在小幅波動(dòng)變化，但總體上呈現(xiàn)先降低而后升高的變化趨勢(shì)，這種變化趨勢(shì)與平流層緯向風(fēng)的季節(jié)轉(zhuǎn)換有很好的對(duì)應(yīng)關(guān)系。夏初平流層中下層由西風(fēng)轉(zhuǎn)為東風(fēng)，之后伴隨東風(fēng)逐漸增強(qiáng)準(zhǔn)零風(fēng)層高度逐漸降低，夏末隨著東風(fēng)逐漸減弱，準(zhǔn)零風(fēng)層高度又逐漸恢復(fù)到夏初的位置。新疆地區(qū)夏季準(zhǔn)零風(fēng)層的形成機(jī)制可以用熱成風(fēng)原理來(lái)解釋。根據(jù)熱成風(fēng)原理平流層緯向風(fēng)的分布取決于對(duì)流層上部的緯向風(fēng)和平流層的緯向溫度梯度，北半球夏季在輻射過(guò)程的影響下，平流層大氣平均緯向溫度梯度由低緯指向高緯(為正值)，因此對(duì)流層上部的西風(fēng)會(huì)隨高度逐漸減小，并在某一高度減小為零，再向上轉(zhuǎn)為東風(fēng)，在下層西風(fēng)與上層?xùn)|風(fēng)之間形成準(zhǔn)零風(fēng)層[8]。

圖3 2017年6—8月烏魯木齊站緯向風(fēng)隨高度-時(shí)間的分布

2.3 水平風(fēng)速、風(fēng)向的頻率分布

2.3.1 水平風(fēng)速的頻率分布

圖4為利用FNL再分析資料統(tǒng)計(jì)的夏季新疆地區(qū)30 hPa、50 hPa、70 hPa各層水平風(fēng)速的頻率分布。可以看出，整個(gè)夏季30 hPa水平風(fēng)速概率集中分布在5～10 m/s之間；其次6月和8月水平風(fēng)速落在0～5 m/s區(qū)間的頻率較高，而7月水平風(fēng)速落在10～15 m/s區(qū)間的頻率較高。50 hPa水平風(fēng)速頻率主要集中分布在0～5 m/s，其次是5～10 m/s，各月水平風(fēng)速頻率的分布比較一致。70 hPa水平風(fēng)速頻率主要分布在0～5 m/s和5～10 m/s兩個(gè)區(qū)間，且各月存在明顯的差異，低風(fēng)速出現(xiàn)的頻率隨月份顯著減少，而高風(fēng)速出現(xiàn)的頻率隨月份顯著增加。對(duì)各層次低風(fēng)速頻率分布的對(duì)比分析還發(fā)現(xiàn)，50 hPa準(zhǔn)零風(fēng)層出現(xiàn)的累計(jì)頻率最高，是新疆地區(qū)夏季準(zhǔn)零風(fēng)層存在的主要層次。30 hPa低風(fēng)速出現(xiàn)的頻率存在先減小而后增大的現(xiàn)象，也從側(cè)面說(shuō)明夏季準(zhǔn)零風(fēng)層高度隨時(shí)間的變化存在先降低而后升高的現(xiàn)象。

圖4 新疆夏季水平風(fēng)速的頻率分布:(a)30 hPa,(b)50 hPa,(c)70 hPa

2.3.2 水平風(fēng)向的頻率分布

圖5為利用FNL再分析資料統(tǒng)計(jì)的夏季新疆地區(qū)30 hPa、50 hPa、70 hPa各層水平風(fēng)向的頻率分布?？梢钥吹剑麄€(gè)夏季30 hPa的水平風(fēng)向主要為東南風(fēng)；50 hPa的水平風(fēng)向主要為偏東風(fēng)，其中以東風(fēng)的頻率最大，其次是東北風(fēng)和東南風(fēng)；70 hPa的水平風(fēng)向主要為偏西風(fēng)，其中以西南風(fēng)的頻率最大，其次是西風(fēng)和西北風(fēng)。

圖5 2017年新疆地區(qū)夏季平流層風(fēng)向玫瑰圖

3 平流層低層風(fēng)場(chǎng)的數(shù)值模擬

3.1 模擬方案

本文數(shù)值模擬采用的模式為WRFV4.1.2版本。模式中心點(diǎn)設(shè)為(87.32°E，42.4°N)，水平格點(diǎn)數(shù)為501×501，水平分辨率為1 km，垂直方向84層。經(jīng)過(guò)多組物理方案的對(duì)比試驗(yàn)，最終采用的物理方案包括：NSSL2-moment云微物理過(guò)程方案，MYNN地表過(guò)程方案，RUC陸面過(guò)程方案，MYNN邊界層過(guò)程方案，rrtmg短波輻射和長(zhǎng)波輻射方案。數(shù)值模擬試驗(yàn)時(shí)間為2019年8月1—10日，每日起報(bào)2次，起報(bào)時(shí)間分別為00:00和12:00(世界時(shí)，下同)，模擬時(shí)效為12 h。采用GFS資料作為驅(qū)動(dòng)場(chǎng)，每3 h更新1次。

3.2 模擬結(jié)果

3.2.1 緯向風(fēng)、經(jīng)向風(fēng)廓線的模擬

利用新疆烏魯木齊站的探空資料對(duì)平流層低層各層經(jīng)向風(fēng)和緯向風(fēng)風(fēng)速誤差進(jìn)行檢驗(yàn)。圖6給出了2019年8月1—10日時(shí)段內(nèi)，緯向風(fēng)和經(jīng)向風(fēng)的均方根誤差隨高度分布，在準(zhǔn)零風(fēng)層高度范圍內(nèi)(一般18～25 km)經(jīng)向風(fēng)和緯向風(fēng)的均方根誤差均小于3 m/s。總體來(lái)看，平流層中低層緯向風(fēng)和經(jīng)向風(fēng)的模擬誤差除個(gè)別高度接近或超過(guò)3 m/s，其余高度均方根誤差在2 m/s左右。

3.2.2 準(zhǔn)零風(fēng)層的模擬

圖7給出了烏魯木齊2019年8月1日00:00—12:00緯向風(fēng)的模擬結(jié)果，從圖中可以清晰的看到平流層準(zhǔn)零風(fēng)層位于20 km附近，下層為弱西風(fēng)，上層為弱東風(fēng)。為了更精細(xì)地刻畫準(zhǔn)零風(fēng)層高度、風(fēng)向、風(fēng)速等參數(shù)的特征，本文在高分辨率數(shù)值模擬的基礎(chǔ)上對(duì)模擬結(jié)果在垂直層次上進(jìn)行了加密插值處理，垂直分辨率達(dá)到200 m，并采用如下的判別方法確定準(zhǔn)零風(fēng)層的高度：①對(duì)水平方向上的任意格點(diǎn)，從對(duì)流層頂開(kāi)始向上判別是否存在緯向風(fēng)的轉(zhuǎn)換層(即下層為西風(fēng)且上層為東風(fēng)的層次)，同時(shí)滿足風(fēng)速小于等于8 m/s。如存在緯向風(fēng)轉(zhuǎn)換層則將該層的高度確定為準(zhǔn)零風(fēng)層的起始高度；②自準(zhǔn)零風(fēng)層起始高度向上逐層判別，將同時(shí)滿足風(fēng)速小于等于8 m/s和風(fēng)向?yàn)槠珫|風(fēng)(即U<0)的最大高度，作為準(zhǔn)零風(fēng)層的頂高。

圖7 2019年8月1日00:00—12:00烏魯木齊上空緯向風(fēng)場(chǎng)模擬

依據(jù)準(zhǔn)零風(fēng)層診斷方法，得到2019年8月1—10日準(zhǔn)零風(fēng)層起始高度及其對(duì)應(yīng)風(fēng)向和風(fēng)速的模擬與探空資料對(duì)比結(jié)果(圖8)。需要說(shuō)明的是，2019年8月8日00:00的探空資料因探測(cè)高度偏低，未能分析得到實(shí)際準(zhǔn)零風(fēng)層的相關(guān)參數(shù)，其他時(shí)次探空數(shù)據(jù)均可分析得到準(zhǔn)零風(fēng)層。從準(zhǔn)零風(fēng)層起始高度及其對(duì)應(yīng)的風(fēng)向、風(fēng)速隨時(shí)間演變的模擬結(jié)果與探空資料對(duì)比來(lái)看，兩者隨時(shí)間的演變趨勢(shì)比較一致，但模擬的準(zhǔn)零風(fēng)層風(fēng)速較實(shí)際偏大。整個(gè)模擬時(shí)段內(nèi)準(zhǔn)零風(fēng)層起始高度的平均絕對(duì)誤差為467 m，起始高度對(duì)應(yīng)的風(fēng)速均方根誤差為1.75 m/s。

圖8 2019年8月1—10日烏魯木齊站模擬與觀測(cè)的準(zhǔn)零風(fēng)層參數(shù)隨時(shí)間演變:(a)準(zhǔn)零風(fēng)層起始高度,(b)風(fēng)速,(c)風(fēng)向

4 結(jié)論與討論

本文利用2017年6—8月的FNL再分析風(fēng)場(chǎng)資料，分析了新疆地區(qū)夏季平流層低層風(fēng)場(chǎng)的時(shí)空演變特征，開(kāi)展了平流層低層風(fēng)場(chǎng)的高分辨率數(shù)值模擬并利用烏魯木齊站的探空資料對(duì)模擬結(jié)果進(jìn)行了檢驗(yàn)。

(1)新疆地區(qū)夏季平流層緯向風(fēng)隨時(shí)間的變化存在經(jīng)向差異，這一特征在50 hPa高度最為明顯；同時(shí)準(zhǔn)零風(fēng)層開(kāi)始和結(jié)束的時(shí)間也存在經(jīng)向差異，準(zhǔn)零風(fēng)層的出現(xiàn)時(shí)間南部地區(qū)早于北部地區(qū)，而準(zhǔn)零風(fēng)層的結(jié)束時(shí)間則相反。

(2)新疆地區(qū)夏季準(zhǔn)零風(fēng)層處于70～40 hPa之間，其高度隨時(shí)間呈先降低而后升高的變化趨勢(shì)，這種變化趨勢(shì)與平流層緯向風(fēng)的季節(jié)轉(zhuǎn)換有很好的對(duì)應(yīng)關(guān)系。

(3)對(duì)各層風(fēng)速、風(fēng)向頻率分布的統(tǒng)計(jì)分析表明，50 hPa準(zhǔn)零風(fēng)層出現(xiàn)的累計(jì)頻率最高，是新疆地區(qū)夏季準(zhǔn)零風(fēng)層存在的主要層次。30 hPa低風(fēng)速出現(xiàn)的頻率存在先減小而后增大的現(xiàn)象，從側(cè)面也說(shuō)明夏季準(zhǔn)零風(fēng)層高度隨時(shí)間的變化存在先降低而后升高的現(xiàn)象。

(4)對(duì)平流層低層風(fēng)場(chǎng)的模擬檢驗(yàn)結(jié)果表明，在準(zhǔn)零風(fēng)層高度范圍內(nèi)(一般18～25 km)經(jīng)向風(fēng)和緯向風(fēng)的均方根誤差均小于3 m/s。模擬的準(zhǔn)零風(fēng)層起始高度及其對(duì)應(yīng)的風(fēng)向、風(fēng)速隨時(shí)間演變趨勢(shì)與探空資料分析的結(jié)果比較一致，模擬時(shí)段內(nèi)準(zhǔn)零風(fēng)層起始高度的平均絕對(duì)誤差為467 m，起始高度對(duì)應(yīng)的風(fēng)速均方根誤差為1.75 m/s。

文中以新疆地區(qū)為代表對(duì)中緯度地區(qū)的準(zhǔn)零風(fēng)層開(kāi)展了數(shù)值模擬，并初步驗(yàn)證了高分辨率數(shù)值模式在準(zhǔn)零風(fēng)層預(yù)報(bào)中的適用性，但由于模擬范圍有限，對(duì)可能影響準(zhǔn)零風(fēng)層變化的天氣系統(tǒng)刻畫能力不足，有待后續(xù)開(kāi)展深入研究。