拉薩河下游河谷區(qū)風(fēng)沙災(zāi)害現(xiàn)狀、成因及發(fā)展趨勢探討*

常春平,原立峰

（1.河北師范大學(xué)資源與環(huán)境科學(xué)學(xué)院,石家莊050016;2.南京郵電大學(xué)電子信息科學(xué)系,南京 210003）

風(fēng)沙災(zāi)害以其發(fā)生頻率高、影響范圍廣、危害程度大對經(jīng)濟(jì)建設(shè)、社會發(fā)展、大氣環(huán)境、交通運輸及人民生命財產(chǎn)和健康帶來嚴(yán)重不利影響[1]。拉薩河下游地帶北依念青唐古拉山脈,南達(dá)喜馬拉雅山脈東段北側(cè)。作為西藏自治區(qū)的政治、經(jīng)濟(jì)、文化中心的拉薩市位于拉薩河下游寬谷區(qū)。隨著青藏鐵路的建設(shè)和拉薩市經(jīng)濟(jì)、旅游及城市建設(shè)的發(fā)展,拉薩河下游河谷風(fēng)沙災(zāi)害已成為影響該區(qū)域首要環(huán)境問題。本文在遙感影像解譯和實地調(diào)查基礎(chǔ)上,對該區(qū)域風(fēng)沙災(zāi)害現(xiàn)狀、危害、成因與發(fā)展趨勢進(jìn)行了分析,為進(jìn)一步風(fēng)沙治理提供依據(jù)。

1 拉薩河下游河谷區(qū)域背景

研究區(qū)東自拉薩大橋、西至曲水的拉薩河下游河谷地帶。沿拉薩河河谷全長64.5 km,寬 2～7 km 不等 。地理坐標(biāo)介于東經(jīng) 91°09′07″-90°43′52″,北緯 29°20′05″-29°39′29″,面積約 727.5 km2。

拉薩河下游河谷區(qū)因受不同等級的斷裂(帶)控制及第四紀(jì)以來強(qiáng)烈的構(gòu)造運動影響,河谷兩側(cè)的山體巖石十分破碎,山勢陡峻。山坡中下部以及山前地帶堆積厚度不等的第四紀(jì)砂礫質(zhì)沉積物;河道成多股辮狀,河漫灘和河心洲沙丘發(fā)育。河谷平原為巨厚的第四系松散堆積物,大都為砂卵石層上覆砂土或亞砂土。本區(qū)屬高原溫帶半干旱氣候,夏季溫涼多雨、冬季干冷多風(fēng)。風(fēng)沙天氣也主要發(fā)生在冬春季節(jié)(2-5月),起沙風(fēng)以偏西風(fēng)為主。地表水資源主要是拉薩河及其支流。植被分區(qū)屬于青藏高原中部溫性草原地帶的雅魯藏布江中游谷地灌叢草原區(qū),地帶性植被為溫性草原。

截止2006年,拉薩市城關(guān)區(qū)、堆龍德慶縣和曲水縣三區(qū)縣農(nóng)村總?cè)丝?09 666人,勞動力35 093人;國民生產(chǎn)總值 60 940萬元,農(nóng)村經(jīng)濟(jì)總收入39 564萬元。

2 沙塵源分布特征與危害

2.1 沙塵源分布特征

根據(jù)IKONOS和Quick Bird影像解譯結(jié)果,風(fēng)沙源分布特征為:

(1)沿河谷及河谷兩側(cè)分布。從沙塵源地的物質(zhì)基礎(chǔ)來看,絕大多數(shù)是拉薩河及其一級支流沖積物,即使有部分沙塵源地也已被風(fēng)力侵蝕改造,沙塵物質(zhì)源頭仍是河流沖積物。極少數(shù)沙塵源地的物質(zhì)基礎(chǔ)為洪積物,它們也分布于拉薩河兩側(cè)的谷地前緣。在河流沖積物基礎(chǔ)上成為沙塵源地(耕地除外)的面積占沙塵源地總面積的31.91%。

(2)分布地貌部位多樣。復(fù)雜多樣的河谷地貌類型,決定了沙塵源地分布地貌部位的多樣性。根據(jù)野外調(diào)查,沙塵源地主要分布在洪積扇、邊灘(高河漫灘)和心灘、河流階地面和河谷兩側(cè)山坡,分布在這四種地貌部位的沙塵源地總面積(耕地除外),自大到小依次為洪積扇＞邊灘(高河漫灘)和心灘＞河谷兩側(cè)山坡＞河流階地面。

(3)呈小面積零星分布。研究區(qū)處于拉薩河下游,河床寬闊,河灘亂流,常常將低級階地、河漫灘和心灘侵蝕切割成零星斑塊。河谷兩側(cè)山體高大,拉薩河沿岸突兀的中小山嘴眾多,河谷地帶近地面風(fēng)場在大尺度大氣流場的控制下,受地形影響顯著,山谷風(fēng)和山嘴附近局部回流十分強(qiáng)盛,形成很多面積不大、相對孤立的山坡覆沙體。

其次,人類經(jīng)濟(jì)活動對沙塵源地的形成和分布特征也產(chǎn)生重大影響,主要是天然草地利用不當(dāng)而退化成沙塵源地,以及大面積耕地在冬春季節(jié)成為沙塵源地等等。

2.2 風(fēng)沙危害

(1)污染城市大氣環(huán)境。根據(jù)在拉薩市1990-2001年的環(huán)境空氣質(zhì)量監(jiān)測數(shù)據(jù)和2003年西藏自治區(qū)環(huán)境狀況公報表明,SO2、NO2和其他NOX的監(jiān)測值無論何時都小于100 μ g/m3,SO2日均值濃度 2～ 4 μ g/m3,年均值濃度為 2 μ g/m3;NOX日均值濃度 2～ 73 μ g/m3,年均值濃度為 26 μ g/m3,符合《環(huán)境空氣質(zhì)量標(biāo)準(zhǔn)》一級標(biāo)準(zhǔn)的要求。但TSP或PM10的小時監(jiān)測值在大風(fēng)時超過1 000 μ g/m3,TSP日均值濃度 41～1 437 μ g/m3,年均值濃度為337 μ g/m3;降塵量 3.32 ～ 31.34t/(km2?月),嚴(yán)重超標(biāo)。所以,懸浮顆粒物為首要污染物,且空氣污染程度較重時段均與風(fēng)沙天氣有關(guān)[2-4]。

(2)危害農(nóng)牧業(yè)生產(chǎn)和水利設(shè)施。首先,蠶食可利用土地。研究區(qū)為西藏自治區(qū)的政治、經(jīng)濟(jì)和文化核心地帶,可利用土地資源十分寶貴,但根據(jù)IKONOS和Quick Bird影像解譯結(jié)果表明現(xiàn)有裸露和半裸露沙塵源地面積27 264.74 hm2,占總土地面積37.48%,且仍在繼續(xù)擴(kuò)大,造成土地失去利用價值。

其次,降低土地質(zhì)量。根據(jù)野外調(diào)查,分布于拉薩河階地面上的農(nóng)田,年風(fēng)蝕深度一般可達(dá)0.5～3.0 cm。風(fēng)蝕模擬研究結(jié)果顯示,采集于調(diào)查區(qū)的農(nóng)田沙土和沙壤土原生樣品,在風(fēng)速分別達(dá)到6.8 m/s和7.6 m/s時便發(fā)生土壤風(fēng)蝕起沙現(xiàn)象,并且隨風(fēng)速增大,風(fēng)蝕量迅速增加。表層耕作土的細(xì)顆粒物質(zhì)和營養(yǎng)元素被大量吹蝕。用當(dāng)?shù)爻Ｒ姷?5 m/s風(fēng)速吹蝕85 min后,不論沙土和沙壤土,土壤表層中不僅＜0.063 mm的粉黏粒組分普遍減少4.2%～4.8%,甚至0.063～1.0 mm的極細(xì)紗和粗砂也減少1.78%～1.85%,而＞1.0 mm的極粗砂和礫石組分增加了4.1%～6.69%。因而其平均粒徑都相應(yīng)增加0.03 mm左右。

第三,直接危害種植業(yè)和畜牧業(yè)。調(diào)查中發(fā)現(xiàn),研究區(qū)大面積的草地因沙埋而退化,不僅草地建群種發(fā)生明顯變化,產(chǎn)草量也大幅降低。風(fēng)沙流對春季稚嫩的禾苗具有強(qiáng)烈打擊作用,常導(dǎo)致禾苗植物組織嚴(yán)重受傷而減產(chǎn)甚至枯死。

第四,填埋水利設(shè)施。調(diào)查中發(fā)現(xiàn),在菜納、協(xié)榮等地修建的灌溉系統(tǒng)普遍遭受風(fēng)沙危害,冬春季節(jié)渠道積沙嚴(yán)重,部分渠道不得不采取地下渠的方式。這不僅增加投資,而且明渠部分每年需要清沙,增大維護(hù)費用,影響水利設(shè)施正常的灌溉功能。

(3)威脅交通安全。冬春季節(jié)經(jīng)常出現(xiàn)大風(fēng)天氣,風(fēng)沙活動強(qiáng)度大,造成能見度大幅降低,嚴(yán)重威脅交通安全。當(dāng)前,公路交通受影響最大的是拉薩河南岸沿線。據(jù)查訪,風(fēng)沙天氣嚴(yán)重時的能見度僅50 m左右,有時甚至小于50 m,已經(jīng)到了車輛無法通行的程度。青藏鐵路開通后,風(fēng)沙天氣對鐵路安全運營的威脅是顯而易見的。

3 風(fēng)沙災(zāi)害成因分析

3.1 氣候因素

拉薩河谷冬春季節(jié)干冷多風(fēng)、年內(nèi)降水分配不均。由于地勢海拔高,大部分地面位于大氣對流層1/3的厚度內(nèi),空氣干潔,太陽總輻射量高,但反輻射和長波有效輻射量很大,輻射平衡值相對偏小。由此決定了年平均氣溫較低(7.7℃),而且變化強(qiáng)烈(月均最高與最低氣溫差17.7℃)。受喜馬拉雅山脈阻擋,印度洋西南氣流進(jìn)入拉薩河下游地帶時已成為強(qiáng)弩之末,年均降水量僅有435 mm,且95%以上集中在夏半年(5-9月)的少風(fēng)時段。冬半年降水稀少,大風(fēng)頻發(fā),而同期干燥度卻高達(dá)10以上。在這種氣候條件下,風(fēng)蝕氣候侵蝕力很強(qiáng),年風(fēng)蝕氣候因子指數(shù)一般在50以上,為風(fēng)沙活動提供了極為有利的氣候條件[5]。

3.2 沙源因素

根據(jù)野外調(diào)查和遙感影像判讀結(jié)果,按沙質(zhì)地表流動性和砂礫質(zhì)地表的植被覆蓋狀況,將沙塵源地劃分為8類,再根據(jù)所處地貌部位劃分為24個亞類。統(tǒng)計表明,沙塵源地總面積達(dá)27 264.74 hm2,占土地總面積的37.48%。從9種典型沙塵源地的地表粒度組成可以看出(圖1),每種地表類型沙粒粒徑集中分布在0.25 mm以下,以細(xì)沙、極細(xì)沙和黏粒成分為主,平均含量占60.69%。根據(jù)董治寶等對風(fēng)蝕率與粒度關(guān)系特征的研究表明[6],風(fēng)成沙粒徑0.4～0.075 mm為易蝕性顆粒。研究區(qū)風(fēng)沙源地表物質(zhì)易蝕性顆粒含量顯示,易蝕性顆粒含量達(dá)到69.64%,總體分布趨勢是流動沙地＞半固定沙地＞固定沙地。這些顆粒極易被風(fēng)力吹揚(yáng)成為大氣的沙塵物質(zhì),為拉薩市風(fēng)沙天氣提供了豐富的物質(zhì)來源。

3.3 植被因素

在山坡下部沙黃土母質(zhì)上發(fā)育的冷棕鈣土,地帶性植被有寬苞金背柳灌叢和砂生槐灌叢,主要受到放牧的影響,成為退化灌叢。植被蓋度和高度都受到明顯的影響,寬苞金背柳灌叢的蓋度在45%左右,群落高度10 cm左右。物種可達(dá)到10種以上,常見的伴生種有砂生槐、固沙草、披堿草、針茅和豆科黃芪屬的幾種植物等。砂生槐灌叢的蓋度也在50%左右,群落平均高度20～30 cm,常常伴生有藏沙蒿、固沙草、青藏苔草、西藏黃芪、勁直黃芪、甘遂、馬先蒿等,可達(dá)到10余種。

圖1 典型沙塵源地地表物質(zhì)粒度組成

沙礫地和沙地上的原生植被砂生槐灌叢和固沙草群落,由于生態(tài)環(huán)境脆弱,受人為活動影響較大,處于沙地植被的不同演替階段。砂生槐灌叢的蓋度最低只有30%,在保護(hù)較好的地區(qū),可以達(dá)到75%。而在覆沙的沙礫地上,固沙草常常成為砂生槐的共建種,有時也形成砂生槐-藏龍蒿或砂生槐-藏沙蒿群落;群落物種數(shù)常常為2～6個,常見的有固沙草、藏沙蒿、藏龍蒿、勁直黃芪、小畫眉草等。固沙草群落的蓋度從10%到60%都會出現(xiàn)。固定沙地上的固沙草草原蓋度可以達(dá)到50%～60%,但常形成單優(yōu)勢種群落,群落中的伴生種很少,常常在1 m×1 m范圍內(nèi)僅有1～2個伴生物種出現(xiàn),而在半固定沙地上,蓋度只有10%～40%,但物種組成有所提高,可達(dá)到3～4種,并常與藏龍蒿、藏沙蒿和砂生槐形成斑塊狀鑲嵌分布的半流動-半固定沙地植被景觀。

隨著沙礫地和沙地上的原生植被遭到破壞,有較大面積的沙地處于流動和半固定狀態(tài),流動沙地上次生演替的早期階段主要為毛瓣棘豆群居,蓋度最大可達(dá)20%～30%,常形成單優(yōu)勢種群落;而在沙礫地上,次生演替的早期階段主要為勁直黃芪-藏沙蒿群落,蓋度常常5%～10%,物種數(shù)3～5個。在半固定沙地上以藏沙蒿灌叢和藏龍蒿灌叢為主,常伴生固沙草、勁直黃芪、毛瓣棘豆等。隨著群落的演替進(jìn)程,漸有砂生槐定居,在景觀尺度上呈斑塊狀鑲嵌分布。在河岸邊水分條件相對較好的沙礫地和沙地上生長有三春柳灌叢和白草草原,三春柳灌叢的蓋度一般為30%～40%,常有白草伴生。白草草原作為演替中間階段的群落,常形成單優(yōu)勢群落,蓋度可達(dá)20%～40%。

對沙地植被不同演替階段的植物群落蓋度、物種數(shù)和高度進(jìn)行了對比(圖2)。群落蓋度可很好地反映沙地的固定程度。隨沙地植被受干擾強(qiáng)度的增加,從固定沙地到半固定沙地,再到半流動沙地,植物群落的蓋度和高度都有明顯降低,而蓋度和高度是影響植物群落防沙和固沙效率的主要因素。研究區(qū)砂礫地和沙地植被都受到了不同程度的干擾,以半固定-半流動沙地植被為主。因此,植被退化是風(fēng)沙活動強(qiáng)度增強(qiáng)的主要因素之一。

圖2 拉薩河谷沙地植物群落特征

3.4 人類生產(chǎn)行為

近年來,研究區(qū)的人類工程行為規(guī)模越來越大,涉及的空間范圍很廣,對原本脆弱的生態(tài)環(huán)境造成的局部損害逐漸加劇。人類工程行為主要有交通建設(shè)、采挖砂石和建筑工程三種類型,其中前兩者造成的沙塵源地面積最大,危害最嚴(yán)重。盡管從絕對面積來看,這三種類型的沙塵源地都不大,但工程建設(shè)過程中對原生地表產(chǎn)生的破壞程度極大,沙質(zhì)或砂礫地表的物質(zhì)內(nèi)聚力幾乎全部喪失。工程建設(shè)場地在缺乏保護(hù)性措施的情況下,冬春季節(jié)的大風(fēng)極易使干燥疏松的地表物質(zhì)風(fēng)蝕,并形成高濃度的沙塵氣流。

圖3 典型地表類型起沙強(qiáng)度

通過對流沙區(qū)、砂礫質(zhì)地、翻耕農(nóng)田、半固定沙地(沙質(zhì)退化灌草地)起沙強(qiáng)度計算分析表明(圖3),隨著風(fēng)速增大,地表起沙強(qiáng)度迅速增強(qiáng),尤其是砂礫質(zhì)地表遭人為破壞后,起沙強(qiáng)度呈現(xiàn)數(shù)倍甚至數(shù)十倍增加。就原生地表而言,流沙地的起沙強(qiáng)度最大;風(fēng)速在23 m/s以下時,起沙強(qiáng)度從大到小依次為翻耕農(nóng)田、半固定沙地(沙質(zhì)退化灌草地)和砂礫質(zhì)地表。農(nóng)業(yè)生產(chǎn)活動場所的耕地盡管起沙起塵強(qiáng)度沒有流沙地大,但其顯著特點是面積廣、可提供的沙塵量大,臨界起沙風(fēng)速相對較小,使耕地對大氣粉塵的貢獻(xiàn)率較高,成為沙塵源的主體之一。牧業(yè)生產(chǎn)利用的土地類型多樣,包括沙質(zhì)、砂礫質(zhì)地表,這些地類因過度利用而退化,是目前主要的沙塵源地。

4 沙害發(fā)展趨勢分析

4.1 自然影響因素的發(fā)展趨勢

氣候因素是最活躍、多變的,是制約植被、土壤、水文等因子的主導(dǎo)因素。因此,對沙害發(fā)展趨勢進(jìn)行分析,主要著重討論氣候因素(重點是溫度、降水和大風(fēng))的可能變化趨勢。據(jù)吳祥定、林振耀等人的研究成果[5-9],氣溫(圖4A)在17世紀(jì)80年代前為寒冷期,17世紀(jì)80年代至18世紀(jì)20年代為溫暖期,隨后到80年代為較寒冷期,從18世紀(jì)80年代至今為相對溫暖期。降水量的變化略有不同(圖4B),19世紀(jì)以前為干濕波動期,19世紀(jì)中葉至20世紀(jì)初為干旱少雨期,20世紀(jì)初至60年代為相對多雨期,70年代又處于相對干旱少雨期。20世紀(jì)70年代以來,我們根據(jù)現(xiàn)代氣象觀測數(shù)據(jù)進(jìn)行了補(bǔ)充研究,結(jié)果表明,氣溫自20世紀(jì)70年代末至今,仍處于顯著的繼續(xù)升高階段(圖5A),而降水也略有增加(圖5B),大風(fēng)日數(shù)略有減少。

在以上中尺度氣候變化中,也明顯存在小尺度(10 a尺度)的波動。從拉薩市氣象臺的記錄數(shù)據(jù)看,近期氣候波動較大,分析其變化并根據(jù)其冷暖、干濕和大風(fēng)期的演替規(guī)律和演替期時間長短,可以粗略地判定其小尺度氣候變化周期。根據(jù)小波分析結(jié)果,溫度變化主要有33年周期,降水以8 a為主要周期,大風(fēng)日數(shù)則以24 a周期最明顯(表1)。按照這一周期推算,未來20 a內(nèi)的氣溫可能會維持在相對較高的水平;降水量在未來10 a內(nèi)可能有所增加,但總體增加幅度不大;大風(fēng)日數(shù)在未來15 a內(nèi)總體上不會減少?；谶@樣的推算結(jié)果,在未來15 a內(nèi),氣候變化將對風(fēng)沙災(zāi)害起到繼續(xù)強(qiáng)化的作用。

表1 氣候要素周期分析

圖4 西藏1978年以前氣溫(A)和降水(B)變化

圖5 拉蒴市自20世紀(jì)50年代以來的氣溫(A)和降水(B)變化

4.2 人文影響因素的發(fā)展趨勢

由于研究區(qū)內(nèi)居民環(huán)護(hù)意識的不斷增強(qiáng),農(nóng)牧業(yè)生產(chǎn)對沙塵天氣發(fā)生的貢獻(xiàn)率沒有明顯上升。根據(jù)西藏自治區(qū)統(tǒng)計資料,從研究區(qū)的人口變化、農(nóng)作物播種面積和牲畜年末存欄數(shù)分析,1992-2002年間研究區(qū)人口由188 662人增長到217 828人,增加了29 166人;農(nóng)作物播種面積由11 786.7 hm2增加到12 667 hm2,增加了880.3 hm2;牲畜總頭數(shù)由77.66萬羊單位下降為71.52萬羊單位,下降了6.14萬羊單位。所以,近10 a來上述三項指標(biāo)變化微小,加之農(nóng)田基本建設(shè)力度逐年加大,農(nóng)田風(fēng)蝕基本維持在原有的水平。牧業(yè)發(fā)展基本上處于多年下降的趨勢,草地過載現(xiàn)象將會減緩;預(yù)計在國家人口政策不發(fā)生改變的情況下,人口增長率將維持目前的水平或有所下降;另外,隨著第三產(chǎn)業(yè)的興起,產(chǎn)業(yè)結(jié)構(gòu)將得到有效調(diào)整,對土地壓力增加的幅度不會過大?？傮w來看,人類活動將繼續(xù)起到增強(qiáng)風(fēng)沙災(zāi)害的作用,但不會有顯著增強(qiáng)趨勢。

總之,在10 a尺度上的氣候變化相對于人文因素變化,其總體趨勢是相對緩慢的,但其影響范圍十分廣泛和深刻。在前述氣候變化的背景下,干燥度將繼續(xù)上升,由此決定了植被的自然恢復(fù)能力不能得到提高,依靠植被自然演替而固定沙塵源地非常困難。因此,氣候變化趨勢在未來15 a內(nèi)不利于削弱風(fēng)沙災(zāi)害強(qiáng)度。人文因素的發(fā)展趨勢在未來15 a內(nèi),盡管對風(fēng)沙災(zāi)害強(qiáng)度的加強(qiáng)作用會逐漸降低,但總體上仍起強(qiáng)化作用。因此,拉薩河下游河谷地帶的風(fēng)沙災(zāi)害如果沒有人為的強(qiáng)迫干擾(即有效治理)活動,風(fēng)沙災(zāi)害將會進(jìn)一步發(fā)展,拉薩市的人居環(huán)境不可能得到顯著改善。

[1]徐海量,陳亞寧.塔里木盆地風(fēng)沙災(zāi)害危險性評價[J].自然災(zāi)害學(xué)報,2003,12(2):35-39.

[2]倪雪.拉薩貢嘎機(jī)場干季風(fēng)沙、浮塵天氣特征及預(yù)報初探[J].四川氣象,2001,21(4):46-48.

[3]周順武,黃瓊中.拉薩市空氣污染物濃度的統(tǒng)計分析及其可能的天氣原因[J].西藏科技,2002(5):39-42.

[4]黃瓊中.拉薩市環(huán)境空氣質(zhì)量與氣象特征分析[J].中國環(huán)境監(jiān)測,2001,17(6):50-53.

[5]繆啟龍,向毓意,顧顯躍.氣候變化對西藏環(huán)境經(jīng)濟(jì)的可能影響[J].應(yīng)用氣象學(xué)報,1998,9(2):225-230.

[6]董治寶,李振山.風(fēng)成沙粒度特征對其風(fēng)蝕可蝕性的影響[J].土壤侵蝕與水土保持學(xué)報,1998,4(4):1-5,12.

[7]代加冼.青藏高原氣候?qū)W[M].北京:氣象出版社,1990.

[8]吳祥定,林振耀.歷史時期青藏高原高原氣候變化特征的初步分析[J].氣象學(xué)報,1981,39(1):90-97.

[9]吳祥定,林振耀.橫斷山區(qū)近代氣候變化的研究[J].地理研究,1987,6(2):48-55.

[10]林振耀,趙昕奕.青藏高原氣溫降水變化的空間特征[J].中國科學(xué):D輯,1996,26(4):354-358.