羅 漢,韓永翔*,李巖瑛(.南京信息工程大學(xué)氣象災(zāi)害預(yù)報預(yù)警與評估協(xié)同創(chuàng)新中心,中國氣象局氣溶膠-云-降水重點開放實驗室,江蘇 南京 0044；.甘肅省武威市氣象局,甘肅 武威 733000)

邊界層高度時間演變及塵卷風(fēng)對總沙塵量的貢獻

羅 漢1,韓永翔1*,李巖瑛2(1.南京信息工程大學(xué)氣象災(zāi)害預(yù)報預(yù)警與評估協(xié)同創(chuàng)新中心,中國氣象局氣溶膠-云-降水重點開放實驗室,江蘇 南京 210044；2.甘肅省武威市氣象局,甘肅 武威 733000)

采用干絕熱曲線法,計算了敦煌地區(qū)每日最大對流邊界層高度(簡稱:CBLmax)并分析了可能的影響因子,在此基礎(chǔ)上計算出塵卷風(fēng)對大氣年沙塵氣溶膠的貢獻.結(jié)果表明: CBLmax具有非常顯著的年變化特征,呈單峰分布,12月最低,5月最高.年平均高度為2.2km,極端時出現(xiàn)接近6km.熱力因素對CBLmax的貢獻具有決定性,云量的多寡可影響CBLmax的變化.塵卷風(fēng)的起沙量對大氣年均沙塵氣溶膠總量的貢獻至少在54.4%以上.

對流邊界層厚度；塵卷風(fēng)；起沙量；貢獻率

沙塵氣溶膠是對流層大氣氣溶膠的重要組成成分,約占全球自然氣溶膠總量的三分之一[1].它通過起沙-傳輸-沉降成為了全球氣候變化的關(guān)鍵影響因素之一[2-5].沙塵暴和揚沙被認為是沙塵氣溶膠的最主要來源[6],但是衛(wèi)星觀測的沙塵氣溶膠含量的時間變化與沙塵暴發(fā)生的時間并不完全匹配[7],研究認為存在塵卷風(fēng)的起沙機制,其在供給大氣沙塵氣溶膠的總量上可能扮演著重要的角色[8].塵卷風(fēng)是一種旋轉(zhuǎn)上升的對流渦,能夠卷起并攜帶地面沙塵上升到對流層中上部,沙塵羽是沒有形成標準渦旋的塵卷風(fēng),本文將塵卷風(fēng)和沙塵羽統(tǒng)稱為塵卷風(fēng).研究表明,塵卷風(fēng)貢獻了大約 35%的全球沙塵氣溶膠[9-10],在中國塔克拉瑪干沙漠中其貢獻達到了53%[11].

Renno等[12]將塵卷風(fēng)看做一個熱力發(fā)動機,提出了其熱力學(xué)理論,塵卷風(fēng)的強度取決于其熱力學(xué)效率,而熱力學(xué)效率與地表溫度和對流邊界層高度關(guān)系密切.利用塵卷風(fēng)的熱力學(xué)理論計算公式[10]得到塔克拉瑪干沙漠的塵卷風(fēng)起沙量[11],但其計算的對流邊界層高度是平均值,缺乏年際間的變化,同時對其他影響對流邊界層高度的因子如云量、降水等沒有考慮,日照長度并不是實際日照時數(shù),因此對塵卷風(fēng)起沙量的計算帶來較大的誤差.

1 研究區(qū)域和研究方法

1.1 研究區(qū)域與資料來源

敦煌地區(qū)總面積為 3.12萬 km2,地處東經(jīng)94°41′,北緯 40°09′,海拔 1140m,位于中國西北干旱區(qū)腹地.該區(qū)的地理景觀為沙漠和戈壁,大部分地方植被覆蓋不足 10%.日照充足,干旱少雨,年降水量在40mm左右,而蒸發(fā)潛力高達3400mm,是沙塵暴和塵卷風(fēng)的多發(fā)區(qū)之一. 2006～2015年地面資料、探空資料以及揚沙、沙塵暴等數(shù)據(jù)來自敦煌國家基準站,氣溶膠指數(shù)來自 http: //www.nasa.gov/.

1.2 計算CBLmax

干絕熱曲線法又稱 T-LogP圖法,能夠較為準確地計算出CBLmax[13-14].其適用于有探空資料的地區(qū)[15],本文通過2006～2015年10a每日08:00的探空資料,利用T-LogP圖法計算得到CBLmax.具體計算方法見文獻[15].

1.3 計算塵卷風(fēng)的起沙量

Rennó等提出了塵卷風(fēng)的熱力學(xué)理論[12],塵卷風(fēng)的強度取決于熱力學(xué)效率η[10].計算式:

式中: ZCBL為每日平均對流邊界層高度,單位為m; Th為地表氣溫,單位為 K; Γad為絕熱遞減率(Γad=10K/km).塵卷風(fēng)的覆蓋區(qū)域面積比為 σ,計算式:

我國公共服務(wù)的投入總量大，但國情決定了這仍難以滿足公眾不斷增加的公共服務(wù)需求，與投入緊密相連的財政機制還沒有制度化，缺乏法律制約，常常出現(xiàn)財政收支失衡、結(jié)構(gòu)紊亂、責(zé)任混亂等現(xiàn)象，完善的財政制度是公共服務(wù)提供的堅實基礎(chǔ)。

式中:無量綱機械能摩擦損耗系數(shù) μ≈12～24;地面到對流邊界層頂?shù)膲簭姴?Δp=gρairZCBL;空氣密度ρair=1kg/m3;重力加速度 g=9.8g/m2;驅(qū)動塵卷風(fēng)的熱量通量 Fin≈11±5kW/m2;對流邊界層有效太陽輻射時間尺度 TR≈9×105s[10].塵卷風(fēng)總起沙量DAEtot的計算式:

式中：Fd為塵卷風(fēng)的起沙通量,g/m2s; Dtime為日照時長,s; S為能夠容易揚起松散顆粒物進入大氣的土壤面積或區(qū)域[16-18],m2.

2 研究結(jié)果

2.1 CBLmax的變化

利用T-LogP圖法計算的2006～2015年10a CBLmax(圖 1)顯示,其呈現(xiàn)出非常顯著的年變化特征,在年內(nèi)呈現(xiàn)峰谷交替,12月到次年1月份為波谷, 5月左右到達波峰.年平均CBLmax變化雖然不大(10a平均2200m),但10a總趨勢緩慢降低,平均逐年降低0.01km.

圖1 CBLmax和14:00溫度變化Fig.1 The daily variation of CBLmaxand air temperature at 14:00

每月平均的 CBLmax變化(圖 2)顯示其呈現(xiàn)出單峰分布,從1月的700m迅速增加到5月的3500m左右,達到峰值,然后緩慢降低, 12月到達最低值600m左右.10年中CBLmax出現(xiàn)3500m以上的天數(shù)達到770d,極端時甚至出現(xiàn)接近6000m高度,這與其他沙漠地區(qū)短期加密觀測的結(jié)果一致[19],出現(xiàn)3500m以上深厚的對流邊界層的原因可能與北半球劇烈的太陽輻射等氣候背景和極端干燥的地表環(huán)境有關(guān)[20].

2.2 CBLmax高度與氣象要素的關(guān)系

2.2.1 CBLmax與14:00的溫度和14:00風(fēng)速的關(guān)系 地表強烈加熱引起的熱力不穩(wěn)定以及熱力湍流是影響大氣邊界層高度和對流邊界層形成的主要熱力原因[21],風(fēng)速則是其動力因素.熱力不穩(wěn)定以及熱力湍流同大氣溫度密切相關(guān).敦煌地區(qū)的CBLmax在14:00最大,因此選用14:00的溫度和 14:00風(fēng)速對 CBLmax的影響進行研究.2006～2015年CBLmax和14:00的氣溫曲線(圖2)顯示CBLmax與14:00溫度具有高度的相關(guān)性,相關(guān)系數(shù)達到0.73,通過了99.999%置信度檢驗,這說明溫度對CBLmax的貢獻具有決定性.但是它們峰值出現(xiàn)的時間有所不同, CBLmax提前溫度2個月到達峰值(圖2).同時,圖2顯示風(fēng)速在4月達到年中最大,早于 CBLmax峰值 1個月,暗示了CBLmax同時受到溫度和風(fēng)速的共同影響.

圖2 CBLmax與14:00時溫度和風(fēng)速的月平均變化Fig.2 The mean monthly variation of CBLmaxwith air temperature and wind speed at 14:00

2.2.2 CBLmax與日均云量的關(guān)系 敦煌地區(qū)年平均降水量僅為 40mm左右,降水過程雖然對CBLmax高度有影響,但太少的降水不足以影響總的趨勢,而日均云量的多少可以影響地面所接受的太陽輻射而影響湍流的發(fā)展,進而影響CBLmax的發(fā)展.

將2006～2015年晴天、多云和陰天的CBLmax進行對比(圖 3),結(jié)果發(fā)現(xiàn)晴天、多云和陰天的CBLmax變化趨勢基本一致,但是晴天 CBLmax高于多云天,而多云天也明顯高于陰天.晴天和多云天6月的CBLmax為年中最高,而陰天5月最高,說明陰天可能也是影響CBLmax提前溫度到達年內(nèi)峰值的因素.將晴天與陰天的CBLmax進行對比,其中夏季8月高度差最大,大于1400m,陰天高度減少率為43.8%; 9月的陰天高度相比晴天少了1300m,高度減少率為年內(nèi)最高47.1%;冬季12月高度差最小,其值小于 100m,陰天高度減少率為年內(nèi)最低為2%;陰天相對晴天CBLmax高度平均減少率為 27.1%.熱力因素如溫度在干旱區(qū)是影響湍流強弱最主要的因素,地表接受的太陽輻射晴天最多、多云天次之、而陰天最少,所以晴天CBLmax高度高于多云天,而多云天高于陰天.因此,云量的多寡可影響CBLmax的變化.

圖3 2006～2015年晴天、多云和陰天CBLmax對比Fig.3 A comparison of CBLmaxon sunny, cloudy and overcast days from 2006 to 2015

2.3 估算塵卷風(fēng)的起沙量

根據(jù)敦煌探空氣球進行的邊界層高度的短期加密觀測,日內(nèi)對流邊界層高度的變化先增加后減小,基本呈現(xiàn)出正態(tài)分布,日落后其高度雖有降低,但不會降低至 0m,而是轉(zhuǎn)為邊界層剩余層,高度約占CBLmax的38.6%[22],因此每日平均的對流邊界層高度與 CBLmax之間的關(guān)系為:日平均對流邊界層高度=CBLmax×61.4%.將地表氣溫、每日平均對流邊界層高度代入式(1)、式(2),可分別獲得熱力學(xué)效率η與塵卷風(fēng)的覆蓋區(qū)域面積比σ.敦煌地區(qū)的總面積為3.12萬km2,能夠容易揚起松散顆粒物進入大氣的土壤面積或區(qū)域 S為1.97×1010m2[23].塵卷風(fēng)和沙塵羽的起沙通量 Fd分別為(0.7±0.3)g/(m2·s)和(0.1±0.03)g/(m2·s)并都可利用式(3)計算各自起沙量[10],由于沙塵羽是塵卷風(fēng)的一種特殊形態(tài),但起沙通量有較大差異,因此本文分開計算兩者的起沙量,但在討論中合并兩者的起沙量并統(tǒng)稱為塵卷風(fēng)的起沙量. Dtime取每日地面觀測資料中的實際日照時間,將η、σ、S、Fd和 Dtime代入式(3),計算得到敦煌地區(qū) 2006～2015年塵卷風(fēng)的每日起沙量.

在有沙塵暴、揚沙、降水、陰天等天氣環(huán)境下很難發(fā)生塵卷風(fēng),因此本文將沙塵暴和揚沙發(fā)生日、日降水大于4mm和日平均總云量大于70%以上的數(shù)據(jù)剔除,獲得了訂正以后的塵卷風(fēng)的每日起沙量.沙漠和戈壁地區(qū)大氣中的沙塵氣溶膠主要來自沙塵暴、揚沙和塵卷風(fēng)的貢獻[11],那么剔除沙塵暴和揚沙后, Aura衛(wèi)星觀測的敦煌氣溶膠指數(shù)AI則反映的是塵卷風(fēng)的貢獻.剔除沙塵暴、揚沙、降水、陰天等天氣后,圖4顯示了每日塵卷風(fēng)的起沙量與Aura-AI之間的變化(共2314d),它們均呈典型的年單峰變化,年初與年末兩者同時到達年內(nèi)低值,年中兩者同時到達年內(nèi)最高,兩者的趨勢具有高度的一致性,相關(guān)系數(shù)為 0.38,并通過了99.999%置信度檢驗,這表明本文用熱力學(xué)公式[11]計算的塵卷風(fēng)的起沙量是可信的.

圖4 敦煌地區(qū)每日塵卷風(fēng)的起沙量與Aura-AI變化(剔除沙塵暴、揚沙、降水和陰天等天氣現(xiàn)象)Fig.4 Daily dust emission of dust devils and Aura-AI change in Dunhuang from 2006 to 2015 (excluding the weather phenomena of dust storm, blowing sand,precipitation and overcast sky)

表1 計算的塵卷風(fēng)各月平均起沙量Table 1 The mean of monthly dust emissions by dust devils

計算出敦煌地區(qū)塵卷風(fēng)各月平均的起沙量在年內(nèi)呈單峰分布(表 1),從 1月開始逐漸增大,至7月達到年內(nèi)最大值,而后逐漸減少.7月平均每天塵卷風(fēng)的起沙量是1月的8倍.塵卷風(fēng)年均起沙量(最大、平均、最小)為4.28×106、1.85×106、和 7.38×105t.塵卷風(fēng)各季節(jié)起沙量分別為:春天5.71×105、夏天7.22×105、秋天4.21×105和冬天1.38×105t,夏天的起沙量是冬天的5.23倍.

2.4 塵卷風(fēng)的起沙量對大氣沙塵氣溶膠的貢獻

沙塵暴和揚沙能夠抬升大量沙塵氣溶膠粒子進入大氣[8].為了估計其潛在的起沙量,需要知道沙塵暴和揚沙的各自垂直起沙通量和其持續(xù)時間.沙塵暴的起沙通量在 40.07×10-7kg/ (m2·s)至 1.54×10-8kg/(m2·s)之間[24-26],揚沙的起沙通量在 4.14×10-8kg/(m2·s)至 0.99×10-8kg/ (m2·s)之間[25-26].敦煌地區(qū)在 2006～2015年期間共發(fā)生35次共119.1h沙塵暴和242次共698.9h揚沙天氣,因此平均每年沙塵暴和揚沙的持續(xù)時間分別為11.91和69.89h.假設(shè)沙塵暴和揚沙可以橫掃整個敦煌地區(qū)內(nèi)能夠容易揚起松散顆粒物進入大氣的土壤面積為 1.97×1010m2[24],則根據(jù)之前所提的沙塵暴和揚沙起沙通量和持續(xù)時間,計算其年均起沙量.沙塵暴的年均最大最小起沙量分別為 3.3845×106和 1.3007×104t;揚沙的年均最大最小起沙量分別為2.0520×105和4.9070×104t.

由于沙漠戈壁區(qū)大氣中沙塵氣溶膠主要來自沙塵暴、揚沙和塵卷風(fēng)的貢獻[12],因此敦煌地區(qū)的年均總起沙量為:1)最大:7.8716×106t(塵卷風(fēng):4.2819×106t,沙塵暴和揚沙:3.5897×106t);2)最小:8.0069×105t(塵卷風(fēng):7.3861×105t,沙塵暴和揚沙:6.2077×104t);平均:3.1283×106t(塵卷風(fēng): 1.8518×106t,沙塵暴和揚沙:1.2765×106t).因此,塵卷風(fēng)年均起沙對于總起沙量的貢獻至少為54.4%,這與Han等[11]估計的塔克拉瑪干沙漠中塵卷風(fēng)所占總起沙量比率為 52.8%的結(jié)果基本一致.

由于目前對塵卷風(fēng)的了解十分有限,雖然它對年起沙總量的貢獻很大,但它對區(qū)域或全球氣候變化的貢獻目前仍難以估計.

3 結(jié)論

3.1 敦煌地區(qū) CBLmax具有非常顯著的年變化特征,呈單峰分布,12月最低,5月最高.年平均高度為2.2km,極端時出現(xiàn)接近6km的超厚對流邊界層.

3.2 熱力因素對 CBLmax的貢獻具有決定性.陰天相對晴天的CBLmax平均減少率為27.1%,云量的多寡可影響CBLmax的變化

3.3 敦煌干旱區(qū)塵卷風(fēng)的起沙量對于年平均總起沙量的貢獻至少為 54.4%,最高可能達到 90%以上.

[1] 欒兆鵬,趙天良,韓永翔,等.干旱半干旱地區(qū)塵卷風(fēng)研究進展[J]. 沙漠與綠洲氣象, 2016,10(2):1-8.

[2] 沈凡卉,王體健,莊炳亮,等.中國沙塵氣溶膠的間接輻射強迫與氣候效應(yīng) [J]. 中國環(huán)境科學(xué), 2011,31(7):1057-1063.

[3] Bishop J K B, Sherman J T. Robotic Observations of Dust Storm Enhancement of Carbon Biomass in the North Pacific [J]. Science, 2002,298(5594):817-21.

[4] 韓永翔,宋連春,趙天良,等.北太平洋地區(qū)沙塵沉降與海洋生物興衰的關(guān)系 [J]. 中國環(huán)境科學(xué), 2006,26(2):157-160.

[5] 王民俊,韓永翔,鄧祖琴,等.全球主要沙源區(qū)沙塵氣溶膠與太陽輻射的關(guān)系 [J]. 中國環(huán)境科學(xué), 2012,32(4):577-583.

[6] Overpeck J, Lacis A, Healy R, et al. Possible role of dust-induced regional warming in abrupt climate change during the last glacial period [J]. Nature, 1996,384(6608):447-449.

[7] 鄧祖琴,韓永翔,白虎志,等.中國北方沙漠戈壁區(qū)沙塵氣溶膠與太陽輻射的關(guān)系 [J]. 中國環(huán)境科學(xué), 2011,31(11):1761-1767.

[8] Han Y, Dai X, Fang X, et al. Dust aerosols: A possible accelerant for an increasingly arid climate in North China [J]. Journal of Arid Environments, 2008,72(8):1476–1489.

[9] 段佳鵬,韓永翔,趙天良,等.塵卷風(fēng)對沙塵氣溶膠的貢獻及其與太陽輻射的關(guān)系 [J]. 中國環(huán)境科學(xué), 2013,33(1):43-48.

[10] Koch J, Renno N O. The role of convective plumes and vortices on the global aerosol budget [J]. Geophysical Research Letters, 2005,32(18):109-127.

[11] Han Y, Wang K, Liu F, et al. The contribution of dust devils and dusty plumes to the aerosol budget in western China [J]. Atmospheric Environment, 2016,126:21-27.

[12] Rennó N O, Burkett M L, Larkin M P. A Simple Thermodynamical Theory for Dust Devils. [J]. Journal of the Atmospheric Sciences, 1998,55(21):3244-3252.

[13] 牛生杰,呂晶晶,岳 平.半干旱荒漠化草原春季邊界層特征的一次綜合探測 [J]. 中國沙漠, 2013,33(6):1858-1865.

[14] 李巖瑛,張 強,胡興才,等.西北干旱區(qū)和黃土高原大氣邊界層特征對比及其對氣候干濕變化的響應(yīng) [J]. 冰川凍土, 2012,34(5):1047-1058.

[15] 廖國蓮.大氣混合層厚度的計算方法及影響因子 [J]. 中山大學(xué)研究生學(xué)刊(自然科學(xué),醫(yī)學(xué)版), 2005,26(4):66-73.

[16] Sinclair P C. General characteristics of dust devils [J]. Journal of Applied Meteorology, 1969,8(1):32-45.

[17] Balme M, Greeley R. Dust devils on Earth and Mars [J]. Reviews of Geophysics, 2006,225(1/2):3-11.

[18] Oke A M C, Tapper N J, Dunkerley D. Willy-willies in the Australian landscape: The role of key meteorological variables and surface conditions in defining frequency and spatial characteristics [J]. Journal of Arid Environments, 2007,71(2): 201-215.

[19] Wang M, Wei W, He Q, et al. Summer atmospheric boundary layer structure in the hinterland of Taklimakan Desert, China [J]. Journal of Arid Land, 2016,8(6):846-860.

[20] 李巖瑛,錢正安,薛新玲,等.西北干旱區(qū)夏半年深厚的混合層與干旱氣候形成 [J]. 高原氣象, 2009,28(1):46-54.

[21] 張 杰,張 強,唐從國.極端干旱區(qū)大氣邊界層厚度時間演變及其與地表能量平衡的關(guān)系 [J]. 生態(tài)學(xué)報, 2013,33(8):2545-2555.

[22] 張 強,趙映東,王 勝,等.極端干旱荒漠區(qū)典型晴天大氣熱力邊界層結(jié)構(gòu)分析 [J]. 地球科學(xué)進展, 2007,22(11):1150-1159.

[23] 尚立照,陳翔舜,王小軍,等.基于RS和GIS的敦煌市沙漠化動態(tài)監(jiān)測 [J]. 水土保持通報, 2016,36(2):125-128.

[24] 楊興華,何 清,艾力·買買提明.塔中地區(qū)一次沙塵暴過程的輸沙通量估算 [J]. 干旱區(qū)研究, 2010,27(6):969-974.

[25] 沈建國,孫照渤,章秋英,等.干旱草原地區(qū)起沙通量的初步研究[J]. 中國沙漠, 2008,28(6):1045-1049.

[26] 沈志寶,申彥波,杜明遠,等.沙塵暴期間戈壁沙地起沙率的觀測結(jié)果 [J]. 高原氣象, 2003,22(6):545-550.

致謝：本文英文部分由南京信息工程大學(xué)趙天良教授幫忙潤色,在此表示衷心的感謝.

Temporal evolution of the boundary layer height and contribution of dust devils to dust aerosols.

L
UO Han1, HAN Yong-xiang1*, LI Yan-ying2(1.Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters, Key Laboratory for Aerosol-Cloud-Precipitation of China Meteorological Administration, Nanjing University of Information Science and Technology, Nanjing 210044, China；2.Wuwei Meteorological Bureau in Gansu Province, Wuwei 733000, China). China Environmental Science, 2017,37(7)：2438～2442

By using the method of dry adiabatic curve, the daily maximum convective boundary layer height (CBLmax) in Dunhuang was calculated and the influence factors were also analyzed. Then the dust devil contribution to dust aerosols were estimated. The results showed that the CBLmaxhad obviously the annual variations with the single-peak distribution from a minimum in December to a maximum in May. The annual average of CBLmaxwas 2.2km with an extreme height of nearly 6km. Thermal factors played a dominant role in the CBLmaxvariation, and cloud amount could affect the variation of CBLmax. More than 54.4% amount of annual dust aerosols was contributed by dust devils.

convective boundary layer thickness；dust devils；dust emission amount；contribution rate

X513

A

1000-6923(2017)07-2438-05

羅 漢(1988-),男,甘肅敦煌人,碩士,主要從事大氣物理學(xué)及大氣環(huán)境研究.發(fā)表論文1篇.

2016-12-02

國家自然科學(xué)基金資助項目(41375158)

* 責(zé)任作者, 教授, han-yx66@126.com