西南典型喀斯特地區(qū)石漠化時(shí)空演變特征
——以貴州省普定縣為例

習(xí)慧鵬,王世杰,白曉永,*,唐 紅,吳路華,陳 飛,2,肖建勇,2,王明明,李匯文,操 玥,陳 歡,2,冉 晨,2,羅旭玲,2

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

中國(guó)南方喀斯特地區(qū)存在著非常嚴(yán)重的石漠化問題[1- 3],在此地區(qū)的自然災(zāi)害和貧窮發(fā)生率極高,嚴(yán)重制約著該地區(qū)的發(fā)展。目前,研究其空間演變規(guī)律已成為治理石漠化的重要前提,以往學(xué)者利用不同方法分析了石漠化的空間演變規(guī)律：熊康寧等基于2000年Landsat TM影像和實(shí)際調(diào)查數(shù)據(jù),將貴州省石漠化分為6級(jí)并探討了貴州省石漠化的現(xiàn)狀及其空間演變趨勢(shì)[4]；胡寶清等基于2000年TM影像提取了廣西省石漠化信息,并利用數(shù)理統(tǒng)計(jì)等方法探討了石漠化的分布演變規(guī)律,得出了人類活動(dòng)是石漠化產(chǎn)生的主要原因[5]；白曉永等基于1986、1995、2000年3期貴州省石漠化數(shù)據(jù),利用轉(zhuǎn)移矩陣等方法探討了貴州省石漠化的時(shí)空演變規(guī)律[6]。可以看出,以往的研究多基于實(shí)際數(shù)據(jù)和數(shù)理統(tǒng)計(jì)的方法綜合分析,所取得的成果對(duì)早期的石漠化治理提供了巨大的科技支撐,然而研究方向多集中于省域級(jí)大尺度低空間分辨率上,且研究年限多在20世紀(jì)80年代—21世紀(jì)初,所取得的研究成果和經(jīng)驗(yàn)方法對(duì)已治理了20余年的石漠化指導(dǎo)意義不足。

此外,如何快速高效地獲取石漠化各因子的貢獻(xiàn)率一直是眾多學(xué)者的研究重點(diǎn),在以往研究中學(xué)者們通過數(shù)學(xué)分析的手段總結(jié)出了很多經(jīng)典的方法。藍(lán)安軍等基于主成分分析、相關(guān)性分析,并建立了有關(guān)石漠化的動(dòng)力指數(shù)公式,對(duì)貴州省的石漠化眾多影響因子分析得出了包括植被覆蓋度、喀斯特面積、河網(wǎng)密度、耕地、未利用地、草地、海拔、土地開墾率和25°以上的土地面積等9個(gè)因子是導(dǎo)致石漠化產(chǎn)生的驅(qū)動(dòng)因素[7]；胡寶清等基于主成分分析,主要分析了巖性、地貌、水熱條件、土壤和NDVI等5個(gè)指標(biāo),計(jì)算得到石漠化地質(zhì)背景區(qū)的環(huán)境生態(tài)指數(shù)回歸模型[8]；皇甫江云等基于專家打分的方法,對(duì)各等級(jí)石漠化對(duì)應(yīng)的影響因子歸類打分,得出草地石漠化中各因子的貢獻(xiàn)率[9]?？梢钥闯?以往的經(jīng)典方法受人為影響較大,且由于喀斯特環(huán)境的復(fù)雜性,石漠化影響因子回歸模型無法推廣。

因此,針對(duì)以往研究中存在的兩大問題,本文選取貴州省中部典型喀斯特地區(qū)普定縣為研究區(qū),基于Landsat、DEM 等30 m分辨率數(shù)據(jù),提取了植被覆蓋度、巖石裸露率和坡度等3個(gè)最具代表性的指標(biāo)來精確反演石漠化[10- 12],旨在解決以下3個(gè)問題：(1)尋找更加高效快速的石漠化模擬反演方法。(2)引入地理探測(cè)器,采用更科學(xué)的方法確定各因子的貢獻(xiàn)率。(3)研究各等級(jí)石漠化分布規(guī)律,分析新形勢(shì)下貴州典型喀斯特地區(qū)的石漠化演變速率,揭示石漠化的時(shí)空演變特征(技術(shù)路線見圖4)。研究結(jié)果中獲得的各因子貢獻(xiàn)率為反演石漠化提供了參考,發(fā)現(xiàn)的石漠化分布演變規(guī)律為石漠化治理提供了方向。

1 研究區(qū)與數(shù)據(jù)源

1.1 研究區(qū)概況

位于貴州省中西部地區(qū)的普定縣(圖1),全縣國(guó)土面積約1081.8 km2,喀斯特面積為869.64 km2,年均降水量約1200 mm,年均溫15.3℃,平均海拔為1100 m,地勢(shì)南北高,中間低。截止2013年,普定縣人口密度達(dá)到440人/km2,遠(yuǎn)遠(yuǎn)高于我國(guó)西部其他地區(qū)。由于巖溶地貌的廣泛發(fā)育,石灰土比重占到了63%,連年的高強(qiáng)度降水和人類活動(dòng),使得普定縣石漠化極為嚴(yán)重,2008年,普定被確定為全國(guó)100個(gè)、全省55個(gè)石漠化綜合治理試點(diǎn)縣,眾多科研單位在此地區(qū)開展石漠化治理工作,普定縣成為了石漠化防治研究的典型地區(qū)。

圖1 研究區(qū)位置圖,左圖為研究區(qū)在貴州省的位置,右圖為生成的普定縣TIN(triangular irregular network)Fig.1 Location map of the study area

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

本研究使用的實(shí)驗(yàn)數(shù)據(jù)主要包括：普定縣巖性圖、普定縣石漠化分布圖、Landsat TM/ETM+、DEM等30 m空間分辨率數(shù)據(jù),具體見表1。

表1 主要數(shù)據(jù)源

2 研究方法

2.1 石漠化指標(biāo)

2.1.1 巖石裸露率

巖石裸露作為石漠化地區(qū)最顯著的特征,在石漠化信息的提取中有著至關(guān)重要的作用,本文擬采用NDRI來獲取研究區(qū)巖石裸露率,NDRI(Normalized Difference Rock Index)是指歸一化巖石指數(shù),基于喀斯特碳酸鹽巖的光譜信息在短波紅外波段具有較大差異[13-14],基于NDVI思想,在ENVI 5.2軟件下進(jìn)行波段運(yùn)算,根據(jù)公式(1)得到研究區(qū)的NDRI：

(1)

根據(jù)像元二分模型原理,假設(shè)像元由巖石裸露和非裸露兩部分組成,根據(jù)NDRI,計(jì)算得到基巖裸露率(Fr)：

(2)

式中, NDRIr為全由巖石裸露得到的像元值,NDRIo為全無巖石裸露得到的像元值。在統(tǒng)計(jì)得到的所有像元值中,取頻率1%的像元值為NDRIo,取累計(jì)頻率為99%的像元值為NDRIr。

2.1.2 FVC(Fractional Vegetation Cover)

FVC植被覆蓋度指的是植被冠層在地面的垂直投影面積與所處地總面積之比[15-16],同樣根據(jù)像元二分模型原理,本文擬采用NDVI來獲取研究區(qū)植被覆蓋度,公式如下：

(3)

在計(jì)算得到像元的NDVI值后,統(tǒng)計(jì)其百分比,取1%的像元值為NDVIo,作為沒有植被存在的像元值,取99%的像元值為NDVIr,作為完全由植被覆蓋的像元值。

2.1.3 Slope

本文選擇的DEM數(shù)據(jù)為GDEMDEM 30 m分辨率數(shù)字高程數(shù)據(jù),數(shù)據(jù)標(biāo)識(shí)為ASTGTM_N26E105,在ArcGIS下對(duì)其進(jìn)行投影轉(zhuǎn)換為ALBERS等面積投影,坐標(biāo)系為Krasovsky_1940_Albers,掩膜提取后制作得到坡度圖。

2.2 貢獻(xiàn)率獲取

地理探測(cè)器是探測(cè)和利用空間分異性的工具[17],在空間一致的情況下,通過采集樣本點(diǎn)上各自變量和因變量的值,可以利用公式計(jì)算自變量對(duì)因變量的解釋能力大小,公式表達(dá)式為：

(4)

圖2 2000年石漠化及其布點(diǎn)圖Fig.2 Karst rocky desertification and its site map in 2000

圖3 各因子十級(jí)分布圖Fig.3 Ten level distribution map of each factor

2.3 石漠化分類定級(jí)標(biāo)準(zhǔn)

基于前人工作建立石漠化分級(jí)分類標(biāo)準(zhǔn)[10, 18- 21](表2),通過地理探測(cè)器得到各因子貢獻(xiàn)率之后,在ArcGIS 10.2軟件平臺(tái)下按照下表進(jìn)行重分類,并疊置分析,基于普定縣巖性圖和土地分類結(jié)果剔除非喀斯特地區(qū)和水體、城鄉(xiāng)建筑等不可能發(fā)生石漠化的地方,得到石漠化空間分布圖。

表2 石漠化分類定級(jí)標(biāo)準(zhǔn)

2.4 石漠化空間動(dòng)態(tài)變化分析

石漠化動(dòng)態(tài)變化無論是在時(shí)間上還是在空間上都是一個(gè)復(fù)雜的過程,傳統(tǒng)的數(shù)字分析表達(dá)能力欠佳,因此本文在探討石漠化時(shí)空演變時(shí)將數(shù)學(xué)模型與制圖綜合在一起,旨在探討以下科學(xué)問題：第一是明確石漠化因子的變化特征,以探討因變量對(duì)自變量的影響大?。坏诙敲鞔_石漠化的轉(zhuǎn)移方向；第三是明確各等級(jí)石漠化的轉(zhuǎn)移快慢。因此,有必要對(duì)石漠化的因子變化和時(shí)空轉(zhuǎn)變分別構(gòu)建模型討論。

2.4.1 石漠化因子變化

本文中石漠化的影響因子主要由植被覆蓋度、巖石裸露率、坡度等組成,而坡度由于變化較小無法捕捉,故不在此討論。在石漠化演變過程中,像元中的植被覆蓋度和巖石裸露率等級(jí)必然也在發(fā)生著變化,因此,兩者的變化必然與石漠化演變具有密切關(guān)系,基于此,構(gòu)建下述公式：

Q=q1-q2

(5)

式中,Q指兩個(gè)不同年份的因子在柵格計(jì)算器下相減后的值,q1指一段時(shí)間內(nèi)該像元變化前的等級(jí),q2指變化后的等級(jí)。

2.4.2 石漠化轉(zhuǎn)移矩陣

在某一時(shí)間段內(nèi),石漠化的演變方向與演變規(guī)模大多較復(fù)雜,且時(shí)常因地而異,因此需要定量化地討論以明確其大致的演變方向和趨勢(shì),構(gòu)建轉(zhuǎn)移變化矩陣。

(6)

式中,S為面積,單位為km2；m,n為一段時(shí)期內(nèi)前期和末期石漠化的等級(jí)；Sk,k為在該時(shí)期內(nèi)石漠化的轉(zhuǎn)移面積,單位為km2。

2.4.3 石漠化轉(zhuǎn)移變化速率

僅僅探討石漠化的轉(zhuǎn)移方向和面積對(duì)于治理石漠化的指導(dǎo)意義是遠(yuǎn)遠(yuǎn)不夠的,石漠化的演變快慢反映了一種類型石漠化在該地區(qū)的變化發(fā)展模式,在橫向比較上更具有科學(xué)價(jià)值,也更能說明石漠化的演變情況。因此需要建立模型以探討石漠化的演變速率。

(7)

式中,Pi為某一類型的石漠化演變速率,單位為km2/a；S2為一定時(shí)間段內(nèi)末期的面積,S1為前期的面積,單位為km；T為該時(shí)間段,單位為a。

2.5 技術(shù)路線

圖4 技術(shù)路線圖Fig.4 Technical roadmap

3 結(jié)果與分析

3.1 權(quán)重分析

在地理探測(cè)器中根據(jù)石漠化和各因子等級(jí)計(jì)算得到各自所占的權(quán)重,結(jié)果表明,三個(gè)指標(biāo)中植被覆蓋度和巖石裸露率主要控制了石漠化的分布：FVC(0.47)>NDRI(0.45)>Slope(0.15),由于三者相加大于1,故各因子權(quán)重所占的百分比即為本文中得出的石漠化貢獻(xiàn)率：FVC(44%)、NDRI(42%)、Slope(14%)(表3)。

表3 石漠化各因子貢獻(xiàn)率

Q：百分比權(quán)重,quantile；P：置信度,power

3.2 石漠化分布及各因子等級(jí)圖

按照表2已劃分的石漠化等級(jí)標(biāo)準(zhǔn),在ENVI 5.2、ArcGIS 10.2軟件下制圖分類,得到普定縣1990—2015年FVC、NDRI、Slope等級(jí)圖,并在地理探測(cè)器中獲得其各自的貢獻(xiàn)率,在ArcGIS下疊置分析,并結(jié)合Landsat影像,在ENVI下得到土地利用分類圖,剔除了水體、城鄉(xiāng)建筑地等不可能發(fā)生石漠化的地區(qū),得到各年份石漠化空間分布圖。

3.2.1 NDVI提取植被覆蓋度

首先在ENVI 5.2下波段運(yùn)算得到的1990—2015不同年份的NDVI,然后再根據(jù)公式(3)計(jì)算得到植被覆蓋度,并在ArcGIS10.2軟件下利用Reclassify工具進(jìn)行重分類,按照表2將其分為>71%、51%—71%、31%—50%、21%—30%、11%—20%、<10%共六類并分別賦值為1、2、3、4、5、6,得到1990—2015年植被覆蓋度分布圖(圖5)。由圖中可以看出,越靠近河流和城鎮(zhèn)地區(qū),植被覆蓋度的等級(jí)越低。

圖5 植被覆蓋度空間分布Fig.5 The spatial distribution of fractional vegetation cover

3.2.2 NDRI提取巖石裸露率

首先根據(jù)公式(1),在ENVI 5.2下波段運(yùn)算得到的1990—2015不同年份的NDRI, 再利用公式(2)計(jì)算得到巖石裸露率,并在ArcGIS下重分類,將其按照表2分為<20%、21%—30%、31%—50%、51%—70%、71%—90%、>90%共六類并分別賦值為1、2、3、4、5、6,得到1990—2015年巖石裸露率分布圖(圖6)。由圖中可以看出,越靠近河流和城鎮(zhèn)地區(qū),巖石裸露率的等級(jí)越高,且1997年的巖石裸露面積明顯最小。

圖6 巖石裸露率空間分布Fig.6 The spatial distribution of rock exposedness

3.2.3 DEM提取坡度

在投影變換后,利用ArcGIS軟件中的Slope工具根據(jù)DEM數(shù)據(jù)提取坡度,并按照表2的標(biāo)準(zhǔn),重分類得到普定縣坡度空間分布圖(圖7),由圖中可以看出,坡度較陡的區(qū)域多沿著普定縣SW—NE方向分布,總體呈現(xiàn)出南北坡度大,中部較平坦。

圖7 坡度分布圖Fig.7 The spatial distribution of Slope

3.2.4 石漠化分布圖

在得到各年段的植被覆蓋度、巖石裸露率和坡度分布后,利用ArcGIS 10.2中的疊加分析工具,按照地理探測(cè)器計(jì)算得出的各因子貢獻(xiàn)率疊加制圖,并剔除非喀斯特地區(qū)和水體、城鄉(xiāng)建筑等不可能發(fā)生石漠化的區(qū)域, 得到1990、1997、2006和2015年的石漠化空間分布圖(圖8)。

圖8 石漠化空間分布圖Fig.8 The spatial distribution of karst rocky desertification

3.2.5 精度驗(yàn)證

為檢驗(yàn)反演效果,本文選取了喀斯特科學(xué)數(shù)據(jù)中心提供的2005年普定縣石漠化實(shí)際分布圖作為標(biāo)準(zhǔn),將制作得到的2006年數(shù)據(jù)進(jìn)行了對(duì)比,由圖9和表4可以看出,石漠化在總體上反演精度較高,達(dá)到了80%,該石漠化反演方法可信度較高。

3.3 時(shí)空演變分析

3.3.1 石漠化面積統(tǒng)計(jì)分析

根據(jù)石漠化空間分布圖,利用柵格計(jì)算器計(jì)算得到普定縣不同年份的各等級(jí)石漠化面積,統(tǒng)計(jì)如下表5,可以看出,無石漠化面積占比最大最穩(wěn)定,而極重度石漠化由于面積小,變化面積不明顯,而潛在石漠化、輕度石漠化和中度以上石漠化的面積在不同年段發(fā)生著劇烈地變化。1990—1997年間,石漠化面積減少了70.95 km2；1997—2006年間,石漠化面積增加了216.45 km2；2006—2015年間,石漠化面積又減少了111.05 km2；1990—2015年間,石漠化面積增加了34.50 km2。 綜上,普定縣輕度以下石漠化分布面積在喀斯特地區(qū)均有大量分布,而中度以上多分布在三岔河流域沿岸附近；在1990—2015年間,普定縣石漠化演變復(fù)雜劇烈,經(jīng)歷了由好轉(zhuǎn)到惡化再到好轉(zhuǎn)的階段,演變總體趨勢(shì)在好轉(zhuǎn),但局部在惡化,說明治理石漠化不可掉以輕心,還需要加大治理力度。

圖9 實(shí)際與模擬石漠化對(duì)比Fig.9 Comparison between actual and simulated karst rocky desertification

年份Year無石漠化NKRD潛在石漠化PKRD輕度石漠化LKRD中度石漠化MKRD重度石漠化SKRD已石漠化KRD2005年/km2379.10164.79225.5792.527.60325.692006年/km2306.99171.23271.96100.6518.80391.36精度 Precision/%80.9896.2482.9491.9240.5383

表中石漠化類型均為英文簡(jiǎn)寫,KRD:石漠化,Karst rocky desertification；NKRD:無石漠化,No karst rocky desertification；PKRD:潛在石漠化,Potential karst rocky desertification; LKRD:輕度石漠化,Light karst rocky desertification；MKRD:中度石漠化,Moderate karst rocky desertification；SKRD:重度石漠化,Severe karst rocky desertification

表5 各級(jí)石漠化在喀斯特區(qū)域的面積和比例

ESKRD：極重度石漠化,Extremely severe karst rocky desertification

3.3.2 石漠化因子時(shí)空演變分析

根據(jù)圖5以及公式(5),在ArcGIS 10.2軟件下運(yùn)用柵格計(jì)算器進(jìn)行減法運(yùn)算,將各年份植被覆蓋度相減,得到1990—1997、1997—2006、2006—2015、1990—2015年間的植被覆蓋度變化圖(圖10),其中,正值表示該地區(qū)的植被覆蓋度等級(jí)在升高,負(fù)值則表示在下降。可以看出,1990—1997年植被覆蓋度大都在減少,而在2006—2015年卻大都在增加,在1997—2006年間,植被覆蓋度的變化則具有明顯的地理位置特征。

圖10 植被覆蓋度時(shí)空演變Fig.10 Fractional vegetation cover changes in each period

再次運(yùn)用柵格計(jì)算器工具,將圖6中的各年份巖石裸露率相減,得到1990—1997、1997—2006、2006—2015、1990—2015年間的變化圖(圖11),同樣,正值表示該地區(qū)的巖石裸露率等級(jí)在升高,負(fù)值則表示在下降?？梢钥闯?1997—2006年巖石裸露率等級(jí)大都在升高,而在其他年段則大都在降低,結(jié)合普定縣行政區(qū)分析,升高的地區(qū)多位于普定縣中東部城關(guān)鎮(zhèn)一帶。

圖11 巖石裸露率時(shí)空演變Fig.11 Rocky exposedness changes in each period

3.3.3 石漠化轉(zhuǎn)移變化時(shí)空演變分析

根據(jù)公式(6),在柵格計(jì)算器下通過柵格疊加得到某一年份對(duì)應(yīng)另一年份的石漠化面積,統(tǒng)計(jì)結(jié)果見圖12(NKRD—無石漠化；PKRD—潛在石漠化；LKRD—輕度石漠化；MKRD—中度石漠化；SKRD—重度石漠化；ESKRD—極重度石漠化；Non—karst：非喀斯特區(qū)域)和圖13,可以看出：1990—1997年,石漠化好轉(zhuǎn)的比例(25.06%)>惡化(16.76%)；1997—2006年,石漠化好轉(zhuǎn)(9.73%)<惡化(37.83%)；2006—2015年,石漠化好轉(zhuǎn)(28.26%)>惡化(14.20%)；1990—2015年,石漠化好轉(zhuǎn)(24.65%)<惡化(25.47%)；此外,在圖中可以看出由無、潛在石漠化轉(zhuǎn)變?yōu)橹囟纫陨鲜拿娣e雖然不大,但占比極高,結(jié)合影像分析(圖14),發(fā)現(xiàn)此種演變方式集中于三岔河流域沿岸附近,由于河道拓寬,導(dǎo)致沿岸一帶植被減少、基巖裸露率增大,且此地區(qū)坡度均較大,是石漠化發(fā)生急劇變化的主要原因。而重度以上石漠化轉(zhuǎn)變?yōu)闊o、潛在石漠化面積較少,結(jié)合影像分析(圖14),發(fā)現(xiàn)此種演變方式多發(fā)生于城關(guān)鎮(zhèn)一帶,此處靠近普定縣城,且海拔較低,有利于土壤積累,極少部分較為嚴(yán)重的石漠化地區(qū)突變?yōu)榱藷o石漠化或潛在石漠化,筆者認(rèn)為是居民的開墾荒地導(dǎo)致附近的石漠化土地演變?yōu)榱烁鼗蛄值?因此可以看到人類活動(dòng)的兩面性。

綜上,在石漠化演變過程中,輕度以下石漠化演變最為劇烈,尤其是潛在石漠化和輕度石漠化,演變面積大,方向難測(cè)；石漠化的演變方式因地而異,在靠近河流沿岸地區(qū),石漠化常發(fā)生突變而惡化,而在城鄉(xiāng)居民地附近,常發(fā)生重度石漠化突變?yōu)闊o石漠化。

圖12 石漠化轉(zhuǎn)移變化分布圖Fig.12 Change distribution map of karst rocky desertification圖中石漠化類型均為英文簡(jiǎn)寫,KRD：石漠化,Karst rocky desertification；NKRD：無石漠化,No karst rocky desertification；PKRD：潛在石漠化,Potential karst rocky desertification；LKRD：輕度石漠化,Light karst rocky desertification；MKRD：中度石漠化,Moderate karst rocky desertification；SKRD：重度石漠化,Severe karst rocky desertification；ESKRD：極重度石漠化,Extremely severe karst rocky desertification

圖13 1990—2015年石漠化轉(zhuǎn)移矩陣Fig.13 Karst rocky desertification transfer matrix in 1990—2015

圖14 石漠化解譯與實(shí)際影像和衛(wèi)星圖對(duì)比Fig.14 Comparison of karst rocky desertification with actual image and satellite map

3.3.4 石漠化轉(zhuǎn)移變化速率

再次運(yùn)用柵格計(jì)算器計(jì)算得到普定縣石漠化轉(zhuǎn)移速率變化圖(圖15),并統(tǒng)計(jì)其轉(zhuǎn)移矩陣如圖16,可以看出重度以上石漠化演變速率極低,幾乎等于0,這也驗(yàn)證了前人的工作基礎(chǔ)[6, 22]。在1990—1997年間,輕度-潛在石漠化變化速率最大,為11.2 km2/a；1997—2006年間,潛在-輕度石漠化變化速率最大,為16.5 km2/a,2006—2015年間,輕度-潛在石漠化變化速率最大,為10.4 km2/a,1990—2015年間,潛在-輕度石漠化變化速率最大,為3.22 km2/a。計(jì)算各等級(jí)石漠化平均變化速率得出(圖17),潛在石漠化惡化的速率最大,為2.75 km2/a,其中,在1997—2006年間中速率達(dá)到了16.5 km2/a,演變十分迅速。此外,在各年段上,中度以上石漠化整體均在好轉(zhuǎn),而輕度以下石漠化則反復(fù)無常,尤其是潛在石漠化,好轉(zhuǎn)的速率極小,而惡化則十分迅速,因此需要加強(qiáng)對(duì)輕度以下石漠化的監(jiān)察與防護(hù)。

圖15 石漠化轉(zhuǎn)移速率分布圖Fig.15 Velocity distribution map of karst rocky desertificationX—X為石漠化等級(jí)無變化區(qū)域,即石漠化演變速率為0；NK為非喀斯特地區(qū),非喀斯特地區(qū)不涉石漠化問題

圖16 1990—2015年石漠化轉(zhuǎn)移速率矩陣Fig.16 Transfer rate matrix of karst rocky desertification in 1990—2015

圖17 各級(jí)石漠化平均變化速率 Fig.17 Average change rate of karst rocky desertification at all grades

4 討論

4.1 石漠化時(shí)空演變分析

1)1990年至1997年間,石漠化面積減少了71.0 km2,實(shí)際上是在好轉(zhuǎn),而在1997年至2006年間,石漠化面積又快速增加了217.0 km2,尤其是潛在向輕度石漠化的轉(zhuǎn)變面積最大,速度最快,這段時(shí)間正處于經(jīng)濟(jì)發(fā)展的黃金期,人類不合理活動(dòng)強(qiáng)度的增大是導(dǎo)致惡化的主要原因,直到2005年后,政府開始實(shí)施退耕還林政策,大量石漠化土地被禁止開發(fā),石漠化面積才又開始逐漸減少[23-24],從2006年至2015年,石漠化面積又減少了110.0 km2。整體來看,從1990年到2015年,石漠化面積增加了34.5 km2,普定縣石漠化經(jīng)歷了由好轉(zhuǎn)到惡化再到好轉(zhuǎn)的演變,總體趨勢(shì)雖然在好轉(zhuǎn),但局部地區(qū)仍在惡化,石漠化的治理不可掉以輕心。

2)石漠化轉(zhuǎn)移變化平均速率最快的是潛在石漠化,平均速率為2.75 km2/a,其中在1997—2006年間,速率達(dá)到了16.5 km2/a；變化速率最慢的是極重度石漠化,平均速率為0.0002 km2/a,這也驗(yàn)證了前人的工作成果,石漠化變化速率與等級(jí)高低成反比。此外,好轉(zhuǎn)速率最快的是中度石漠化,平均速率為0.80 km2/a,惡化速率最快的是無石漠化,平均速率為3.10 km2/a,由圖17看出,輕度以下石漠化變化方向不易確定,而中度以上石漠化則幾乎都在好轉(zhuǎn),因此,在治理中度以上石漠化的同時(shí),更要注重輕度以下尤其是潛在石漠化的防治工作。

3)結(jié)合石漠化因子轉(zhuǎn)移變化圖與石漠化空間變化圖可以看出,植被覆蓋度與巖石裸露率的變化成相反關(guān)系,在植被覆蓋度增加高的區(qū)域,石漠化演變多以好轉(zhuǎn)為主,多集中在普定縣中部、東南一帶；而在巖石裸露率增加高的區(qū)域,石漠化演變多以惡化為主,多集中在三岔河一帶。

4)由圖8分析得出,在石漠化空間分布上,本研究發(fā)現(xiàn)普定縣中度以上石漠化多集中分布在三岔河沿岸附近,而輕度以下石漠化則分布較散亂,無太大規(guī)律。

5)由圖13、圖14看出,在三岔河附近,常發(fā)生著無或潛在石漠化突變?yōu)橹囟仁默F(xiàn)象,這往往是由河道拓寬,雨水沖刷,導(dǎo)致無或潛在石漠化地區(qū)植被劇減,巖石裸露率暴增導(dǎo)致。此外,在城關(guān)鎮(zhèn)靠近縣城一帶,也會(huì)出現(xiàn)重度石漠化突變?yōu)闊o石漠化的現(xiàn)象,筆者認(rèn)為由于此處地勢(shì)較低,易于土層堆積,人類的開墾荒地等城鎮(zhèn)化活動(dòng)使得重度或極重度石漠化也變?yōu)榱藷o或潛在石漠化,從而這也體現(xiàn)了人類活動(dòng)對(duì)石漠化演變的積極一面。

4.2 展望與不足

本文使用的數(shù)據(jù)由于依賴于遙感影像,因此難免受到衛(wèi)星過境時(shí)間和天氣云量等的影響,加之獲得的植被覆蓋度、巖石裸露率等常受到季節(jié)變化的干擾,很難去計(jì)算全年的平均值,此為最大的局限性。本文采取的方法在總體上反演石漠化精度較高,因子權(quán)重較為合理,在喀斯特中小尺度地區(qū)內(nèi)反演石漠化具有一定的優(yōu)勢(shì)。

5 結(jié)論

本文基于Landsat、DEM等數(shù)據(jù)提取得到的普定縣1990—2015年4期石漠化影像展開分析,重點(diǎn)討論了普定縣的石漠化空間演變規(guī)律,主要認(rèn)識(shí)如下：1)在1990—2015年間,普定縣石漠化演變復(fù)雜劇烈,經(jīng)歷了由好轉(zhuǎn)到惡化再到好轉(zhuǎn)的階段,總體演變趨勢(shì)雖然好轉(zhuǎn),但局部地區(qū)仍舊在惡化,石漠化治理力度需要加大。2)普定縣中度以上石漠化集中分布在三岔河沿岸附近,輕度以下石漠化則無明顯地理分布特征。3)潛在和輕度石漠化演變方向復(fù)雜,潛在石漠化變化速率最大(2.75 km2/a),最大時(shí)在1997—2006年間達(dá)到了16.5 km2/a。4)植被覆蓋度與巖石裸露率的變化呈負(fù)相關(guān),前者主導(dǎo)了石漠化的好轉(zhuǎn),后者控制了石漠化的惡化。5)石漠化演變存在突變,在三岔河沿岸附近,存在無/潛在石漠化有突變?yōu)橹囟仁默F(xiàn)象；在靠近城鄉(xiāng)居民地附近,也存在重度以上石漠化突變?yōu)闊o石漠化的現(xiàn)象。6)植被覆蓋度、巖石裸露率、坡度等能較準(zhǔn)確反演石漠化,貢獻(xiàn)率分別為44%、42%、14%。