鄭新倩,彭冬梅,吳 燁,周成龍,產(chǎn)秀媚,買買提艾力.買買提依明*

(1.新疆農(nóng)業(yè)氣象臺,新疆烏魯木齊830002;2.中國氣象局烏魯木齊沙漠氣象研究所,新疆烏魯木齊830002;3.塔克拉瑪干沙漠大氣環(huán)境觀測試驗(yàn)站,新疆塔中841000;4.布爾津縣氣象局,新疆布爾津836600)

巴丹吉林沙漠北緣拐子湖近地層湍流能譜特征分析

鄭新倩1,彭冬梅1,吳 燁2,3,周成龍2,3,產(chǎn)秀媚4,買買提艾力.買買提依明2,3*

(1.新疆農(nóng)業(yè)氣象臺,新疆烏魯木齊830002;2.中國氣象局烏魯木齊沙漠氣象研究所,新疆烏魯木齊830002;3.塔克拉瑪干沙漠大氣環(huán)境觀測試驗(yàn)站,新疆塔中841000;4.布爾津縣氣象局,新疆布爾津836600)

利用2013年6月巴丹吉林沙漠拐子湖地區(qū)流沙下墊面的陸氣通量觀測資料,計算并分析了該研究區(qū)不同大氣穩(wěn)定層結(jié)條件下的湍流速度各分量譜和溫度譜及湍流的局地各向同性特征.結(jié)果表明:該研究區(qū)風(fēng)速各分量的湍流強(qiáng)度隨平均風(fēng)速的增加而逐漸減小,當(dāng)平均風(fēng)速gt;2 m/s時,風(fēng)速各分量的湍流強(qiáng)度逐漸穩(wěn)定且基本限定在0.5以內(nèi).在不同大氣穩(wěn)定度下,湍流速度和溫度能譜曲線在慣性副區(qū)內(nèi)逐漸有合并趨勢且遵循Kolmogorov提出的-2/3定律,速度各分量譜在高頻段均滿足各向同性且符合低頻限制理論.隨著穩(wěn)定度的增加,風(fēng)速分量的能譜曲線逐漸降低且向高頻端移動,風(fēng)速分量和溫度能譜對應(yīng)的譜峰長度則逐漸減小.該研究區(qū)水平湍流尺度范圍為9.0～600 m,垂直湍流譜峰波長為10.79～75 m.該結(jié)果介于草地和森林下墊面之間,與塔克拉瑪干沙漠腹地的試驗(yàn)結(jié)果較為接近.

巴丹吉林沙漠;湍流;大氣穩(wěn)定度;能譜分析;局地各向同性

大氣湍流是一種高頻微尺度運(yùn)動,由許多不同大小的湍渦相互疊加而成[1,2],是大氣和地表之間的能量、物質(zhì)轉(zhuǎn)化與交換的最重要途徑,直接與大氣邊界層中的動量、熱量、水汽、污染物的垂直和水平輸送密切相關(guān)[3-5].大氣湍流強(qiáng)度用湍流能譜密度表示,可描述湍渦的無序分布,不僅可以反應(yīng)湍流總能量和不同尺度湍渦對整體湍能貢獻(xiàn)的相對大小[6],還可用于確定渦度相關(guān)系統(tǒng)對高頻和低頻率信號的響應(yīng)能力,而且對于探討近地面層物質(zhì)、能量的輸送機(jī)制也有重要意義[4,7].

20世紀(jì)40年代,Kolmogorov提出了局地均勻各向同性理論,同時首次通過量綱分析得到湍流能譜與頻率之間的-5/3指數(shù)關(guān)系[8].在此之后,Monin和Obukhov給出了湍流相似理論[9].該階段由于監(jiān)測儀器響應(yīng)較慢,不能滿足相關(guān)要素湍流脈動量的測量要求,很大程度上限制了湍流理論的進(jìn)一步發(fā)展.20世紀(jì)60年代以來,隨著大氣探測手段和電子計算能力的不斷提高,使得對于大氣湍流特征的分析越來越深入.1968年,AFCRL(Air Force Cambridge Research Laboratories)開展的Kansas試驗(yàn)是第一次大規(guī)模使用快速響應(yīng)探測儀器,給出了平坦均一下墊面風(fēng)速和溫度序列,并得到了很多標(biāo)志性的結(jié)果[10].Kaimal等利用熱線風(fēng)速計給出了湍流動能與穩(wěn)定度之間的關(guān)系,同時指出地形對湍流譜的低頻擾動可使其峰值產(chǎn)生位移[11,12].通過HEIFE和HUBEX試驗(yàn)分別得出了中國西北戈壁和淮河流域地區(qū)的湍流能譜特征[13,14].劉樹華等對長白山森林冠層上下兩個高度的湍流譜結(jié)構(gòu)、局地各向同性和耗散率進(jìn)行分析,揭示了森林粗糙表面在湍流輸送過程中的動力擾動和對大渦的破碎作用[15].此后通過EBEX-2000試驗(yàn)分析了美國加州棉田湍流速度、溫度和濕度在不同穩(wěn)定度下的特征,并與Kansas和長白山森林實(shí)驗(yàn)結(jié)果進(jìn)行比較.倪攀等針對科爾沁半干旱風(fēng)沙草原地區(qū)進(jìn)行了湍流譜和協(xié)譜的分析,發(fā)現(xiàn)半干旱風(fēng)沙草原下墊面相對森林,其湍渦在輸送過程中的破碎作用較小,物質(zhì)、能量主要以大渦形式輸送[4].溫雅婷等發(fā)現(xiàn)塔克拉瑪干沙漠腹地水平湍流尺度范圍介于森林和草地之間且不穩(wěn)定層結(jié)比穩(wěn)定層結(jié)的譜峰更偏向低頻端,峰值尺度更大[6].總體來說,針對沙漠下墊面湍流譜特征的研究相對較少,沙漠下墊面是否滿足湍流能譜與頻率之間的指數(shù)關(guān)系及渦度相關(guān)系統(tǒng)在該下墊面能否滿足對高頻和低頻率信號的響應(yīng)仍需開展進(jìn)一步分析.

拐子湖位于巴丹吉林沙漠北緣,屬內(nèi)蒙古自治區(qū)阿拉善盟額濟(jì)納旗地區(qū).這里自然條件異常惡劣,干旱、高溫、大風(fēng)、沙塵暴頻繁是這里的主要?dú)夂蛱攸c(diǎn).整體為典型的內(nèi)陸溫帶荒漠氣候,處于中國夏季風(fēng)邊緣區(qū),對全球氣候變化的響應(yīng)十分敏感.干旱少雨,蒸發(fā)量較高,年均氣溫為9.2℃,最高溫度44.8℃,最低溫度-30.7℃.年均降水量僅42.9 mm,潛在蒸發(fā)量高達(dá)4 217.9 mm[16,17].該區(qū)域風(fēng)沙災(zāi)害頻發(fā),是內(nèi)蒙古強(qiáng)和特強(qiáng)沙塵暴的主要多發(fā)區(qū)之一,同時也是我國沙塵暴北方路徑的必經(jīng)之地[18].主風(fēng)向?yàn)槲黠L(fēng),年均8級以上大風(fēng)約58次,沙塵暴30 d左右,歷史最高風(fēng)速為38 m/s.地表以大面積流沙覆蓋為主,平均粒徑為250～500 μm,與塔克拉瑪干、古爾班通古特、騰格里沙漠相比,其沙粒粒徑相對較粗.目前,針對巴丹吉林沙漠的研究多集中于高大沙山的形成和演化、湖泊-地下水關(guān)系、地表風(fēng)蝕起沙、沙漠氣候變化等問題,而關(guān)于該地區(qū)陸氣相互作用的研究,特別是陸氣湍流能譜分析卻相對較少.本文利用位于巴丹吉林沙漠北緣拐子湖通量觀測站2013年6月的湍流通量資料,分析該地區(qū)流沙下墊面不同穩(wěn)定度條件下,湍流速度各分量譜和溫度譜及湍流的局地各向同性特征.這對進(jìn)一步了解沙漠地區(qū)湍流特征及物質(zhì)、能量的輸送機(jī)制有重要意義,并且為分析和評價沙漠下墊面不同大氣穩(wěn)定度下的渦度相關(guān)實(shí)測數(shù)據(jù)質(zhì)量提供參考依據(jù).

1 研究設(shè)備及方法

1.1 試驗(yàn)設(shè)備

所用資料來自于拐子湖氣象站以南1 km處的陸氣湍流通量觀測系統(tǒng)(41°22′N,102°22′E,海拔960 m)(圖1).該觀測系統(tǒng)由中國氣象局烏魯木齊沙漠氣象研究所于2012年6月建成,整體位于平坦的流動沙地之上.主要設(shè)備包括:10 m梯度風(fēng)、溫、濕探測系統(tǒng),渦度相關(guān)探測系統(tǒng),輻射及土壤溫、濕度觀測系統(tǒng).其中,渦度相關(guān)探測系統(tǒng)安裝在距地面3 m高度處,由三維超聲風(fēng)速儀(R3-50,GILL,UK)和紅外CO2/H2O分析儀(LI-COR7500A,Campbell Scientific,USA)組成.數(shù)據(jù)采集頻率為10 Hz,存儲于美國Campbell Scientific公司生產(chǎn)的CR1000型數(shù)據(jù)采集器中.

1.2 湍流譜分析理論基礎(chǔ)

在大氣邊界層中觀測到的湍渦是時間和空間上的廣闊結(jié)構(gòu),對其理想的分析方法是采用譜分析[2].根據(jù)Taylor湍流冰凍理論,利用傅里葉變換,可將時間序列的信號轉(zhuǎn)化成空間分布,可發(fā)現(xiàn)不同時間和空間尺度的湍渦在整個湍流狀態(tài)中所起的作用,確定渦度相關(guān)系統(tǒng)對高頻湍流信號的響應(yīng)能力,以此來檢測評估湍流數(shù)據(jù)質(zhì)量的好壞.湍流能譜一般分為3個部分:含能渦區(qū)、慣性副區(qū)和耗散區(qū).Kolmogorov提出在慣性副區(qū)中風(fēng)速的湍流能譜密度表達(dá)式[8]為:

圖1 拐子湖地理位置及通量觀測站示意圖

依據(jù)泰勒假設(shè),對湍流能譜密度函數(shù)積分,將波數(shù)空間轉(zhuǎn)化為頻率空間,即得到湍流能譜與頻率的關(guān)系.湍流的熱力和動力生成主要發(fā)生在靠近譜峰的低頻區(qū)域,低頻湍渦所含能量逐漸向較高頻率的小尺度湍渦傳遞,最終完全被分子熱耗散.高頻部分存在一段慣性副區(qū).慣性副區(qū)屬于符合局地均勻和各向同性的小尺度湍流中尺度稍小的子區(qū),沒有顯著的湍流能量產(chǎn)生或耗散,且遵循-2/3定律向更小尺度傳遞,湍流能量只是通過慣性力從低波數(shù)(較大尺度湍渦)向高波數(shù)(較小尺度湍渦)輸送,并在壓力脈動的作用下在空間均勻分布[11,19].這里以風(fēng)速分量α和溫度T為例:

式中,α分別為u、v、w分量;aα為風(fēng)速分量的Kolmogorov常數(shù),Kaimal給出au=0.51,av=aw=4au/3,數(shù)值為0.5～0.6;ε為湍流耗散率;為無因次頻率,f是周期頻率,z為觀測高度為平均風(fēng)速.

1.3 湍流譜的計算和處理

將2013年6月的10 Hz原始湍流通量數(shù)據(jù)分割成每30 min為一組數(shù)據(jù).通過Eddypro 4.0軟件對分割好的數(shù)據(jù)進(jìn)行相應(yīng)的低頻濾波、各分量脈動計算、快速傅里葉變換(FFT)等處理后得到相應(yīng)各分量的功率譜.在此基礎(chǔ)上,對得到的功率譜高頻區(qū)進(jìn)行低通濾波、滑動平均及湍流譜歸一化等處理后得到最終的湍流譜,相應(yīng)處理的詳細(xì)步驟見劉樹華和溫雅婷等[2,6]分別對EBXB-2000實(shí)驗(yàn)和塔克拉瑪干沙漠腹地的湍流譜分析.

2 結(jié)果與分析

2.1 湍流強(qiáng)度分析

湍流強(qiáng)度指被水平風(fēng)速標(biāo)準(zhǔn)化后的標(biāo)準(zhǔn)差[19],定義為Iα= σα/U,其中Iα為各風(fēng)速分量(u、v、w)的湍流強(qiáng)度,σα為相應(yīng)風(fēng)速分量的標(biāo)準(zhǔn)差,U為平均風(fēng)速.從2013年6月的通量觀測資料中選取天氣系統(tǒng)較為穩(wěn)定的14-17日和19-20日作為巴丹吉林沙漠拐子湖地區(qū)湍流強(qiáng)度的分析對象,分析結(jié)果如圖2所示.從圖中可以看出,整體上風(fēng)速各分量的湍流強(qiáng)度隨平均風(fēng)速的增加而逐漸減小,其中平均風(fēng)速在1 m/s左右時,各風(fēng)速分量的湍流強(qiáng)度快速下降,當(dāng)平均風(fēng)速gt;2 m/s時,風(fēng)速各分量的湍流強(qiáng)度逐漸穩(wěn)定,基本限定在0.5以內(nèi).風(fēng)速各分量的湍流強(qiáng)度平均值依次為0.31、0.30和0.09,Iu和Iv大小相似且均大于Iw.整體上符合泰勒假說的條件.

圖2 風(fēng)速各分量的湍流強(qiáng)度與平均風(fēng)速的關(guān)系

2.2 湍流譜分析及譜峰頻率和波長

2.2.1 湍流譜分析

利用2013年6月的渦度相關(guān)通量數(shù)據(jù),從中選取天氣晴朗且湍流充分發(fā)展的6月20日作為重點(diǎn)研究對象,對其全天的大氣穩(wěn)定度((z-d)/L)進(jìn)行計算(圖3).大氣穩(wěn)定度((z-d)/L)的定義為:

其中,z為傳感器的安裝高度,此處為3.0 m;d為零平面位移,研究區(qū)為裸露的流沙下墊面,此處為0 cm;L為奧布霍夫長度;為平均溫度的特征值;u為摩擦風(fēng)速;k為卡曼常數(shù),此處取0.4;g為重力加速度;為熱量通量.當(dāng)(時,大氣處于不穩(wěn)定層結(jié);當(dāng)(z-d)/L=0時,大氣處于近中性層結(jié);當(dāng)(z-d)/Lgt;0時,大氣處于穩(wěn)定層結(jié).6月20日全天的大氣穩(wěn)定度變化如圖3所示,從中特選取11:30-12:00和12:00-12:30代表大氣層結(jié)的不穩(wěn)定狀態(tài),平均大氣穩(wěn)定度分別為-1.810和-1.402;5:00-5:30和9:00-9:30代表大氣層結(jié)處于近中性狀態(tài),平均穩(wěn)定度分別為0.011和0.010;23:00-23:30和00:00-00:30代表大氣層結(jié)處于穩(wěn)定狀態(tài),平均穩(wěn)定度分別為0.975和0.638,相應(yīng)的湍流通量數(shù)據(jù)進(jìn)行湍流能譜分析.

圖3 2013年6月20日拐子湖地區(qū)流沙下墊面大氣穩(wěn)定度日變化情況

在不同大氣穩(wěn)定度條件下,對巴丹吉林沙漠北緣拐子湖地區(qū)流沙下墊面3 m高度處的風(fēng)速分量和溫度進(jìn)行能譜分析.結(jié)果如圖4所示,其中橫坐標(biāo)為無因次頻率,縱坐標(biāo)為歸一化的湍流能譜在雙對數(shù)坐標(biāo)系下,速度能譜曲線在低頻區(qū)相對分散,在慣性副區(qū)內(nèi)逐漸有合并趨勢且歸一化湍流能譜在該區(qū)域內(nèi)大體遵循Kolmogorov提出的-2/3定律(圖4a、4b、4c).整體來看,在不同的大氣穩(wěn)定度下各風(fēng)速分量的低頻區(qū)均呈現(xiàn)能譜曲線隨大氣穩(wěn)定度逐漸增加而降低且向高頻端移動的排列趨勢.這一結(jié)果與Kansas實(shí)驗(yàn)得到的能譜曲線分布規(guī)律是一致的,但與之相比本研究中風(fēng)速能譜曲線依次排列的規(guī)律并不十分明顯.這是由于隨著穩(wěn)定度的增加,各風(fēng)速分量相應(yīng)的湍流能量逐漸減小,進(jìn)而造成能譜曲線隨之降低且湍流尺度減小,向高頻端移動.大氣不穩(wěn)定時,譜峰能量偏向低頻端,那里存在更多能量且較大尺度的湍渦,同時也顯示了沙漠地區(qū)對湍流相對于其他下墊面較小的破碎作用.此外,與溫雅婷[6]和倪攀[4]分別對塔克拉瑪干沙漠和科爾沁草地湍流能譜的研究結(jié)果類似,發(fā)現(xiàn)本研究中各風(fēng)速分量的能譜曲線在穩(wěn)定條件下,能譜的高頻端出現(xiàn)類似的上翹現(xiàn)象,但相對微弱.這說明此時的信號被噪聲干擾,需要對數(shù)據(jù)進(jìn)行相應(yīng)的濾波處理.同時,值得注意的是u譜和v譜在部分穩(wěn)定度下的低頻端具有上翹現(xiàn)象,以至于能譜曲線沒有表現(xiàn)出明顯的譜峰.這種能譜曲線的低頻端上翹現(xiàn)象在Cava等利用南極Nansen Ice Sheet觀測數(shù)據(jù)進(jìn)行湍流譜分析過程中也有發(fā)現(xiàn)[20],并且在倪攀針對科爾沁草地湍流能譜分析中也有所提及[4],但究其產(chǎn)生的原因還不是十分清楚,還有待進(jìn)一步研究.

對于溫度能譜來說,許多研究認(rèn)為,在不穩(wěn)定層結(jié)下,溫度能譜曲線在慣性副區(qū)遵從-2/3冪次率,無因次頻率大于0.3;在穩(wěn)定層結(jié)下,無因次頻率大于0.8時,能譜曲線服從-2/3定律.同時,在低頻的含能區(qū)域,溫度能譜曲線隨穩(wěn)定度參數(shù)呈現(xiàn)規(guī)則的依次變化.對照本研究中的溫度能譜,整體上其分布的規(guī)律性較弱.在不同穩(wěn)定度下,各溫度能譜曲線基本遵從-2/3冪次率.相比各風(fēng)速譜,溫度能譜在z/L為0.011和0.010兩個近中性層結(jié)下,分布整體較為分散.同時在各大氣穩(wěn)定度下,溫度能譜也沒有表現(xiàn)出如風(fēng)速能譜按照穩(wěn)定度不同依次排列的現(xiàn)象.

2.2.2 譜峰頻率和波長

圖4 拐子湖地區(qū)流沙下墊面3 m高度處風(fēng)速分量和溫度在不同大氣穩(wěn)定度下的歸一化譜

譜峰頻率是湍流譜能量最大值對應(yīng)的頻率,譜峰波長反映的是對湍流能量貢獻(xiàn)最大的渦的尺度,譜峰波長越大,對應(yīng)的湍渦尺度也越大.根據(jù)Taylor假設(shè):λmα=(z-d)/nmα,其中nmα是風(fēng)速分量(u、v、w)和溫度(T)譜對應(yīng)的譜峰頻率;λmα為譜峰波長.從表1拐子湖地區(qū)流沙下墊面風(fēng)速分量和溫度在不同大氣穩(wěn)定度下,湍流能譜的譜峰頻率和譜峰波長可見,隨著大氣穩(wěn)定度的逐漸增加,風(fēng)速分量和溫度能譜對應(yīng)的譜峰頻率均有不同程度的增加,對應(yīng)的譜峰長度則逐漸減小,其中在不穩(wěn)定層結(jié)下各分量譜峰波長達(dá)到最大值.這說明隨大氣穩(wěn)定度的增加,各分量對應(yīng)的湍渦尺度逐漸減小.同時,通過橫向?qū)Ρ瓤砂l(fā)現(xiàn)w譜峰波長整體上小于u譜和v譜.這符合一般垂直湍渦尺度相對較小的認(rèn)識.

此外,與其他下墊面湍流尺度相比,Kansas實(shí)驗(yàn)得出典型平坦草原下墊面湍渦的水平尺度范圍為22～2200 m,垂直尺度為5～200 m;劉樹華等對長白山原始森林進(jìn)行湍流觀測研究發(fā)現(xiàn),其水平尺度范圍為3～160 m,垂直尺度為1～50 m[15];EBEX-2000實(shí)驗(yàn)得出水平湍流譜峰波長范圍為4～810 m,垂直湍流譜峰波長為1.2～81 m[16];溫雅婷等對與本研究下墊面同為流動沙漠的塔克拉瑪干沙漠腹地的湍流尺度進(jìn)行分析發(fā)現(xiàn),水平湍流譜峰波長范圍為15～909 m,垂直湍流譜峰波長為5～76 m.本研究得出巴丹吉林沙漠北緣拐子湖地區(qū)水平湍流尺度范圍為9.01～600 m,垂直湍流譜峰波長為10.79～75 m.該結(jié)果明顯小于Kansas實(shí)驗(yàn)結(jié)果,大于長白山原始森林,與塔克拉瑪干沙漠和美國加州棉花下墊面結(jié)果較為接近.其中巴丹吉林沙漠拐子湖地區(qū)的水平湍流尺度的最高值小于塔克拉瑪干沙漠和美國加州棉花下墊面,垂直湍流譜峰波長的最低值高于其二者.對湍流尺度的差異比較,證明了不同下墊面對湍渦的破碎作用,即下墊面越粗糙,對湍渦能量貢獻(xiàn)最大的渦的水平尺度越小.同時,太陽輻射及地表強(qiáng)迫的不同引起了垂直湍流尺度的差異.

2.3 局地各向同性

根據(jù)Kolmogorov湍流理論,在慣性副區(qū)各風(fēng)速分量(u、v、w)能譜滿足Sv(f)/Su(f)=Sw(f)/Su(f)=4/3條件,則說明湍流符合各向同性.圖5為研究區(qū)3 m高度處不同大氣穩(wěn)定度條件下,橫向速度譜、垂直速度譜分別與縱向速度譜比值隨頻率的變化情況.整體來看,Sv(f)/Su(f)和Sw(f)/Su(f)在不同大氣穩(wěn)定條件下大體均處于4/3直線以下且低頻段曲線相對分散,高頻段集中收斂在4/3直線附近.隨著大氣穩(wěn)定度的增加,其中Sv(f)/Su(f)譜線穩(wěn)定達(dá)到各向同性條件對應(yīng)的頻率分別為0.040、0.173、0.249;對于Sw(f)/Su(f)譜線來說,對應(yīng)的頻率分別為0.084、0.619、0.622.可見在低頻段湍流譜各向同性較差,高頻段各向同性較好.這符合Kaimal提出的湍流各向同性的低頻限制理論[12].其中,Sv(f)/Su(f)的譜線隨著大氣穩(wěn)定度的增加相對更為遠(yuǎn)離4/3直線且到達(dá)4/3直線的頻率也逐漸增加.與其他下墊面的結(jié)果相比,本研究得到的各向同性數(shù)據(jù)在低頻段比Kansas實(shí)驗(yàn)、EBEX-2000實(shí)驗(yàn)及塔克拉瑪干腹地的結(jié)果分散的多.對比橫向速度譜和垂直速度譜分別與縱向速度譜的比值隨頻率的變化情況,其中Sw(f)/Su(f)譜線在低頻段更為遠(yuǎn)離4/3直線且達(dá)到4/3直線時對應(yīng)的頻率相對更大,這種現(xiàn)象在大氣穩(wěn)定和近中性條件下表現(xiàn)的要比不穩(wěn)定條件更為明顯,這說明Sw(f)/Su(f)譜線的低頻限制相對較強(qiáng),Sw(f)在低頻處能量更小.

3 結(jié)論

通過對巴丹吉林沙漠拐子湖地區(qū)流沙下墊面3 m高度的陸氣通量資料進(jìn)行分析,得到了研究區(qū)近地層湍流速度譜、溫度譜及局地各向同性特征,結(jié)果分析表明:

表1 拐子湖地區(qū)流沙下墊面3 m高度處風(fēng)速分量和溫度在不同大氣穩(wěn)定度下湍流能譜的譜峰頻率和譜峰波長

圖5 拐子湖地區(qū)流沙下墊面3 m高度處不同大氣穩(wěn)定度條件下風(fēng)速譜的湍流各向同性特征

(1)在巴丹吉林沙漠北緣拐子湖地區(qū)流動沙漠下墊面,風(fēng)速各分量的湍流強(qiáng)度隨平均風(fēng)速的增加而逐漸減小.平均風(fēng)速gt;2 m/s時,風(fēng)速各分量的湍流強(qiáng)度逐漸穩(wěn)定,基本限定在0.5以內(nèi),Iu和Iv大小相似且均大于Iw.整體上符合泰勒假設(shè).

(2)在不同大氣穩(wěn)定度條件下,速度和溫度能譜曲線在低頻區(qū)相對分散,在慣性副區(qū)內(nèi)逐漸有合并趨勢且大體遵循Kolmogorov提出的-2/3定律.風(fēng)速分量的能譜曲線隨著穩(wěn)定度增加逐漸降低且向高頻端移動,湍流能量和尺度均逐漸減小.大氣不穩(wěn)定時,譜峰能量偏向低頻端,那里存在更多能量且較大尺度的湍渦,同時也顯示了沙漠地區(qū)對湍流相對于其他下墊面較小的破碎作用.溫度譜線整體較為分散,隨穩(wěn)定度參數(shù)的變化相應(yīng)的排列規(guī)律不明顯.

(3)隨著大氣穩(wěn)定度的逐漸增加,風(fēng)速分量和溫度能譜對應(yīng)的譜峰長度則逐漸減小,即湍渦尺度逐漸減小.本研究得出研究區(qū)水平湍流尺度范圍為9.01～600 m,垂直湍流尺度范圍為10.79～75 m.該結(jié)果與其他下墊面的對比,進(jìn)一步證明了粗糙下墊面對湍渦的破碎作用.

(4)通過各向同性分析結(jié)果顯示,在不同大氣穩(wěn)定條件下,速度各分量譜在高頻段均能滿足各向同性,符合Kaimal提出的湍流各向同性的低頻限制理論.

[1] 魏偉,張宏升.希爾伯特-黃變換技術(shù)及在邊界層湍流研究中的應(yīng)用[J].氣象學(xué)報,2013,71(6):1183-1193.

[2] 劉樹華,李潔,劉和平,等.在EBEX-2000實(shí)驗(yàn)資料中湍流譜和局地各向同性特征[J].大氣科學(xué),2005,29(2):213-224.

[3] 李富余,張宏升,李詩明,等.青藏高原當(dāng)雄地區(qū)湍流脈動耗散率特征研究 [J].北京大學(xué)學(xué)報(自然科學(xué)版),2004,40(5):781-787.

[4] 倪攀,金昌杰,王安志,等.科爾沁草地不同大氣穩(wěn)定度下湍流特征譜分析[J].生態(tài)學(xué)雜志,2009,28(12):2495-2502.

[5] Dyer A J,Garratt J R,Francey R J,et al.An international turbulence comparison experiment (ITCE 1976)[J].Boundary-Layer Meteorology,1982,24(2):181-209.

[6] 溫雅婷,焦冰,繆啟龍,等.塔克拉瑪干沙漠腹地近地層湍流能譜特征分析[J].中國沙漠,2012,32(6):1716-1722.

[7] 張明達(dá),張一平,孫曉敏,等.西雙版納熱帶季節(jié)雨林林冠上方湍流譜特征[J].生態(tài)學(xué)雜志,2009,28(9):1779-1786.

[8] Kolmogorov A N.The localstructure ofyurbulence in incompressible viscous fluid for very large reynolds numbers[J].Proceedings of the Royal Society of London,1968,434(1890):9-13.

[9] Monin A S,Obukhov A M.Basic laws of turbulent mixing in the atmosphere near the ground[J].Basic Laws of Turbulent Mixing in the Ground Layer of the Atmosphere,1959,64:1963-1987.

[10] Haugen D A,Kaimal J C,Bradley E F.An experimental study of Reynolds stress and heat flux in the atmospheric surface layer[J].Quarterly Journal of the Royal Meteorological Society,2010,97(412):168-180.

[11] KaimalJ C,Wyngaard J C,IzumiY,etal.Spectral characteristics of surface-layer turbulence[J].Quarterly Journal of the Royal Meteorological Society,1972,98(417):563-589.

[12] Kaimal J C,Finnigan J J.Atmospheric boundary layer flows:theirstructure and measurement [M].New York:Oxford University Press,1994:213-214.

[13] 劉曉虎,王介民.戈壁地區(qū)湍流能量耗散率和溫度脈動耗散率[J].高原氣象,1992,11(4):344-352.

[14] 徐敏,蔣維楣,胡非,等.湍流濕度脈動特性研究--根據(jù)淮河流域?qū)嶒?yàn)資料[J].地球物理學(xué)報,2002,45(1):17-25.

[15] 劉樹華,劉和平,X Mei,等.森林冠層上下湍流譜結(jié)構(gòu)和耗散率[J].中國科學(xué):地球科學(xué),1998(5):469-480.

[16] 胡文峰,何清,楊興華,等.一次沙塵暴過程中的沙地面微氣象要素及起沙參數(shù)分析[J].干旱區(qū)資源與環(huán)境,2012,26(5):73-78.

[17] 楊帆,鄭新倩,努爾阿米娜.依明,等.巴丹吉林沙漠北緣沙塵天氣過程中近地面氣象要素變化及風(fēng)沙流結(jié)構(gòu)分析 [J].沙漠與綠洲氣象,2015,9(4):67-74.

[18] 劉景濤,鄭明倩.內(nèi)蒙古中西部和特強(qiáng)沙塵暴的氣候?qū)W特征[J].高原氣象,2003,22(2):54-62.

[19] 張翔,李云波,張學(xué)宏,等.大連海域近海面湍流結(jié)構(gòu)及譜特征[J].海洋通報,2012,31(1):9-14.

[20] Cava D,Giostra U,Tagliazucca M.Spectral Maxima In A Perturbed Stable Boundary Layer [J]. Boundary-Layer Meteorology,2001,100(3):421-437.

Characteristics of Surface-layer Turbulence Spectra under Different Atmospheric Stability in Guaizi Lake Region of Badain Jaran Desert

ZHENG Xinqian1,PENG Dongmei1,WU Ye2,3,ZHOU Chenglong2,3,CHAN Xiumei4,Ali Mamtimin2,3*
(1 Xinjiang Agro-Meteorolical Observatory,Urumqi 830002,China;2 Institute of Desert Meteorology,China Meteorological Administration,Urumqi 830002,China;3 Taklimakan Desert Atmosphere and Environment Station,Tazhong 841000,China;4 Burqin Meteorological Bureau,Burqin 836600,China)

In this paper,we analyzed the turbulence spectral characteristics of wind speed(u,v,w)and temperature and the turbulence spectral local isotropy of the Guaizi lake shifting sandy land in different atmospheric stable stratification condition by use of land-atmosphere turbulent flux data in June 2013.The result indicated that the turbulence intensity of wind speed(u,v,w)decreased with the increase of wind speed,which was gradually stabilized and less than 0.5 when average wind speedgt;2 m/s.Under different atmospheric stable stratification condition,the wind speed and temperature spectra had concentrated tendency and they followed the-2/3 law in the higher frequency.The wind speed (u,v,w)spectra fit lower frequency limit for isotropy.With the increase of atmospheric stable stratification condition,the curve of wind speed (u,v,w)spectra decreased gradually and then moved toward the region of higher frequency,and the peak wavelength of the wind speed and temperature spectra decreased gradually.In Guaizi lake shifting sandy land,the horizontal peak wavelength of the wind speed ranged from 9.01 to 600 m,and the vertical peak wavelength of the wind speed ranged from 10.79 to 75 m.The values of the peak wavelength were between the value in forest and in grassland,and they were similar to that of Taklimakan Desert.

Badain Jaran desert;turbulence;atmospheric stability;spectra analysis;local isotropy

P425.2

A

1002-0799(2017)05-0063-07

鄭新倩,彭冬梅,吳燁,等.巴丹吉林沙漠北緣拐子湖近地層湍流能譜特征分析[J].沙漠與綠洲氣象,2017,11(5):63-69.

10.12057/j.issn.1002-0799.2017.05.009

2017-03-01;

2017-04-01

國家自然科學(xué)基金(41505008、41575008、41305107).

鄭新倩(1987-),女,工程師,主要從事為農(nóng)業(yè)氣象服務(wù).

E-mail:zhengxinqian1987@126.com

買買提艾力.買買提依明(1978-),男,維吾爾族,副研究員,主要從事沙漠陸面過程研究.E-mail:ali@idm.cn