霍 文，何 清，楊興華，艾力·買(mǎi)買(mǎi)提明，黃有志

（1.中國(guó)氣象局 烏魯木齊沙漠氣象研究所，烏魯木齊830002；2.塔克拉瑪干沙漠大氣環(huán)境觀測(cè)試驗(yàn)站，新疆 塔中841000；3.布爾津縣氣象局，新疆 布爾津836600）

眾所周知，沙漠作為一種特殊的下墊面形式，通過(guò)大氣動(dòng)力的作用以及不同沙漠的局地微氣候或者局地天氣孕育了沙塵暴、濃浮塵、揚(yáng)沙等災(zāi)害性天氣，不斷影響著周邊以及下游地區(qū)。事實(shí)上，沙漠氣候和沙漠天氣對(duì)非沙漠地區(qū)以及全球的的綜合影響是十分顯著的。關(guān)于沙漠的研究國(guó)外有很多的先驅(qū)者，Aufrere于1928—1934年間利用當(dāng)時(shí)可以收集到的氣象資料，試圖建立沙丘形態(tài)、走向與風(fēng)況的關(guān)系，這是有關(guān)沙漠氣象研究由定性向定量轉(zhuǎn)變的初步嘗試，雖然并不完善，卻是重要的研究思想突破［1］。1929—1930年，Bagnold在利比亞沙漠中進(jìn)行了人類(lèi)最早的汽車(chē)沙漠考察，并從沙粒的運(yùn)動(dòng)特征著手進(jìn)行研究［2］，他于1935年在英國(guó)劍橋大學(xué)建立了第一個(gè)用于風(fēng)沙研究的風(fēng)洞，并創(chuàng)造性地應(yīng)用和推廣了流體力學(xué)原理，建立了風(fēng)沙物理學(xué)，并著《風(fēng)沙和荒漠沙丘物理學(xué)》（The Physics of Blown Sandand Desert Dunes）一書(shū)。隨著國(guó)內(nèi)外沙塵暴的廣泛研究以及預(yù)報(bào)業(yè)務(wù)的客觀需要，更多的預(yù)報(bào)模型得以開(kāi)發(fā)，如前蘇聯(lián)科學(xué)家提出的波查羅夫（Bocharov）模型，德克薩斯理工大學(xué)提出的侵蝕分析模型（TEAM），澳大利亞科學(xué)家提出的風(fēng)蝕評(píng)價(jià)模型（WEAM）、美國(guó)農(nóng)業(yè)部的修正風(fēng)蝕方程（RWEQ）和風(fēng)蝕預(yù)報(bào)系統(tǒng)（WEPS）等。關(guān)于風(fēng)蝕模擬的代表成果集中反映在《風(fēng)蝕物理學(xué)及其模擬》（Physics and Modeling of Wind Erosion）一書(shū)中［3］。沙塵災(zāi)害預(yù)報(bào)質(zhì)量和準(zhǔn)確率的提高在于模式中起沙計(jì)算模塊的不斷改進(jìn)，即物理化起沙參數(shù)的不斷優(yōu)化，起沙參數(shù)在不同的沙漠區(qū)域各異，需要大量的野外試驗(yàn)與理論相結(jié)合不斷的更新和量化確定。目前，仍有許多學(xué)者致力于風(fēng)蝕、地表起沙參數(shù)等方面的深入研究［4－15］，但是地表起沙是一個(gè)極為復(fù)雜的物理過(guò)程，加之世界沙漠分布廣泛，沙漠地區(qū)自然條件苛刻，不利于科學(xué)家深入研究，很多的研究局限于沙漠的外圍或者個(gè)別沙漠腹地，因此地表起沙物理特性的研究還任重而道遠(yuǎn)。

同時(shí)，沙害問(wèn)題依然是人們和全社會(huì)關(guān)注的重點(diǎn)，研究沙塵災(zāi)害的環(huán)境意義和社會(huì)效益不再贅述，實(shí)際上災(zāi)害性沙塵暴天氣的預(yù)報(bào)工作才是氣象科研工作者的重中之重，還需要有關(guān)精密的、確切的資料。沙塵天氣過(guò)程中沙塵通量的準(zhǔn)確計(jì)算是一個(gè)極為關(guān)鍵的地表起沙參數(shù)問(wèn)題，目前的很多沙塵天氣數(shù)值預(yù)報(bào)模式及起沙模式的計(jì)算值與實(shí)際值存在很大的差異，有很多空?qǐng)?bào)、漏報(bào)情況。需要許多學(xué)者為數(shù)值模式的改進(jìn)提供較為準(zhǔn)確的代表性區(qū)域的參數(shù)計(jì)算方法或者實(shí)測(cè)值及計(jì)算值。本文基于此目的，結(jié)合實(shí)驗(yàn)觀測(cè)給出塔里木盆地多種沙源類(lèi)型沙通量的計(jì)算值，以期為相關(guān)數(shù)值模式預(yù)報(bào)參數(shù)化方案的改進(jìn)、為沙塵天氣預(yù)報(bào)準(zhǔn)確率的提高盡微薄之力。

1 資料與方法

1.1 觀測(cè)場(chǎng)概況

在沙漠地區(qū)研究沙塵天氣及近地層沙塵物質(zhì)的運(yùn)動(dòng)特征、物理過(guò)程、起沙參數(shù)需要長(zhǎng)期的、實(shí)時(shí)的觀測(cè)數(shù)據(jù)，基礎(chǔ)性科研成果的獲取必須基于長(zhǎng)期系統(tǒng)的基礎(chǔ)性觀測(cè)研究。中國(guó)氣象局烏魯木齊沙漠氣象研究所先后于2008年7月在塔中沙漠腹地、2008年9月在肖塘沙漠北緣過(guò)渡帶、2009年6月在若羌綠洲邊緣沙地、2010年5月在策勒綠洲荒漠過(guò)渡帶、塔里木盆地不同沙塵源地建立了4個(gè)風(fēng)沙觀測(cè)試驗(yàn)場(chǎng)，獲取了寶貴的塔里木盆地多種沙源類(lèi)型的實(shí)時(shí)觀測(cè)數(shù)據(jù)，用于科研基礎(chǔ)研究。其中，塔中風(fēng)沙觀測(cè)場(chǎng)（38°58′N(xiāo)，83°39′E，海拔1 082m）深入塔克拉瑪干沙漠腹地220km，該地區(qū)地表為流沙覆蓋，植被覆蓋度極低。地貌以沙丘地貌為主，沙丘地貌表現(xiàn)為一系列線狀的高大復(fù)合型縱向沙壟與壟間地相間分布，沙壟走向?yàn)镹NE—SSW或EN—SW方向，相對(duì)高度為40～50m。壟間平坦沙地寬1～3km，長(zhǎng)2～5km。高大沙壟的前緣分布有低矮的新月形沙丘和沙丘鏈。觀測(cè)場(chǎng)設(shè)在兩大沙壟間的平坦沙地上。該地區(qū)地表沙源的平均粒徑為136μm［16］，且分選情況較好。肖塘風(fēng)沙觀測(cè)場(chǎng)（40°49′N(xiāo)，84°17′E，海拔1 099m）地處塔克拉瑪干沙漠北緣荒漠過(guò)渡帶上，下墊面為平坦的風(fēng)蝕裸露古河床。往北2km為胡楊林，西距沙漠公路約1km，周?chē)行略滦紊城鸷蛷?fù)合型新月形沙丘鏈，一般寬6～15m，高1～10m。復(fù)合型沙丘單體走向?yàn)镹W—SE，其間廣泛分布有土平地、風(fēng)蝕殘丘等地貌類(lèi)型。該區(qū)域地表沙源平均粒徑為159 μm［16］，由于地處沙漠邊緣過(guò)渡帶和下伏河床沉積土，沙樣分選較差。若羌風(fēng)沙觀測(cè)場(chǎng)位于若羌綠洲邊緣的沙地上（39°24′N(xiāo)，88°15′E，海拔798m），觀測(cè)場(chǎng)周邊地勢(shì)平坦開(kāi)闊，地表裸露。地表沙粒較細(xì)且均勻，平均粒徑僅103μm［16］，分選比較好。策勒風(fēng)沙觀測(cè)場(chǎng)（37°02′37″N，80°40′52″E，海拔1 328m），地勢(shì)平坦開(kāi)闊，地表狀況為自然平沙地，周邊自然分布灌叢沙包，平均粒徑僅115μm［16］，分選比較好。

1.2 儀器與數(shù)據(jù)獲取

塔中、肖塘、策勒和若羌風(fēng)向、風(fēng)速數(shù)據(jù)來(lái)自MetOne公司生產(chǎn)的010C／020C風(fēng)速／風(fēng)向傳感器。觀測(cè)場(chǎng)實(shí)地安裝高度為5，10，20，50，100，200cm，風(fēng)向安裝高度為20，200cm，在200cm高度安裝有空氣溫濕度傳感器，用于監(jiān)測(cè)沙塵天氣過(guò)程中貼地層風(fēng)速廓線特征。H11—LIN型風(fēng)蝕傳感器由美國(guó)Sensit公司生產(chǎn)，可以更加精細(xì)化地監(jiān)測(cè)、研究地表沙?；顒?dòng)。觀測(cè)場(chǎng)實(shí)地傳感器的安裝高度為5，10cm，可測(cè)量沙塵天氣過(guò)程的沙塵撞擊顆粒數(shù)、撞擊動(dòng)能及沙粒運(yùn)動(dòng)速度，風(fēng)蝕傳感器上端有一個(gè)環(huán)狀壓電式晶片（直徑為25mm，高為13mm），當(dāng)沙粒撞擊在晶片上時(shí)，會(huì)輸出脈沖信號(hào)給數(shù)據(jù)采集器，記錄下撞擊顆粒的數(shù)目，并記為一次躍移活動(dòng)，數(shù)據(jù)采集頻率設(shè)為1 Hz。經(jīng)過(guò)實(shí)驗(yàn)室標(biāo)定Sensit對(duì)直徑大于50μm沙粒的撞擊具有很好的反應(yīng)。Sensit傳感器經(jīng)過(guò)長(zhǎng)期野外試驗(yàn)驗(yàn)證，對(duì)地表沙粒躍移運(yùn)動(dòng)的測(cè)量效果良好，結(jié)合貼地層風(fēng)速梯度儀，可求出沙塵天氣過(guò)程的起沙臨界摩擦速度。塔中、若羌、肖塘和策勒平均粒徑的計(jì)算使用英國(guó)Malvern公司生產(chǎn)的Mastersizer 2000型激光粒度儀進(jìn)行測(cè)量，粒度參數(shù)計(jì)算方法如下：

粒度參數(shù)指從累積曲線上求出來(lái)的能表示樣品粒度分布特征的數(shù)值，例如橫坐標(biāo)為粒徑大小，縱坐標(biāo)為粒度頻率累積分布即為累積百分含量。樣品累積曲線中的5%處的粒徑大小并換算為φ值即可用φ5表示，其他以此類(lèi)推。

由于各個(gè)風(fēng)沙觀測(cè)場(chǎng)地域以及儀器損壞等原因，本文以實(shí)際觀測(cè)已有的資料為基礎(chǔ)，通過(guò)撞擊顆粒數(shù)連續(xù)性強(qiáng)，數(shù)據(jù)穩(wěn)定等特點(diǎn)分別遴選出4個(gè)風(fēng)沙觀測(cè)場(chǎng)代表性沙塵暴過(guò)程，塔中（2009年5月26日）、肖塘（2009年7月17日）、若羌（2010年5月8日）、策勒（2010年8月24日）的一個(gè)對(duì)時(shí)數(shù)據(jù)作為樣本分析。本文考慮到由于地域和天氣的差異，下墊面條件的不同，做垂直沙通量與水平沙通量的橫向比較的可比性不強(qiáng)，因此給出風(fēng)速參考標(biāo)準(zhǔn)，可以界定在貼地面同等高度，同等風(fēng)速條件下各個(gè)實(shí)驗(yàn)區(qū)域沙塵通量方向上的差別。

1.3 沙通量計(jì)算方法與臨界摩擦速度

沙塵通量可以分為水平沙塵通量，即輸沙率與垂直輸沙通量，分別用Q與F表示，Q的物理意義是與水平風(fēng)向垂直從地表到觀測(cè)高度范圍內(nèi)單位寬度單位時(shí)間通過(guò)的沙粒質(zhì)量，Q是沙粒漂移強(qiáng)度q（z）沿高度的積分［17］；F的物理意義是某一高度處與地表面平行的單位面積內(nèi)單位時(shí)間地表向大氣中輸送的塵粒質(zhì)量［15］。關(guān)于沙通量的量綱國(guó)內(nèi)外研究尚未統(tǒng)一，本文為了方便，均采用國(guó)際通用單位kg／（m·s）表示。關(guān)于水平沙塵通量Q（輸沙率）的計(jì)算方法，國(guó)內(nèi)外很多學(xué)者通過(guò)大量的野外觀測(cè)研究與風(fēng)洞模擬研究總結(jié)出各自的計(jì)算公式（表1），根據(jù)塔中地區(qū)實(shí)測(cè)沙塵通量與理論沙塵通量的比較，Bagnold與Zingg公式的結(jié)果都比實(shí)測(cè)的要偏小，Zingg的結(jié)果最小，而Kawamura公式的結(jié)果則偏大，只有Lettau公式的計(jì)算結(jié)果與實(shí)測(cè)最接近，風(fēng)速20.0m／s以下的時(shí)候，兩條曲線幾乎重合，當(dāng)風(fēng)速大于20.0m／s的時(shí)候，實(shí)測(cè)結(jié)果又比Lettau的結(jié)果稍微偏大，不過(guò)仍然很接近［18］，因此本文計(jì)算水平輸沙量Q時(shí)采用Lettau公式。關(guān)于垂直輸沙通量F的計(jì)算，Lu等［19］、Marticorena等［20］、Shao［21－22］、Gillette等［23］都進(jìn)行了研究，并建立了自己的計(jì)算模型。其中Lu、Marticorena和Shao的輸沙通量模式理論較為完善，但是模式形式比較復(fù)雜，計(jì)算中需要的參數(shù)較多且有的參數(shù)很難確定，實(shí)際運(yùn)用中有很大難度，而Gillette等的計(jì)算方案在形式上非常簡(jiǎn)單，只需要u＊和u＊t兩個(gè)參數(shù)，比較適合簡(jiǎn)單或近似的估算，因此本文計(jì)算垂直輸沙量F時(shí)采用Gillette等公式，其公式具體如下：

根據(jù)風(fēng)蝕傳感器中的計(jì)算公式本文將摩擦速度的計(jì)算公式做簡(jiǎn)單變換，具體如下：

式中：k——卡曼常數(shù)取值0.4；z——高度取值2m；u——z處即2m高的風(fēng)速；z0——地表粗糙度。鑒于策勒、若羌風(fēng)沙觀測(cè)場(chǎng)所獲取的數(shù)據(jù)無(wú)法計(jì)算z0，本文考慮策勒、若羌與塔中自然平緩沙地上（B站）標(biāo)校自動(dòng)站地表狀況相似，因此利用西側(cè)2.2km處的兩層風(fēng)、溫同步資料，對(duì)塔中地區(qū)的地表粗糙度長(zhǎng)度z0進(jìn)行確定。計(jì)算得出年平均粗糙長(zhǎng)度z0為1.175×10－4m。塔克拉瑪干沙漠北部肖塘地區(qū)年平均粗糙度為6.05×10－5m［24］，C2是經(jīng)驗(yàn)常數(shù)，近似取1.4×10－15g·cm－6·s3。

表1 代表性輸沙率模型

塔中、肖塘、若羌3地的臨界摩擦速度分別在0.25，0.27，0.21m／s左右時(shí)，地表有沙粒起動(dòng)。因此可以近似地認(rèn)為塔中、肖塘、若羌3地的臨界摩擦速度分別為0.25，0.27，0.21m／s［25］。如圖1所示，根據(jù)H11—LIN型風(fēng)蝕傳感器所測(cè)得的撞擊顆粒數(shù)與摩擦速度做相關(guān)計(jì)算得出策勒的臨界摩擦速度為0.22m／s。

圖1 策勒風(fēng)蝕傳感器所測(cè)的摩擦速度與顆粒撞擊數(shù)的關(guān)系

2 結(jié)果與分析

在以下的表述過(guò)程中，參數(shù)單位統(tǒng)一為國(guó)際標(biāo)準(zhǔn)單位的量綱kg／（m·s），F(xiàn)量級(jí)為10－8kg／（m2·s），Q的量級(jí)為10－4kg／（m·s），為了方便比較分析，下述的F和Q省略單位和量級(jí)。由圖2—5可知，垂直沙通量的分布時(shí)間段多集中于白天及午后，這與沙塵暴多發(fā)于這一時(shí)間段有關(guān)。若羌的垂直沙通量F要明顯大于策勒、肖塘、塔中等地區(qū)。在地面2m出風(fēng)速達(dá)到2m／s時(shí)，若羌F值接近于50，最大，策勒在7.5左右次之，塔中和肖塘比較接近于3左右；在貼地面2m處的風(fēng)速達(dá)到6m／s時(shí)，若羌F值達(dá)到150最大，策勒F為20，塔中F為10，肖塘F為10。由于起沙量的大小主要取決于風(fēng)速和沙粒粒徑，肖塘本地的沙塵粒徑要小于塔中，在風(fēng)速大小基本一致的情況下，垂直沙通量F卻大致相同，這很好的印證了垂直沙通量F還受到下墊面、地形以及地表土壤含水率等因素的影響。4個(gè)站垂直沙通量F與風(fēng)速都表現(xiàn)出較好的線性相關(guān)，風(fēng)速遞增，垂直沙通量遞增；風(fēng)速遞減，垂直沙通量遞減。

圖2 2010年5月8日若羌垂直沙通量F與風(fēng)速日變化

圖3 2010年8月24日策勒垂直沙通量F與風(fēng)速日變化

圖4 2010年7月17日肖塘垂直沙通量F與風(fēng)速日變化

圖5 2010年5月26日塔中垂直沙通量F與風(fēng)速日變化

從圖6—9可知，水平沙塵通量Q與實(shí)測(cè)2m處風(fēng)速具有一致性，水平沙塵通量Q值隨風(fēng)速的變化而變化，水平沙塵通量Q集中在天氣過(guò)境時(shí)，時(shí)間節(jié)點(diǎn)多在午后，為沙塵暴的高發(fā)時(shí)段。整個(gè)水平沙塵通量的量級(jí)要大于垂直沙塵通量。當(dāng)貼地面風(fēng)速在2 m／s左右時(shí)，若羌Q值在100左右，策勒Q值在40左右，塔中和肖塘在20左右，塔中Q值稍大；當(dāng)貼地面的風(fēng)速在4m／s時(shí)，若羌Q值在300左右，策勒Q值在80左右，塔中和肖塘在40左右，塔中Q值較大；當(dāng)貼地面的風(fēng)速在6m／s時(shí)，若羌Q值在400左右，策勒Q值在125左右，塔中Q值在70左右和肖塘Q值在65左右?？梢?jiàn)在相同風(fēng)速下，若羌Q值最大，策勒次之，塔中第三，肖塘最小。在風(fēng)速遞增幅度相同條件下，若羌Q值遞增幅度最大，策勒次之，塔中、肖塘基本相同。當(dāng)貼地面風(fēng)速遞增2m／s，若羌的Q值遞增速率約為100，策勒約為40，塔中、肖塘約為20。臨界摩擦速度、粒徑、風(fēng)速、土壤含水量都是輸沙量的決定性因素，在本文4個(gè)實(shí)驗(yàn)區(qū)域中風(fēng)速等同條件下，沙地的土壤含水率基本相同，值趨于零，臨界摩擦速度及粒徑大小成為主導(dǎo)因素，在若羌臨界摩擦速度最小，策勒較大，塔中第三，肖塘最大的條件下，沙塵通量值出現(xiàn)了上述結(jié)果。但是平均粒徑只代表這一區(qū)域的平均狀況，并不能較好地詮釋相對(duì)細(xì)顆粒和相對(duì)粗顆粒的百分比含量，每次過(guò)程采樣也不能完全界定其大小，存在一定難度，所以平均粒徑大小的影響程度只能作為定性參考，定量化計(jì)算尚需努力。

圖6 2010年5月8日若羌水平沙通量Q與風(fēng)速日變化

圖7 2010年5月24日策勒水平沙通量Q與風(fēng)速日變化

圖8 2010年7月17日肖塘水平沙通量Q與風(fēng)速日變化

圖9 2010年5月26日塔中水平沙通量Q與風(fēng)速日變化

3 結(jié)論

塔里木盆地多種沙源類(lèi)型的沙塵通量（F與Q）值均與實(shí)測(cè)2m處風(fēng)速變化趨勢(shì)一致，水平沙塵通量的量級(jí)要大于垂直沙塵通量。在相同風(fēng)速下，若羌沙塵通量最大，策勒次之，塔中第三，肖塘最小。在風(fēng)速遞增幅度相同條件下，若羌水平沙塵通量遞增幅度最大，策勒次之，塔中、肖塘基本相同。本文沙通量的計(jì)算值與楊興華［15］相較，F(xiàn)與Q 的量級(jí)處于同一水平，F(xiàn)值相近，Q值稍?。慌c成天濤［7］相較，塔里木盆地輸沙量要大于渾善達(dá)克沙地，這是由臨界起沙風(fēng)速較小，粒徑較細(xì)引起的；與申彥波［12］相較，塔里木盆地4個(gè)代表站點(diǎn)的F與Q平均值的量級(jí)分別為10－8kg／（m2·s），10－4kg／（m·s）與敦煌戈壁計(jì)算的順風(fēng)向沙粒通量和垂直塵粒通量平均值的量級(jí)一致。

4 討論

就本文采用的計(jì)算輸沙量的方法而言，臨界摩擦速度、粒徑、風(fēng)速都是輸沙量的決定性影響因素，風(fēng)的動(dòng)力作用不僅僅表現(xiàn)為風(fēng)速、風(fēng)向數(shù)值的變化，還有很多間接作用。風(fēng)可以塑造沙丘的形態(tài)，風(fēng)是大自然的搬運(yùn)工，沙塵粒度的組份，與風(fēng)的起動(dòng)與遷移有著重要的聯(lián)系，沙塵顆粒能否被搬運(yùn)，很大程度上取決于風(fēng)速的大小。沙塵粗顆粒不能長(zhǎng)期在空氣中懸浮搬運(yùn)，大部分在地面跳躍遷移，沙塵細(xì)顆?；蛘邩O細(xì)顆粒可以隨風(fēng)做長(zhǎng)途懸浮搬運(yùn)，假設(shè)研究區(qū)域的盛行風(fēng)向具有較好的一致性，那么沙地長(zhǎng)期風(fēng)化的結(jié)果就是平均粒徑增大，導(dǎo)致輸沙量減小。臨界摩擦速度與輸沙量成反比關(guān)系，粒徑與輸沙量成反比關(guān)系，而風(fēng)速與輸沙量為正比關(guān)系，但是臨界摩擦速度取值與粒徑值與風(fēng)速又有著較為復(fù)雜的關(guān)系，還有其它影響要素，如液態(tài)或者固態(tài)降水改變了表層土壤含水率從而導(dǎo)致臨界摩擦速度增大，輸沙量減小。不同于臨界起沙風(fēng)速，地面起沙率和垂直起沙障礙物高度的增大而減少，另外植物株間距或障礙物間距也是一個(gè)重要的影響因素［7］。因此，輸沙量的很多計(jì)算公式都是經(jīng)驗(yàn)性的總結(jié)，鑒于各個(gè)研究區(qū)域下墊面環(huán)境的差異性及地表形態(tài)微變化，所有的經(jīng)驗(yàn)性公式都有著自己的適用性和局限性。風(fēng)蝕起沙的臨界摩擦速度、起沙風(fēng)速以及沙塵通量，也隨著天氣系統(tǒng)動(dòng)力、熱力條件改變而改變、也隨著地表狀況變化而變化。因此，沙通量變化不確定性會(huì)直接影響到沙塵數(shù)值預(yù)報(bào)模式的參數(shù)化方案的改進(jìn)，期望能通過(guò)大量的野外試驗(yàn)來(lái)確定和驗(yàn)證沙通量的變化閾值。

在沙塵暴數(shù)值預(yù)報(bào)模式中，u＊t是判定起沙過(guò)程的重要參量，一般通過(guò)參數(shù)化方案實(shí)時(shí)計(jì)算得出。本文所計(jì)算的u＊t是典型地點(diǎn)典型過(guò)程的臨界摩擦速度。通過(guò)對(duì)同一次過(guò)程野外試驗(yàn)獲取的u＊t和模式計(jì)算出的u＊t進(jìn)行比較，找出其中的差異，探索其原因，從而可為模式相關(guān)參數(shù)化方案的優(yōu)化提供依據(jù)。水平沙塵通量Q計(jì)算方案的選取本身就是一個(gè)與實(shí)驗(yàn)結(jié)合的檢驗(yàn)工作，本文利用Lettau水平沙塵通量公式正是我們研究團(tuán)隊(duì)在典型沙塵暴過(guò)程中實(shí)測(cè)收集的結(jié)果與多種計(jì)算方案擬合檢驗(yàn)后的最優(yōu)選取，可為改進(jìn)沙塵暴數(shù)值模式沙塵通量預(yù)報(bào)參數(shù)化方案提供參考。

［1］ 魏文壽，崔彩霞，尚華明，等，譯.沙漠氣象學(xué)［M］.北京：氣象出版社，2008.

［2］ Bagnold R A.The Physics of Blown Sand and Desert Dunes［M］.London：Methuen，1941.

［3］ Shao Yaping.Physics and Modeling of Wind Erosion［M］.Sydney：Kluwer Academic Publishers，2000.

［4］ Dong Zhibao，Liu Xiaoping，Wang Hongtao，et al.The flux profile of a blowing sand cloud：a wind tunnel investigation［J］.Geomorphology，2002，49（3／4）：219－230.

［5］ 李萬(wàn)源，沈志寶，呂世華，等.風(fēng)蝕影響因子的敏感性試驗(yàn)［J］.中國(guó)沙漠，2007，27（6）：985－993.

［6］ 沈志寶，申彥波，杜明遠(yuǎn)，等.沙塵暴期間戈壁沙地起沙率的觀測(cè)結(jié)果［J］.高原氣象，2003，22（6）：545－550.

［7］ 成天濤，呂達(dá)仁，徐永福.渾善達(dá)克沙地起沙率和起沙量的估計(jì)［J］.高原氣象，2006，25（2）：236－241.

［8］ 陳渭南，董治寶，楊佐濤，等.塔克拉瑪干沙漠的起沙風(fēng)速［J］.地理學(xué)報(bào)，1995，50（4）：360－367.

［9］ Gillette A D，Ranjit P.Modeling dust emission caused by wind erosion［J］.Journal of Geophysical Research，1988，93（11）：14234－14242.

［10］ 梅凡民，張小曳，鹿化煜，等.若干風(fēng)蝕粉塵釋放模型述評(píng)［J］.中國(guó)沙漠，2004，24（6）：791－796.

［11］ 王鵬祥，孫蘭東，岳平，等.鹽池地區(qū)沙塵暴期間風(fēng)沙運(yùn)動(dòng)若干特征研究［J］.中國(guó)沙漠，2007，27（6）：1077－1079.

［12］ 申彥波，沈志寶.敦煌戈壁地表風(fēng)蝕起沙量的計(jì)算［J］.高原氣象，2004，23（5）：648－653.

［13］ 申彥波，沈志寶，杜明遠(yuǎn)，等.風(fēng)蝕起沙的影響因子及其變化特征［J］.高原氣象，2005，24（4）：613－616.

［14］ 李江風(fēng).沙漠氣候［M］.北京：氣象出版社，2002.

［15］ 楊興華，何清，艾力·買(mǎi)買(mǎi)提明.塔中地區(qū)一次沙塵暴過(guò)程的輸沙通量估算［J］.干旱區(qū)研究，2010，27（6）：971－974.

［16］ 霍文，何清，劉新春，等.塔克拉瑪干沙漠及其大氣降塵點(diǎn)粒度特征分析［J］.水土保持研究，2009，16（4）：207－211.

［17］ 楊興華，何清，艾力·買(mǎi)買(mǎi)提明.塔克拉瑪干沙漠塔中地區(qū)春夏季風(fēng)蝕起沙研究［J］.中國(guó)沙漠，2010，30（4）：770－776.

［18］ 楊興華.塔克拉瑪干沙漠塔中地區(qū)春夏季沙粒運(yùn)移及動(dòng)力參數(shù)研究［D］.烏魯木齊：新疆師范大學(xué)，2009.

［19］ Lu Hua，Shao Yaping.A new model for dust emission by saltation bombardment［J］.Journal of Geophysical Research，1999，104（14）：16827－16841.

［20］ Marticorena B，Bergametti G，Aumont B，et al.Modeling the atmospheric dust cycle.Part 2：simulation of Saharan dust sources［J］.Journal of Geophysical Research，1997，102（4）：4387－4404.

［21］ Shao Y，Raupach M R，F(xiàn)indlater P A.Effect of saltation bombardment on the entrainment of dust by wind［J］.Journal of Geophysical Research，1993，98（7）：12719－12726.

［22］ Shao Yaping，Lu Hua.A simple expression for wind erosion threshold friction velocity［J］.Journal of Geophysical Research，2000，105（17）：22437－22443.

［23］ 何清，繆啟龍，張瑞軍，等.塔克拉瑪干沙漠肖塘地區(qū)空氣動(dòng)力學(xué)粗糙度分析［J］.中國(guó)沙漠，2008，28（6）：1011－1016.

［24］ 何清，楊興華，艾力·買(mǎi)買(mǎi)提明.塔克拉瑪干沙漠風(fēng)蝕起沙觀測(cè)研究［J］.中國(guó)沙漠，2011，31（2）：315－322.