1980—2015年典型巖溶區(qū)土地利用動(dòng)態(tài)變化研究

劉南健 安佑志 李松

摘 要：基于1980—2015年7期的Landsat TM遙感影像解譯數(shù)據(jù)，結(jié)合地學(xué)統(tǒng)計(jì)方法對(duì)典型巖溶區(qū)縣域尺度下的土地利用時(shí)空演變格局進(jìn)行定量分析。結(jié)果表明，35年來(lái)，研究區(qū)的用地類(lèi)型以林地、耕地、草地為絕對(duì)主導(dǎo)用地，3者面積比例之和穩(wěn)定在98%以上，整體變幅較小。研究區(qū)用地變化劇烈程度依次是建設(shè)用地>水域>林地>未利用地>草地>耕地;從綜合土地利用角度看，1995—2005是研究區(qū)35年來(lái)用地結(jié)構(gòu)變化最快的時(shí)段，2005—2015年間變化速度減慢。從用地的轉(zhuǎn)移情況可知，1980—1995年期間各用地的轉(zhuǎn)移量較少;1995—2005年間最顯著的特征是耕地轉(zhuǎn)給其他地類(lèi)的量最多，達(dá)22.36km2，2005—2015年各用地轉(zhuǎn)移趨向縮減，轉(zhuǎn)移不顯著。35年來(lái)研究區(qū)土地的轉(zhuǎn)出量耕地>草地>林地>水域>建設(shè)用地，其中，耕地的轉(zhuǎn)出量遠(yuǎn)多余其他用地，達(dá)33.21km2。

關(guān)鍵詞：土地利用;時(shí)空變化;GIS;轉(zhuǎn)移矩陣

中圖分類(lèi)號(hào) F301.2文獻(xiàn)標(biāo)識(shí)碼 A文章編號(hào) 1007-7731（2019）24-0112-05

Abstract：Based on Landsat TM remote sensing image interpretation data of 7 periods from 1980 to 2015，the spatial and temporal evolution pattern of land use at county scale in typical karst areas was quantitatively analyzed with geostatistical method. The results show that：（1）In the past 35 years，woodland、 cultivated land and grassland are the absolute dominant land use types in the study area，and the total proportion of the three areas is more than 98%，and the overall change is small. （2）The intensity of land use change in the study area is construction land>water area>woodland>unused land>grassland>cultivated land in turn;from the perspective of comprehensive land use，the 10-year period 1995-2005 is the fastest period of land use structure change in the study area in 35 years，and the change rate slows down in 2005-2015. （3）From the situation of land transfer，it can be seen that the amount of land transfer from 1980 to 1995 is relatively small;the most significant feature in the 10 years from 1995 to 2005 is that the amount of cultivated land transferred to other land types is the largest，reaching 22.36 km2. The trend of land transfer from 2005 to 2015 is decreasing，but the transfer is not very significant. In the past 35 years，the amount of land transferred from the study area was cultivated land>grassland>woodland>water area>construction land. The amount of arable land transferred is far more than that of other land，reaching 33.21 km2.

Key words：Land use;Spatial and temporal change;GIS;Transfer matrix

土地利用/覆被變化（LUCC）被視為人-地時(shí)空耦合復(fù)雜系統(tǒng)下非常重要的科學(xué)問(wèn)題。作為人類(lèi)廣泛生產(chǎn)活動(dòng)與天然環(huán)境作用下的產(chǎn)物，土地利用變化越來(lái)越成為人類(lèi)與環(huán)境協(xié)調(diào)發(fā)展的重點(diǎn)研究?jī)?nèi)容。與此同時(shí)，土地利用時(shí)空格局的動(dòng)態(tài)監(jiān)測(cè)與模擬分析對(duì)土地生態(tài)系統(tǒng)以及各子系統(tǒng)的良性循環(huán)具有突出和長(zhǎng)遠(yuǎn)的實(shí)際運(yùn)用價(jià)值[1-3]。

特定區(qū)域內(nèi)土地屬性的變化與一定范圍內(nèi)的社會(huì)實(shí)踐活動(dòng)有密切關(guān)聯(lián)。土地景觀格局演變的眾多研究中，以大中區(qū)域、典型地區(qū)、流域和地理交界地帶等地區(qū)居多[4-7]。在研究方法方面，有學(xué)者通過(guò)運(yùn)用分形理論和地理信息圖譜理論探討了區(qū)域的用地結(jié)構(gòu)及時(shí)空演變模式[8-9]。在技術(shù)硬件支持的基礎(chǔ)上，國(guó)內(nèi)學(xué)者通過(guò)計(jì)量地理和計(jì)算機(jī)模擬方法對(duì)土地利用變化也做了大量研究。通過(guò)運(yùn)用多種土地利用量化指數(shù)能夠更加客觀反映一個(gè)區(qū)域地段的景觀復(fù)雜變化[10]。通過(guò)解譯歷史時(shí)期影像數(shù)據(jù)的方法，利用計(jì)算機(jī)模擬未來(lái)的土地利用情景變化對(duì)區(qū)域土地資源配置具有一定的用地規(guī)劃參考作用[11]?？λ固氐貐^(qū)作為一個(gè)典型生態(tài)脆弱的水土保持與土地覆被變化研究區(qū)域，利用長(zhǎng)時(shí)段的遙感影像數(shù)據(jù)對(duì)其進(jìn)行土地利用時(shí)序-空間格局演化進(jìn)程分析，無(wú)論對(duì)區(qū)域的生態(tài)安全評(píng)估還是經(jīng)濟(jì)發(fā)展都具有深刻意義[12-13]。生態(tài)脆弱區(qū)的土地利用與變化有它獨(dú)特的演化特性[14]，而喀斯特山區(qū)因石漠化和水土流失以及土地退化等生態(tài)問(wèn)題也受到了一些學(xué)者的關(guān)注[15-16]。

土地形態(tài)、屬性和性質(zhì)等特征的變化和區(qū)域內(nèi)的生境具有較強(qiáng)的互動(dòng)反饋?zhàn)饔?，土地作為生態(tài)系統(tǒng)的一部分也反映著生境的動(dòng)態(tài)狀況。前人對(duì)喀斯特的研究多數(shù)集中于石漠化和水土保持[17-20]。本研究以鎮(zhèn)寧和關(guān)嶺縣域聯(lián)合尺度為整體研究單元，使用1980—2015年共7個(gè)時(shí)間節(jié)點(diǎn)的遙感影像解譯數(shù)據(jù)，通過(guò)ArcGIS10.6軟件對(duì)典型巖溶區(qū)土地利用演變進(jìn)行定量分析，旨在為生態(tài)脆弱區(qū)土地的科學(xué)合理規(guī)劃利用與管理提供理論依據(jù)。

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

1.1 研究區(qū)域 以典型喀斯特巖溶縣鎮(zhèn)寧和關(guān)嶺為整體研究單元。鎮(zhèn)寧縣和關(guān)嶺縣位于貴州省西南典型喀斯特地貌區(qū)（105°15′～106°1′E，25°25′～26°11′N(xiāo)），隸屬于安順市，巖溶地貌分布廣。其中鎮(zhèn)寧地處貴州高原向南傾斜的斜坡地帶，位于世界上最大的巖溶區(qū)-華南巖溶的中心位置，屬黔中丘原中山類(lèi)型，海拔自北向南逐漸降低，地形南北狹長(zhǎng)且跨度大。關(guān)嶺地處云貴高原東部的一段斜坡地帶，貴州省中部偏西南，地勢(shì)上西北面較高，東南側(cè)較低，東北面與鎮(zhèn)寧縣相銜接。關(guān)嶺喀斯特地貌分布廣泛，氣候垂直分異明顯。關(guān)嶺巖溶景觀大面積發(fā)育，石灰?guī)r分布比較突出。

1.2 研究數(shù)據(jù) 本研究所使用的Landsat TM遙感影像解譯數(shù)據(jù)集來(lái)源于中科院資源環(huán)境數(shù)據(jù)中心（http：//www.resdc.cn/），分別是1980、1990、1995、2000、2005、2010和2015年共7期數(shù)據(jù)，分辨率為30m，該精度符合本研究區(qū)域要求。根據(jù)實(shí)際情況，本研究采用一級(jí)類(lèi)地類(lèi)分類(lèi)，分別為耕地、林地、草地、水域、建設(shè)用地和未利用地。利用現(xiàn)有的研究區(qū)矢量邊界數(shù)據(jù)對(duì)獲取的柵格執(zhí)行裁剪，從而得出本文的土地利用分類(lèi)柵格數(shù)據(jù)集。使用ArcGIS10.6加載7期的柵格數(shù)據(jù)集，通過(guò)轉(zhuǎn)換工具得到研究區(qū)的矢量數(shù)據(jù)集。

2 研究方法

2.1 土地利用動(dòng)態(tài)度

2.1.1 單一土地利用動(dòng)態(tài)度 單一土地利用動(dòng)態(tài)度用來(lái)定量描述特定時(shí)間區(qū)間內(nèi)變化中的土地狀態(tài)[21]，公式為：

式中，Xi即是地類(lèi)i的單一土地利用動(dòng)態(tài)度，Xi值越小，表明其轉(zhuǎn)出的量越少;Lit2和Lit1表示的地類(lèi)i在t1和t2時(shí)間節(jié)點(diǎn)的面積;T為起始時(shí)間與終止時(shí)間之差。

2.1.2 綜合土地利用動(dòng)態(tài)度 綜合土地利用度用來(lái)反映土地對(duì)象在開(kāi)始時(shí)段與結(jié)束時(shí)段之差內(nèi)轉(zhuǎn)移的整體動(dòng)態(tài)過(guò)程[22]，公式為：

式中，LC表示對(duì)象區(qū)域t時(shí)間范圍內(nèi)的綜合土地利用動(dòng)態(tài)度，計(jì)算獲得的LC越大則用地的演化狀態(tài)越強(qiáng)烈;LUi-j代表開(kāi)始時(shí)間點(diǎn)至終止時(shí)間點(diǎn)某i型用地轉(zhuǎn)化到其他類(lèi)別用地的量;LUi表示研究時(shí)段內(nèi)初期某i類(lèi)土地的數(shù)量;t則是研究的時(shí)間跨度。

2.2 轉(zhuǎn)移矩陣 土地動(dòng)態(tài)轉(zhuǎn)移矩陣通過(guò)對(duì)一定時(shí)間范圍內(nèi)的用地土地轉(zhuǎn)移進(jìn)行定量研究，可以宏觀量化研究區(qū)的各地類(lèi)的轉(zhuǎn)移情況，還能預(yù)測(cè)未來(lái)的變化趨勢(shì)，能夠較好地反映土地的轉(zhuǎn)移數(shù)量和方向[23]。

3 結(jié)果與分析

3.1 鎮(zhèn)寧關(guān)嶺1980—2015年用地格局 通過(guò)ArcGIS10.6軟件平臺(tái)對(duì)1980、1990、1995、2000、2005、2010和2015年7個(gè)時(shí)間節(jié)點(diǎn)的遙感解譯數(shù)據(jù)進(jìn)行幾何統(tǒng)計(jì)和可視化分析，得到研究對(duì)象在各時(shí)間節(jié)點(diǎn)的用地分布格局（圖1）。從用地的區(qū)域結(jié)構(gòu)來(lái)看，鎮(zhèn)寧關(guān)嶺的林地、耕地和草地比重顯著大于其他類(lèi)型，這3種地類(lèi)廣泛散布在研究區(qū)內(nèi)，并且面積長(zhǎng)期穩(wěn)定居于主導(dǎo)位置，林地在全區(qū)域廣泛分布，又以東南部部居多，其次是耕地和林地幾乎交錯(cuò)分布，耕地分布的區(qū)域幾乎都有林地。耕地則大片分布于海拔相差較小的北部區(qū)域，西面和東北面較高海拔區(qū)也分布較多。耕地總體分散在西北側(cè)。草地分布在研究區(qū)的中心位置分布較多，西側(cè)海拔較高區(qū)域和南端分布明顯少于其他方位。35年來(lái)草地的分布是較為穩(wěn)定的。其次與社會(huì)經(jīng)濟(jì)活動(dòng)密切的建設(shè)用地也團(tuán)聚在北部海拔差距較小的位置，以及靠近耕地和水域的中西部位。2010年作為建設(shè)用地發(fā)展的北方建設(shè)用地呈成長(zhǎng)條形迅速擴(kuò)張，同時(shí)向周邊分散發(fā)展，主要向耕地地區(qū)擴(kuò)張。鎮(zhèn)寧和關(guān)嶺相比較而言，社會(huì)經(jīng)濟(jì)活動(dòng)的區(qū)域也是集中于鎮(zhèn)寧的北部地區(qū)，其次是關(guān)嶺的中西側(cè)零星散布。水域主體分布在偏北部地區(qū)，中西部也有少量獨(dú)立的水域區(qū)，從研究區(qū)整體來(lái)看，水域面積十分稀少，空間上的分布也不均勻。在時(shí)間演進(jìn)基礎(chǔ)上來(lái)看，據(jù)表1可知，近35年來(lái)，鎮(zhèn)寧-關(guān)嶺的土地利用結(jié)構(gòu)以林地、耕地和草地為主，面積從大小依次是：林地>耕地>草地>建設(shè)用地>水域>未利用地，這3種地類(lèi)所占的比重穩(wěn)定達(dá)到了90%以上，林地占比在50%左右，變動(dòng)幅度十分微小。面積最小的分別是未利用地，水域和建設(shè)用地，35年來(lái)這3種用地的總面積和占比分別為87.7km2和2.76%。耕地和草地的數(shù)量總體表現(xiàn)為下降趨勢(shì)，1980—2015年兩者的數(shù)量分別縮減29.06km2和4.53km2。耕地面積每年平均減少0.83km2，草地每年平均減少0.13km2。耕地?cái)?shù)量的減少最多。林地總體為上升趨勢(shì)，1980—2015年增加了22.4797km2，平均每年增加0.64km2。水域數(shù)量變化比較特殊，1980—2005年總體數(shù)量在縮減，1980年的數(shù)量是6.38km2，而2010年銳減到3.65km2，減幅達(dá)到50%左右。值得注意的是，2010年北部建設(shè)用地區(qū)域出現(xiàn)大塊水域面積，這可能是人工湖，水域在2005—2010年有所增加并保持穩(wěn)定。從時(shí)空演變格局來(lái)看，研究區(qū)近35年來(lái)用地結(jié)構(gòu)變化最突出的是耕地，這與城區(qū)擴(kuò)建、交通運(yùn)輸用地的建設(shè)、水土流失和荒漠化有密切聯(lián)系。因此加強(qiáng)水土保持和荒漠化治理工作極為重要，這直接關(guān)系到當(dāng)?shù)厣鷳B(tài)的可持續(xù)發(fā)展和人地協(xié)調(diào)。

3.2 土地利用動(dòng)態(tài)度

3.2.1 單一土地利用動(dòng)態(tài)度 根據(jù)公式（1）可算出研究范圍內(nèi)1980—2015年分時(shí)段的單一土地動(dòng)態(tài)結(jié)果（表2）。近35年來(lái)，研究區(qū)的土地利用變化穩(wěn)定程度耕地>草地>未利用地>林地>水域>建設(shè)用地。變化最明顯的是2005—2015年的水域和建設(shè)用地。

3.2.2 綜合土地利用動(dòng)態(tài)度（LC） 根據(jù)公式可計(jì)算出研究區(qū)內(nèi)1980—2015年不同時(shí)段的LC值。1980—1995的LC值僅0.02%，表明該時(shí)期用地結(jié)構(gòu)變化較慢。1995—2005的10年間，用地變化最快。2005—2015年變化速度減緩。35年間的LC為0.09%，說(shuō)明鎮(zhèn)寧關(guān)嶺縣整體的土地利用處于緩慢演變階段。

3.3 土地利用轉(zhuǎn)移矩陣分析 根據(jù)7期遙感解譯數(shù)據(jù)制成1980—1995、1995—2005、2005—2015和1980—2015的轉(zhuǎn)移矩陣，從細(xì)節(jié)到整體分析研究區(qū)35年來(lái)用地的轉(zhuǎn)移情況。從表3能夠得知，1980—1995的15年間，研究區(qū)地類(lèi)轉(zhuǎn)移以草地最為顯著。15年間草地共轉(zhuǎn)出2.38km2，其中轉(zhuǎn)為林地最多，為2.36km2。水域轉(zhuǎn)移量也較顯著，每轉(zhuǎn)出2.72km2，就有1km2轉(zhuǎn)為林地。轉(zhuǎn)變量最少的是耕地、林地和建設(shè)用地，3種用地分別的總轉(zhuǎn)變量不到1km2?？傮w來(lái)看，1980—1995年的土地轉(zhuǎn)變量不明顯。本時(shí)期變化最顯著的用地是草地。

從表4可以看出，耕地僅僅在這10年的轉(zhuǎn)移量最大，達(dá)22.36km2，其中轉(zhuǎn)移最顯著的是由耕地轉(zhuǎn)出19.30km2到林地，而林地轉(zhuǎn)移給耕地的僅有1.34km2，這也是本階段林地面積增加的一個(gè)重要原因。較上一個(gè)時(shí)期，草地轉(zhuǎn)給耕地的面積明顯增加，達(dá)到1.77km2，而草地總轉(zhuǎn)移量也達(dá)到了4.54km2。此時(shí)水域的轉(zhuǎn)移最少，轉(zhuǎn)移最顯著的是耕地。

從表5可知，2005—2015這10年間用地轉(zhuǎn)化的動(dòng)態(tài)進(jìn)程整體較小。耕地是轉(zhuǎn)出量最多的一種地類(lèi)，轉(zhuǎn)移量共計(jì)11km2，耕地轉(zhuǎn)出量從小到大依次是建設(shè)用地>水域>草地>林地。這10年是建設(shè)用地快速發(fā)展的階段，每一種用地或多或少地都有轉(zhuǎn)進(jìn)。當(dāng)然，以耕地的轉(zhuǎn)進(jìn)最多。

從表6可以看出，研究區(qū)1980—2015年，耕地轉(zhuǎn)出面咱們最大，共計(jì)達(dá)到33.21km2，其中由耕地轉(zhuǎn)移到林地的20.44km2最為顯著。林地的轉(zhuǎn)移量次之，共計(jì)4.59km2，大部分林地轉(zhuǎn)為建設(shè)用地。草地的轉(zhuǎn)移也比較顯著，共轉(zhuǎn)出10.12km2，其中2.14km2轉(zhuǎn)成建設(shè)用地，2015年有0.27km2的草地轉(zhuǎn)成了未利用地。35年間建設(shè)用地的轉(zhuǎn)進(jìn)量主要來(lái)自于耕地、林地和草地，3者轉(zhuǎn)為建設(shè)用地共10.47km2。建設(shè)用地的轉(zhuǎn)移比較少，這可能是由于該用地主要是人居生活場(chǎng)所用地，較為固定，很難發(fā)生轉(zhuǎn)移。水域在35年間的轉(zhuǎn)移方向主要是林地、草地和耕地。耕地轉(zhuǎn)出量最顯著，轉(zhuǎn)給林地和建設(shè)用地的最多，分別達(dá)到20.44km2和5.60km2。

4 結(jié)論

本研究將典型巖溶區(qū)鎮(zhèn)寧縣與關(guān)嶺縣作為整體分析主體，通過(guò)遙感解譯數(shù)據(jù)結(jié)合數(shù)理統(tǒng)計(jì)方法與用地轉(zhuǎn)移指標(biāo)量化了鎮(zhèn)寧關(guān)嶺35年來(lái)土地演變的動(dòng)態(tài)過(guò)程。結(jié)果表明，研究區(qū)以林地、耕地、草地3類(lèi)土地為主體，3者比例明顯多余其他用地類(lèi)型比例;35年來(lái)，鎮(zhèn)寧和關(guān)嶺2縣耕地面積的減少最為突出，達(dá)29.06km2。建設(shè)用地以增長(zhǎng)數(shù)量10.21km2成為增量最大的地類(lèi);從不同用地的動(dòng)態(tài)度來(lái)看，研究區(qū)的用地穩(wěn)定程度依次是耕地>草地>未利用地>林地>水域>建設(shè)用地;用地結(jié)構(gòu)變化最顯著的時(shí)段是1995—2005年。1980—1995的用地變化十分微弱，2005—2015年用地變化程度則逐漸減緩。通過(guò)土地利用轉(zhuǎn)移矩陣可得知，研究區(qū)在1995—2005的10年間土地轉(zhuǎn)移量最大，為28.37km2。耕地在1995—2005年轉(zhuǎn)為林地最顯著。就研究區(qū)1980—2015年整體來(lái)看，耕地的轉(zhuǎn)出量最大，達(dá)33.21km2;而轉(zhuǎn)進(jìn)量較少，其中有20.44km2轉(zhuǎn)為林地，這也是林地面積增多的主要用地來(lái)源。

參考文獻(xiàn)

[1]Turner II B L，Lambin E F，Reenberg A.The emergence of land change science for global environmental change and sustainability.Proceedings of the National Academy of Sciences of the UnitedStates of America，2007，104（52）：20666-20671.

[2]劉紀(jì)遠(yuǎn)，張?jiān)鱿?，徐新良，?21世紀(jì)初中國(guó)土地利用變化的空間格局與驅(qū)動(dòng)力分析[J].地理學(xué)報(bào)，2010，64（12）：1411-1420.

[3]Sanneke van Asselen，Peter H.Verburg.Land cover change or land-use intensification： simulating land system change with a global-scale land change model[J].Global Change Biology 2013，19（12）： 3648-3667.

[4]劉紀(jì)遠(yuǎn)，匡文慧，張?jiān)鱿?，?20世紀(jì)80年代末以來(lái)中國(guó)土地利用變化的基本特征與空間格局[J].地理學(xué)報(bào)，2014，69（1）：3-14.

[5]宋開(kāi)山，劉殿偉，王宗明，等.1954年以來(lái)三江平原土地利用變化及驅(qū)動(dòng)力[J].地理學(xué)報(bào)，2008，63（1）：93-104.

[6]程建，程久苗，吳九興，等.2000～2010年長(zhǎng)江流域土地利用變化與生態(tài)系統(tǒng)服務(wù)功能變化[J].長(zhǎng)江流域資源與環(huán)境，2017，26（6）：894-901.

[7]蒲羅曼，張樹(shù)文，李飛，等.近40年?yáng)|北農(nóng)牧交錯(cuò)帶土地利用變化對(duì)人類(lèi)活動(dòng)的響應(yīng)——以吉林省西部地區(qū)為例[J].江蘇農(nóng)業(yè)科學(xué)，2016，44（6）：522-525.

[8]龍靜，胡寶清，劉占明.基于分形理論的區(qū)域喀斯特地區(qū)土地利用動(dòng)態(tài)分析——以廣西都安瑤族自治縣為例[J].廣西師范學(xué)院學(xué)報(bào)（自然科學(xué)版），2011，28（3）：50-56.

[9]龔文峰，袁力.土地利用及生態(tài)服務(wù)價(jià)值變化的圖形信息特征分析——以松花江流域哈爾濱段為例[J].環(huán)境科學(xué)與技術(shù)，2010，33（8）：200-205.

[10]羅婭，董國(guó)濤，宋文龍，等.黃河河口鎮(zhèn)—潼關(guān)區(qū)間1998-2010年土地利用變化特征[J].地理學(xué)報(bào)，2014，69（1）：42-53.

[11]李志，劉文兆，鄭粉莉.基于CA-Markov模型的黃土塬區(qū)黑河流域土地利用變化[J].農(nóng)業(yè)工程學(xué)報(bào)，2010，26（1）：346-352.

[12]許月卿，羅鼎，馮艷，等.西南喀斯特山區(qū)土地利用/覆被變化研究——以貴州省貓?zhí)恿饔驗(yàn)槔齕J].資源科學(xué)，2010，32（9）：1752-1760.

[13]冉晨，白曉永，譚秋，等.典型喀斯特地區(qū)石漠化景觀格局對(duì)土地利用變化的響應(yīng)[J].生態(tài)學(xué)報(bào)，2018，38（24）：240-249.

[14]Yan X ，Cai Y L . Multi-Scale Anthropogenic Driving Forces of Karst Rocky Desertification in Southwest China[J]. Land Degradation & Development，2015，26（2）：193-200.

[15]胡馨月.喀斯特石漠化貧困地區(qū)文化生態(tài)發(fā)展研究——以貴州省平壩縣桃花村為例[J].湖北農(nóng)業(yè)科學(xué)，2015，54（2）：474-479.

[16]劉愿理，廖和平，巫芯宇，等.西南喀斯特地區(qū)耕地破碎與貧困的空間耦合關(guān)系研究[J].西南大學(xué)學(xué)報(bào)（自然科學(xué)版），2019，41（01）：10-20.

[17]習(xí)慧鵬，王世杰，白曉永，等.西南典型喀斯特地區(qū)石漠化時(shí)空演變特征——以貴州省普定縣為例[J].生態(tài)學(xué)報(bào)，2018，38（24）：258-272.

[18]張琦，蔡雄飛，汪發(fā)勇，等.喀斯特石漠化與水土流失強(qiáng)度耦合關(guān)系分析—以六盤(pán)水市為例[J].亞熱帶水土保持，2016，28（3）：7-11.

[19]熊康寧，李晉，龍明忠.典型喀斯特石漠化治理區(qū)水土流失特征與關(guān)鍵問(wèn)題[J].地理學(xué)報(bào)，2012，67（7）：878-888.

[20]王世杰，李陽(yáng)兵，李瑞玲.喀斯特石漠化的形成背景、演化與治理[J].第四紀(jì)研究，2003，23（6）：657-666.

[21]童小容，楊慶媛，畢國(guó)華.重慶市2000～2015年土地利用變化時(shí)空特征分析[J].長(zhǎng)江流域資源與環(huán)境，2018，27（11）：89-103.

[22]馬禮，蘇海霞.河北沽源縣土地利用動(dòng)態(tài)變化研究[J].干旱區(qū)資源與環(huán)境，2010，24（1）：131-136.

[23]張超，王琳，張秋霞，等.布哈河流域土地利用轉(zhuǎn)移矩陣及空間變化研究[J].水利水電技術(shù)，2016，47（5）：6-11.

（責(zé)編：張 麗）