塔克拉干沙漠腹地沙塵暴天氣下風(fēng)速脈動(dòng)特征

鄭先念, 楊興華, 何 清,4, 金莉莉, 楊 帆,艾力·買(mǎi)買(mǎi)提依明, 周成龍, 張建濤

(1.新疆氣候中心, 烏魯木齊 830002; 2.中國(guó)氣象局烏魯木齊沙漠氣象研究所, 烏魯木齊 830002;3.中國(guó)氣象局塔克拉瑪干沙漠氣象野外科學(xué)試驗(yàn)基地, 新疆 塔中 841000; 4.新疆氣象臺(tái), 烏魯木齊 830002)

瞬時(shí)風(fēng)速和平均風(fēng)速隨時(shí)間的變化總趨勢(shì)是一致的，風(fēng)沙研究中通常將平均風(fēng)速代替瞬時(shí)風(fēng)速來(lái)研究風(fēng)對(duì)輸沙量的影響，這種處理易于掌握風(fēng)沙運(yùn)移的宏觀特征[1]。但也有局限性，因?yàn)樽匀粻顟B(tài)下引起沙粒運(yùn)動(dòng)的風(fēng)幾乎是湍流的，各點(diǎn)的流速大小和方向隨時(shí)間脈動(dòng)表現(xiàn)出一定的陣性特征[2]。如果不考慮風(fēng)速脈動(dòng)，根據(jù)定常風(fēng)速計(jì)算得到的地表輸沙過(guò)程是平穩(wěn)的，而野外實(shí)際輸沙過(guò)程非平穩(wěn)且具有一定結(jié)構(gòu)特征[3]。并且風(fēng)速變化的周期、頻率對(duì)輸沙量有很大影響[4]，因此需要考慮風(fēng)速的脈動(dòng)性。前人研究還發(fā)現(xiàn)風(fēng)速脈動(dòng)使床面沙粒受到躍移上升力和剪切力，引發(fā)沙粒振動(dòng)，使沙粒易于起運(yùn)[5]，影響沙粒的躍移軌跡，且風(fēng)速脈動(dòng)與輸沙率之間具有良好的相關(guān)關(guān)系，而沙粒的運(yùn)動(dòng)又可以影響風(fēng)速脈動(dòng)，進(jìn)而影響沙粒自身的輸移[6]。所以湍流對(duì)于風(fēng)沙輸移有重要的意義。

脈動(dòng)風(fēng)速在沙漠一綠洲過(guò)渡帶典型下墊面(裸樣沙地、檉柳樣地和駱駝刺樣地)上都近似服從高斯分布，風(fēng)速脈動(dòng)規(guī)律沒(méi)有因植被的存在而發(fā)生改變[7]，在晴天和沙塵暴天氣下，近地表風(fēng)速與風(fēng)速脈動(dòng)強(qiáng)度呈正相關(guān)，并受到沙塵一定的抑制作用[8]。風(fēng)速脈動(dòng)強(qiáng)度和風(fēng)速脈動(dòng)的相對(duì)值整體上都隨著植被覆蓋面積的增大呈增大趨勢(shì)[9]。而在戈壁地表風(fēng)速脈動(dòng)具有非平穩(wěn)性，隨著風(fēng)速增強(qiáng)而增大，尤其是在大風(fēng)天氣下研究風(fēng)蝕不能忽略風(fēng)速的脈動(dòng)性[10]。同時(shí)，通過(guò)對(duì)北方兩次沙塵天氣的研究發(fā)現(xiàn)，沙塵暴期間大氣中可吸入顆粒物的濃度與風(fēng)速存在顯著的正相關(guān)關(guān)系[11]。能見(jiàn)度、PM10，PM2.5也隨風(fēng)速脈動(dòng)而產(chǎn)生波動(dòng)變化，地面風(fēng)速、PM10，PM2.5之間相互作用，呈明顯脈動(dòng)特征[12]。

作為世界第二大流動(dòng)沙漠，塔克拉瑪干沙漠不僅是我國(guó)沙塵暴頻發(fā)地區(qū)之一，也是亞洲沙塵暴產(chǎn)塵的重要源地，該地區(qū)大風(fēng)條件下?lián)P起的沙塵會(huì)隨著高空西風(fēng)環(huán)流影響整個(gè)東亞地區(qū)。然而，前人針對(duì)野外沙塵暴天氣條件下的沙漠腹地近地表80 m高度范圍的風(fēng)速脈動(dòng)特征及影響的研究較少。鑒于此，通過(guò)野外實(shí)測(cè)塔克拉瑪干沙漠腹地沙塵暴天氣條件下的風(fēng)速資料分析不同高度風(fēng)速脈動(dòng)特征，完善了沙漠腹地風(fēng)與沙粒運(yùn)動(dòng)的關(guān)系、防沙治沙以及風(fēng)沙互饋機(jī)制等方面的研究[13]，為進(jìn)一步研究湍流情況下風(fēng)蝕和風(fēng)揚(yáng)粉塵釋放奠定基礎(chǔ)。

1 試驗(yàn)材料與方法

1.1 研究區(qū)概況

研究區(qū)位于塔克拉瑪干沙漠腹地塔中地區(qū)(83°39′E，38°58′N(xiāo))。該地區(qū)氣候極端干旱，植被稀少，地表為流沙覆蓋，只有塔中作業(yè)區(qū)所在的壟間地有人工建造面積約為2×2 km2地形呈NE—WS走向的人工綠地，種植有梭梭、檉柳、駱駝刺和一些野生蘆葦?shù)裙躺衬秃的望}堿的植物，生長(zhǎng)高度夏秋兩季高度可達(dá)2 m左右，覆蓋度較高。試驗(yàn)站位于塔中作業(yè)區(qū)內(nèi)東面為高大縱向沙隴，西面為坡度相對(duì)平緩的縱向沙隴，北面距沙漠公路1.2 km，沙漠公路兩側(cè)也是人工種植的固沙植物，其余多為流沙地表；該區(qū)年平均氣溫為13.6℃，極端最高氣溫達(dá)46.0℃，極端最低氣溫達(dá)-32.0℃；平均年降水量25.9 mm左右，平均年蒸發(fā)量高達(dá)3 812.3 mm；平均年風(fēng)速為2.3 m/s，平均年揚(yáng)沙日數(shù)為59.6 d，平均年沙塵暴日數(shù)為15.8 d，主導(dǎo)風(fēng)向以NE，NNE主[14]。

1.2 數(shù)據(jù)資料與研究方法

采用風(fēng)速和風(fēng)向資料來(lái)自于安裝在沙壟間谷地的80 m鐵塔梯度探測(cè)系統(tǒng)，鐵塔共有10個(gè)呈對(duì)數(shù)間距分布的測(cè)量梯度，高度分別為0.5 m，1 m，2 m，4 m，10 m，20 m，32 m，47 m，63 m和80 m，因?yàn)轱L(fēng)速在近地層呈對(duì)數(shù)分布。探測(cè)的基本要素為：溫度、濕度、氣壓、風(fēng)速和水平風(fēng)向等。其中風(fēng)速儀采用英國(guó)Gill公司生產(chǎn)的WindObserver Ⅱ型二維風(fēng)速風(fēng)向傳感器，啟動(dòng)風(fēng)速0.01 m/s,風(fēng)速量程為0～75 m/s；風(fēng)向測(cè)量范圍0°～359°，分辨率為1°，觀測(cè)的風(fēng)速和風(fēng)向數(shù)據(jù)有3 s數(shù)據(jù)、5 s數(shù)據(jù)、10 s數(shù)據(jù)和分鐘平均數(shù)據(jù)及小時(shí)平均數(shù)據(jù)。本文計(jì)算風(fēng)速脈動(dòng)特征采用3 s間隔觀測(cè)的數(shù)據(jù)，分別選取了北京時(shí)間2018年5月20日15：26—15：35風(fēng)向?yàn)镹NE，5月24日16：39—16：48風(fēng)向?yàn)镋NE的兩次典型沙塵暴天氣的數(shù)據(jù)；計(jì)算風(fēng)速脈動(dòng)與能見(jiàn)度，PM10，PM2.5關(guān)系采用的是1 min間隔觀測(cè)的風(fēng)速數(shù)據(jù)(為了與能見(jiàn)度、PM10，PM2.5采集間隔時(shí)間一致)，PM10，PM2.5采用的是Grimm180觀測(cè)的1 min間隔的數(shù)據(jù)，選取的是2018年5月24日15：36—18：15(北京時(shí))一次沙塵暴過(guò)程的數(shù)據(jù)。

脈動(dòng)風(fēng)速值以瞬時(shí)風(fēng)速與相應(yīng)時(shí)段10 min平均風(fēng)速之差表示，用以反映風(fēng)速波動(dòng)幅度；以脈動(dòng)風(fēng)速的均方根表示風(fēng)速脈動(dòng)強(qiáng)度(uv)：

(1)

2 結(jié)果與分析

2.1 風(fēng)速脈動(dòng)

圖1，2可看出，不同高度風(fēng)速脈動(dòng)值在時(shí)間序列上波動(dòng)性均具有一致性，即同步增大或減小，脈動(dòng)頻率不穩(wěn)定，無(wú)明顯周期性，說(shuō)明自然風(fēng)除具有脈動(dòng)性外，還有間歇性特點(diǎn)。2次沙塵暴天氣過(guò)程的風(fēng)速脈動(dòng)最大波動(dòng)幅度分別為13.59 m/s和10.09 m/s，最小值波動(dòng)幅度分別也達(dá)到了7.53 m/s，6.58 m/s，這么大的脈動(dòng)風(fēng)速對(duì)輸沙過(guò)程定會(huì)產(chǎn)生很大的影響，因?yàn)楦鶕?jù)定常風(fēng)速計(jì)算得到的地表輸沙過(guò)程是平穩(wěn)的，而野外實(shí)際輸沙過(guò)程非平穩(wěn)且具有一定結(jié)構(gòu)特征[3]。同時(shí),相鄰高度間均具有顯著相關(guān)性，經(jīng)計(jì)算相關(guān)系數(shù)均大于0.49、5月20日相鄰高度相關(guān)性在0.49以上、5月24日相鄰高度相關(guān)性在0.62以上，不同高度間相關(guān)性不是很顯著，高度間隔越大相關(guān)性越不顯著。

圖1 5月20日沙塵暴天氣條件下風(fēng)速脈動(dòng)變化

圖2 5月24日沙塵暴天氣條件下風(fēng)速脈動(dòng)變化

各高度層脈動(dòng)風(fēng)速概率分布近似符合高斯函數(shù)，且高度越高，概率分布擬合曲線越陡，脈動(dòng)風(fēng)速分布越集中。2次沙塵暴天氣下風(fēng)速脈動(dòng)的波動(dòng)范圍在0.5～10 m高度層內(nèi)均隨高度的增高而增大，以0.5 m高度的風(fēng)速脈動(dòng)強(qiáng)度與其他高度的表現(xiàn)最顯著，風(fēng)速脈動(dòng)值在0.5 m高度層分別為-3.21～4.32 m/s，-2.77～3.81 m/s，到了10 m高度層上為-5.91～5.16 m/s，-4.78～5.49 m/s；而在32～80 m高度層上脈動(dòng)風(fēng)速的波動(dòng)范圍隨著高度越高而減小，20 m高度層的波動(dòng)范圍5月20日為-6.42～4.99 m/s，5月24日為-4.50～4.42 m/s，80 m高度層的波動(dòng)范圍5月20日為-4.04～3.97 m/s，5月24日為-4.76～3.16 m/s。由此進(jìn)一步驗(yàn)證脈動(dòng)風(fēng)速的分布規(guī)律與高度無(wú)關(guān)，高度僅對(duì)脈動(dòng)風(fēng)速的波動(dòng)范圍和概率分布的集中性產(chǎn)生影響。

2.2 風(fēng)速脈動(dòng)強(qiáng)度

從表1，2可知，2次沙塵暴天氣條件下，隨著高度的降低和風(fēng)速平均值減小，風(fēng)速脈動(dòng)強(qiáng)度呈先增長(zhǎng)后降低的趨勢(shì)，側(cè)面驗(yàn)證了風(fēng)速脈動(dòng)范圍隨高度的增大而先增后減小的規(guī)律。具體表現(xiàn)為，5月20日的沙塵暴天氣下，除2 min和10 min時(shí)段的外，其他各時(shí)段在20 m高度層上風(fēng)速脈動(dòng)強(qiáng)度達(dá)到最大值，在0.5～20 m高度范圍內(nèi)均隨高度的增高而增大、在20～80 m范圍內(nèi)整體呈減小后增大的趨勢(shì)，2 min和10 min時(shí)段則是在10 m高度上增長(zhǎng)到最大值隨后呈曲線增大的變化趨勢(shì)。5月24日的沙塵暴天氣下，除2 min，4 min和5 min時(shí)段外，同樣在其他各時(shí)段風(fēng)速脈動(dòng)強(qiáng)度在20 m高度達(dá)到最大值，在0.5～20 m高度范圍內(nèi)呈隨高度的增高而增大、而20～80 m高度范圍內(nèi)脈動(dòng)強(qiáng)度呈先降低后增的規(guī)律。2 min，4 min和5 min時(shí)段則是在10 m高度上增長(zhǎng)到最大值隨后呈曲線增大的趨勢(shì)。因?yàn)轱L(fēng)速脈動(dòng)強(qiáng)度隨著不同高度間風(fēng)速增加的幅度而變化。

在同一高度層和時(shí)距上，風(fēng)速脈動(dòng)強(qiáng)度與平均風(fēng)速正相關(guān)，即風(fēng)速脈動(dòng)強(qiáng)度值隨平均風(fēng)速值的增大而增大，如10 m高度層10 min時(shí)距上風(fēng)速平均值分別為5月20日11.97 m/s，5月24日9.24 m/s，相對(duì)應(yīng)風(fēng)速脈動(dòng)強(qiáng)度值分別為2.41 m/s，2.09 m/s，風(fēng)速脈動(dòng)、風(fēng)速脈動(dòng)強(qiáng)度均與風(fēng)速相關(guān)，風(fēng)速越大，風(fēng)速脈動(dòng)和風(fēng)速脈動(dòng)強(qiáng)度越大。相同高度不同時(shí)距上，風(fēng)速脈動(dòng)強(qiáng)度隨平均風(fēng)速的變化上比較復(fù)雜，總體呈現(xiàn)出時(shí)距越長(zhǎng)風(fēng)速脈動(dòng)強(qiáng)越大，如10 m高度層風(fēng)速脈動(dòng)強(qiáng)度值，在1 min時(shí)距為5月20日是1.52 m/s，5月24日是0.90 m/s，而增10 min時(shí)距為5月20日增到1.60 m/s，5月24日增到1.21 m/s，同時(shí)，風(fēng)速平均值隨時(shí)距增長(zhǎng)則呈現(xiàn)出減小的變化規(guī)律。因?yàn)椴煌瑫r(shí)距的風(fēng)況資料反應(yīng)真實(shí)情況的準(zhǔn)確度不同[17]平均風(fēng)速時(shí)距越小，越接近瞬時(shí)風(fēng)速狀態(tài)，所以計(jì)算平均風(fēng)速的時(shí)段越小，算出的輸沙率的值更精確。

表1 5月20日沙塵暴天氣下不同高度平均風(fēng)速、脈動(dòng)強(qiáng)度和湍流度隨時(shí)間的變化

同樣，2次沙塵暴天氣下的湍流度0.04～0.35之間變化，屬于高湍流強(qiáng)度，80 m高度范圍內(nèi)均與高度成反比關(guān)系，隨著高度的增高而減小；底層0.5 m和1 m高度湍流度較大，原因是檉柳和梭梭等植物冠層對(duì)該高度層平均風(fēng)速消減的結(jié)果。

2.3 風(fēng)速脈動(dòng)與水平能見(jiàn)度、PM10，PM2.5濃度變化

圖3看出，沙塵暴發(fā)生前，能見(jiàn)度都保持在1 100 m以上，PM10濃度維持在4 000 μg/m3以下波動(dòng)，PM2.5濃度在800 μg/m3以下。沙塵暴發(fā)生時(shí)，PM10濃度在瞬間上升到4 500 μg/m3，PM2.5濃度也上升到800 μg/m3以上，沙塵暴過(guò)程中能見(jiàn)度降到700 m以下。沙塵暴后，轉(zhuǎn)化成揚(yáng)沙天氣，空氣中含有大量的氣溶膠，短時(shí)間內(nèi)能見(jiàn)度還比較低，但開(kāi)始呈現(xiàn)上升趨勢(shì)，PM10，PM2.5濃度最大值區(qū)域?qū)?yīng)能見(jiàn)度最低區(qū)域，沙塵暴過(guò)程中，能見(jiàn)度都維持在1 000 m以下，而PM10濃度值多次達(dá)到峰值6 500 μg/m3，PM2.5濃度值也多次達(dá)到1 000 μg/m3以上，峰值到達(dá)了1 355.9 μg/m3。沙塵暴過(guò)后能見(jiàn)度呈上升趨勢(shì)，對(duì)應(yīng)PM10，PM2.5濃度值呈波動(dòng)下降趨勢(shì)，因?yàn)樯硥m暴結(jié)束后天氣轉(zhuǎn)為揚(yáng)沙天氣。

分析圖4，5可知，風(fēng)速波動(dòng)與PM10，PM2.5/能見(jiàn)度分別呈現(xiàn)明顯的正/反向變化，風(fēng)速波動(dòng)與PM10，PM2.5濃度增大，能見(jiàn)度減小，其中PM10達(dá)到6 500 μg/m3時(shí)不隨風(fēng)速的變化而變化。沙塵暴發(fā)生時(shí)，瞬時(shí)風(fēng)速峰值達(dá)到10.50 m/s時(shí)，PM10，PM2.5的濃度值最高處于6 500 μg/m3，1 355.9 μg/m3，能見(jiàn)度也達(dá)到最低值。地面產(chǎn)生強(qiáng)風(fēng)，使得沙塵粒子懸浮在空氣中，一旦風(fēng)速減小，受到重力沉降作用，懸浮在近地面的PM10，PM2.5迅速增加。此后伴隨著風(fēng)速的周期波動(dòng)，PM10，PM2.5的濃度也相應(yīng)作出調(diào)整，直到沙塵結(jié)束減小到3 000 μg/m3，6 009 μg/m3，能見(jiàn)度又升回1 000 m以上。由于沙塵暴結(jié)束后天氣轉(zhuǎn)為揚(yáng)沙天氣，所以沙塵暴過(guò)后，風(fēng)速?zèng)]有驟減，而是隨時(shí)間的變化呈波動(dòng)減小。

表2 5月24日沙塵暴天氣下不同高度平均風(fēng)速、脈動(dòng)強(qiáng)度和湍流度隨時(shí)間的變化

圖3 能見(jiàn)度與PM10，PM2.5變化

圖4 地面風(fēng)速與PM10，PM2.5隨時(shí)間的變化

圖5 地面風(fēng)速與能見(jiàn)度隨時(shí)間的變化

圖6 風(fēng)速脈動(dòng)與PM10，PM2.5隨時(shí)間變化

為了更加詳細(xì)了解風(fēng)速脈動(dòng)對(duì)PM10，PM2.5濃度值的影響，圖6給出他們之間的關(guān)系。從圖6可見(jiàn)，本次沙塵暴出現(xiàn)時(shí)，風(fēng)速脈動(dòng)最大值達(dá)到2.93 m/s，對(duì)應(yīng)PM10，PM2.5濃度值也達(dá)到了第一次的峰值，隨后風(fēng)速脈動(dòng)仍然存在，其值基本維持在-2.75～3.45 m/s，PM10，PM2.5濃度值隨之出現(xiàn)躍增趨勢(shì)，而負(fù)風(fēng)速脈動(dòng)除少數(shù)幾個(gè)點(diǎn)外均未出現(xiàn)PM10，PM2.5濃度值相應(yīng)增加現(xiàn)象。根據(jù)風(fēng)速脈動(dòng)定義，風(fēng)速脈動(dòng)為正值，意味著短時(shí)間內(nèi)的瞬時(shí)風(fēng)速大于平均風(fēng)速，因此能搬運(yùn)沙粒的數(shù)量較多，相應(yīng)空氣中PM10，PM2.5濃度值增大。

3 討 論

通過(guò)對(duì)塔克拉瑪干沙漠塔中地區(qū)人工綠地的風(fēng)速脈特征分析發(fā)現(xiàn)，在風(fēng)速脈動(dòng)強(qiáng)度隨高度的變化方面與前人[7-8]在塔克拉瑪干沙漠南緣分析得出的隨高度的降低而降低不同,筆者認(rèn)為，一方面這是由于所選的高度層不同，本文所選的是近地層80 m高度范圍。另一方面是本文研究區(qū)的下墊面為植被覆蓋的壟間谷地地形。至于不同風(fēng)速下(5月20日的風(fēng)速明顯比5月24日的大)風(fēng)速脈動(dòng)強(qiáng)度隨高度和的變化在1 min、和2 min，4 min，5 min，10 min時(shí)段與其他時(shí)段變化上不一致，這主要是瞬時(shí)風(fēng)速具陣性特征所致。

PM10和PM2.5在沙塵暴發(fā)生和結(jié)束時(shí)的變化與沈潔等在民勤研究的結(jié)果不同，筆者認(rèn)為，一是二者研究時(shí)所用時(shí)間尺度不同，本文采用的是分鐘數(shù)據(jù)，而沈潔等用的是小時(shí)數(shù)據(jù)；二是本文研究?jī)H研究沙塵暴過(guò)程的變化情況，沙塵暴開(kāi)始前和結(jié)束后都是揚(yáng)沙天氣，因此沒(méi)出現(xiàn)突增突降現(xiàn)象。

4 結(jié) 論

(1) 脈動(dòng)風(fēng)速相鄰高度間均具有顯著相關(guān)性，高度間隔越大相關(guān)性越不顯著。脈動(dòng)風(fēng)速近似符合高斯函數(shù)。風(fēng)速脈動(dòng)隨高度的變化上，呈現(xiàn)出在0.5～32 m高度層脈動(dòng)風(fēng)速的波動(dòng)范圍整體上隨高度的增高而增大，而在32～80 m高度層脈動(dòng)風(fēng)速的波動(dòng)范圍與高度成反比，高度越高，波動(dòng)范圍越小。

(2) 風(fēng)速脈動(dòng)強(qiáng)度在0.5～20 m高度范圍內(nèi)呈隨高度的增高而增大、而20～80 m高度范圍內(nèi)脈動(dòng)強(qiáng)度呈先降低后增的規(guī)律。同一高度層上風(fēng)速平均值隨時(shí)距增長(zhǎng)則呈現(xiàn)出減小的變化規(guī)律。2次沙塵暴天氣下的湍流度在0.04～0.35之間變化，屬于高湍流強(qiáng)度，80 m高度范圍內(nèi)均與高度成反比關(guān)系，并隨著高度的增高而減小。

(3) 沙塵暴過(guò)程中地面風(fēng)速的大小對(duì)水平能見(jiàn)度、PM10，PM2.5具有直接影響。因?yàn)轱L(fēng)對(duì)沙塵粒子的搬運(yùn)具有直接動(dòng)力作用，所以風(fēng)速脈動(dòng)變化與水平能見(jiàn)度呈現(xiàn)反比關(guān)系，與PM10，PM2.5呈正比關(guān)系。