王 軍,嚴有龍,王金滿,應凌霄,唐 倩

1 自然資源部國土整治中心,自然資源部土地整治重點實驗室, 北京 100035 2 中國地質(zhì)大學(北京)土地科學技術(shù)學院, 北京 100083

流域作為自然地理和經(jīng)濟發(fā)展的復合型區(qū)域,兼有重要的資源、生態(tài)、經(jīng)濟和文化功能[1]。流域生態(tài)系統(tǒng)服務囊括了人類從流域系統(tǒng)獲得的供給、調(diào)節(jié)、文化與支持服務等所有惠益[2]。流域生境是生物用于生長、發(fā)育、繁殖與分布的資源與條件,是區(qū)域生物多樣性最重要的體現(xiàn)[3- 4],生境質(zhì)量是基于生存資源可獲得性,生態(tài)系統(tǒng)提供適合于生物生長、發(fā)育、繁殖與分布條件的能力[5-7]。進入工業(yè)社會以來,人類活動導致了生境破碎、退化甚至消失[7-8],而土地利用變化程度體現(xiàn)了人類活動的強度,成為生境質(zhì)量最重要的威脅因子[9]。因此,基于土地利用變化研究區(qū)域生境質(zhì)量格局,對區(qū)域可持續(xù)發(fā)展具有重要意義。

由于人們對土地利用變化的持續(xù)關(guān)注,土地利用模擬模型研究隨之興起。元胞自動機(Cellular Automata, CA)作為諸多模型中最主要且應用較廣泛的方法之一,是諸多模型研究的基礎(chǔ)[10],譬如Logistic-CA模型[11]、ANN-CA(Artificial Neural Network-Cellular Automata)模型[12]、CLUE-S(the Conversion of Land Use and its Effects at Small regional extent)模型[13]與FLUS(Future Land Use Simulation)模型[14]。而CA-Markov模型是其衍生的較為成功的模擬方法,綜合了CA模型模擬復雜系統(tǒng)空間變化的能力和Markov模型的長期預測優(yōu)勢,能有效預測土地利用類型相互轉(zhuǎn)化概率并模擬其空間轉(zhuǎn)化模式[15-17]。土地利用模擬也逐漸從單一用地擴展為多地類模擬[16-17]。例如Han等采用CA-Markov模型較好的模擬了佛山地區(qū)不同發(fā)展情景下土地利用空間格局[16]。目前國內(nèi)外有關(guān)流域生境質(zhì)量研究主要集中在兩方面：對單點、樣地尺度生境質(zhì)量評價[18-19]和區(qū)域整體生境質(zhì)量評價[20-21]。研究方法較多采用水文水力法、河流地貌法、生境預測模擬法和綜合評估法[18-27]。例如Kail等使用土壤和水文評估工具(Soil and Water Assessment Tool, SWAT),結(jié)合河流水力學與2D水動力學模型探索生境質(zhì)量空間分布信息[26-27]；陳淼等對三峽庫區(qū)36條重要支流的多個樣點進行了綜合評價等[28]。諸多學者將生態(tài)系統(tǒng)服務評估與權(quán)衡模型(Integrated Valuation of Ecosystem Services and Tradeoffs, InVEST)應用于區(qū)域生境質(zhì)量評估,取得了明顯成效[23-24,28- 29]。隨著土地利用模型技術(shù)的成熟,對歷史時期與未來生境質(zhì)量的探索逐漸成為研究熱點[7,30]。但目前利用這些技術(shù)預測不同情景下生境質(zhì)量變化特征的研究尚不多見。

福建省是我國東南沿海的天然屏障,閩江是省內(nèi)第一大河,生態(tài)功能十分重要,為多種動植物提供繁殖、生活的棲息生境,也是我國流域生態(tài)補償重要試點之一[31]。本文選用InVEST模型,并結(jié)合CA-Markov模型模擬未來不同情景下土地利用狀況,分析2000—2040年閩江流域生境質(zhì)量的時空演變規(guī)律,為閩江流域土地資源可持續(xù)利用與生物多樣性保護提供科學依據(jù)和決策參考。

1 數(shù)據(jù)源與研究方法

1.1 研究區(qū)概況

閩江發(fā)源于福建省建寧縣均口鎮(zhèn),總長2872 km,干流全長577 km,年徑流量達621億m3,流域面積6萬多km2,90%以上位于福建省內(nèi),占全省陸域面積48.87%,地理位置約為116°—120°E,25°—29°N。考慮到研究區(qū)域完整性及特殊性,結(jié)合閩江流域自然邊界最終確定上游(浦城縣、武夷山市、光澤縣、松溪縣、建陽市、邵武市、順昌縣、建甌市、政和縣、南平市、連城縣),中游(寧化縣、清流縣、永安市、大田縣、建寧縣、明溪縣、三明市、泰寧縣、將樂縣、沙縣、尤溪縣)與下游(德化縣、古田縣、閩清縣、閩侯縣、永泰縣、福州市、長樂市)為研究范圍(圖1)。流域內(nèi)地形復雜,盆地、山地丘陵與峽谷相互交錯,主要土壤類型為紅壤與黃壤；該區(qū)為亞熱帶海洋性季風氣候,降雨量及產(chǎn)水量豐富,植被覆蓋率高,主要植被類型為常綠闊葉林、常綠針葉林與竹林,區(qū)內(nèi)武夷山具有我國東南現(xiàn)存面積最大、保留最完整的中亞熱帶常綠闊葉林森林生態(tài)系統(tǒng)。流域內(nèi)動植物資源豐富；閩江口濕地和三明市沙溪流域明溪段等地還是全球性的候鳥遷徙通道[31]。

圖1 閩江流域位置圖Fig.1 Location map of Min River Basin

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

研究數(shù)據(jù)源包括InVEST模型“Habitat Quality”模塊所需數(shù)據(jù)與IDRISI模型數(shù)據(jù),本文中CA模型的元胞大小為1 km2,統(tǒng)一調(diào)整圖像柵格為1 km×1 km。土地利用與土地覆被數(shù)據(jù)為2000年、2005年、2010年與2015年地理空間數(shù)據(jù)云平臺下載Landsat系列遙感圖通過鑲嵌、校正、裁剪、解譯等處理獲得,解譯精度達95%以上,劃分為13個地類(圖2)。在地理空間數(shù)據(jù)云下載研究區(qū)30 m分辨率數(shù)字高程數(shù)據(jù),通過鑲嵌、裁剪與填洼處理獲得高程、坡度與坡向等數(shù)據(jù)。距水域、道路與居民點距離利用GIS提取的研究區(qū)水域、道路與居民點進行歐氏距離分析獲取。土壤質(zhì)地來源于世界土壤數(shù)據(jù)庫(HWSD)。人口密度與國內(nèi)生產(chǎn)總值來源于資源環(huán)境科學與數(shù)據(jù)中心已發(fā)布的公里網(wǎng)格GDP和人口密度數(shù)據(jù)。夜間的燈光指數(shù)來源于資源環(huán)境科學與數(shù)據(jù)中心與NOAA網(wǎng)站下載獲得1 km柵格數(shù)據(jù)。半飽和參數(shù)K參照模型參數(shù)與學者研究獲取[32-33]。

圖2 閩江流域2015年土地利用現(xiàn)狀數(shù)據(jù)Fig.2 Land use data of Min River Basin in 2015

1.3 研究方法

1.3.1土地利用模擬預測

CA模型被廣泛應用于城市空間擴展和土地利用模擬,模擬結(jié)果為合理利用土地和區(qū)域經(jīng)濟發(fā)展規(guī)劃提供理論依據(jù)；Markov模型是將馬爾科夫鏈理論模型和方法研究應用于概率論中,以分析隨機事件的變化并預測未來趨勢[15-17]。

(1)Logistic模型

Logistic模型屬于非線性模型,主要應用于對多因素綜合影響下的事件概率進行預測,近年來廣泛運用于土地利用變化的因子分析中[15-17]。每個柵格可能出現(xiàn)某種土地利用類型的概率可以通過將因變量與自變量的數(shù)據(jù)采取Logistic回歸分析得出。ROC(Relative Operating Characteristic)是驗證Logistic回歸方程擬合度的一個精度指標,一個完整的隨意模型確定的ROC值為0.5,而滿意適合結(jié)果的ROC值為1。

(2)CA-Markov模型

CA是一種基于不連續(xù)的時空動態(tài)模擬模型。其特點是時間、空間和狀態(tài)都是離散狀態(tài),復雜系統(tǒng)可以由一些很簡單的局部規(guī)則來產(chǎn)生[15-17]。通常包括單元、狀態(tài)、鄰近范圍和轉(zhuǎn)換規(guī)則4個要素,可用如下表達式表示：

S(t+1)=f(St,N)

(1)

式中,S為元胞離散、有限的狀態(tài)集合,N為元胞的鄰域,t、t+1分別表示2個不同的時刻,f為元胞狀態(tài)轉(zhuǎn)化規(guī)則。

Markov模型可實現(xiàn)系統(tǒng)由一種狀態(tài)轉(zhuǎn)移至另一種狀態(tài),預測結(jié)果由初始狀態(tài)向量和轉(zhuǎn)移概率矩陣決定,其特點為無后效性和穩(wěn)定性,土地利用變化滿足該模型使用條件[5,26,30]。運用Markov模型進行土地利用預測的關(guān)鍵在于確定在n時刻的轉(zhuǎn)移概率Pij。其表達式為：

(2)

利用此概率建立土地利用Markov預測模型為：

S(k+1)=S(k)Pij=S(0)Pij(k+1)

(3)

式中,S(k+1)為預測地類在t=k+1時的狀態(tài)向量,即預測結(jié)果；S(0)為初始狀態(tài)量,即預測初期各土地利用類型的面積；S(k)為預測地類在t=k時刻的狀態(tài)向量；Pij為研究時段內(nèi)土地利用類型i轉(zhuǎn)換為j的概率。

(3)模擬精度檢驗

模擬精度檢驗包括數(shù)量與空間精度檢驗[15-17]。數(shù)量精度檢驗公式定義如下：

(4)

式中,Ei為第i類土地利用類型的數(shù)量精度；miy為第i類土地利用類型的模擬面積；mix為第i類土地利用類型的實際面積。Ei的值越小說明模擬精度越高。

IDRISI模型中Kappa系數(shù)是一個空間精度檢驗系數(shù),在評價空間模擬精度時,計算得出的Kappa系數(shù)>0.80時,精度通過檢驗；當0.4

1.3.2生境質(zhì)量測算

InVEST模型由美國斯坦福大學聯(lián)合世界自然基金會等開發(fā),初衷是為了權(quán)衡區(qū)域發(fā)展與保護間的關(guān)系,尋求最優(yōu)自然資源管理和經(jīng)濟發(fā)展模式。模型中“Habitat Quality”可作為生境質(zhì)量的反映,結(jié)合土地利用和土地覆被與生境質(zhì)量威脅因素的信息生成生境質(zhì)量地圖,將不同的土地利用和土地覆被類型視為相應的生態(tài)系統(tǒng)類型或人類活動的干擾因子,依據(jù)各生態(tài)系統(tǒng)類型對動植物的生境適宜度和人類干擾因子的威脅強度來模擬生境質(zhì)量的空間分布[3,29,32,34]。生境質(zhì)量計算公式如下：

(5)

(6)

式中,Qij指土地利用類型j中柵格i的生境質(zhì)量,處于0—1區(qū)間內(nèi),值越高生境質(zhì)量越好；Dij是柵格i的生境脅迫水平；Hj為土地利用類型j的生境適合度；k是半飽和常數(shù)；Z為模型默認參數(shù)的歸一化常量；R表示脅迫因子；y表示脅迫因子r柵格圖層的柵格數(shù)；r表示圖層的柵格數(shù)；Yi表示脅迫因子所占柵格數(shù)；Wi表示脅迫因子的權(quán)重,介于0—1之間；ry表示柵格y的脅迫因子值(0或1)；iriy表示柵格y的脅迫因子值ry對生境柵格i的脅迫水平；Sjr表示生境類型j對脅迫因子r的敏感度。

閩江流域作為福建省工農(nóng)業(yè)的重要生產(chǎn)基地,人類活動與工農(nóng)業(yè)生產(chǎn)等為區(qū)內(nèi)生境的主要脅迫因子。研究選取受到人類干擾較大的建設(shè)用地、交通用地、采礦場、裸地、水田與旱地作為生態(tài)環(huán)境的威脅源頭。模型通過空間距離這一變量計算生態(tài)威脅因子對各個土地利用類型的影響程度,鑒于空間關(guān)系的復雜性,提供了線性衰退與指數(shù)衰退兩種威脅因子對土地利用類型影響距離的方法。參照模型參數(shù)與相關(guān)學者研究進行賦值(表1與表2)[3,7,20,29,32-34]。

表1 閩江流域脅迫因子的最大影響距離及其權(quán)重

表2 閩江流域不同生境類型對不同脅迫因子的敏感度

1.3.3土地利用預測情景設(shè)置

基于2000—2015年的土地利用變化數(shù)據(jù),設(shè)立2025—2040年現(xiàn)有發(fā)展與生態(tài)保護兩種情景?，F(xiàn)有發(fā)展情景下根據(jù)研究區(qū)土地利用特點選擇高程、坡度、坡向、到水域距離、土壤質(zhì)地、人口密度、夜間燈光指數(shù)、國內(nèi)生產(chǎn)總值、到公路距離、到居民點距離10個驅(qū)動因子為自變量,因變量為各地類分布圖,通過Logistic回歸分析制作土地適宜性圖集。

生態(tài)保護情景下遵從全國國土規(guī)劃綱要設(shè)立的目標：堅持集聚開發(fā)與均衡發(fā)展相協(xié)調(diào)、提高土地利用效率、加強重點生態(tài)功能區(qū)保護、強化耕地資源保護等,嚴格限制未利用地增長,控制建設(shè)用地增長,保護耕地、林地、水域等,基本限制其他地類轉(zhuǎn)為未利用地,水田、水域和濕地限制其轉(zhuǎn)為建設(shè)用地,利用MCE模型加入限制因子更新土地利用適宜性圖集,模擬土地利用變化。

2 結(jié)果與分析

2.1 閩江流域土地利用變化

2.1.1現(xiàn)狀土地利用變化分析

由圖2與圖3可知,2000—2015年間研究區(qū)土地利用有如下規(guī)律：耕地面積呈下降趨勢,降速有所擴大,其中耕地利用以水田為主,沿山腳平原與河漫灘地勢平坦區(qū)域零散分布；林地是主要用地類型,占流域總面積的87%左右,面積保持在相對穩(wěn)定的狀態(tài),以針葉林與闊葉林為主,面積有所增長,生態(tài)用地保護的較好；針闊混交林面積占比相對較小且不斷縮減,零散分布在上游和下游區(qū)域,表明研究區(qū)森林類型相對單一；灌木林的占比較大且面積較穩(wěn)定,林質(zhì)林相有很大提升空間；草地占比較小且面積呈不斷縮減趨勢,主要集中在上游地勢相對平坦的區(qū)域與現(xiàn)有城市綠化區(qū)域,現(xiàn)狀條件下減幅最大,達44.64%；流域用地主要為水域,同時有相對規(guī)模的濕地分布于水域附近,面積有所擴展,其中以水域增長為主,生態(tài)建設(shè)至關(guān)重要,構(gòu)成了閩江流域的生命線；工礦、交通與建設(shè)用地占比相對較小,但面積在不斷擴展,建設(shè)用地在流域下游福州、長樂等區(qū)域分布集中,交通用地是連接各大城市的樞紐,采礦場零散分布于流域上中游區(qū)域,現(xiàn)狀條件下建設(shè)用地增長面積最大,但交通用地漲幅最大,達到227.27%；裸地規(guī)模相對較大,現(xiàn)狀條件下有緩慢增長趨勢。

圖3 閩江流域2000—2040年土地利用面積統(tǒng)計圖Fig.3 Land use area statistics in Min River Basin from 2000 to 2040

2.1.2閩江流域土地利用變化模擬分析

(一)模擬精度檢驗

研究利用Logistic模型制作土地利用適宜性圖集,模擬單項土地利用變化,ROC檢驗結(jié)果均大于等于0.83,表明所選驅(qū)動因子具有較好的模擬精度。為驗證CA-Markov模型模擬精度,首先基于2005、2010年兩期土地利用數(shù)據(jù)和2010年土地利用驅(qū)動因子預測2015年土地利用圖,再與實測數(shù)據(jù)進行對比分析。面積檢驗結(jié)果見表3,除草地與采礦場外其余各地類模擬的數(shù)量精度均小于15%,具有較好的模擬精度。

表3 2015年土地利用模擬面積精度檢驗表

(二)土地利用變化模擬分析

據(jù)圖3可知,在現(xiàn)有發(fā)展情景下,耕地面積總體表現(xiàn)出波動下降的趨勢；林地總面積基本能保持穩(wěn)定,針葉林面積減至相對穩(wěn)定的數(shù)值,針闊混交林與灌木林面積在波動狀態(tài)下緩慢上升；草地面積則下降到相對穩(wěn)定的水平；水域總面積基本能保持穩(wěn)定,存在部分水域向濕地及建設(shè)用地轉(zhuǎn)換的現(xiàn)象；受經(jīng)濟建設(shè)影響,工礦、交通與建設(shè)用地面積都表現(xiàn)出一定程度增加,其中建設(shè)用地增速最高；裸地面積表現(xiàn)為先增加后減少的趨勢,主要分布于流域上中游區(qū)域。

據(jù)圖3與圖4可知,在生態(tài)保護情景下,實現(xiàn)集聚開發(fā)與均衡發(fā)展相協(xié)調(diào),耕地面積總體表現(xiàn)出波動下降的趨勢,較現(xiàn)有發(fā)展情景降幅??；林地總面積基本能保持穩(wěn)定,其中針葉林面積減少至相對穩(wěn)定的數(shù)值,闊葉林、針闊混交林與灌木林面積在波動狀況下緩慢上升,較現(xiàn)有發(fā)展情景下林地受到較好保護；草地面積則下降到相對穩(wěn)定的水平,在未來保護實施中草地應作為重點關(guān)注對象；水域總面積基本能保持穩(wěn)定,存在部分水域向濕地及建設(shè)用地轉(zhuǎn)換的現(xiàn)象,向建設(shè)用地轉(zhuǎn)換的幅度較現(xiàn)有發(fā)展情景要??；工礦、交通與建設(shè)用地面積都表現(xiàn)出一定程度的增加,其中建設(shè)用地增速最高,生態(tài)保護情景并沒有限制建設(shè)用地增長；裸地面積逐漸減小,分布也較為規(guī)律,生態(tài)保護情景能有效控制裸地增加。

圖4 閩江流域土地利用預測圖Fig.4 Land use forecast map of Min River Basin

2.2 閩江流域生境質(zhì)量時空動態(tài)分析

根據(jù)圖5可知,工礦建設(shè)用地區(qū)域為研究區(qū)生境質(zhì)量低值的主要分布區(qū),下游的福州與長樂市分布有最大的生境質(zhì)量低值聚集區(qū),同時生境質(zhì)量分布與地形條件存在極大相關(guān)性,低值區(qū)多為海拔較低的區(qū)域,沿流域附近的河谷盆地零散分布,高海拔地區(qū)多為生境質(zhì)量良好的區(qū)域；土地利用規(guī)劃及有效的建設(shè)開發(fā)與保護行為導致未來生境破碎度減小,生境質(zhì)量的聚集效應更為明顯,而且生態(tài)保護情景下的生境質(zhì)量分布最為規(guī)律,也體現(xiàn)出其較好保護了生境質(zhì)量。

圖5 2000—2040年閩江流域生境質(zhì)量時空分布與變化Fig.5 Temporal-spatial distribution and change of habitat quality in Min River Basin from 2000 to 2040

據(jù)統(tǒng)計結(jié)果圖6可知,2000、2005、2010和2015年閩江流域生境質(zhì)量的均值分別為0.83、0.83、0.82、0.82,生境狀況良好,2005—2015年3期的生境質(zhì)量降幅分別為0.13%、0.23%和0.24%,其中2010—2015階段下降最快,與生境質(zhì)量現(xiàn)實結(jié)果變化規(guī)律剛好相反,生境質(zhì)量總體呈現(xiàn)衰退趨勢,且衰退變化率逐期升高。2015—2025年在現(xiàn)有發(fā)展情景下降幅達0.81%,生態(tài)保護情景下降幅僅為0.38%；同時,現(xiàn)有發(fā)展情景下,2025—2040年閩江流域生境質(zhì)量的變化率為0.01%,生境質(zhì)量下降的趨勢未得到有效緩解；生態(tài)保護情景下,2025—2040年閩江流域生境質(zhì)量的變化率為0.03%,生境質(zhì)量下降的趨勢得到有效緩解,呈V字型波動增長。參照已有研究的劃分標準[5],按照生境質(zhì)量每0.2分劃分一個級別,劃分為差、較差、一般、良好、優(yōu)等五個級別,2000、2005、2010與2015年閩江流域生境質(zhì)量為優(yōu)等等級的柵格個數(shù)占比分別為86.87%、86.89%、86.69%、87.06%,呈現(xiàn)出波動上升的態(tài)勢；現(xiàn)有發(fā)展情景下,2025、2030、2035與2040年閩江流域生境質(zhì)量為優(yōu)等等級的柵格個數(shù)占比分別為86.93%、86.98%、87.06%、87.00%,呈現(xiàn)倒V型波動態(tài)勢,較2015年略有降低；生態(tài)保護情景下,2025、2030、2035與2040年閩江流域生境質(zhì)量為優(yōu)等等級的柵格個數(shù)占比分別為87.36%、87.39%、87.41%、87.42%,呈現(xiàn)緩慢上升的態(tài)勢。

圖6 2000—2040年閩江流域生境質(zhì)量變化 Fig.6 Habitat quality changes in the Min River Basin from 2000 to 2040

通過與研究區(qū)行政區(qū)劃圖疊加(圖7)分析可知,2000—2015年上中游區(qū)域生境水平整體較高,均高于0.8的水平；下游區(qū)域整體水平較低,尤其是靠近閩江口的福州與長樂市,但生境質(zhì)量水平與區(qū)位分布并不是簡單的線性關(guān)系,比如下游也存在永泰縣生境質(zhì)量達0.84的水平,該區(qū)域分布有藤山自然保護區(qū)?，F(xiàn)有發(fā)展情景下,發(fā)現(xiàn)蒲城、光澤、松溪、建寧、福州與長樂等生境質(zhì)量下降最為明顯,區(qū)域生境質(zhì)量差距呈現(xiàn)擴張趨勢；生態(tài)保護情景下基本維持2015年的生境質(zhì)量水平,值得注意的是三明市生境質(zhì)量明顯降低,其原因可能是三明市建設(shè)用地在未來的擴張潛力較大；研究區(qū)整體水平較現(xiàn)有發(fā)展情景略有提高。

圖7 閩江流域2000—2040年各縣域生境質(zhì)量時空變化Fig.7 Temporal-spatial change in habitat quality of counties in Min River Basin from 2000 to 2040

3 討論

研究區(qū)土地利用格局受到地形、區(qū)位、氣候、社會經(jīng)濟、人口與政策等多重因素影響。根據(jù)土地利用現(xiàn)狀分析結(jié)果可知,林地在閩江流域占據(jù)了絕對的優(yōu)勢,以針葉林及闊葉林為主,南方丘陵山地帶土壤肥沃,降雨充沛,適宜林木生長,生態(tài)用地保存較為完整,但林地結(jié)構(gòu)有待提升,對殘次林的改造是下一步保護的重點；自建國以來福建省在擴大耕地面積方面做了很多努力,但耕地為補充建設(shè)用地潛力的主要來源,新增耕地面積依然不及被占用的耕地面積,耕地面積減小的態(tài)勢明顯,此外區(qū)內(nèi)低海拔的耕地以水田為主；濕地與水域為研究區(qū)提供了良好的生境,維持在較高的水平,是流域生態(tài)補償重點關(guān)注對象,構(gòu)成了閩江流域的生命線；草地也是高生境質(zhì)量分布的區(qū)域,但規(guī)模不斷減小,原因在于其轉(zhuǎn)變的彈性較大；建設(shè)用地與裸地是低生境質(zhì)量分布的主要區(qū)域,均存在一定程度擴展,其中以建設(shè)用地為主,這與閩江流域依托東南沿海經(jīng)濟帶優(yōu)勢大力發(fā)展經(jīng)濟及人口擴展密切相關(guān),這些結(jié)論與白健、鐘莉娜等研究相似[10,20, 35-37]?，F(xiàn)有發(fā)展情景下耕地與林地等由于未受到嚴格保護,縮減的態(tài)勢明顯,同時,隨著經(jīng)濟社會發(fā)展與人口增長,建設(shè)用地與裸地等增長明顯,基于規(guī)劃設(shè)立的生態(tài)保護情景,通過設(shè)置限制條件有效的緩解了該類問題,保護了耕地、林地與水域等,建設(shè)用地向著更加有序的方向擴展,裸地得到有效遏制,土地利用效率提高,滿足社會發(fā)展與生態(tài)保護的需求,但設(shè)計的保護策略也存在對草地關(guān)注度不足的問題,草地占比依然保持在較低水平。

研究區(qū)生境質(zhì)量與土地利用格局密切相關(guān)。生境質(zhì)量總體良好得益于研究區(qū)高植被覆蓋度、水域、濕地及相應的生態(tài)保護措施,維持生境是流域生態(tài)補償政策設(shè)計的重要目的[38-39]。流域上中游是經(jīng)濟欠發(fā)達的南平、三明等地區(qū),山地面積占比大,是重要的生態(tài)功能區(qū),而流域下游河流入?？诘貐^(qū)地勢相對平坦,作為海峽西岸的重要口岸,是經(jīng)濟較為發(fā)達的福州、長樂市,生態(tài)環(huán)境相對較差,同時低海拔區(qū)域人為擾動更為劇烈,導致空間上表現(xiàn)出上中游區(qū)域生境質(zhì)量較下游區(qū)域高,呈現(xiàn)出高海拔地區(qū)總體生境質(zhì)量比低海拔地區(qū)高的變化規(guī)律[1]。經(jīng)濟發(fā)展與城市擴展導致區(qū)域土地利用類型轉(zhuǎn)變,區(qū)域自然生境被進一步分割,低生境區(qū)域擴展上升,整體來看閩江流域生境質(zhì)量有降低的趨勢,這種趨勢在經(jīng)濟發(fā)達區(qū)域尤為明顯,原因是這些區(qū)域社會經(jīng)濟與人為活動更為激烈,需要采取必要措施,防止生境的進一步惡化。現(xiàn)有發(fā)展情景下由于保護力度較弱致使生境質(zhì)量進一步降低,生態(tài)保護情景下強化了相應保護措施,生態(tài)用地得到有效保護,建設(shè)用地有序擴張,生境質(zhì)量得到良好維持。

本研究綜合了CA-Markov模型與InVEST模型分析閩江流域過去與未來長時間序列的生境質(zhì)量格局,為閩江流域生境質(zhì)量的維護與優(yōu)化提供科學依據(jù)和決策參考,探索了回溯歷史時期和預測未來長時間序列的生境質(zhì)量研究思路。研究佐證了生態(tài)保護的重要性,區(qū)域上中下游存在明顯的生態(tài)質(zhì)量差異,在未來發(fā)展中應繼續(xù)保持生態(tài)補償策略,提升林質(zhì)林相,嚴格保護流域范圍內(nèi)的耕地、草地與水域等,特別是關(guān)注有限的草地,在下游區(qū)域加強對建設(shè)用地與未利用地的管控,促進流域上中下游共建共治共享流域生態(tài)環(huán)境?；贑A-Markov模型模擬的土地利用精度經(jīng)檢驗其Kappa系數(shù)為89.09%,精度較高,滿足了較大尺度的研究需要,小尺度上則需要更高的精度,同時對柵格尺度的探索也值得深化研究。

4 結(jié)論

作為海峽西岸的重要生態(tài)屏障,閩江流域森林覆蓋率較高,生態(tài)用地保存較為完整,但伴隨著經(jīng)濟發(fā)展與人口增長,建設(shè)用地占用耕地、草地與濕地等現(xiàn)象依然存在,現(xiàn)狀條件下水田、旱地、針闊混交林、灌木林、草地與濕地面積呈下降的趨勢,其中草地減幅最大,達44.64%,針葉林、闊葉林、水域、建設(shè)用地、交通用地、采礦場與裸地面積呈上升的趨勢,其中交通用地漲幅最大,達到227.27%。本文構(gòu)建的CA-Markov模型模擬土地利用變化精度較高,設(shè)立現(xiàn)有發(fā)展與生態(tài)保護兩種情景對照分析,基于生態(tài)保護情景在未來土地利用格局優(yōu)化上表現(xiàn)良好,可為研究區(qū)未來發(fā)展提供決策參考。

本研究利用InVEST模型中的“Habitat Quality”模塊模擬閩江流域生境質(zhì)量空間格局,探究研究區(qū)從2000—2040年由柵格到整體再到縣域的生境質(zhì)量空間格局。閩江流域總體呈現(xiàn)較高的生境質(zhì)量,基本維持在0.82的水平,工礦、交通與建設(shè)用地區(qū)域為研究區(qū)生境質(zhì)量低值的主要分布區(qū),下游的福州與長樂分布有最大的生境質(zhì)量低值聚集區(qū)。2000—2015年閩江流域生境質(zhì)量總體呈現(xiàn)衰退趨勢,且降幅逐期升高,最高達0.24%；現(xiàn)有發(fā)展情景下生境質(zhì)量下降的趨勢未得到有效緩解,生態(tài)保護情景下生境質(zhì)量下降的趨勢得到有效緩解,同時生境質(zhì)量為優(yōu)等級的柵格個數(shù)占比最大,最高達87.42%。總體看上中游區(qū)域生境水平較下游高,2000、2005、2010與2015年上中游區(qū)域平均生境質(zhì)量均高于0.8的水平,下游區(qū)域整體水平較低,尤其是靠近閩江口的福州與長樂市,但下游也存在永泰縣生境質(zhì)量達0.84的水平?，F(xiàn)有發(fā)展情景下區(qū)域生境質(zhì)量差距出現(xiàn)擴張趨勢,生態(tài)保護情景下各區(qū)域則基本維持了2015年的生境質(zhì)量水平,但也存在例如三明市生境質(zhì)量明顯降低的情況出現(xiàn),整體水平較現(xiàn)有發(fā)展情景有明顯提高。