孫曉莉，鄭 毅，趙 然，沈金祥，田淑靜，費(fèi)麗娜

(1.云南國(guó)土資源職業(yè)學(xué)院，昆明 652501；2.徐州工程學(xué)院管理學(xué)院，江蘇 徐州 221000)

【研究意義】土地是人類生產(chǎn)生活和社會(huì)發(fā)展過(guò)程中的重要自然資源和物質(zhì)保障[1]，不同的土地利用類型在空間分布上的變化直接影響景觀格局變動(dòng)[2]。通過(guò)對(duì)流域景觀格局研究有助于了解景觀格局演變的機(jī)制與規(guī)律，進(jìn)一步探索自然、人類活動(dòng)等與景觀變化的關(guān)系，實(shí)現(xiàn)土地資源可持續(xù)利用[3]。【前人研究進(jìn)展】目前，流域土地利用景觀格局研究已引起國(guó)內(nèi)外廣泛關(guān)注，主要集中在流域景觀格局及其時(shí)空變化[4-7],流域土地利用景觀格局及其穩(wěn)定性分析[8]，流域景觀格局變化對(duì)流域植被[9]、生態(tài)風(fēng)險(xiǎn)[10]、水質(zhì)變化[11]影響等方面。流域土地利用景觀格局變化研究對(duì)流域綜合治理和景觀調(diào)控等具有非常重要的理論和現(xiàn)實(shí)意義[12]?！颈狙芯壳腥朦c(diǎn)】學(xué)者對(duì)滇池、劍湖、洱海等流域研究較多，對(duì)陽(yáng)宗海流域開(kāi)展相關(guān)研究的較少。因此，以陽(yáng)宗海流域作為研究區(qū)域，開(kāi)展流域土地利用景觀格局時(shí)空變化分析對(duì)流域內(nèi)土地利用、湖泊保護(hù)等具有重要意義?！緮M解決的關(guān)鍵問(wèn)題】文中選擇2005年、2010年、2015年、2020年4期遙感影像數(shù)據(jù)，采用目視解譯和外業(yè)核查方式獲得研究區(qū)矢量數(shù)據(jù)圖斑，從土地利用變化、地學(xué)信息圖譜、景觀格局變化不同角度分析研究區(qū)土地利用景觀格局變化特征；從景觀空間角度探討陽(yáng)宗海湖泊面積及質(zhì)心變化，以期為陽(yáng)宗海湖泊及流域內(nèi)土地利用、生態(tài)安全體系構(gòu)建等提供數(shù)據(jù)支撐和決策依據(jù)，進(jìn)一步促進(jìn)湖泊流域生態(tài)環(huán)境的有效保護(hù)、合理利用與可持續(xù)發(fā)展。

1 材料與方法

1.1 研究區(qū)概況

陽(yáng)宗海(102°58′～103°01′E，24°51′～24°58′N)位于云南省昆明市以東30 km(圖1)，海拔1770 m，地跨澄江、呈貢、宜良3個(gè)縣(區(qū))，是云南九大高原湖泊之一[13]。陽(yáng)宗海湖區(qū)屬亞熱帶高原季風(fēng)氣候，年均氣溫14.5 ℃，年均降水量963.5 mm。

圖1 研究區(qū)位置示意圖

陽(yáng)宗海流域?qū)僦榻饔蚰媳P江水系，是典型的高原構(gòu)造型亞熱帶石灰?guī)r淡水深水湖泊，流域面積為34 050.45 hm2。陽(yáng)宗海湖泊入湖河流主要有陽(yáng)宗河、七里河、湯池河以及其它臨時(shí)性溪流等。地勢(shì)總體上為北高南低，其四面環(huán)山，山系多呈南北走向[14]。

1.2 數(shù)據(jù)來(lái)源與處理

數(shù)據(jù)來(lái)源：2005年、2010年、2015年、2020年4期陽(yáng)宗海流域土地利用數(shù)據(jù)來(lái)源于遙感影像，圖像空間分辨率為0.5 m；行政區(qū)劃等數(shù)據(jù)來(lái)源于云南省第三次全國(guó)國(guó)土調(diào)查數(shù)據(jù)。

數(shù)據(jù)處理：通過(guò)對(duì)獲取的遙感影像進(jìn)行幾何校正、圖像融合、裁剪等處理，依據(jù)GB/T21010—2007《土地利用現(xiàn)狀分類》，將陽(yáng)宗海流域的土地利用類型劃分為耕地、園地、林地、草地、建設(shè)用地、水域及水利設(shè)施用地、未利用地7類。按照上述分類系統(tǒng)采用目視解譯及外業(yè)核查的方式獲取2005年、2010年、2015年、2020年4期陽(yáng)宗海流域景觀類型的矢量數(shù)據(jù)圖斑(圖2)。

圖2 研究區(qū)2005—2020年土地利用類型

主要基于ArcGIS 10.0和Fragstats 4.2等軟件對(duì)數(shù)據(jù)進(jìn)行處理,為下文研究提供數(shù)據(jù)基礎(chǔ)。

1.3 研究方法

1.3.1 土地利用轉(zhuǎn)移矩陣 土地利用轉(zhuǎn)移矩陣可以描述各地類之間變化的來(lái)源和流向，計(jì)算出不同土地利用類型之間相互轉(zhuǎn)移的方向和數(shù)量，通過(guò)對(duì)某一時(shí)期內(nèi)各種土地利用類型的演變方向和程度進(jìn)行相應(yīng)描述[15-16]，清晰的表達(dá)土地利用狀態(tài)的時(shí)空演化過(guò)程[1]。其表達(dá)式為：

(1)

式中：S為土地總面積；n為土地利用的類型數(shù)；i，j分別為研究期初期和末期的土地利用類型。

1.3.2 地學(xué)信息圖譜 地學(xué)信息圖譜具有同時(shí)表達(dá)地理空間結(jié)構(gòu)特征和時(shí)間動(dòng)態(tài)變化的功能[17]。利用地學(xué)信息圖譜方法研究陽(yáng)宗海流域土地利用時(shí)空格局變化，其本質(zhì)是在GIS和遙感(Remote Sensing，RS)的技術(shù)支持下，對(duì)研究區(qū)的土地利用地類圖斑數(shù)據(jù)進(jìn)行相應(yīng)的提取[18]，并采用空間分析等方法生成一系列具有聯(lián)系的土地利用地類圖斑變化信息。文中主要采用地學(xué)信息圖譜對(duì)研究區(qū)各土地利用類型時(shí)空變化特征進(jìn)行分析。

1.3.3 景觀格局指數(shù) 景觀格局指數(shù)是能夠高度濃縮景觀格局信息并反映其結(jié)構(gòu)組成和空間分布特征的定量指標(biāo)[19]。為了既全面反映研究區(qū)景觀格局時(shí)空變化特征，又有效的減少信息冗余。文中借助于ArcGIS 10.0和Fragstats 4.2軟件對(duì)各地類圖斑進(jìn)行相應(yīng)處理，主要選擇斑塊個(gè)數(shù)(Number of patches，NP)、平均斑塊面積(Mean patch size，MPS)、形狀指數(shù)(Landscape shape index，LSI)、聚合度(Aggregation index，AI)、蔓延度(CONTAG)、香農(nóng)均勻度(Shannon’s eveness index，SHEI)、香農(nóng)多樣性指數(shù)(Shannon’s diversity，SHDI)7個(gè)指標(biāo)[20]，通過(guò)景觀格局信息圖譜分析研究區(qū)2005—2020年土地利用景觀及陽(yáng)宗海湖面景觀時(shí)空格局變化。

1.3.4 景觀質(zhì)心量測(cè) 景觀質(zhì)心量測(cè)主要用于研究各類景觀的動(dòng)態(tài)變化，通過(guò)計(jì)算景觀質(zhì)心分析景觀格局的空間變化規(guī)律和趨勢(shì)[12]。其計(jì)算公式如下[21]。

(2)

式中：xt、yt分別第t年表示湖泊分布質(zhì)心的經(jīng)緯度坐標(biāo)；Cti表示第t年湖泊斑塊的面積;xi、yi分別表示第i個(gè)湖泊班塊的經(jīng)緯度坐標(biāo)。

2 結(jié)果與分析

2.1 陽(yáng)宗海流域土地利用變化

2005—2020年建設(shè)用地和林地面積增加，耕地和未利用地面積減少，其他地類面積變化不大。受“退耕還林”政策實(shí)施影響，2010—2015年耕地面積出現(xiàn)大幅度減少，達(dá)1186.24 hm2；林地面積出現(xiàn)顯著增加，達(dá)970.47 hm2。2010—2020年，隨著學(xué)校、華僑城溫泉水公園以及周邊配套設(shè)施的建設(shè)，研究區(qū)內(nèi)建設(shè)用地增幅明顯，達(dá)683.42 hm2；2015—2020年隨著陽(yáng)宗海南岸“一帶三區(qū)”濕地項(xiàng)目打造，實(shí)施“退耕還湖還濕還林”政策，研究區(qū)內(nèi)水域及水利設(shè)施用地出現(xiàn)小幅增加，達(dá)98.46 hm2。

從表1可以看出，陽(yáng)宗海流域土地利用景觀類型以林地和耕地為主，二者景觀類型面積占總面積的比例達(dá)65%以上，其中又以林地面積占比最大，2005—2020年間林地面積占研究區(qū)總面積的比例均在40%以上，其次為耕地，占研究區(qū)總面積的比例均在20%以上。此外，未利用地面積占研究區(qū)總面積的比例在15%左右，充分說(shuō)明研究區(qū)以農(nóng)林生產(chǎn)為主，水域用地充沛，隨著周邊學(xué)校、溫泉水上公園、陽(yáng)宗海南岸濱海生態(tài)濕地等的建設(shè)，研究區(qū)建設(shè)用地呈現(xiàn)逐年增加趨勢(shì)。

表1 研究區(qū)土地利用景觀類型面積

2.2 陽(yáng)宗海流域土地利用地學(xué)信息圖譜分析

陽(yáng)宗海流域土地利用景觀時(shí)空格局變化可以通過(guò)地學(xué)信息圖譜進(jìn)行分析，一般包括增勢(shì)圖譜和減勢(shì)圖譜[19]。增勢(shì)圖譜是指由其他土地利用資源轉(zhuǎn)入各土地利用資源的部分，稱為新增地類來(lái)源；減勢(shì)圖譜是指各土地利用資源轉(zhuǎn)為其他土地利用類型的部分，稱為各地類資源減少流向[22]。

從表2可看出，2005—2020年，各土地利用類型面積均發(fā)生變化。從總體變化來(lái)看：耕地、園地和未利用地面積出現(xiàn)減少，其中耕地和未利用地面積減少較大，分別為816.45、707.92 hm2，園地減少規(guī)模較小，為33.00 hm2；林地、草地、建設(shè)用地和水域及水利設(shè)施用地出現(xiàn)增加，其中建設(shè)用地和林地面積增加較大，分別為809.31、627.23 hm2，草地和水域及水利設(shè)施用地增加規(guī)模較小，分別為67.35、53.48 hm2。

從表2和圖3可知，耕地增加來(lái)源地類主要為林地和未利用地，二者占比達(dá)85.15%，除陽(yáng)宗海東岸耕地增加較少之外，其他區(qū)域均有分布，說(shuō)明隨著人口增加，主要通過(guò)開(kāi)墾未利用地等促使耕地增加；園地增加來(lái)源地類主要為耕地和未利用地，二者占比達(dá)74.84%，主要分布于陽(yáng)宗海西南側(cè)，少量分布于陽(yáng)宗海流域最北部等區(qū)域，馬郎櫻桃文化旅游節(jié)帶動(dòng)了當(dāng)?shù)厝罕姺N植櫻桃等當(dāng)?shù)厮?，使得?yáng)宗海西南側(cè)等部分耕地流向園地；林地增加地類來(lái)源主要為耕地和未利用地，二者占比達(dá)87.16%，整個(gè)研究區(qū)域均有分布，主要是受到當(dāng)?shù)亍巴烁€林”政策實(shí)施影響，通過(guò)退耕和開(kāi)墾質(zhì)量較好的未利用地進(jìn)行林地種植；草地增加來(lái)源地類主要為未利用地和林地，二者占比達(dá)61.45%，其次為耕地，主要分布于陽(yáng)宗海周邊區(qū)域，較為分散。距離居民地較遠(yuǎn)、且耕種條件差的耕地容易轉(zhuǎn)化為草地，此外，判讀時(shí)容易受到“林下草”的影響，即在影像判讀過(guò)程中由于林地種植具有周期性，當(dāng)樹(shù)木較小時(shí)，影像判讀為草地，使得新增草地來(lái)源于林地，而當(dāng)樹(shù)木生長(zhǎng)茂密時(shí)，出現(xiàn)遮蓋草地的情況，使得影像判讀為林地，出現(xiàn)草地流向林地的情況；建設(shè)用地增加來(lái)源地類主要為耕地、未利用地和林地，三者占比達(dá)95.84%，建設(shè)用地增加除了新修建的道路等基礎(chǔ)配套設(shè)施外，主要為陽(yáng)宗海西南側(cè)的云南國(guó)土資源職業(yè)學(xué)院、陽(yáng)宗海北部的華僑城溫泉公園，建設(shè)用地的增加占用較多的耕地、未利用地和林地等；水域及水利設(shè)施用地增加來(lái)源地類主要為耕地和林地，二者占比達(dá)66.76%，主要表現(xiàn)為坑塘水面及陽(yáng)宗海沿岸部分“退耕退林還湖”增加；未利用地增加來(lái)源地類主要為林地和耕地，二者占比達(dá)87.40%，整個(gè)研究區(qū)域均有分布，主要包括耕種條件較差、距居民地較遠(yuǎn)等導(dǎo)致耕種荒廢的耕地以及受到天氣等因素影響退化為未利用地的灌木林地等。

圖3 2005—2020年土地利用增勢(shì)圖譜

表2 研究區(qū)2005—2020年土地利用轉(zhuǎn)移矩陣

從表2和圖4可知，耕地流向的地類主要為林地和建設(shè)用地，二者占比為64.27%，主要原因?yàn)樾略鼋ㄔO(shè)用地占用耕地及退耕還林政策的影響；園地流向的地類主要為耕地和林地，二者占比為74.74%，主要受到當(dāng)?shù)厝罕姺N植果園和農(nóng)作物的影響；林地流向的地類主要為未利用地和耕地，二者占比為81.10%，受天氣等影響,一部分灌木林轉(zhuǎn)化為未利用地，雖當(dāng)?shù)貙?shí)施退耕還林政策，但仍存在將林地開(kāi)墾為耕地的情況；草地流向的地類主要為未利用地和林地，二者占比為58.33%，一部分草地受到天氣等因素影響轉(zhuǎn)化為未利用地，此外受到上述“林下草”因素影響；建設(shè)用地流向的地類主要為耕地和林地，二者占比為66.13%，主要原因?yàn)椴糠洲r(nóng)村建設(shè)用地復(fù)墾為耕地，采礦用地復(fù)墾為林地；水域及水利設(shè)施用地流向的地類主要為耕地和建設(shè)用地，二者占比為56.57%，主要原因?yàn)閲焯锸沟貌糠炙蜣D(zhuǎn)化為耕地，農(nóng)村居民點(diǎn)用地占用部分坑塘水面等；未利用地流向的地類包含所有地類，主要為林地、耕地和建設(shè)用地，三者占比為94.56%，存在新增建設(shè)用地占用未利用地以及將未利用地開(kāi)墾為耕地和林地的情況。

圖4 2005—2020年土地利用減勢(shì)圖譜

2.3 陽(yáng)宗海流域空間格局景觀特征分析

將研究區(qū)4期矢量數(shù)據(jù)圖斑轉(zhuǎn)化為柵格數(shù)據(jù)(柵格單元大小為10 m×10 m)。通過(guò)選取7個(gè)景觀格局指數(shù)，從景觀尺度分析陽(yáng)宗海流域2005—2020年景觀格局變化特征并揭示其所蘊(yùn)含的景觀效應(yīng)[20]。

由圖5可見(jiàn)，陽(yáng)宗海流域景觀斑塊個(gè)數(shù)呈現(xiàn)先增加后減少趨勢(shì)，2005—2010年斑塊個(gè)數(shù)出現(xiàn)上升，2010—2020年斑塊個(gè)數(shù)出現(xiàn)減少，階段性特征明顯，其中2005—2010年斑塊個(gè)數(shù)增加，平均斑塊面積減少，表明研究區(qū)景觀格局受人類社會(huì)經(jīng)濟(jì)活動(dòng)干擾，景觀破碎化程度增加；2010—2020年斑塊個(gè)數(shù)減少，平均斑塊面積增加，表明研究區(qū)景觀格局逐步向著聚集方向發(fā)展。15年間研究區(qū)景觀蔓延度指數(shù)雖有小幅波動(dòng)，但整體變化不大，說(shuō)明研究區(qū)內(nèi)優(yōu)勢(shì)斑塊類型連接性較好，未出現(xiàn)明顯破壞情況，且隨著研究區(qū)道路等基礎(chǔ)設(shè)施不斷完善，進(jìn)一步增強(qiáng)了景觀連通性。聚合度呈現(xiàn)下降趨勢(shì)，但變化較緩，說(shuō)明研究區(qū)相同斑塊類型的聚集程度降低，有景觀破碎化趨勢(shì)，但不明顯。香農(nóng)多樣性指數(shù)和香農(nóng)均勻度雖出現(xiàn)一定的波動(dòng)變化，呈現(xiàn)先上升、再下降、然后上升的趨勢(shì)，但變化幅度不大，說(shuō)明研究區(qū)景觀多樣性和景觀優(yōu)勢(shì)度相對(duì)較為穩(wěn)定。總體而言，陽(yáng)宗海流域2005—2020年空間格局景觀特征具有一定的空間分異性，2010年以后研究區(qū)逐步進(jìn)行治理和保護(hù)，景觀聚集度、景觀多樣性和景觀優(yōu)勢(shì)度等均保持相對(duì)穩(wěn)定，人為因素干擾得到緩解，生態(tài)環(huán)境得到保護(hù)。

圖5 陽(yáng)宗海流域空間格局景觀指數(shù)變化

2.4 陽(yáng)宗海湖泊質(zhì)心量測(cè)

2.4.1 陽(yáng)宗海湖泊面積變化 利用研究區(qū)2005—2020年的陽(yáng)宗海湖泊矢量圖斑數(shù)據(jù)進(jìn)行統(tǒng)計(jì)得到陽(yáng)宗海湖泊面積，如表3所示。2005—2015年陽(yáng)宗海面積呈逐年減少趨勢(shì)，其中2005—2010年減少50.44 hm2，2010—2015年減少32.94 hm2。2010—2015年與2005—2010年相比縮減速度減緩，2015—2020年陽(yáng)宗海面積出現(xiàn)增加，面積達(dá)70.24 hm2?？傮w來(lái)看，2005—2020年呈現(xiàn)先減少后增加的趨勢(shì)。這與當(dāng)?shù)貙?shí)施《陽(yáng)宗海保護(hù)條例》、“退耕還湖”以及入湖河道治理等一系列舉措密不可分。

表3 陽(yáng)宗海湖泊景觀面積變化情況

面積加權(quán)平均形狀指數(shù)和分維數(shù)用來(lái)度量斑塊形狀的復(fù)雜程度，其值越高說(shuō)明形狀越復(fù)雜。從表3可得出，形狀指數(shù)和分維數(shù)2005—2010年出現(xiàn)減少，2010—2020年出現(xiàn)連續(xù)增加，但2015—2020年增幅較緩。陽(yáng)宗海2008年發(fā)生污染以來(lái)，引起當(dāng)?shù)卣叨戎匾?，?shí)施一系列舉措治理和保護(hù)陽(yáng)宗海。2010年相較于2005年、2015年出現(xiàn)較為重要的轉(zhuǎn)折點(diǎn)，尤其是陽(yáng)宗海南岸和北岸入湖口形狀變化較為明顯，這與形狀指數(shù)和分維數(shù)值反映的情況是吻合的。2015—2020年隨著陽(yáng)宗海濕地項(xiàng)目落地，陽(yáng)宗海湖泊形狀和分維數(shù)出現(xiàn)小幅增加。

2.4.2 陽(yáng)宗海湖泊質(zhì)心變化 質(zhì)心是描述地理空間對(duì)象分布的重要指標(biāo)，其可以用來(lái)表示地理空間對(duì)象分布所發(fā)生的變化情況。從圖6可得出，陽(yáng)宗海質(zhì)心2005—2010年總體向東北方向偏移3.51 m，2010—2015年總體向東北方向偏移15.92 m，2015—2020年總體向西南方向偏移22.74 m?？蓪㈥?yáng)宗海湖泊質(zhì)心偏移分為2個(gè)時(shí)期：2005—2015年整體向東北方向偏移，且2005—2010年期間偏移較為緩慢；2015—2020年質(zhì)心偏移方向發(fā)生變化，由2005—2015年的東北方向偏移向西南方向發(fā)生偏移，且偏移速度較快。陽(yáng)宗海南岸及湖泊周邊“泥沙淤積”到“圍湖造田”“圍湖建塘”再到“退耕還湖”的各時(shí)期轉(zhuǎn)變直接影響著陽(yáng)宗海湖泊質(zhì)心發(fā)生變化。

圖6 2005—2020年陽(yáng)宗海湖泊質(zhì)心變化

3 討 論

許多學(xué)者從陽(yáng)宗海流域重金屬污染[23]、生態(tài)環(huán)境質(zhì)量變化[24]等方面進(jìn)行分析研究，而從土地利用景觀格局角度對(duì)其開(kāi)展相關(guān)研究的較少。文中選擇2005年、2010年、2015年、2020年4期數(shù)據(jù)，采用土地利用轉(zhuǎn)移矩陣、景觀格局指數(shù)等方法對(duì)研究區(qū)土地利用數(shù)量、景觀格局指數(shù)、湖泊質(zhì)心變化進(jìn)行分析，研究結(jié)果與樊凱等[8]、朱泓等[24]的研究結(jié)論基本一致，2010—2020年，隨著人類活動(dòng)強(qiáng)度降低，斑塊數(shù)量、景觀多樣性指數(shù)和優(yōu)勢(shì)度等雖有小幅變化，但景觀結(jié)構(gòu)整體趨于穩(wěn)定，生態(tài)環(huán)境質(zhì)量持續(xù)變好。若能獲取更高分辨率、長(zhǎng)時(shí)間序列的影像數(shù)據(jù)進(jìn)行分析，文中研究結(jié)果可能會(huì)更具有系統(tǒng)性和針對(duì)性。流域土地利用景觀格局變化是一個(gè)長(zhǎng)期、動(dòng)態(tài)的過(guò)程，陽(yáng)宗海流域景觀格局變化是否受到氣溫、降雨量、植被等自然因素以及人類活動(dòng)因素、政策因素等驅(qū)動(dòng)機(jī)制影響，陽(yáng)宗海湖泊的水質(zhì)變化是否與流域內(nèi)土地利用景觀格局變化存在直接關(guān)聯(lián)等內(nèi)容將是今后進(jìn)一步探索和研究的方向。

4 結(jié) 論

從整體來(lái)看，陽(yáng)宗海流域在2005—2020年各地類既有增加，也有減少。隨著社會(huì)經(jīng)濟(jì)發(fā)展、人口增加及當(dāng)?shù)卣畬?shí)施“退耕還林還草”政策等因素影響，建設(shè)用地、林地呈現(xiàn)增加趨勢(shì)，未利用地、耕地面積減少，園地等其余地類面積變化幅度不大。2005—2010年，受人為等干擾因素影響，流域內(nèi)土地利用景觀破碎度增加，湖泊面積出現(xiàn)減少。2010—2020年，研究區(qū)景觀格局逐步向著聚集方向發(fā)展，景觀多樣性和景觀優(yōu)勢(shì)度相對(duì)較為穩(wěn)定，湖泊面積由減向增變化，生態(tài)環(huán)境治理和保護(hù)成效得以體現(xiàn)。此外，隨著陽(yáng)宗海湖泊周邊“退耕還湖”“河道治理”等措施實(shí)施，陽(yáng)宗海湖泊形狀和分維數(shù)越來(lái)越復(fù)雜，質(zhì)心發(fā)生偏移。本研究結(jié)果可為構(gòu)建湖泊流域生態(tài)安全體系提供數(shù)據(jù)支撐和決策依據(jù)，對(duì)促進(jìn)湖泊流域生態(tài)環(huán)境的有效保護(hù)、合理利用與可持續(xù)發(fā)展具有重要的意義[25]。