陳起偉,熊康寧,蘭安軍
( 1.貴州師范學院 地理與旅游學院,貴州 貴陽 550003;2.貴州師范大學 中國南方喀斯特研究院,貴州 貴陽 550001 )
中國南方地區連片分布著世界最大面積的碳酸鹽巖[1-3],石漠化問題已經到了很嚴重的程度。中國南方巖溶地區的石漠化問題得到了國家的高度重視,國家投入了大量的資金,實施了“長治工程”、“珠治試點工程”、石漠化綜合治理和鞏固退耕還林成果專項資金建設項目等一系列工程。但由于缺乏對巖溶地區特殊水土流失及石漠化過程和機理的深入了解,沒有考慮到流域處于河段的位置、地貌、基巖等環境差異性,千篇一律采用相同的治理模式治理,許多地區所采取的治理措施、方法針對性不強,效果并不十分明顯[4]。由于不同地貌類型、不同河段的流域的石漠化發生背景不同,石漠化特征和演變規律存在差異[5][6],所以研究流域內不同河段石漠化及演變特征的差異性對于更好地有針對性地采取不同的石漠化治理模式具有重要意義。沿河淇灘小流域、遵義龍壩小流域和畢節鴨池小流域分別代表了貴州烏江流域下游、中游和上游三種不同的流域河段和地貌類型,分析其石漠化的特征和變化規律,可反映貴州典型流域不同河段巖溶石漠化的特征和演變規律。
貴州沿河淇灘小流域、遵義龍壩小流域和畢節鴨池小流域分別代表了貴州烏江流域下游、中游和上游三種不同的流域河段及中低山峽谷、中低山溝谷盆地和巖溶高原山地三種不同的地貌類型。綜合分析當前各種石漠化調查方法,結合研究需要,石漠化調查與分析利用遙感技術進行。為了科學地反映烏江流域不同河段石漠化特征和變化情況,選取IRS-P6融合數據(5.8m分辨率)和ALOS(10m分辨率)開展該地區石漠化類型的對比調查。其中,遵義龍壩小流域使用影像為IRS-P6融合數據,獲取時間為2005年10月和2007年11月;沿河淇灘和畢節鴨池小流域使用ALOS影像數據,時間為2005年10月和2007年11月。石漠化分類參照采用喀斯特石漠化的六級分類體系[7],結合研究區域各種石漠化特點和遙感信息的可判性分為非喀斯特、非石漠化和石漠化等 6種類型。對遙感圖像進行幾何校正和地理校正與配準,鑲嵌和裁切,并進行適當的增強處理,收集流域土地利用圖(1:10000)、數字高程模型(1:10000)、地質圖、植被類型圖和野外考察數據和資料等進行輔助分類,采用計算機自動監督分類與目視解譯相結合的方式進行圖像解譯,分類后再對圖像進行野外實地核查校正。分類精度采用隨機抽樣的方法,將分類數據與原始圖像、實地考察數據對比,2005年和2007年二期數據的總分類精度分別為89.2%和90.3%。根據三個不同河段兩期石漠化圖數據庫,對比分析烏江流域不同河段石漠化差異及其變化特征。
沿河淇灘小流域地處沿河縣淇灘鎮境內,屬烏江流域下游段,出露地層主要為寒武系及奧陶系的白云巖、白云質灰巖,局部地帶分布著志留系碎屑巖。沿河淇灘小流域巖溶地貌發育,受烏江水系切割影響,境內海拔高差懸殊,地貌多以中低山峽谷、溝谷為主,地表破碎,山高坡陡。流域屬亞熱帶季風濕潤氣候區,多年平均降雨量為1136mm,多年平均氣溫為17.5℃。土壤以黃壤為主,植被為天然亞熱帶針葉林和闊葉林混交植物,樹種為松、杉、柏、楓樹等。
遵義龍壩小流域位于遵義縣龍坪鎮境內,處于烏江流域中游區,流域地形起伏較大,山地、丘陵、河谷縱橫排列,相間交錯,平均海拔為920m。巖性以灰巖、泥灰巖、白云質灰巖為主。遵義龍壩小流域氣候溫和,雨量充沛,冬無嚴寒,夏無酷暑,年平均降雨量為1036mm,年均氣溫為14.6℃。土壤主要分布有石灰土、水稻土、黃壤等。遵義龍壩小流域屬闊葉林和針葉林混交地帶,針葉林居多,現有天然植被較少,主要有馬尾松、柏樹、灌木、刺叢、櫟樹等。
畢節鴨池小流域位于畢節市東南部鴨池鎮,處于烏江流域上游段,屬巖溶高原山地地貌類型。流域地勢起伏大,海拔在 1400m~1750m之間,相對高差達350m,區內巖石以碳酸巖類的石灰巖及侏羅紀紫色砂頁巖為主。受成土母巖影響,土壤主要為黃壤和紫砂土,區域氣候溫涼濕潤,年均氣溫為14.0℃,年降水量為863mm。植被為亞熱帶常綠闊葉林,原生植被多被破壞,現為次生林,以馬尾松、柏樹、灌木為主。沿河淇灘小流域、遵義龍壩小流域和畢節鴨池小流域在烏江流域的位置如圖1所示。

圖1 淇灘、龍壩和鴨池小流域在烏江流域的位置Fig.1 Location of Yachi, Longba and Qitan research area in Wujiang River Basin
根據石漠化分級標準,利用研究區兩期遙感影像,經遙感解譯和野外實地調查,繪制淇灘、龍壩和鴨池小流域2005年和2007年石漠化分布圖,統計生成各種等級石漠化面積及變化,見表1。因各小流域面積大小不同,不能直接對比上述小流域各種石漠化的面積和變化面積,于是從石漠化景觀差異、發生率、各種石漠化占巖溶面積的比例和石漠化年變化率等方面對比分析三個小流域石漠化的特征和變化規律,以揭示烏江流域不同河段石漠化差異性。
采用圖斑平均面積和景觀破碎度對小流域景觀結構進行定量化描述。石漠化圖斑數量和斑塊平均面積反映了各種石漠化的內部構成,經統計分析,畢節鴨池小流域各種石漠化平均斑塊面積較大,平均為3.6畝,而遵義龍壩和沿河淇灘小流域平均斑塊面積分別為1.7畝和1.9畝,相差不大。2005年~2007年間,畢節鴨池小流域平均斑塊變化較大,遵義龍壩小流域石漠化斑塊變化面積最小。畢節鴨池小流域石漠化分布較連續,變化也較快。景觀破碎度表征石漠化景觀被分割的破碎程度,反映了各種石漠化景觀空間結構的復雜性,在一定程度上反映了人類對景觀的干擾程度。畢節鴨池小流域強度和中度石漠化圖斑的破碎度最小,而遵義小流域輕度石漠化圖斑破碎度最大。遵義龍壩和沿河淇灘小流域各種石漠化類型的破碎度指數相差不大,而畢節鴨池小流域破碎度指數相差較大,除畢節鴨池小流域的中度和強度石漠化區外,三個小流域破碎度指數變化曲線較和緩,見圖2。

表1 研究區2005年~2007年石漠化動態變化表(單位:km2)Tab.1 2005~2007 dynamic change of rocky desertification in research area

圖2 畢節、遵義和沿河小流域圖斑破碎度圖Fig.2 Fragmentation of rocky desertification in Yachi, Longba and Qitan research area
石漠化發生率反映了區域石漠化發生的概率,石漠化發生率受各種石漠化影響因子的作用,研究石漠化發生率可以分析不同河段發生石漠化的危險性。為此,統計分析三個小流域石漠化面積和巖溶面積的比例。經統計計算,沿河淇灘小流域石漠化發生率最高,達 40.3%;其次是畢節鴨池小流域,石漠化發生率為 37.1%;遵義巖溶高原盆地石漠化發生率最低,只有 22.4%。流域的上游和下游石漠化發生率明顯高于中游地區,特別要注意流域上游和下游的保護問題,防治石漠化的發生。
流域不同河段各種等級石漠化的面積比例存在一定差異。分別對三個小流域各種石漠化所占的巖溶面積比例進行統計分析(見圖 3),研究不同河段各種等級石漠化的構成特征。

圖3 沿河、遵義和畢節鴨池小流域各種石漠化比例Fig.3 rocky desertification scale in Yachi,Longba and Qitan research area
對流域不同位置各等級石漠化構成比例進行分析,發現流域不同河段無石漠化比例相差較大,流域中游無石漠化比例最高,其次是下游地區,流域上游無石漠化比例最低;流域不同地段潛在石漠化比例相差較小,其中流域上游潛在石漠化比例最高,中游和下游地區潛在石漠化比例相當;上游區畢節鴨池小流域輕度石漠化比例最高,其次是沿河烏江流域下游地區,遵義烏江流域中游地區輕度石漠化比例最小,流域不同地段輕度石漠化比例相差較大;烏江流域下游地區中度石漠化比例最高,其次是中游地區,烏江流域上游地區中度石漠化比例最小;烏江流域下游強度石漠化比例最高,而流域上游區強度石漠化比例最小。中度和強度石漠化都是烏江流域下游淇灘小流域比例最高,其石漠化等級程度在流域內最為嚴重,中游和上游小流域中、強度石漠化比例逐漸降低,石漠化等級變低。
石漠化的發生率和各種石漠化的面積比例只能反映出不同地貌類型石漠化的構成。為了研究不同地貌類型區石漠化的變化情況,對各小流域的石漠化變化面積和石漠化年變化率進行統計分析,研究不同地貌類型區石漠化的變化特征。分別對代表不同河段的三個小流域2005年~2007年間的石漠化面積變化進行統計,研究各小流域各種等級石漠化變化量。
對流域不同河段石漠化變化面積(表1)進行分析,發現各小流域無石漠化面積都在增加,其中烏江流域中游面積增加最大;和無石漠化一樣,潛在石漠化面積在流域不同河段同樣都是增加的,但是流域下游段淇灘小流域增加的面積很小,在流域上游地區增加的面積最大;輕度石漠化面積在流域不同河段都在減少,流域中游變化量最大,其次是流域的上游和下游;中度石漠化面積只有流域的下游在增加,烏江流域上游和中游地區中度石漠化面積都在減少,其中流域中游減少面積最大;流域內強度石漠化面積都有所降低,從下游到上游地區變化面積逐漸減低。
石漠化程度綜合指數主要反映了小流域石漠化發生的廣度和深度,反映石漠化屬性,同時反映了人類因素和自然環境的綜合效應。小流域的石漠化程度變化是各種石漠化類型變化的綜合結果,石漠化程度變化量可定量表達該地區石漠化的綜合水平和變化趨勢,石漠化程度綜合變化率反映各小流域石漠化變化快慢。石漠化程度綜合變化率表達為:

式中R為石漠化程度綜合變化率,A為第i級石漠化分級數,Cib和Cia分別為b時間和a時間第i等級的石漠化程度所占區域面積的百分比。當R>0 時,表示該區域石漠化處于下降階段;R<0則相反。
從流域不同河段的石漠化綜合指數統計數據可以看出,在2005年~2007年間,沿河淇灘小流域石漠化程度變化量最大,達 6.32,而烏江流域上游的畢節鴨池小流域變化量較小,只有 1.98,兩個小流域之間石漠化程度變化量相差較大。從石漠化程度的年變化率分析,年變化率最高的是沿河小流域3.0%,其次是遵義2.4%,畢節小流域的年變化率最低,只有1.4%,石漠化綜合指數年變化率從流域的上游到下游逐漸升高。
根據各種石漠化的年變化率可以更好地分析不同河段各種石漠化的變化快慢。從圖 4可以看出,沿河烏江流域下游區無石漠化年變化率最大,其次是遵義烏江流域中游區和畢節上游區;潛在石漠化年變化率上游和中游區變化較大,而流域下游區變化率最小;流域中游區輕度石漠化的年變化率明顯高于流域上游和下游區,上游和下游地區年變化率相當;中度石漠化只有流域下游年變化率為正值,遵義和畢節鴨池小流域都為負值,其中,遵義烏江流域中游區變化率最大;強度石漠化各種地貌類型區年變化率都為負值,其中,流域下游區年變化率最大,畢節流域上游區和遵義流域中游區變化率也較大。從流域不同河段分析,潛在、輕度和中度石漠化都是遵義烏江流域中游區年變化率最大,說明其石漠化治理相對比較容易,石漠化治理效果要比上游和下游地區明顯。而沿河烏江流域下游段無石漠化和強度石漠化年變化率最高,其他幾種石漠化類型年變化率最低,說明其石漠化治理趨于兩級分化,中、輕度石漠化治理效果不明顯。

圖4 畢節鴨池、遵義龍壩和沿河淇灘小流域石漠化年變化率Fig.4 annual change rate of rocky desertification in Yachi, Longba and Qitan research area
因遙感影像限制,不能全面分析烏江流域不同河段石漠化特征,選擇代表烏江流域不同河段的三個小流域進行石漠化特征及演變規律對比分析。受影像時段限制,動態變化選取的時間不長,最好選用時間間隔更長的影像進行對比分析,更能反映出烏江流域石漠化變化規律。
通過對烏江不同河段三個小流域石漠化對比分析發現,烏江流域不同河段石漠化現狀及演變具有一定規律性,烏江流域上游石漠化平均斑塊面積較大,石漠化圖斑的破碎度最小,中游及下游段石漠化圖斑破碎度最大。烏江流域下游和上游河段石漠化發生率明顯高于中游地區,下游段石漠化發生率最高;從石漠化等級強度差異分析,中度和強度石漠化都是烏江流域下游比例最高,其石漠化等級程度在流域內最為嚴重,中游和上游小流域中強度石漠化比例逐漸降低,石漠化等級降低。
2005年~2007年間,烏江流域各河段無石漠化和潛在石漠化面積都增加,流域下游沿河淇灘小流域中度石漠化面積增加,其余各等級石漠化面積都減少,烏江流域石漠化處于好轉趨勢。其中,流域中段龍壩小流域石漠化減少面積最大,其次是上游鴨池小流域,淇灘小流域石漠化面積減少最小。石漠化綜合年變化率最高的是沿河小流域,畢節小流域的年變化率最低,石漠化綜合年變化率從流域的上游到下游逐漸升高。烏江流域中段龍壩小流域石漠化治理效果要比上游和下游地區明顯。下游段沿河淇灘小流域石漠化治理趨于兩級分化,中、輕度石漠化治理效果不明顯。
[1] 孫凡,徐勝旺,馬生麗,等.典型喀斯特地區季節性石漠化與生態環境建設[J].西南農業大學學報:社會科學版,2011,9(5):1-6.
[2] PRENTICEIC.The Carbon Cycle and Atmospheric Carbon Dioxide[M].Climate Change 2001: The Scientific Basis(IPCC).Cambridge: Cambridge University Press, 2001:184-237.
[3] Yuan Dao-xian.Rock Desertification in the Subtropical Karst of South China.In Webmaster.Copyrights by Dynamics Laboratory and Network Center of Guangxi Normal University, 1995-1999.
[4] 王魁,馬祖陸,蔡德所,等.我國巖溶地區水土流失研究進展及發展趨勢[J].中國水土保持,2011,(9):30-32.
[5] 張信寶,王世杰,白曉永,等.貴州石漠化空間分布與喀斯特地貌、巖性、降水和人口密度的關系[J].地球與環境,2013,41(1):1-6.
[6] 黃秋昊,蔡運龍,王秀春.我國西南部喀斯特地區石漠化研究進展[J].自然災害學報,2007,16(2):106-111.
[7] 熊康寧,黎平,周忠發,等.喀斯特石漠化的遙感—GIS典型研究[M].北京:地質出版社,2002.
[8] 蘇維詞,楊華,李晴,等.我國西南喀斯特山區土地石漠化成因及防治[J].土壤通報,2006,37(3):446-450.