羅旭玲,白曉永,*,譚 秋,陳 歡,冉 晨,習(xí)慧鵬

1 貴州師范大學(xué) 地理與環(huán)境科學(xué)學(xué)院, 貴陽(yáng) 550025 2 中國(guó)科學(xué)院地球化學(xué)研究所, 貴陽(yáng) 550081

西南典型喀斯特地區(qū)是一種脆弱的生態(tài)環(huán)境,土壤侵蝕和石漠化已成為制約西南喀斯特山區(qū)社會(huì)經(jīng)濟(jì)可持續(xù)發(fā)展的問(wèn)題[1-7]。喀斯特山地土壤流失具有地上地下重疊、地表產(chǎn)流產(chǎn)沙少等特點(diǎn)[8],研究喀斯特地區(qū)石漠化與土壤侵蝕變得極其重要。

為了更好的防治石漠化與土壤侵蝕的危害,當(dāng)前對(duì)特定喀斯特地貌背景下的研究有很多。例如：何永彬等[9]對(duì)西南喀斯特山地土壤侵蝕研究探討；鄧艷等[10]西南巖溶石漠化綜合治理水-土-植被關(guān)鍵技術(shù)進(jìn)展與建議；戴全厚等[11]對(duì)西南喀斯特石漠化與水土流失研究進(jìn)展的分析；李瑞玲等[12]對(duì)貴州巖溶地區(qū)巖性與土地石漠化的相關(guān)分析,等研究表明石漠化治理與土壤侵蝕息息相關(guān)。張信寶等[13]對(duì)貴州石漠化空間分布與喀斯特、巖性、降水和人口密度的關(guān)系研究；白曉永等[14]對(duì)貴州碳酸鹽巖巖性基地對(duì)土地石漠化時(shí)空演變的控制；李陽(yáng)兵等[1]對(duì)中國(guó)典型石漠化地區(qū)土地利用與石漠化的關(guān)系的研究,以上研究均取得較好的成果。喀斯特地區(qū)土壤侵蝕和石漠化受坡度影響[15],由石漠化[16]與土壤侵蝕[17]的概念可知,影響因子不同對(duì)兩者的作用存在差異。

因此,本文在借鑒前人的研究基礎(chǔ)上,以沿河縣為例,利用ArcGIS軟件對(duì)土壤侵蝕、石漠化、巖性三者進(jìn)行疊加分析,對(duì)石漠化與土壤侵蝕空間分布、不同巖性下的石漠化與土壤侵蝕分布狀況、土壤侵蝕與石漠化相關(guān)性進(jìn)行分析,以期獲得較好的分析結(jié)果,為研究區(qū)石漠化與土壤侵蝕防治工作提供指導(dǎo)依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

沿河土家族自治縣位于貴州省黔東北(圖1),隸屬銅仁市,108°03′49″—108°37′5″E,28°12′45″—29°05′23″N。縣內(nèi)屬于中亞熱帶季風(fēng)溫潤(rùn)氣候區(qū),年均溫13—18℃,年降水量1050—1220 mm,年日照1100—1400 h。巖性以碳酸鹽巖類(石灰?guī)r、白云質(zhì)灰?guī)r、白云巖等)為主,碳酸鹽巖類與碎屑巖類交錯(cuò)分布,巖性組合較為復(fù)雜[18]。在全縣面積中,山地占69.9%,丘陵占27%,槽壩占3.1%,屬于典型的巖溶山區(qū)[19-20]。沿河縣水土流失得到有效控制[21]。位于烏江流域下游的沿河縣比流域上游石漠化等級(jí)較高,但近幾年來(lái)石漠化面積有減少,石漠化處于好轉(zhuǎn)趨勢(shì)[21]。

圖1 研究區(qū)位置圖Fig.1 Location map of the study area

1.2 研究方法

1.2.1 數(shù)據(jù)源

DEM用于土壤侵蝕LS因子計(jì)算,Landsat OLI 2015年遙感影像用于土地利用解譯,NDVI用于C因子計(jì)算,均來(lái)自地理空間數(shù)據(jù)云。用于R因子計(jì)算的降雨數(shù)據(jù),來(lái)自中國(guó)氣象數(shù)據(jù)共享網(wǎng),用于K因子計(jì)算的沙粒、粉粒、粘粒、有機(jī)質(zhì)成分,來(lái)自HWSD數(shù)據(jù)庫(kù)。土壤侵蝕所有因子具有統(tǒng)一的投影(WGS_1984_UTM_Zone_48N)和分辨率(30 m)。沿河縣巖性圖和石漠化分布圖源于喀斯特科學(xué)數(shù)據(jù)中心。數(shù)據(jù)詳細(xì)來(lái)源及網(wǎng)站見(jiàn)表1。

表1 主要數(shù)據(jù)來(lái)源

1.2.2 土壤侵蝕計(jì)算

研究計(jì)算公式為Wischemier和Smith于1958年提出[22],經(jīng)美國(guó)農(nóng)業(yè)部修正的通用土壤流失方程(USLE)。表達(dá)式如下：

A=R×K×LS×C×P

(1)

式中,A為年土壤侵蝕量,R為降雨侵蝕力因子,K為土壤侵蝕力因子,LS為地形因子,C植被與作物管理因子,P土壤保持措施因子。

(1) 降雨侵蝕力因子

研究采用周伏建、黃炎和[23-24]等人根據(jù)南方實(shí)測(cè)數(shù)據(jù)提出的R值計(jì)算式：

(2)

式中,Pi為月降水量,由氣象數(shù)據(jù)共享網(wǎng)下載的貴州省19個(gè)氣象站點(diǎn)的月值數(shù)據(jù),利用EXCEL表格統(tǒng)計(jì)獲得,借助ArcGIS地圖代數(shù)幾何計(jì)算工具,按照公式(2)計(jì)算得每個(gè)站點(diǎn)的年降雨侵蝕力R,進(jìn)行Kriging插值得到R值的空間分布圖,對(duì)R重采樣(30 m)和裁剪。

(2) 土壤可蝕性因子

土壤可蝕性因子K值是經(jīng)過(guò)試驗(yàn)獲得的定量數(shù)值,通常以標(biāo)準(zhǔn)小區(qū)單位降水侵蝕力所造成的土壤流失量獲得,但在該區(qū)域大規(guī)模布設(shè)天然小區(qū)很難做到,因此本研究以Williams等[25]發(fā)展的僅由土壤顆粒和土壤有機(jī)碳組成數(shù)據(jù)來(lái)估算K值,其表達(dá)式為：

(3)

(3) 地形坡長(zhǎng)坡度因子

地形地貌對(duì)土壤侵蝕有著重要的影響,區(qū)域不同,USLE模式中斜坡長(zhǎng)度因子(L)和坡度因子(S)有變化。本文根據(jù)水流功率理論和物理過(guò)程,以及充分考慮了徑流和水文過(guò)程的影響,選擇以下公式[26-27]進(jìn)行計(jì)算LS。

L=(λ/22.13)m

(5)

λ=flowacc×cellsize

(6)

m=n/(1+n)

(7)

n=(sinθ/0.0896)/(3.0×sin0.8θ+0.56)

(8)

(9)

式中,L為坡長(zhǎng)因子；m為可變坡長(zhǎng)指數(shù)；22.13為標(biāo)準(zhǔn)小區(qū)坡長(zhǎng)；S為坡度因子；θ為坡度。

(4) 植被與作物管理因子

采用Van der Knijff M提出的利用NDVI確定C因子的方法[28-29]。兩者關(guān)系為：

(10)

式中,α=2,β=1,NDVI為歸一化植被指數(shù)。

(5) 土壤保持措施因子

水土保持措施因子的值在0—1之間。貴州喀斯特山區(qū),海拔落差較大,坡度較陡,水土保持措施常與土地利用類型緊密相關(guān)。所以,把水體、建設(shè)用地和裸地賦值為0,水田0.15,旱地0.5,其他用地0.8,有林地、灌木用地、草地為1。

2 結(jié)果與分析

2.1 石漠化與土壤侵蝕空間分布

研究區(qū)石漠化和土壤侵蝕空間分布情況,通過(guò)ArcGIS空間分析得出。喀斯特面積為623.04 km2,占總面積的25.24%；非石漠化面積為1845.76 km2,占總面積的74.76%(表2)。沿河縣石漠化面積達(dá)284.44 km2,占全縣國(guó)土面積的11.52%。已發(fā)生石漠化現(xiàn)象區(qū)域,中度石漠化面積最大,為135.87 km2,占石漠化面積的47.77%；輕度石漠化與重度石漠化次之,面積分別為126.33,21.58 km2,分別占石漠化面積的44.41%,7.59%；極重度石漠化分布面積最小,為0.66 km2,僅占0.23%。分析得出,沿河縣已發(fā)生石漠化等級(jí)以輕度、中度和重度為主。石漠化變化面積整體呈現(xiàn)隨等級(jí)的增加而減少。

表2 沿河縣喀斯特與非喀斯特區(qū)石漠化分布情況

非石漠化=無(wú)石漠化+非喀斯特;PKRD,潛在石漠化 Potential Rock Desertification；LKRD,輕度石漠化 Light Rock Desertification；MKRA,中度石漠化 Medium Rock Desertification；SKRD,重度石漠化 Severe Rock Desertification；ESKRD,極重度石漠化 Extremely Severe Rock Desertification；NKRD,非石漠化區(qū) Non-Rock Desertification Area

土壤侵蝕面積1838.97 km2,占全縣總面積的74.48%,微度侵蝕面積為629.83 km2,占全縣總面積的25.51%(表3)。已經(jīng)發(fā)生土壤侵蝕現(xiàn)象區(qū)域,輕度侵蝕面積分布最廣,1637.69 km2,占土壤侵蝕面積的89.05%；中度侵蝕和強(qiáng)度侵蝕面積次之,分別為154.15,28.51 km2,占土壤侵蝕面積的8.38%,1.55%,極強(qiáng)度侵蝕面積為13.63 km2,占侵蝕面積0.74%,劇烈侵蝕面積最少,為4.99 km2,僅占土壤侵蝕面積的0.27%。分析數(shù)據(jù)可知,沿河縣境內(nèi)土壤侵蝕較嚴(yán)重,以輕度和中度侵蝕為主,隨侵蝕等級(jí)的升高,其面積越小。

表3 沿河縣水土流失情況

PSE,微度侵蝕 Potential soil erosion；LSE,輕度侵蝕 Light soil erosion；MSE,中度侵蝕 Medium soil erosion；SSE,強(qiáng)度侵蝕 Strength soil erosion；ESSE,極強(qiáng)度侵蝕 Extremely strength soil erosion；FSE,劇烈侵蝕 Fierce soil erosion

2.2 不同巖性下的石漠化與土壤侵蝕分布狀況

對(duì)巖性、石漠化和土壤侵蝕在ArcGIS中進(jìn)行疊加分析,統(tǒng)計(jì)得到不同巖性所占全縣、石漠化和土壤侵蝕面積比例(表4)。非碳酸鹽巖分別占石漠化與土壤侵蝕面積的5.08%,14.75%,土壤侵蝕在碳酸鹽巖上分布較廣,但石漠化分布卻是最少?；?guī)r與碎屑巖占石漠化面積7.50%,占土壤侵蝕面積9.38%?；?guī)r與白云巖混合面積最少,占沿河縣面積8.18%,占土壤侵蝕面積8.69%,是土壤侵蝕分布最少的巖性,占石漠化面積11.56%,該區(qū)人類活動(dòng)較強(qiáng),導(dǎo)致石漠化面積較大,但因地表較緩,所以土壤侵蝕較小。灰?guī)r夾碎屑巖互層面積最大,占沿河縣面積30.66%,分別占石漠化與土壤侵蝕面積的36.53%和30.90%,巖性基底大,加上人類活動(dòng)與自然環(huán)境的綜合影響,因此,石漠化和土壤侵蝕面積分布均為最大。連續(xù)性石灰?guī)r占石漠化面積21.46%,占土壤侵蝕面積19.55%,連續(xù)性石灰?guī)r分別占石漠化與土壤侵蝕面積的17.86%和16.74%。由分析可知,同一巖性上石漠化與土壤侵蝕發(fā)生不具有單一的相關(guān)性,石漠化發(fā)生比例小的地方土壤侵蝕并不是最小,說(shuō)明石漠化與土壤侵蝕的發(fā)生有著一定的滯后性。

2.3 石漠化與土壤侵蝕相關(guān)性分析

2.3.1 不同巖性下各石漠化等級(jí)區(qū)的土壤侵蝕

依據(jù)沿河縣巖性、石漠化和土壤侵蝕疊加分析得到的數(shù)據(jù),在EXCLE表格中統(tǒng)計(jì)得到表5和圖2。土壤侵蝕主要發(fā)生在潛在、輕度和中度石漠化區(qū)(表5,圖2),在非碳酸鹽巖、灰?guī)r夾碎屑巖、灰?guī)r與碎屑巖互層和連續(xù)性石灰?guī)r中,潛在石漠化區(qū)侵蝕面積最大,分別為18.40,18.53,112.08,56.01 km2。因?yàn)橹卸仁瘏^(qū)植被較發(fā)育,土被基本連續(xù),土壤侵蝕明顯,所以,中度石漠化區(qū)侵蝕面積比輕度石漠化區(qū)要大。極重度石漠化區(qū)最少,其中非碳酸鹽上極重度石漠化區(qū)侵蝕面積0.06 km2最小,土壤侵蝕與石漠化等級(jí)在中度至極重度石漠化區(qū)間呈負(fù)相關(guān),因?yàn)槭燃?jí)越高,巖石裸露率越大,表層土壤越薄,會(huì)出現(xiàn)少有或無(wú)土可流。灰?guī)r與白云巖混合和連續(xù)性白云巖上,在潛在石漠化區(qū)面積最大,分別為38.91 km2與94.44 km2,在極重度石漠化區(qū)最小,分別為0.17 km2與0.07 km2。土壤侵蝕面積隨著石漠化等級(jí)的增加而減少,存在負(fù)相關(guān)。

表4 不同巖性下土壤侵蝕與石漠化面積比例/%

NCRC,非碳酸鹽巖組合Non-carbonate rock combination；LCRC,灰?guī)r夾碎屑巖組合 Limestone clastic rock combination；MCLD,灰?guī)r與白云巖混合組合 Mixed combination of limestone and dolomite；LLRI,灰?guī)r與碎屑巖互層 Limestone and lastic rock interbed；CD,連續(xù)性白云巖 Continuous dolomite；CL,連續(xù)性石灰?guī)r Continuous limestone

表5 不同巖性下各石漠化等級(jí)區(qū)已發(fā)生土壤侵蝕面積/km2

圖2 不同巖性下各石漠化等級(jí)區(qū)域的土壤侵蝕 Fig.2 Soil erosion in different rocky desertification areas under different lithologiesNCRC,非碳酸鹽巖組合Non-carbonate rock combination；LCRC,灰?guī)r夾碎屑巖組合 Limestone clastic rock combination；MCLD,灰?guī)r與白云巖混合組合 Mixed combination of limestone and dolomite；LLRI,灰?guī)r與碎屑巖互層 Limestone and lastic rock interbed；CD,連續(xù)性白云巖 Continuous dolomite；CL,連續(xù)性石灰?guī)r Continuous limestone,PSE,微度侵蝕 Potential soil erosion；LSE,輕度侵蝕 Light soil erosion；MSE,中度侵蝕 Medium soil erosion；SSE,強(qiáng)度侵蝕 Strength soil erosion；ESSE,極強(qiáng)度侵蝕 Extremely strength soil erosion；FSE,劇烈侵蝕 Fierce soil erosion.潛在石漠化 Potential rock desertification；LKRD,輕度石漠化 Light rock desertification；MKRA,中度石漠化 Medium rock desertification；SKRD,重度石漠化 Severe rock desertification；ESKRD,極重度石漠化 Extremely severe rock desertification；NKRD,非石漠化區(qū) Non-rock desertification area

2.3.2 不同巖性下各土壤侵蝕等級(jí)區(qū)的石漠化

通過(guò)ArcGIS疊置分析,導(dǎo)出面積進(jìn)行Excel統(tǒng)計(jì)得到表6和圖3。沿河縣石漠化主要發(fā)生在微度和輕度侵蝕區(qū)。各巖性上不同土壤侵蝕等級(jí)區(qū)石漠化主要發(fā)生在微度、輕度和中度。石漠化發(fā)生主要分布在灰?guī)r與碎屑巖互層、連續(xù)性白云巖、連續(xù)性石灰?guī)r上,其面積發(fā)生最大值在輕度侵蝕區(qū),分別為72.20、42.88 km2和30.96 km2,非碳酸鹽巖上分布最少。各巖性中土壤侵蝕等級(jí)與石漠化在微度侵蝕與輕度侵蝕內(nèi)呈正相關(guān),在輕度與劇烈侵蝕內(nèi)呈負(fù)相關(guān)。石漠化治理是預(yù)防潛在石漠化向石漠化的轉(zhuǎn)移,以及對(duì)已發(fā)生石漠化的改善。石漠化的改善有利于固土持水,所以輕度侵蝕區(qū)石漠化面積最大。由表6可知,各土壤侵蝕等級(jí)區(qū),石漠化面積最大的巖性是灰?guī)r與碎屑巖互層,灰?guī)r與碎屑巖互層主要是人類居住聚集的巖層,地表人類活動(dòng)(耕種,開(kāi)采,放牧等)強(qiáng)烈,再加上灰?guī)r與碎屑巖互層本地較大,所以,這可能是石漠化較高的原因。石漠化面積最小的巖性為非碳酸鹽巖。研究區(qū)非碳酸鹽巖分布不是最少,但石漠化分布是最少的,說(shuō)明石漠化主要發(fā)生在碳酸鹽巖上。

3 結(jié)論

(1)沿河縣石漠化面積為284.44 km2,占全縣國(guó)土面積的11.52%。已發(fā)生石漠化現(xiàn)象區(qū)域,以輕度、中度、強(qiáng)度較顯著,分別占石漠化面積的44.41%,47.77%,7.59%,石漠化現(xiàn)象較輕。土壤侵蝕面積為1838.97 km2,占全縣總面積的74.48%,微度侵蝕面積為629.83 km2,占全縣總面積的25.51%。已經(jīng)發(fā)生土壤侵蝕現(xiàn)象區(qū)域,以輕度、中度、強(qiáng)度侵蝕為主,分別占沿河縣土壤侵蝕面積的89.05%,8.38%,1.55%,土壤侵蝕程度較嚴(yán)重。

表6 不同巖性下各土壤侵蝕等級(jí)區(qū)已發(fā)生石漠化面積/km2

圖3 不同巖性下各土壤侵蝕等級(jí)區(qū)域的石漠化Fig.3 Rocky desertification in different soil erosion grades under different lithologies

(2)在不同石漠化等級(jí)區(qū),潛在、輕度和中度石漠化區(qū)土壤侵蝕發(fā)生面積最大,重度和極重度石漠化區(qū),土壤侵蝕發(fā)生面積最小,因?yàn)樵缙趧×彝寥狼治g發(fā)生,導(dǎo)致石漠化嚴(yán)重區(qū)域,地表可流物質(zhì)減少。在不同的土壤侵蝕區(qū),石漠化主要發(fā)生在微度和輕度侵蝕區(qū),輕度侵蝕區(qū)石漠化發(fā)生面積最大。當(dāng)土壤侵蝕嚴(yán)重區(qū)域石漠化發(fā)生更加劇烈時(shí),地表土層薄、巖石裸露程度大、植被覆蓋率低,土壤侵蝕會(huì)逐漸減少。所以,石漠化與土壤侵蝕的發(fā)生是受多種因素綜合作用的結(jié)果。

(3)各巖性上,石漠化與土壤侵蝕的發(fā)生是相互影響的,兩者之間沒(méi)有完全單一正相關(guān)和負(fù)相關(guān)關(guān)系。非碳酸鹽巖、灰?guī)r夾碎屑巖、灰?guī)r與碎屑巖互層和連續(xù)性石灰?guī)r上,石漠化等級(jí)與土壤侵蝕在潛在和輕度、中度和極重度石漠化內(nèi)呈負(fù)相關(guān),在輕度和中度石漠化內(nèi)呈正相關(guān),土壤侵蝕等級(jí)與石漠化在微度和輕度侵蝕內(nèi)呈正相關(guān),在輕度和劇烈侵蝕內(nèi)呈負(fù)相關(guān)?；?guī)r與白云巖混合和連續(xù)性白云巖上,石漠化等級(jí)與土壤侵蝕存在單一的負(fù)相關(guān)關(guān)系,土壤侵蝕等級(jí)與石漠化在土壤侵蝕等級(jí)與石漠化在微度和輕度侵蝕內(nèi)呈正相關(guān),在輕度和劇烈侵蝕內(nèi)呈負(fù)相關(guān)。