廖雨，冉丹陽，張麗芳，李佳玉，楊存建，

1. 四川師范大學地理與資源科學學院，成都 610101

2. 四川師范大學西南土地資源評價與監(jiān)測教育部重點實驗室，成都 610068

0 前言

國土空間是人類生產(chǎn)、生活的載體，也是我們進行生態(tài)文明建設(shè)的介質(zhì)[1]。城市發(fā)展、工業(yè)進步、鄉(xiāng)村振興，使國土空間格局發(fā)生了巨大的轉(zhuǎn)變，導致國土空間的利用與保護之間的矛盾日益突出，國土空間的持續(xù)開發(fā)迎來了新的挑戰(zhàn)[2-3]。隨著國土空間格局的不斷轉(zhuǎn)變，導致生物生存的環(huán)境趨于破碎，生境功能退化甚至喪失，降低了生物的多樣性，對生態(tài)環(huán)境造成巨大破壞[4-5]。因此，將國土空間格局變化與生境退化的研究結(jié)合起來對自然資源合理利用、生態(tài)系統(tǒng)的保護及區(qū)域可持續(xù)發(fā)展至關(guān)重要。

國外的專家學者認為將土地利用轉(zhuǎn)型與生態(tài)效應(yīng)結(jié)合起來對統(tǒng)籌城鄉(xiāng)發(fā)展及生態(tài)文明建設(shè)具有重要的指導意義[6-7]。而對土地利用分類進一步劃分為“三生”空間是研究國土空間的新趨勢，關(guān)于“三生”空間的研究，我國學者做了大量的研究并取得豐碩成果，如張紅旗[8]、劉繼來[9]、朱鳳武[10]、劉長青[11]、于莉[12]等?；谏鲜鰧＜业姆诸愃枷肱c研究，將不同的土地利用類型劃分為生產(chǎn)、生態(tài)與生活用地[13-14]。關(guān)于國土空間格局及生境退化的研究主要通過構(gòu)建生境評價體系及基于模型的評價方法兩種形式[15]。前者需要較高數(shù)據(jù)量，成本較大；后者基于模型評價方法，操作簡單，對數(shù)據(jù)需求量較低，較為直觀，空間可視化程度高，便于從景觀格局的角度出發(fā)來揭示生境退化程度[16-17]。所以本文主要通過INVEST模型中的生境質(zhì)量模塊來揭示不同景觀格局下的生境退化程度。

山地丘陵區(qū)是我國國土空間中重要的組成要素之一，約占我國陸地面積的 65%。由于其本身生態(tài)環(huán)境的脆弱性、敏感性較大，以及人地關(guān)系矛盾突出等條件約束，迫使我國山區(qū)國土空間開發(fā)面臨的問題更加復雜，人地關(guān)系更加突出[1]。林州市位于太行山的東麓，其境內(nèi)山地、丘陵占其總面積的 86%，近些年隨著其經(jīng)濟快速發(fā)展，精準扶貧政策的實施，使得國土空間的開發(fā)及生態(tài)環(huán)境面臨著巨大的挑戰(zhàn)。因此本文選取林州市為研究區(qū)域，探索其國土空間格局時空演變規(guī)律，揭示國土空間格局的變化導致景觀破碎化對其生境退化程度產(chǎn)生的影響。以期能為林州市今后的規(guī)劃發(fā)展提供科學的建議。

1 研究區(qū)域概況

林州市位于河南省的北部、太行山脈東麓，是河南、河北、山西三省交匯之處，主要地形為山地、丘陵，占林州市總面積的 86%。境內(nèi)海拔最高達1600 米左右，最低 200 米左右，海拔高差較大，其西北的地勢高于東南。地處半濕潤地區(qū)，屬于暖溫帶大陸性季風氣候，降水量季節(jié)分配不均勻，夏季高溫多雨，冬季寒冷干燥，溫度變化較大。轄15個鎮(zhèn)、1個鄉(xiāng)及林街道、開元街道、龍山街道、桂園街道等 4個街道辦事處，本文依據(jù)林州實際情況，將其 4個街道辦事處統(tǒng)稱為中心城區(qū)。截止2015年，林州市戶籍人口 112.63萬人，人口密度每平方千米 550人，是我國人口密度較高的縣級市之一。

2 數(shù)據(jù)來源與方法

2.1 數(shù)據(jù)來源

研究利用Landsat 5TM、 Landsat 8 OLI遙感影像數(shù)據(jù)進行圖像預處理后，依據(jù)中國土地利用類型分類體系[18]，通過人工目視解譯獲得林州市(1990年、2000年、2010年、2015年)四期土地利用數(shù)據(jù)，然后借助Google Earth平臺對目視解譯的數(shù)據(jù)進行了二次修正，數(shù)據(jù)驗證精度達到 85%以上，解譯精度可靠。其中土地利用類型分為兩級，6個一級地類(耕地、林地、草地、建設(shè)用地、水域、未利用地)，25個二級地類，基礎(chǔ)遙感影像從地理空間數(shù)據(jù)云網(wǎng)站(http: //www.gscloud.cn/)下載獲取。

2.2 “三生”空間劃分

國土空間是人與自然、社會經(jīng)濟相互作用而表現(xiàn)出的一種形態(tài)，它是一種由生態(tài)功能、生產(chǎn)功能、生活功能等多種功能構(gòu)成的綜合體[3，13]。目前國土空間的劃分依據(jù)也多種多樣，很多學者提出依據(jù)國土空間功能來劃分，但是一種國土空間可能兼有多種功能，以單一功能劃分存在一定的難度[3，11，13]，因此多數(shù)專家們主要依據(jù)國土空間的主體功能對國土空間進行劃分。本文主要借鑒張紅旗等[8]、朱鳳武等[10]的分類思想，依據(jù)國土空間的生產(chǎn)功能、生活功能、生態(tài)功能以及林州市的實際情況將其國土空間劃分為3個大類及7個小類別，如表1。

表1 林州市三生空間分類體系Table 1 Linzhou ecological-production-living spaces classifycation system

2.3 土地利用轉(zhuǎn)移矩陣

土地利用轉(zhuǎn)移矩陣表示系統(tǒng)分析中對系統(tǒng)狀態(tài)與狀態(tài)轉(zhuǎn)移的定量描述。它將土地轉(zhuǎn)移的面積利用矩陣的形式呈現(xiàn)出來，不僅反映研究區(qū)內(nèi)部的土地類型結(jié)構(gòu)，還追溯其地類變化的來源及去向[3，20]。矩陣中行表示k時期的i種土地利用類型，列表示k+1時期的土地利用類型，其表達式為:

如式(1)中所示，L為土地面積，n表示研究區(qū)內(nèi)土地利用類型數(shù)。i、j分別為研究期初與研究期末的土地利用類型。運用 Arcgis10.2軟件對林州市1990年、2000年、2010年、2015年4期土地利用數(shù)據(jù)進行分別疊加，然后用Excel進行數(shù)據(jù)透視表處理，得到 1990—2000年、2000—2010年、2010—2015年、1990—2015年土地轉(zhuǎn)移矩陣，以此分析林州市近25年來國土空間的結(jié)構(gòu)及用地功能類型轉(zhuǎn)化特征。

2.4 景觀指數(shù)法

景觀格局是對一定區(qū)域范圍內(nèi)的生態(tài)環(huán)境體系的綜合反映，利用景觀指數(shù)能夠從多維角度來反映其結(jié)構(gòu)組成和空間配置等方面的特征[19]。本文利用FRAGSTATS4.2軟件對選取的斑塊密度(PD)、斑塊平均面積(MPS)、蔓延度指數(shù)(CONTAG)、香濃多樣性指數(shù)(SHDI)、香濃均勻度指數(shù)(SHEI)、景觀形狀指數(shù)(LSI)等多個景觀指數(shù)進行計算，各個景觀指數(shù)的生態(tài)學含義參考文獻[22-24]。本文利用類型水平及景觀水平2個維度來分析林州市25年間國土空間景觀格局特征及其演變規(guī)律，并能清楚、直觀的反映林州市的景觀異質(zhì)性、連通性、破碎化程度及其景觀的多樣性。

2.5 生境退化度

本文主要運用 InVEST模型中的生境質(zhì)量模塊來評估林州市景觀水平上的生境退化特征[4，25]。其主要是利用土地利用類型與威脅源建立聯(lián)系，根據(jù)不同生境對威脅源的響應(yīng)程度，評估不同景觀格局下的生境分布和退化情況，計算得到生境質(zhì)量和稀缺程度可以反映該區(qū)域的生物多樣性，生境退化特征。

生境退化度的表達式如下:

如式(2)所示，其中R為威脅因子個數(shù)；wr為威脅因子的權(quán)重；yr為威脅層在地類圖層上的柵格個數(shù)；ry為威脅因子的強度；irxy為威脅因子對生境的威脅水平；βx為法律保護程度；Sjr為j類土地類型對威脅因子的敏感性大?。籨xy為柵格x(生境)與柵格y(威脅因子)的距離；drmax為威脅因子r的最大影響范圍。

參考前人研究經(jīng)驗[4，15，17，25]并結(jié)合研究區(qū)的區(qū)域特征，選取 5個威脅因子，分別是城鎮(zhèn)居民點、農(nóng)村居民點、工礦用地、道路、耕地。并確定了各個威脅因子的最大影響距離、權(quán)重以及不同土地利用類型的生境適宜度、各土地利用類型對不同威脅因子的相對敏感程度。其參數(shù)設(shè)置詳見表2、表3。

表2 威脅因子、最大影響距離及權(quán)重Table 2 The maximum influence distance and weight of the threat factor

表3 土地利用類型生境適宜度及威脅因子的敏感程度Table 3 Habitat suitability of land use type and the sensitivity to threat factors

3 結(jié)果與分析

3.1 土地空間利用變化

從數(shù)量結(jié)構(gòu)上的變化來看，1990年—2015年間，生態(tài)空間、生產(chǎn)空間始終為林州市的主要國土空間類型。25年間林州市生活空間持續(xù)增加，增加了2.89%；生態(tài)空間發(fā)生小幅度的波動變化；生產(chǎn)空間從 1990—2010年間持續(xù)減少，減少了 2.87%，但從2010—2015年間，生產(chǎn)空間又有小幅度的上漲，增加了0.64%。林州市三生空間如表4所示。

從鄉(xiāng)鎮(zhèn)分布特征來看，可以得出林州市各鄉(xiāng)鎮(zhèn)空間內(nèi)部分布差異較大，陵陽鎮(zhèn)是除林州市中心城區(qū)之外城鎮(zhèn)生活空間占比最大的鄉(xiāng)鎮(zhèn)，此外城郊鄉(xiāng)、姚村鎮(zhèn)、河順鎮(zhèn)與其他鄉(xiāng)鎮(zhèn)相比城鎮(zhèn)生活空間占比分布面積較大。林州市的農(nóng)業(yè)生產(chǎn)空間主要集中在陵陽鎮(zhèn)、姚村鎮(zhèn)、橫水鎮(zhèn)、采桑鎮(zhèn)，相對應(yīng)的這些鄉(xiāng)鎮(zhèn)的綠被生態(tài)空間面積較少。綠被生態(tài)空間主要集中在石板巖鎮(zhèn)、任村鎮(zhèn)、原康鎮(zhèn)、石板巖鎮(zhèn)、五龍鎮(zhèn)，其中石板巖鎮(zhèn)的綠被面積最大，達到 95%，該鄉(xiāng)鎮(zhèn)的其它空間面積占比很小。陵陽鎮(zhèn)、城郊鄉(xiāng)、姚村鎮(zhèn)的農(nóng)村生活空間占比較大。水域生態(tài)空間主要分布在石板巖鎮(zhèn)、任村鎮(zhèn)。如圖1所示。

表4 林州市25年間土地利用分布變化情況Table 4 Distribution change of land use types in the 25 years in Linzhou

從各個空間類型的轉(zhuǎn)化情況來看，25年間林州市生產(chǎn)空間、生活空間的變化幅度較大，生產(chǎn)空間中農(nóng)業(yè)生產(chǎn)空間大幅度減少，變化面積為 91.11km2，主要與綠被生態(tài)空間相互轉(zhuǎn)化數(shù)量最大；工礦生產(chǎn)空間面積明顯增加，其增加的面積主要源于生態(tài)生產(chǎn)空間以及農(nóng)業(yè)生產(chǎn)空間；農(nóng)村生活空間、城鎮(zhèn)生活空間面積增加，其主要的來源為農(nóng)業(yè)生產(chǎn)空間；生態(tài)空間變化幅度不大，其中綠被生態(tài)空間略有增加，主要由農(nóng)業(yè)生產(chǎn)空間和水域生態(tài)空間轉(zhuǎn)入；水域生態(tài)空間略微減少，其主要轉(zhuǎn)化為綠被生態(tài)空間，如表5、表6、表7所示。

圖1 林州市各鄉(xiāng)鎮(zhèn)土地類分布差異Figure 1 Distribution difference of land use types Townships in Linzhou

表5 1990—2000年林州市土地利用變化轉(zhuǎn)移矩陣(km2)Table 5 Transition matrix of land use types in the Linzhou during 1990-2000 (km2 )

表6 2000—2010林州市國土空間面積轉(zhuǎn)移矩陣(km2)Table 6 Transition matrix of land use types in the Linzhou during 2000-2010 (km2 )

表7 2010—2015林州市國土空間面積轉(zhuǎn)移矩陣(km2)Table 7 Transition matrix of land use types in the Linzhou during 2010-2015(km2 )

3.2 景觀格局特征及其演變

通過Frgstats4.2軟件對選取多個景觀指數(shù)進行計算，得到林州市25年間類型水平及景觀水平上的景觀指數(shù)，如表8、圖2所示。

綠被生態(tài)空間的斑塊面積所占景觀面積比例及平均斑塊面積遠大于其他類型，其次是耕地，這說明它們主要以大面積的斑塊存在，分布比較集中，聚集程度較高，但是綠被生態(tài)空間、農(nóng)業(yè)生態(tài)空間斑塊數(shù)量增多，平均斑塊面積逐年下降，表明其景觀日益趨于破碎化，聚集程度下降。與之相反，城鎮(zhèn)生活空間不斷的向城市外圍擴張，其PLAND、MPS逐年增加，說明城鎮(zhèn)生活空間面積不斷的增加，聚集程度逐年提高。

表8 林州市25年類型水平上景觀指數(shù)變化Table 8 Landscape indexes changes at type level in the 25 years in Linzhou

林州市 25年間景觀蔓延度指數(shù)偏高，其不同斑塊類型的團聚程度較高，景觀連通性較強，1990-2015年間蔓延度指數(shù)(CONTAG)總體呈下降趨勢；景觀形狀指數(shù)(LSI)、斑塊密度呈上升趨勢；香農(nóng)多樣性指數(shù)及香農(nóng)均勻度指數(shù)呈略微上升趨勢。蔓延度指數(shù)的下降，以及景觀形狀指數(shù)、斑塊密度的增加，表明林州市景觀破碎程度加大，景觀連通性減弱；而香農(nóng)多樣性指數(shù)及香農(nóng)均勻度指數(shù)有略微波動，表明各類型在景觀中呈均衡化趨勢分布?？傮w而言，隨著林州市社會經(jīng)濟的快速發(fā)展和人口規(guī)模的擴大，景觀格局破碎化加劇，景觀連通性減弱和多樣性增強表明國土空間利用更加系統(tǒng)化及精細化。

圖2 林州市25年景觀水平上景觀指數(shù)變化Figure 2 Landscape indexes changes at Landscape level in the 25 years in Linzhou

3.3 生境退化度時空變化特征

通過 InVEST 模型生境質(zhì)量(Habitat Quality)模塊來計算生境退化度。如圖3所示。

1990—2015年間林州市土地退化程度逐年增加，并且表現(xiàn)出明顯的空間差異。其主要退變方向是圍繞著城鎮(zhèn)生活空間向四周擴散展開，生活空間周圍土地退化程度最高，依次蔓延到生產(chǎn)空間，其次是生產(chǎn)空間與生態(tài)空間的交界，其交界處土地退化度大幅度降低，最后蔓延到生態(tài)空間，生態(tài)空間周圍的土地退化程度最低。

利用2015年的生境退化度與1990年的生境退化度差值，得到林州市 25年間的生境退化變化情況。其中，“+”值代表2015年林州市生境退化程度高于1990年的生境退化程度，“-”值代表2015年的生境退化程度低于1990年。如圖4所示。

通過 1990年—2015年林州市生境退化度變化可以看出，林州市生境退化度主要是生活空間向四周退化。中心城區(qū)生境退化程度最高；其次為陵陽鎮(zhèn)、姚村鎮(zhèn)、城郊鄉(xiāng)；原康鎮(zhèn)與合澗鎮(zhèn)雖2015年生境退化度高于 1990年，但是其退化程度低于陵陽鎮(zhèn)、姚村鎮(zhèn)、城郊鄉(xiāng)，影響其生境退化的主要因素是農(nóng)業(yè)生產(chǎn)空間的擴大；任村鎮(zhèn)、桂林鎮(zhèn)部分區(qū)域生境退化度低于 1990年，生態(tài)質(zhì)量有變好的趨勢，其生態(tài)質(zhì)量好轉(zhuǎn)的原因主要由于綠被生態(tài)空間的增加；石板沿鎮(zhèn)生境退化程度趨于穩(wěn)定。

4 討論

國土空間格局的演變對景觀格局的改變直接影響生態(tài)環(huán)境的質(zhì)量，生態(tài)空間是城市系統(tǒng)的重要組成部分，對經(jīng)濟的可持續(xù)發(fā)展及居民生活質(zhì)量的提高具有重要的生態(tài)服務(wù)價值。從1990年到2015年林州市生態(tài)空間用地面積雖有所增加，但是其生境退化度卻逐年上升，生態(tài)環(huán)境惡化，這主要是由于生活空間的不斷擴張，各種威脅因子的相互影響，導致綠被生態(tài)空間平均斑塊的破碎度加大，破壞生態(tài)服務(wù)功能的連通性，這是引起生境退化的原因之一，為了避免林州市生態(tài)用地的進一步流失，生態(tài)環(huán)境的進一步惡化，需要對其進行合理的配置，不僅要注意其空間服務(wù)價值更要注重其空間結(jié)構(gòu)合理配置。在精準扶貧的大背景下，要把帶領(lǐng)人民走上脫貧致富的道路與生態(tài)環(huán)境保護相結(jié)合。

圖3 1990—2015年林州市生境退化度Figure 3 Habitat degradation in Linzhou from 1990-2015

圖4 1990—2015年林州市生境退化度變化Figure 4 The change of habitat degradation in Linzhou from 1990-2015

林州市不同區(qū)域的生境退化度存在著差異，生境退化度上升最明顯的區(qū)域分別為中心城區(qū)、姚村鎮(zhèn)、陵陽鎮(zhèn)、城郊鄉(xiāng)，主要是因為這些區(qū)域分布著大量的人工表面，并且近 25年來，林州市經(jīng)濟快速發(fā)展，建設(shè)用地大量擴張，耕地面積持續(xù)減少，導致了生態(tài)用地景觀破碎加劇，生境退化度高，使得生態(tài)環(huán)境下降明顯；而任村鎮(zhèn)及桂林鎮(zhèn)，人工表面較少，人為因素干擾不大，且綠被生態(tài)空間面積較大，生態(tài)功能維持較好。由于城市不同的地區(qū)，其生態(tài)質(zhì)量不同，在今后的發(fā)展建設(shè)中，需因地制宜，結(jié)合不同區(qū)域的生態(tài)特征、人文特征，對城市進行規(guī)劃，可以更好的體現(xiàn)城市功能用地的服務(wù)價值。林州市主城區(qū)既要保護其本身的生態(tài)空間、還要規(guī)劃增加小范圍的綠地等，來提升生態(tài)服務(wù)功能，提高居民生活質(zhì)量。而生態(tài)環(huán)境較好的任村鎮(zhèn)及桂林鎮(zhèn)，可以很好的利用其自然資源，適度開發(fā)旅游項目。但是林州市屬于多山地區(qū)，山地丘陵面積占總面積的 86%，由于山區(qū)具有特殊的屬性特質(zhì)，其自然資源呈梯度分布、具有潛在地質(zhì)災(zāi)害，過度開發(fā)容易造成地質(zhì)災(zāi)害頻發(fā)，生態(tài)惡劣。所以要嚴格把控林州市山地丘陵區(qū)的開發(fā)程度，要避免低效率的生產(chǎn)、生活空間的擴張，多方面考慮，總體調(diào)控，做到國土空間的合理開發(fā)及利用。文章主要通過分析林州市由國土空間格局變化引起的景觀格局改變對生態(tài)環(huán)境的影響，但是怎么在城鎮(zhèn)化及精準扶貧開發(fā)的背景下避免城鎮(zhèn)景觀趨于破碎化，怎么建立起生境退化與景觀破碎指數(shù)之間的聯(lián)系，這個問題值得繼續(xù)研究和探討。

5 結(jié)論

本文通過對林州市25年的“三生”空間格局變化及生境退化度分析，得出如下結(jié)論:

(1) 林州市的國土空間主要分為生態(tài)空間、生產(chǎn)空間、生活空間。其中生態(tài)空間與生產(chǎn)空間是其主要的空間類型。1990—2015年期間林州市生活空間逐年增加，城鎮(zhèn)化程度逐年加大，中心城區(qū)的城鎮(zhèn)快速向四周擴張，農(nóng)業(yè)生產(chǎn)空間大幅度減少，其主要轉(zhuǎn)化成綠被生態(tài)空間及生活空間。生態(tài)空間基本保持原有水平，并未發(fā)生較大波動。經(jīng)濟發(fā)展及人類活動對林州市國土空間格局的變化起到關(guān)鍵性作用。

(2) 1990—2015年間林州市蔓延度指數(shù)下降，景觀形狀指數(shù)、斑塊密度增加，香農(nóng)多樣性指數(shù)及香農(nóng)均勻度指數(shù)有略微波動。其中城鎮(zhèn)生活空間不斷的向外圍擴張，聚集程度逐年提高，綠被生態(tài)空間、農(nóng)業(yè)生態(tài)空間的景觀破碎化程度不斷的加大?？偟膩砜矗捎谌祟惢顒蛹敖?jīng)濟的發(fā)展，林州市景觀格局破碎化加劇，景觀連通性減弱和景觀多樣性增強。

(3) 林州市土地退化程度逐年增加，并且表現(xiàn)出明顯的空間差異。其中西北部山區(qū)及東南部地區(qū)，生境退化程度下降，生態(tài)環(huán)境好轉(zhuǎn)；中心城區(qū)的生境退化程度明顯上升，其生境退化程度最高，其次是與其接壤的三個鄉(xiāng)鎮(zhèn)，分別為陵陽鎮(zhèn)、姚村鎮(zhèn)、城郊鄉(xiāng)。中心城區(qū)、陵陽鎮(zhèn)、姚村鎮(zhèn)、城郊鄉(xiāng)是林州市規(guī)劃發(fā)展的敏感區(qū)域，在未來的城鄉(xiāng)規(guī)劃中，要重點關(guān)注這些鄉(xiāng)鎮(zhèn)，發(fā)展規(guī)劃建設(shè)要與生態(tài)保護相結(jié)合。