杜家穎， 崔麗峰， 涂成龍， 盛茂銀， 劉瑞祿， 李 勇， 孫權(quán)忠

(1.貴州師范大學(xué)喀斯特研究院，貴州貴陽(yáng) 550001； 2.中國(guó)科學(xué)院地球化學(xué)研究所環(huán)境地球化學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室，貴州貴陽(yáng) 550081；3.貴州省水土保持監(jiān)測(cè)站，貴州貴陽(yáng) 550002)

我國(guó)的喀斯特地貌主要分布在以貴州省為中心的西南地區(qū)。由于脆弱的生態(tài)環(huán)境和復(fù)雜的人地系統(tǒng)，加上歷史上不合理的人為活動(dòng)[1-3]，使喀斯特地區(qū)生態(tài)較為脆弱，在經(jīng)歷了長(zhǎng)期大規(guī)模的破壞性、掠奪性開(kāi)墾后，地表植被遭到嚴(yán)重破壞，引發(fā)了嚴(yán)重的水土流失問(wèn)題。我國(guó)西南喀斯特區(qū)域處于世界三大集中連片喀斯特發(fā)育區(qū)之一的東亞片區(qū)中心，喀斯特地區(qū)土地貧瘠、水資源短缺，加上脆弱的生態(tài)系統(tǒng)，已經(jīng)成為制約我國(guó)西南喀斯特地區(qū)區(qū)域經(jīng)濟(jì)社會(huì)發(fā)展的一個(gè)重大生態(tài)問(wèn)題[4]。

在喀斯特地區(qū)，水土流失的影響因素眾多，土壤不僅隨地表徑流流失，同時(shí)還沿巖溶孔(裂)隙、漏斗、落水洞等向地下流失[5]。在制約水土流失的眾多因素中，覆被的影響一直是受到普遍關(guān)注的熱點(diǎn)問(wèn)題。不同覆被下產(chǎn)沙、產(chǎn)流隨著降雨量、降雨強(qiáng)度[6-7]以及植被覆蓋度的變化而變化[8-9]。目前，對(duì)于長(zhǎng)江和黃河流域產(chǎn)流和產(chǎn)沙已有較深入的研究[10]，而對(duì)于西南喀斯特地區(qū)覆被變化對(duì)產(chǎn)沙和產(chǎn)流的影響研究較少；研究西南喀斯特地區(qū)黃色石灰土覆被變化對(duì)產(chǎn)沙和產(chǎn)流的影響，將有助于深入了解喀斯特地區(qū)覆被對(duì)地表徑流形成與產(chǎn)沙、產(chǎn)流變化的影響，揭示不同覆被下的降雨對(duì)產(chǎn)沙和產(chǎn)流的影響。針對(duì)喀斯特地區(qū)覆被類型不同、水土流失嚴(yán)重的特點(diǎn)，本研究以關(guān)嶺縣螞蝗田小流域?yàn)槔?，基?種不同覆被類型(水保林、經(jīng)果林和撂荒地)的坡面徑流小區(qū)，分析喀斯特黃壤區(qū)坡地產(chǎn)流產(chǎn)沙特征和規(guī)律，為喀斯特地區(qū)水土保持的生態(tài)系統(tǒng)恢復(fù)重建提供理論支撐。

1 材料與方法

1.1 研究區(qū)概況

研究區(qū)位于貴州省關(guān)嶺布依族苗族自治縣，屬于云貴高原東部脊?fàn)钚逼履蟼?cè)向廣西丘陵傾斜的斜坡地帶。地勢(shì)西北高、東南低。境內(nèi)最高點(diǎn)位于永寧鎮(zhèn)的舊屋基大坡，海拔 1 850 m，最低點(diǎn)在打邦河注入北盤江的三江口處，海拔 370 m。大部分地區(qū)海拔在800～1 500 m?？h城海拔為 1 025 m。境內(nèi)山脈屬烏蒙山系，山體多起伏綿延。地貌具有高低起伏大、類型復(fù)雜多樣的特征，碳酸鹽巖分布廣泛。關(guān)嶺境內(nèi)氣候呈立體狀，跨越南溫帶、北亞熱帶、中亞熱帶，主要以中亞熱帶季風(fēng)濕潤(rùn)氣候?yàn)橹鳎瑹崃砍渥?，水熱同季。境?nèi)12.5%的低熱河谷地區(qū)有天然溫室之稱。

螞蝗田水土保持監(jiān)測(cè)點(diǎn)始建于2008年，地理坐標(biāo)為105°34′36″E、25°49′2″N。該監(jiān)測(cè)點(diǎn)所在的螞蝗田小流域面積2.43 km2。小流域?qū)賮啛釒駶?rùn)季風(fēng)氣候，降水集中于6—9月，多年平均降水量為1 236 mm，多年平均氣溫為19.2 ℃；現(xiàn)有植被類型多為次生喬木林、疏幼林，喬木品種主要有柳杉、杉木、馬尾松、女貞、滇柏等，土壤類型為黃壤土。

1.2 試驗(yàn)設(shè)計(jì)

本研究采用坡面徑流小區(qū)觀測(cè)法。經(jīng)過(guò)實(shí)地考察研究，選擇不同覆被的直行坡上設(shè)立徑流場(chǎng)，坡向東偏南，堰型為無(wú)側(cè)向收縮矩形薄壁堰。在研究區(qū)內(nèi)設(shè)有6個(gè)坡面徑流小區(qū)來(lái)開(kāi)展坡面尺度的地面徑流和泥沙觀測(cè)。如表1所示，1、4、5號(hào)小區(qū)為經(jīng)果林(喬木，黃花梨)，2、3號(hào)小區(qū)為水保林(灌木，女貞)，6號(hào)小區(qū)為撂荒草地。用面積為100 m2或200 m2的徑流集水小區(qū)控制坡面徑流，其長(zhǎng)度為20 m，寬度為5 m或10 m。小區(qū)下方設(shè)有集流槽、分流池和徑流池，分流級(jí)別為一級(jí)，分流孔數(shù)量為9個(gè)，分流系數(shù)為9。每條流域的徑流小區(qū)旁設(shè)置1臺(tái)自記雨量計(jì)，記錄單次降雨量、降雨歷時(shí)和降雨強(qiáng)度。水位計(jì)為SWZ型周記式水位計(jì)。

表1螞蝗田監(jiān)測(cè)點(diǎn)徑流小區(qū)資料統(tǒng)計(jì)

1.3 測(cè)定項(xiàng)目與方法

1.3.1 產(chǎn)沙、產(chǎn)流 每次降雨采集1次，每次每個(gè)小區(qū)分別采集3個(gè)水樣和泥沙樣品，然后對(duì)樣品進(jìn)行實(shí)驗(yàn)室過(guò)濾、烘干、稱質(zhì)量等分析處理，得出每次降雨后產(chǎn)流量和產(chǎn)沙量。通過(guò)對(duì)卡口站自動(dòng)水位計(jì)的觀測(cè)記錄，用自動(dòng)水位取樣器進(jìn)行水樣分析處理，得出卡口站水位與流量關(guān)系和觀測(cè)河道(段)泥沙懸移質(zhì)、推移質(zhì)數(shù)量，以及降雨洪峰流量及洪水過(guò)程，從而可計(jì)算出流域侵蝕模數(shù)、侵蝕量等。地表產(chǎn)流徑流量計(jì)算公式為f(x)=(0.602+0.083×h/0.46)×8.854×(h+0.001 2)3/2；式中：h為水位，m。徑流量(m3)=流量×產(chǎn)流時(shí)間；產(chǎn)沙量(kg)=含沙量×徑流量。

1.3.2 蓄水能力、保土能力 根據(jù)各小區(qū)內(nèi)的徑流系數(shù)、侵蝕模數(shù)進(jìn)行換算[10]。

1.3.3 數(shù)據(jù)處理 應(yīng)用Excel 2010和SPSS 22.0統(tǒng)計(jì)分析軟件對(duì)測(cè)定數(shù)據(jù)進(jìn)行整理，分析降雨、地表產(chǎn)流、地表產(chǎn)沙特征，同時(shí)采用Excel 2010進(jìn)行制圖。

2 結(jié)果與分析

2.1 降雨特征

按照氣象學(xué)對(duì)于降雨等級(jí)的劃分方法[小雨(P<10 mm)、中雨(10≤P<25 mm)、大雨(25≤P<50 mm)、暴雨(50≤P<100 mm)、大暴雨(P>100 mm)]進(jìn)行分級(jí)。關(guān)嶺縣螞蝗田水土保持監(jiān)測(cè)點(diǎn)降雨總量為1 039.7 mm，降雨歷時(shí)為183 d，侵蝕性降雨總量為504.5 mm。單次最大降雨量為 57.4 mm。從圖1可以看出，降雨集中在5—10月，降雨量占全年的87.5%。7月降雨量最高，為298 mm，4月降雨量最低，為2.2 mm。其中，2月、4月和12月降雨量最少，不超過(guò) 20 mm，5月、6月和7月降雨量最多，超過(guò)150 mm，為大暴雨級(jí)別，而其他月份降雨都超過(guò)25 mm。1月降雨時(shí)間最多(24 d)，4月降雨時(shí)間最少(2 d)。單次最大30 min降雨強(qiáng)度 57.00 mm/h 和單次最大60 min降雨強(qiáng)度37.50 mm/h均發(fā)生在5月22日的降雨事件中。

由此可知，全年中最大降雨量主要集中在5月、6月和7月，且占全年降雨量的66%。除了2月、4月和12月外，其他月份降雨量大多數(shù)在大雨級(jí)別以上。因此，大雨頻繁地出現(xiàn)，對(duì)地表產(chǎn)生沖刷，而喀斯特地區(qū)的黃色石灰土容易流失，從而導(dǎo)致地表徑流的產(chǎn)生，還有可能會(huì)通過(guò)巖溶孔隙、漏斗、落水洞等向地下流失。

2.2 降雨對(duì)土壤侵蝕量的影響

土壤顆粒隨著徑流遷移或者損失等過(guò)程是侵蝕泥沙量的主要表現(xiàn)形式，侵蝕泥沙量受到降雨特征、土壤剖面、地形因子等綜合影響[11]。由于研究小區(qū)的坡度基本一致，可以認(rèn)為在本研究中坡度差異對(duì)地表產(chǎn)流產(chǎn)沙的影響不大。從圖2可以看出降雨量與侵蝕泥沙量呈相關(guān)的變化趨勢(shì)。在坡度一致的情況下，降雨量和泥沙侵蝕量對(duì)于不同土地覆被類型的影響不大。泥沙侵蝕量隨著單次降雨量的增加，侵蝕力度也隨著增加，這是由于土壤中非毛管-毛管孔隙是水分入滲的主要通道，當(dāng)降雨強(qiáng)度增大時(shí)，雨滴動(dòng)能逐漸增大，導(dǎo)致雨滴打擊所產(chǎn)生的擠壓力增大，這種擠壓力加快了入滲水流的速度，增加了侵蝕力度[12]。土壤中非毛管孔隙和毛管孔隙達(dá)到最大極限后，土壤含水量保持不變，多余水分則產(chǎn)流，隨著產(chǎn)流量的增加，泥沙也隨之流失。

2.3 地表產(chǎn)流、產(chǎn)沙與降雨特征關(guān)系

2.3.1 地表產(chǎn)流與降雨特征關(guān)系 坡度對(duì)地表及地下產(chǎn)流的影響明顯，研究區(qū)內(nèi)的小區(qū)坡度大致相同。3、4月降雨次數(shù)少且降雨強(qiáng)度較小，加上冬季降水少，土壤含水量低，導(dǎo)致坡面產(chǎn)流產(chǎn)沙量很小。從圖3可以看出，在5—9月的降雨事件中，由不同覆被小區(qū)的徑流系數(shù)分布特征可知，在同一區(qū)域，喀斯特地區(qū)不同覆被下地表的產(chǎn)流差異明顯，且撂荒草地的徑流系數(shù)最大，明顯大于經(jīng)果林和水保林2種覆被條件；經(jīng)果林小區(qū)之間的徑流系數(shù)差異不大；2個(gè)水保林小區(qū)的產(chǎn)流系數(shù)均明顯低于經(jīng)果林和撂荒草地，顯示了很好的減流作用。

為研究不同覆被下降雨與坡面產(chǎn)流關(guān)系，筆者對(duì)不同覆被下降雨量與產(chǎn)流進(jìn)行線性擬合。從圖4 可以看出，總體上，不同覆被下徑流小區(qū)的產(chǎn)流量與降雨量之間有較好的線性正相關(guān)性。相關(guān)系數(shù)均超過(guò)了0.5，經(jīng)果林為0.763、水保林為 0.628、撂荒草地為0.610。另外，降雨量在0～25 mm之間時(shí)，各覆被下徑流小區(qū)的產(chǎn)流與降雨量可能存在負(fù)相關(guān)性，這可能與當(dāng)?shù)囟嚓幱赀B綿，而降雨強(qiáng)度小的因素有關(guān)。

2.3.2 地表產(chǎn)沙與降雨特征關(guān)系 從圖5可以看出，在研究區(qū)5—9月降雨事件中，不同覆被小區(qū)單位面積沖刷量排序?yàn)榱袒牟莸?經(jīng)果林3>經(jīng)果林1>經(jīng)果林2>水保林2>水保林1。經(jīng)果林樹(shù)冠冠幅大，有利于降雨截留和緩沖雨滴，減小降雨對(duì)地表的擊濺侵蝕?？赡苡捎诮?jīng)果林葉片上的雨滴容易形成滴流，所以減沙作用沒(méi)有以灌木為主的水保林明顯。在強(qiáng)降雨中，水保林的單位面積沖刷量只有撂荒草地的1/4左右。因此，以灌木為主的水保林顯示出了很好的減沙作用。

由圖6可以看出，與降雨量和產(chǎn)流關(guān)系不同，降雨量與產(chǎn)沙量的線性關(guān)系不明顯。其他研究也指出，降雨量與產(chǎn)沙量并非簡(jiǎn)單的線性關(guān)系。產(chǎn)沙主要受短期最大降雨強(qiáng)度和植被的郁閉度影響。在5月份，可能由于植被的冠蓋還沒(méi)有完全發(fā)育，而短期的降雨強(qiáng)度較大，所以導(dǎo)致在小降雨量的情況下產(chǎn)生了較大的產(chǎn)沙量。

2.4 地表產(chǎn)流與產(chǎn)沙的關(guān)系

降雨量直接影響坡面徑流量，而徑流是地表水土流失的主要載體。從表2、圖7可以看出，除撂荒草地和水保林2外，其他徑流小區(qū)的產(chǎn)沙量和產(chǎn)流量擬合度較好(r2>0.5)。其中，經(jīng)果林2的產(chǎn)沙產(chǎn)流線型正相關(guān)擬合關(guān)系最好(r2=0.833)，而水保林2的擬合關(guān)系最差(r2=0.107)。隨著植被覆蓋率的增加，泥沙總量就越低，擬合曲線及決定系數(shù)也就越高。但是，在同樣覆蓋率下的產(chǎn)沙總量也有所不同，這有可能與不同土地利用方式及種植的植被有關(guān)。

3 討論

土壤侵蝕變化的衡量指標(biāo)主要包括產(chǎn)流強(qiáng)度、侵蝕泥沙量、徑流系數(shù)等[13]。通過(guò)6種試驗(yàn)小區(qū)、3種不同覆被類型和降雨量的對(duì)比研究，分析了貴州關(guān)嶺布依族苗族自治縣喀斯特地區(qū)不同覆被下徑流小區(qū)的產(chǎn)沙產(chǎn)流特征。研究表明，在不同的降雨量下，土壤的產(chǎn)流時(shí)間隨之增加，在相同坡度下，降雨量與產(chǎn)流量呈正比。研究區(qū)年降雨總量1 039.7 mm，侵蝕性降雨總量占50%左右，降雨歷時(shí)183 d，降雨多以小雨形式產(chǎn)生。相比其他地區(qū)，研究區(qū)的降雨量總量較大，降雨天數(shù)較多，但是降雨強(qiáng)度較小，因此導(dǎo)致產(chǎn)流量相對(duì)較大。研究還表明，降雨強(qiáng)度與產(chǎn)沙量呈正比；在大雨級(jí)別以上，產(chǎn)沙量也會(huì)隨之增加，這與其他研究得出的結(jié)論相同。研究區(qū)在同一坡度下，產(chǎn)流強(qiáng)度和入滲強(qiáng)度都隨著降雨強(qiáng)度的增加而增加，但是在研究區(qū)內(nèi)，土壤為黃色石灰土，容易隨著流水而流失。

表2不同覆被類型累積產(chǎn)沙量與累積產(chǎn)流量的關(guān)系

注：y為累積產(chǎn)沙量，g；x為累積產(chǎn)流量，L。

不同土地覆蓋變化對(duì)徑流的影響研究主要集中在年徑流量上[14]。張淑敏研究長(zhǎng)江流域大面積森林流域表明，森林覆蓋率高的流域比森林覆蓋率低的流域、有林地比無(wú)林地流域河川年徑流量均有所增加[15]。這更有力地證明了本研究的結(jié)論：各徑流小區(qū)坡面產(chǎn)流產(chǎn)沙特征與降雨的年內(nèi)分布密切相關(guān)，相關(guān)性以經(jīng)果林最大，水保林次之，撂荒草地最低。

不同植被的降雨產(chǎn)流產(chǎn)沙特征差異明顯，以灌木為主的水保林減沙減流效果明顯，其坡面產(chǎn)流產(chǎn)沙明顯低于經(jīng)果林和撂荒地，由此可知，植被對(duì)地表徑流的影響是由冠層、枯枝落葉層和土壤層的綜合效能決定的[16]。覆被的不同可能對(duì)林冠截留、凋落物持水性、土壤理化性質(zhì)等產(chǎn)生影響，進(jìn)而直接影響地表產(chǎn)流產(chǎn)沙過(guò)程[17-18]。陳奇伯等研究表明，灌木林和天然次生林均有較好的調(diào)節(jié)徑流和減少土壤流失的作用[19]。張文源等研究表明，高覆蓋的灌草優(yōu)于喬木，可大幅度地抑制土壤侵蝕的發(fā)生[20]。研究表明，不同植被的降雨產(chǎn)流產(chǎn)沙特征差異顯著，以灌木為主的水保林減沙減流效果明顯，其坡面產(chǎn)流產(chǎn)沙明顯低于經(jīng)果林和撂荒地。這與陳奇伯等的研究結(jié)果[19-20]類似。

4 結(jié)論

研究區(qū)年降水總量1 039.7 mm、降雨歷時(shí)183 d，降雨多以小雨形式產(chǎn)生，侵蝕性降雨總量占50%左右；各徑流小區(qū)坡面產(chǎn)流產(chǎn)沙特征與降雨的年內(nèi)分布密切相關(guān)，相關(guān)性以經(jīng)果林最大，水保林次之，撂荒草地最低；不同植被的降雨產(chǎn)流產(chǎn)沙特征差異明顯，以灌木為主的水保林減沙減流效果明顯，其坡面產(chǎn)流產(chǎn)沙明顯低于經(jīng)果林和撂荒地；不同覆被下各小區(qū)產(chǎn)沙與產(chǎn)流大多具有較好的線性相關(guān)性。

參考文獻(xiàn)：

[1]Legrand H E. Hydrological and ecological problems of karst regions：hydrological actions on limestone regions cause distinctive ecological problems[J]. Science，1973，179(4076)：859-864.

[2]Wardle D A，Walker L R，Bardgett R D，et al. Ecosystem properties and forest decline in contrasting long-term chronosequences[J]. Science，2004，305(5683)：509-513.

[3]張 偉，陳洪松，王克林. 典型喀斯特峰叢洼地坡面土壤養(yǎng)分空間變異性研究[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2008，24(1)：69-71.

[4]Xu E Q，Zhang H Q，Li M X. Mining spatial information to investigate the evolution of karst rocky desertification and its human driving forces in Changshun，China[J]. Science of the Total Environment，2013(458/459/460)：419-426.

[5]彭旭東，戴全厚，楊 智，等. 喀斯特山地石漠化過(guò)程中地表地下侵蝕產(chǎn)沙特征[J]. 土壤學(xué)報(bào)，2016，53(5)：1237-1248.

[6]Reaney S M，Bracken L J，Kirkby M J.Use of the connectivity of runoff model(CRUM) to investigate the influence of storm characteristics on runoff generation and connectivity in semi-arid areas[J]. Hydrological Processes，2007，21(7)：894-906.

[7]陳曉安. 黃土丘陵溝壑區(qū)坡面土壤侵蝕規(guī)律與坡面侵蝕經(jīng)驗(yàn)?zāi)Ｐ偷难芯縖D]. 武漢：華中農(nóng)業(yè)大學(xué)，2010.

[8]謝林妤，白玉潔，張風(fēng)寶，等. 沙層厚度和粒徑組成對(duì)覆沙黃土坡面產(chǎn)流產(chǎn)沙的影響[J]. 土壤學(xué)報(bào)，2017，54(1)：60-72.

[9]李 坤，姚文藝，肖培青，等. 植被對(duì)土壤入滲和地表產(chǎn)流過(guò)程的影響研究進(jìn)展[J]. 中國(guó)水土保持，2017(3)：27-30.

[10]涂利華，謝財(cái)永，胡庭興，等. 華西雨屏區(qū)幾種牧草的水土保持能力研究[J]. 水土保持學(xué)報(bào)，2005，19(5)：35-38.

[11]蘇忖安，松同清，王文麗，等. 模擬降雨對(duì)亞熱帶闊葉林土壤坡面產(chǎn)沙產(chǎn)流及養(yǎng)分流失的影響[J]. 水土保持學(xué)報(bào)，2016(4)：25-32.

[12]王 輝，王全九，邵明安. 表層土壤容重對(duì)黃土坡面養(yǎng)分隨徑流遷移的影響[J]. 水土保持學(xué)報(bào)，2007，21(3)：10-13，18.

[13]汪 軍，王德建，張 剛，等. 連續(xù)全量秸稈還田與氮肥用量對(duì)農(nóng)田土壤養(yǎng)分的影響[J]. 水土保持學(xué)報(bào)，2010，24(5)：40-44.

[14]趙 輝，郭索彥，解明曙，等. 南方花崗巖紅壤區(qū)不同土地利用類型坡地產(chǎn)流與侵蝕產(chǎn)沙研究[J]. 水土保持通報(bào)，2008，28(2)：6-10.

[15]張淑敏. 基于森林作用的流域降雨徑流模型研究[D]. 泰安：山東農(nóng)業(yè)大學(xué)，2012.

[16]紀(jì)啟芳，張興奇，張科利，等. 貴州省喀斯特地區(qū)坡面產(chǎn)流產(chǎn)沙特征[J]. 水土保持研究，2012，19(4)：1-5.

[17]稅 偉，白劍平，簡(jiǎn)小枚，等. 若爾蓋沙化草地恢復(fù)過(guò)程中土壤特性及水源涵養(yǎng)功能[J]. 生態(tài)學(xué)報(bào)，2017，37(1)：277-285.

[18]常宗強(qiáng)，王金葉，常學(xué)向，等. 祁連山水源涵養(yǎng)林枯枝落葉層水文生態(tài)功能[J]. 西北林學(xué)院學(xué)報(bào)，2001，16(增刊)：8-13.

[19]陳奇伯，寸玉康，劉芝芹，等. 滇西高原不同地類坡面產(chǎn)流產(chǎn)沙規(guī)律研究[J]. 水土保持研究，2005，12(2)：71-73.

[20]張文源，王百田，楊光檄，等. 喀斯特黃壤區(qū)侵蝕性降雨及產(chǎn)沙特征分析[J]. 生態(tài)環(huán)境學(xué)報(bào)，2014，23(11)：1776-1782.