蘇紅帆, 侯西勇, 邸向紅

?

北部灣沿海土地利用變化時空特征及情景分析

蘇紅帆1, 侯西勇1, 邸向紅2

(1. 中國科學(xué)院煙臺海岸帶研究所, 山東煙臺 264003; 2. 德州學(xué)院資源環(huán)境與規(guī)劃學(xué)院, 山東德州253023)

基于遙感技術(shù)監(jiān)測北部灣沿海2000年~2010年土地利用變化特征, 應(yīng)用多要素數(shù)據(jù)和Logistic- CA-Markov耦合模型模擬土地利用變化并進行多情景分析。結(jié)果表明: 2000年~2010年, 耕地是土地轉(zhuǎn)化的主要來源, 城鎮(zhèn)、農(nóng)村居民點及交通工礦的擴張主要來自于耕地, 城鎮(zhèn)化加劇增加了耕地保護的壓力; 2005年~2010年, 沿海圍填海的速度較快, 導(dǎo)致海岸線形態(tài)發(fā)生顯著的變化以及海岸帶脆弱性的加劇; 基于區(qū)域環(huán)境與生態(tài)保護、經(jīng)濟社會發(fā)展的歷史過程特征和相關(guān)規(guī)劃, 設(shè)計既有趨勢延續(xù)、城鎮(zhèn)化優(yōu)先、生態(tài)保護優(yōu)先3種情景, 模擬得到2020年、2025年的土地利用分布數(shù)據(jù); 綜合土地利用變化的歷史特征和情景分析的結(jié)果, 提出未來土地利用優(yōu)化的策略與措施, 包括: 優(yōu)化城鎮(zhèn)化空間布局及其發(fā)展質(zhì)量, 嚴格保護沿海濕地尤其是紅樹林濕地, 變化熱點區(qū)域合理規(guī)劃未來的土地利用方向等。

土地利用變化; 北部灣經(jīng)濟帶; Logistic回歸; CA-Markov模型; 情景分析

沿海區(qū)域是陸海交接地帶, 生態(tài)脆弱, 具有較為獨特的土地利用結(jié)構(gòu)、格局與過程。土地利用變化與沿海的資源、環(huán)境、生態(tài)緊密相關(guān), 高強度的經(jīng)濟開發(fā)加劇土地利用變化和多種環(huán)境與生態(tài)問題。例如, 城鎮(zhèn)化進程加快濱海地帶的開發(fā), 圍填海規(guī)模擴大, 造成紅樹林等類型濕地破壞、耕地保護壓力加大、工業(yè)用地利用效率提升壓力加劇等問題[1-5]。合理開發(fā)利用沿海土地資源、協(xié)調(diào)經(jīng)濟發(fā)展與環(huán)境保護之間的關(guān)系日益受到重視, 因此, 許多學(xué)者對沿海土地利用及相關(guān)問題進行研究, 例如, 侯西勇等[6]分析了21世紀初中國海岸帶土地利用的數(shù)量、結(jié)構(gòu)和空間格局特征; 傅伯杰等[7]研究了環(huán)渤海土地利用變化的生態(tài)環(huán)境效應(yīng)及可持續(xù)性; 吳泉源等[8]從土地利用類型相互轉(zhuǎn)換及海岸區(qū)位效應(yīng)等方面分析龍口市海岸帶動態(tài)變化特征; 許艷等[9]分析了江蘇省海岸帶1980年以來的土地利用/覆被變化, 揭示了海岸帶土地利用的地域差異與梯度分異特征。分析土地利用類型變化規(guī)律, 模擬未來土地利用變化的趨勢, 可以為決策者科學(xué)規(guī)劃用地提供有益的指導(dǎo)[10-11]。

環(huán)北部灣區(qū)域是我國與東盟國家既有海上通道、又有陸地接壤的交通樞紐區(qū)域, 區(qū)位優(yōu)勢明顯。2007年廣東省湛江市加快建設(shè)粵西城鎮(zhèn)群中心; 2008年國家大力支持北部灣經(jīng)濟區(qū)開放發(fā)展; 2015年一帶一路戰(zhàn)略北部灣區(qū)域成為21世紀“海上絲綢之路”南海對外海上貿(mào)易的重要門戶。北部灣沿海正在進入經(jīng)濟社會快速發(fā)展的新階段, 學(xué)術(shù)界對該區(qū)域土地資源開發(fā)和利用的關(guān)注也在逐漸提升, 例如:何文等[12]分析了廣西欽州市的土地利用格局; 華璀[13]研究了玉林市土地資源的可持續(xù)利用問題; 王倩[14]研究了防城港海岸帶景觀格局變化; 周飛等[15]、陳炳祿等[16]研究了廣東省湛江市的土地利用與生態(tài)系統(tǒng)服務(wù)價值變化; 姚靜林等[17]、高義等[18]分別分析了北部灣經(jīng)濟區(qū)、廣東省海島海岸帶土地利用變化特征及其驅(qū)動力; 劉雙等[19]、蔡曉樺等[20]分析了北部灣土地資源可持續(xù)利用問題及土地利用結(jié)構(gòu)特征。本文旨在分析北部灣沿海近期土地利用變化特征, 并針對至2025年的變化進行多情景分析和模擬, 對比不同目標導(dǎo)向經(jīng)濟社會發(fā)展?jié)撛谟绊懽饔玫牟町? 為北部灣土地資源的優(yōu)化利用決策提供依據(jù)和支持。

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

1.1 研究區(qū)

北部灣地處我國大陸海岸帶西南端、南海西北部, 是一個半封閉海灣。研究區(qū)選擇廣西沿海區(qū)域以及雷州半島, 包括防城港市、欽州市、北海市、玉林市和湛江市。位于東經(jīng)107°27′~110°58′、北緯20°07′~ 23°07′之間, 總面積4.64萬km2(圖1)。西北部以山地、丘陵為主, 東南部以臺地平原為主, 沿海灘涂分布廣泛; 為亞熱帶向熱帶過渡的海洋性季風(fēng)氣候, 光熱資源豐沛; 水資源尚豐, 但時空分布不均, 多暴雨; 沿海港口、旅游資源豐富。

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

根據(jù)北部灣沿海土地資源開發(fā)利用的歷史與現(xiàn)狀特征以及30 m分辨率Landsat TM/ETM+影像的分辨能力, 區(qū)分耕地、林地、草地、濕地、城鎮(zhèn)、農(nóng)村居民點、交通工礦、未利用地和海水9種土地利用類型: 耕地包括水田和坡耕地; 林地包括有林地、疏林地和灌木林; 草地區(qū)分高、中、低不同覆蓋度; 濕地包括內(nèi)陸水體(湖泊、河流、灘地和水庫坑塘等)、濱海濕地(灘涂、河口水域和瀉湖等)和人工濕地(養(yǎng)殖); 城鎮(zhèn)指縣鎮(zhèn)以上的建成區(qū); 農(nóng)村居民點指鎮(zhèn)以下的居民點用地; 交通工礦是城鎮(zhèn)以外的廠礦、工業(yè)區(qū)、倉庫、機場、鹽田等; 未利用地主要是山區(qū)裸巖和沿海鹽堿地; 海水是指未被利用的近海海域。利用2000年、2005年、2010年的Landsat TM/ETM+影像, 結(jié)合解譯標志庫[21]通過目視解譯獲得三個時期的土地利用數(shù)據(jù), 解譯時最小上圖圖斑2~3個像元, 對應(yīng)1︰10萬比例尺制圖標準。

土地利用分布與變化的自然與人文影響因子選取高程、坡度、坡向、溫度、降雨量、到公路距離、到鐵路距離、到河流距離、到海岸線距離、到居民點距離、人口、GDP共12個[22]。選用30 m分辨率的SRTM DEM信息, 據(jù)其生成坡度、坡向數(shù)據(jù); 溫度與降雨量取多年平均值, 基于氣象站監(jiān)測數(shù)據(jù)通過地統(tǒng)計分析方法獲得空間數(shù)據(jù); 多種距離因子數(shù)據(jù)是利用1︰100萬基礎(chǔ)地理信息通過GIS空間分析而得到; 人口、GDP信息為中國2010年分縣普查與統(tǒng)計資料, 基于夜間燈光數(shù)據(jù)等進行空間化[23]。以上數(shù)據(jù)轉(zhuǎn)換為IDRISI軟件識別的柵格數(shù)據(jù), 柵格大小為500 m。

2 研究方法

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

基于土地利用轉(zhuǎn)移矩陣進行土地利用變化特征分析。土地利用轉(zhuǎn)移矩陣定量描述一定時期內(nèi)土地利用狀態(tài)之間的相互轉(zhuǎn)換關(guān)系[22], 可直觀反映土地的流向或來源。

2.2 土地利用變化模型

越來越多的模型強調(diào)基于遙感監(jiān)測信息在時空雙維度分析土地利用變化的動力機制, 綜合評估與預(yù)測土地利用/覆蓋變化[2]。元胞自動機(Cellular Automata, CA)、CLUE-S、CA-Markov、IMAGE等是較有代表性的模型[24-28]。模型應(yīng)用效果受空間位置、空間尺度、數(shù)據(jù)等因素影響, 集成多種模型可以更好的表征和模擬具有多尺度特征和復(fù)雜動力機制的土地利用變化過程。本文選取Logistic-CA-Markov耦合模型: 二分類Logistic回歸分析是一種概率統(tǒng)計分類模型, 用于預(yù)測某一分類的依賴變量(因變量)對一個或多個預(yù)測變量(自變量)的二進制響應(yīng)。通常采用ROC(Receiver operating characteristic)對回歸方程的解釋能力進行檢驗[29-30], ROC值一般介于0.5~1之間, 趨近于0.5表示方程解釋能力較差, 趨近于1表示效果較好。在土地利用變化研究中, 將每個柵格可能出現(xiàn)某種地類的概率進行運算判斷, 模型的目標變量用二分類變量數(shù)據(jù)表示, 模型的解釋變量通過一系列的自然、社會驅(qū)動因素來表示并確定它們之間的定量關(guān)系。其公式為:

式中,P為類型發(fā)生的概率,1,2, …,x為驅(qū)動因素。

元胞自動機(CA)是一種時間、空間、狀態(tài)都離散且具有時空計算特征的動力學(xué)模型, 用來仿真局部規(guī)則和局部聯(lián)系。柵格化的地理特征數(shù)據(jù)具有與元胞自動機相似的時空分布規(guī)則, 因此, 元胞自動機模型適用于地理特征數(shù)據(jù)的動態(tài)模擬[31-32]。在土地利用變化模擬中, 每個元胞對應(yīng)一定的土地利用類型, 土地利用變化的歷史趨勢、土地適宜性以及相關(guān)的政策、經(jīng)濟因素構(gòu)成轉(zhuǎn)換規(guī)則, 當(dāng)可能性超過閾值時, 土地利用類型發(fā)生轉(zhuǎn)化。元胞自動機可表示為:

式中,為元胞有限、離散的狀態(tài)集合,和+1為不同時刻,為元胞鄰域。

馬爾科夫(Markov)過程是一個典型的無后效性隨機過程。馬爾科夫模型應(yīng)用馬爾科夫過程方法分析隨機事件變化規(guī)律, 并預(yù)測未來變化趨勢。設(shè)()是一隨機過程, 當(dāng)時刻所處的狀態(tài)已知時, 在時刻+1的狀態(tài)與在時刻之前的狀態(tài)無關(guān), 各個時刻的狀態(tài)的轉(zhuǎn)變由一個狀態(tài)轉(zhuǎn)移概率矩陣控制。在土地利用變化研究中, 某一時刻的土地利用類型對應(yīng)Markov過程中的“可能狀態(tài)”(), 土地利用類型之間相互轉(zhuǎn)換的面積或比例即為狀態(tài)轉(zhuǎn)移概率P。其公式為:

(3)

式中,()、(+1)分別為、+1時刻的系統(tǒng)狀態(tài);P為狀態(tài)轉(zhuǎn)移概率矩陣。

Markov與CA均為時間離散、狀態(tài)離散的模型, IDRISI軟件將二者集成, 形成了用于土地利用變化模擬的CA-Markov模型; 通過Logistic回歸分析可以獲取轉(zhuǎn)變適宜性圖像集, 即CA轉(zhuǎn)換規(guī)則, 進而形成Logistic-CA-Markov耦合模型; 其綜合了元胞自動機模型強大的模擬復(fù)雜系統(tǒng)空間演變功能和Markov模型長期預(yù)測的優(yōu)勢, 既提高土地利用類型轉(zhuǎn)化的預(yù)測精度, 又可有效模擬土地利用格局的空間變化。研究表明其模擬精度高、模擬結(jié)果可信、有較強的科學(xué)性與實用性[26], 能夠有效地對土地利用的數(shù)量結(jié)構(gòu)和空間分布進行模擬和預(yù)測[22]。

2.3 情景分析

情景分析是假定某種現(xiàn)象或趨勢將持續(xù)到未來的前提下, 對預(yù)測對象可能出現(xiàn)的情況或引起的后果進行預(yù)測的方法。土地利用變化機制復(fù)雜, 發(fā)展態(tài)勢具有不確定性, 適合應(yīng)用情景分析對其未來態(tài)勢進行探測: 通過對土地利用類型加以定性、定量分析, 模擬外部因素可能發(fā)生的多種情景, 對各種可能情景進行預(yù)測, 進而, 結(jié)合土地利用變化模型, 可以模擬各種情景下土地利用變化的時空特征, 幫助人們發(fā)現(xiàn)未來的機會和威脅、優(yōu)勢和劣勢, 實現(xiàn)土地資源的優(yōu)化配置[33-35]。

3 數(shù)據(jù)與結(jié)果

3.1 土地利用結(jié)構(gòu)與格局特征

北部灣沿海多時相土地利用遙感解譯結(jié)果如圖2 (海水未顯示)。土地利用空間格局特征如下: 林地比重最大, 約為50%, 主要分布于西部與北部, 沿海散布有紅樹林; 耕地比重約為30%, 位于坡度較小區(qū)域; 草地較少, 一般分布在山麓地區(qū); 濕地類型多樣且分布廣泛; 農(nóng)村居民點數(shù)量多、分布廣; 城鎮(zhèn)面積較大, 尤其是沿海港口城市; 交通工礦多為獨立工業(yè)區(qū), 部分沿海地區(qū)有鹽田; 未利用地主要為鹽堿地和裸土地, 分布于沿海區(qū)域, 面積較小。

3.2 土地利用面積轉(zhuǎn)移特征

2000年-2005年、2005年-2010年的土地利用轉(zhuǎn)移矩陣分別如表1、表2所示, 反映土地利用的分布數(shù)量及轉(zhuǎn)移特征: 10年間, 耕地減少220 km2, 前5年主要變?yōu)榱值亍竦?、城?zhèn)和交通工礦, 后5年流失速度減緩, 主要變?yōu)闈竦亍⒊擎?zhèn)和農(nóng)村居民點; 草地減少18 km2, 主要變?yōu)榱值睾蜐竦? 未利用地減少51.5 km2, 主要變?yōu)闈竦?、草? 林地增加121.25 km2, 前5年增加較多, 來源主要為耕地, 后5年增加較少, 主要來源為濕地和未利用地, 少量林地轉(zhuǎn)化為濕地和交通工礦。10年間濕地增加了107.25 km2; 前5年未利用地、耕地、林地是其增加的主要來源, 而濕地則主要變?yōu)榻煌üさV(鹽田)和林地; 后5年濕地主要是由海域、耕地轉(zhuǎn)化而來, 水庫坑塘等的面積增大, 而濱海濕地主要變?yōu)榻煌üさV(鹽田), 灘地、灘涂、沿海潟湖面積則顯著減小。10年間城鎮(zhèn)面積增長了65.75 km2, 前、后5年分別以耕地、交通工礦為主要轉(zhuǎn)化來源。10年間農(nóng)村居民點面積略有增加, 耕地為主要來源, 后5年農(nóng)村居民點主要變?yōu)槌擎?zhèn)。10年間交通工礦面積增加49.25 km2, 前5年主要是海水養(yǎng)殖區(qū)變?yōu)榻煌üさV, 后5年圍填海為主要來源。10年間圍填海造陸使得陸域面積約增加61.5 km2, 2005年前圍填海規(guī)模較小, 之后則速度加快、規(guī)模增大, 并導(dǎo)致海岸線的顯著變化; 圍填海主要形成以海水養(yǎng)殖為主的人工濕地、以鹽田為主的交通工礦以及城鎮(zhèn)。

表1 2000年-2005年土地利用類型轉(zhuǎn)換面積(km2)

表2 2005年-2010年土地利用類型轉(zhuǎn)換面積(km2)

上述變化反映了區(qū)域經(jīng)濟社會發(fā)展的特征及其對土地資源的影響: 1)耕地持續(xù)減少、林地不斷增加, 耕地變?yōu)榱值乇砻魍烁€林政策執(zhí)行情況較好, 生態(tài)保護工作得到推進; 但與此同時, 城鎮(zhèn)、交通工礦發(fā)展又占用較多的耕地、林地、濕地, 這表明僅依靠退耕還林(草)政策不足以滿足生態(tài)保護工作的需求。2)圍填海主要用于交通工礦、圍海漁業(yè)養(yǎng)殖和城鎮(zhèn)建設(shè), 侵占了大量的濱海濕地(灘涂、灘地和沿海潟湖), 生物棲息地逐漸被蠶食, 使生態(tài)環(huán)境面臨惡化風(fēng)險。3)生態(tài)保護政策和工作缺乏“陸海統(tǒng)籌”、“陸海聯(lián)動”,因此, 在海岸帶區(qū)域, 土地類型轉(zhuǎn)換加快、經(jīng)濟社會發(fā)展使土地資源和海洋生態(tài)保護的壓力加大。

3.3 Logistic-CA-Markov模型建立

選取2005年土地利用數(shù)據(jù), 進行Logistic回歸分析, ROC檢驗結(jié)果如表3, 其中, 7種地類的ROC值均大于0.9, 表明回歸模型較為可靠。基于回歸方程建立500 m分辨率的土地轉(zhuǎn)移適宜性圖集, 在CA-Markov模型中, 基于2005年土地利用圖, 設(shè)置5×5的鄰域, 迭代次數(shù)為5, 得到2010年土地利用模擬圖(圖3), 將其與對應(yīng)的觀測值(圖2)進行比較, kappa系數(shù)為0.9894, 證明模型模擬效果理想, 能夠支持土地利用變化情景分析與模擬。

3.4 至2025年土地利用變化情景

從區(qū)域土地利用、生態(tài)保護和經(jīng)濟社會發(fā)展的歷史與現(xiàn)狀及政府部門相關(guān)規(guī)劃出發(fā), 設(shè)計三種情景。

表3 Logistic回歸方程系數(shù)及ROC檢驗結(jié)果

注: 1-9分別為耕地、林地、草地、濕地、城鎮(zhèn)、農(nóng)村居民點、交通工礦、未利用地、海水9種土地利用類型。

情景Ⅰ: 既有趨勢延續(xù)(Business as usual)。Markov理論符合既有趨勢延續(xù)情景的特征, 因此, 直接基于上述過程所建立的Logistic-CA-Markov模型, 以2005年和2010年土地利用數(shù)據(jù)為基礎(chǔ), 預(yù)測得到2020年、2025年的土地利用圖。

情景Ⅱ: 城鎮(zhèn)化(經(jīng)濟發(fā)展)優(yōu)先。優(yōu)先滿足城鎮(zhèn)化、經(jīng)濟發(fā)展對土地資源的需求。主要參考《廣西北部灣經(jīng)濟區(qū)發(fā)展規(guī)劃》(2006年-2020年)、《廣東省土地利用總體規(guī)劃》(2006年-2020年)以及北海市、防城港市、欽州市城市總體規(guī)劃(2008年-2025年)等所確定的城鎮(zhèn)化目標, 如: 廣西北部灣經(jīng)濟區(qū)發(fā)展規(guī)劃中北海市、欽州市、防城港市2020年城市建成區(qū)建設(shè)用地分別控制在140 km2、120 km2和70 km2以內(nèi); 北海城市總體規(guī)劃2025年城市建設(shè)200 km2;防城港市城市總體規(guī)劃至2025年城市建設(shè)用地140 km2; 欽州市總體規(guī)劃建設(shè)用地規(guī)模2025年為234.97 km2, 其中填海79 km2; 玉林中心城區(qū)城市建設(shè)用地規(guī)模至2020年為98 km2; 到2020年湛江市中心城區(qū)城市建設(shè)用地規(guī)?？刂圃?11.62 km2以內(nèi)。綜合以上信息, 設(shè)置2020年、2025年城鎮(zhèn)面積分別增長至約650 km2、約800 km2; 通過Markov面積轉(zhuǎn)移矩陣控制總量目標, 修訂Markov轉(zhuǎn)移概率矩陣中各種地類轉(zhuǎn)向城鎮(zhèn)和交通工礦用地的概率以控制類型轉(zhuǎn)化的合理性, 模擬得到2020年、2025年土地利用圖。

情景Ⅲ: 生態(tài)保護優(yōu)先。設(shè)定原則: 限制城鎮(zhèn)發(fā)展速度; 保護山區(qū)的林草地; 嚴格保護主要水源地和坡度15°以上區(qū)域; 保護沿海灘涂和紅樹林濕地, 限制圍填海工程發(fā)展; 嚴格限制自然保護區(qū)內(nèi)的土地利用變化。重點選取北侖河口國家級自然保護區(qū)、山口紅樹林國家級自然保護區(qū)及湛江紅樹林國家級自然保護區(qū)及眾多水源地列為生態(tài)保護的優(yōu)先區(qū)域, 限制其內(nèi)部的土地利用變化, 模擬得到2020年和2025年土地利用圖。

三種情景的模擬結(jié)果如表4所示。人類活動、經(jīng)濟發(fā)展促進海水養(yǎng)殖區(qū)、圍填海用地等擴大, 呈擴張趨勢。與2010年觀測值相比, 情景Ⅰ未來時期耕地、林地、草地及未利用地面積將減少, 城鎮(zhèn)、農(nóng)村居民點、交通工礦面積將擴大, 主要占用耕地; 情景Ⅱ中, 城鎮(zhèn)向周圍擴張速度更迅速, 對耕地和濕地的占用更為嚴重, 并進一步導(dǎo)致海岸線向海推進; 情景Ⅲ中, 耕地、林地、草地面積緩慢增長, 城鎮(zhèn)化速度受到節(jié)制, 城鎮(zhèn)、農(nóng)村居民點與交通工礦的擴張速度均低于既有趨勢延續(xù)情景和城鎮(zhèn)化(經(jīng)濟發(fā)展) 優(yōu)先情景。三種情景模擬結(jié)果在結(jié)構(gòu)、變化態(tài)勢方面與情景設(shè)計初衷較為一致。

表4 北部灣沿海土地利用模擬情景面積(km2)

將三種情景的空間格局進行比較, 差異較顯著的“熱點”區(qū)域包括城鎮(zhèn)和自然保護區(qū)邊緣、海岸線附近等, 例如: 北海市區(qū)(圖4-A), 情景Ⅰ與Ⅱ均表現(xiàn)為城鎮(zhèn)顯著擴張趨勢, 占用耕地和圍填海特征突出; 沿海區(qū)域圍填海形成較大面積的城鎮(zhèn)與交通工礦, 將海岸線向海推進, 內(nèi)陸區(qū)域城鎮(zhèn)擴張占據(jù)大量的農(nóng)田; 在情景Ⅲ中圍填海與耕地的占用均有所減緩, 因此, 沿海地區(qū)應(yīng)著重控制城鎮(zhèn)發(fā)展規(guī)模、保護耕地、合理規(guī)劃圍填海區(qū)域。山口紅樹林國家級自然保護區(qū)(圖4-B), 情景Ⅰ與Ⅱ中, 西部與東部沿海的林地與濕地被破壞變?yōu)榻煌üさV用地, 而情景Ⅲ保護區(qū)內(nèi)的植被和濕地得以保留, 可見, 自然保護區(qū)邊緣的土地利用較易被改變。將2010年廣西濱海的保護區(qū)與圍填海變化“熱點”區(qū)域進行疊加分析, 重合部分達到24.74 km2, 距“熱點”0~1 km、1~2 km、2~3 km范圍內(nèi)的保護區(qū)面積分別為77.34、105.96、113.19 km2, 可見, 圍填海實際上已經(jīng)嚴重威脅濱海保護區(qū); 為保證自然保護區(qū)內(nèi)的土地利用類型得以保持, 優(yōu)化保護區(qū)及其周邊土地利用結(jié)構(gòu), 嚴格的保護措施迫在眉睫。

A. 北海市區(qū); B. 山口紅樹林國家級自然保護區(qū)

A. Beihai city; B. Shankou mangrove ecological national natural reserve

4 結(jié)論

基于遙感技術(shù)監(jiān)測北部灣沿海2000年-2010年的土地利用變化特征, 應(yīng)用Logistic-CA-Markov耦合模型模擬土地利用變化并進行多情景分析和模擬, 主要結(jié)論如下:

1) 土地利用變化反映經(jīng)濟社會發(fā)展的基本特征: 2000年-2010年, 耕地、草地、濕地等類型的面積減少, 林地、城鎮(zhèn)和農(nóng)村居民點的不斷擴張均占用大量耕地, 耕地保護壓力加大; 2005年之后進入海岸帶開發(fā)提速階段, 因而2005年-2010年與2000年-2005年相比, 土地利用類型間轉(zhuǎn)換的規(guī)模與速度均有提高, 尤其是沿海區(qū)域大量灘涂與濕地被圍墾, 主要變?yōu)榻煌üさV、城鎮(zhèn)、漁業(yè)養(yǎng)殖等, 不斷將海岸線向海推進。

2) 情景分析與模擬結(jié)果揭示了未來時期土地利用變化的可能態(tài)勢: 既有趨勢延續(xù)情景, 耕地、林地、草地與未利用土地將繼續(xù)減少, 城鎮(zhèn)、農(nóng)村居民點與交通工礦則持續(xù)增加, 濕地主要表現(xiàn)為水庫坑塘、河口水域、淺海水域面積增大, 而灘地、灘涂、沿海潟湖面積顯著減小; 城鎮(zhèn)化(經(jīng)濟發(fā)展)優(yōu)先情景, 城鎮(zhèn)發(fā)展將消耗大量耕地、林地, 沿海城市圍填海特征突出; 生態(tài)保護優(yōu)先情景, 對保護區(qū)和沿海濕地的保護成效非常顯著, 耕地、林地、草地有所增加, 城鎮(zhèn)建設(shè)規(guī)模過快增加的態(tài)勢得到抑制。

3) 北部灣沿海2000年-2010年土地利用變化的“熱點”區(qū)域集中在城鎮(zhèn)周邊及海岸帶區(qū)域, 與此相應(yīng), 至2025年不同變化情景模擬結(jié)果之間的差異以及未來潛在的熱點變化區(qū)域也主要集中于城鎮(zhèn)邊緣、海岸帶、自然保護區(qū)邊緣等生態(tài)敏感的區(qū)域, 因此, 未來時期土地利用優(yōu)化管理和生態(tài)保護的基本策略和優(yōu)先目標應(yīng)該包括: 生態(tài)保護政策及工作應(yīng)強化“陸海統(tǒng)籌”和“陸海聯(lián)動”; 加強海洋功能區(qū)劃, 明確保護區(qū)、圍填海建設(shè)區(qū)的范圍, 并加強管制; 優(yōu)化城鎮(zhèn)空間布局及發(fā)展質(zhì)量、提高建設(shè)用地集約化水平、協(xié)調(diào)城市化與耕地及濕地保護之間的關(guān)系、協(xié)調(diào)經(jīng)濟發(fā)展與沿海生態(tài)保護之間的關(guān)系等。

參考文獻：

[1] Turner BL, David S, Steven S. Land use and land cover change[J]. Earth Science Frontier (China University of Geosciences, Beijing), 1997, 4(1-2): 26-33.

[2] Brown D G, Walker R, Manson S, et al. Modeling land use and land cover change[C]//Gutman G, Janetos A C, Justice C O et al. Land Change Science: Observing, Monitoring and Understanding Trajectories of Change on the Earth’s Surface. New York: Springer, 2004: 395- 409.

[3] 劉紀遠, 匡文慧, 張增祥, 等. 20世紀80年代末以來中國土地利用變化的基本特征與空間格局[J]. 地理學(xué)報, 2014, 69(1): 3-14. Liu Jiyuan, Kuang Wenhui, Zhang Zengxiang, et al. Spatio-temporal characteristics, patterns and causes of land use changes in China since the late 1980s[J]. Acta Geographica Sinica, 2014, 69(1): 3-14.

[4] 李興鋼, 梁成華, 王延松, 等. 基于CA-Markov模型的遼河三角洲濕地景觀格局預(yù)測[J]. 環(huán)境科學(xué)與技術(shù), 2013, 36(5): 188-192, 200. Li Xinggang, Liang Chenghua, Wang Yansong, et al. Predicting landscape patterns of Liaohe Delta wetland by CA-Markov Model[J]. Environmental Science and Technology, 2013, 36(5): 188-192, 200.

[5] 吳季秋, 俞花美, 肖明, 等. 基于CA-Markov模型的海灣土地利用動態(tài)預(yù)測[J]. 環(huán)境工程技術(shù)學(xué)報, 2012, 2(6): 531-539. Wu Jiqiu, Yu Huamei, Xiao Ming, et al. Dynamic prediction for gulf land use based on CA Markov Model[J]. Journal of Environmental Engineering Technology, 2012, 2(6): 531-539.

[6] 侯西勇, 徐新良. 21世紀初中國海岸帶土地利用空間格局特征[J]. 地理研究, 2011, 30(8): 1370-1379. Hou Xiyong, Xu Xinliang. Spatial patterns of land use in coastal zones of China in the early 21st century[J]. Geographical Research, 2011, 30(8): 1370-1379.

[7] 傅伯杰, 陳利頂, 蔡運龍, 等. 環(huán)渤海地區(qū)土地利用變化及可持續(xù)利用研究[M].北京: 科學(xué)出版社, 2004. Fu Bojie, Chen Liding, Cai Yunlong, et al. Studies on land use change and its sustainable utilization in Circum Bohai Bay[M]. Beijing: Science Press, 2014.

[8] 吳泉源, 侯志華, 于竹洲, 等. 龍口市海岸帶土地利用動態(tài)變化分析[J]. 地理研究, 2006, 25(5): 921-929. Wu Quanyuan, Hou Zhihua, Yu Zhuzhou, et al. Analysis of the dynamic change of land use in Longkou city's coastal zone based on remote sensing technology[J]. Geographical Research, 2006, 25(5): 921-929.

[9] 許艷, 濮勵杰, 張潤森, 等. 近年來江蘇省海岸帶土地利用/覆被變化時空動態(tài)研究[J]. 長江流域資源與環(huán)境, 2012, 21(5): 565-571. Xu Yan, Pu Lijie, Zhang Runsen, et al. Spatiao temporal dynamics of land use and land cover change in the coastal zone of Jiangsu Province[J]. Resources and Environment in the Yangtze Basin, 2012, 21(5): 565-571.

[10] 嚴靜, 王靜, 周小丹, 等. 鹽城市沿海灘涂圍墾區(qū)土地利用時空動態(tài)變化研究[J]. 安徽農(nóng)業(yè)科學(xué), 2015, 43(29): 323-325. Yan Jing, Wang Jing, Zhou Xiaodan, et al. Study on the inning flats land use spatial-temporal dynamic change in Yancheng City[J]. Journal of Anhui Agricultural Sciences, 2015, 43(29): 323-325.

[11] 宋向群, 蔣欣昀, 王文淵, 等. 基于CA-Markov模型的港口陸域演化研究[J]. 港工技術(shù), 2015, 52(3): 5-8. Song Xiangqun, Jiang Xinyun, Wang Wenyuan, et al. Research of port land area evolution by using CA-Markov Model[J]. Port Engineering Technology, 2015, 52(3): 5-8.

[12] 何文, 盧遠, 賀佩琳, 等. 北部灣經(jīng)濟區(qū)土地利用格局模擬及其生態(tài)環(huán)境效益評價—— 以廣西壯族自治區(qū)欽州市為例[J]. 水土保持通報, 2015, 35(1): 281-286. He Wen, Lu Yuan, He Peilin, et al. Dynamic simulation and eco-environmental benefits evaluation of land use pattern in Beibu Gulf Economic Zone——A case study of Qinzhou City in Guangxi Zhuang Autonomous Region[J]. Bulletin of Soil and Water Conservation, 2015, 35(1): 281-286.

[13] 華璀. 廣西玉林市土地資源可持續(xù)利用研究[J]. 經(jīng)濟地理, 1999, 19(4): 85-89. Hua Cui. Research on the sustainable use of land resource in Yulin city, Guangxi[J]. Economic Geography. 1999, 19(4): 85-89.

[14] 王倩. 防城港海岸帶景觀格局變化研究[D]. 廣西大學(xué). 2013. Wang Qian. The research about the changes on landscape pattern in Fang chenggang coastal zone[D]. Guangxi University, 2013.

[15] 周飛, 陳士銀, 鐘來元, 等. 區(qū)域土地利用與生態(tài)系統(tǒng)服務(wù)價值變化研究—— 以廣東省湛江市為例[J]. 土壤, 2008, 40(5): 847-851. Zhou Fei, Chen Shiyin, Zhong Laiyuan, et al. On changes of regional land use and ecosystem services value-A case study of Zhanjiang City, Guangdong Province[J]. Soils, 2008, 40(5): 847-851.

[16] 陳炳祿, 陳新庚, 吳群河. 湛江市土地利用生態(tài)適宜性評價[J]. 中山大學(xué)學(xué)報(自然科學(xué)版), 1998, 37(suppl2): 217-220. Chen Binglu, Chen Xingeng, Wu Qunhe. Ecological evaluation of land use suitability in Zhanjiang[J]. Acta Scientiarum Naturalium Universitatis Sunyatseni, 1998, 37(suppl2): 217-220.

[17] 姚靜林, 楊小雄. 廣西北部灣經(jīng)濟區(qū)土地利用變化及其驅(qū)動力研究[J]. 大眾科技, 2012, 14(5): 224-227. Yao Jinglin, Yang Xiaoxiong. Research of land use change and its driving force in Guangxi Beibu Gulf Economic Zone[J]. Popular Science & Technology, 2012, 14(5): 224-227.

[18] 高義, 蘇奮振, 孫曉宇, 等. 近20年廣東省海島海岸帶土地利用變化及驅(qū)動力分析[J]. 海洋學(xué)報, 2011, 33(4): 95-103. Gao Yi, Su Fenzhen, Sun Xiaoyu, et al. A study on driving forces of land use change of Guangdong Province coastal zone and islands in recent 20a[J]. Acta Oceanologica Sinica, 2011, 33(4): 95-103. (in Chinese with English abstract)

[19] 劉雙, 胡動剛. 北部灣土地資源可持續(xù)利用評價研究[J]. 農(nóng)業(yè)技術(shù)經(jīng)濟, 2013, 11: 57-63. Liu Shuang, Hu Donggang. Evaluation of sustainable utilization of land resource in Beibu Gulf[J]. Agricultural Technical Economy, 2013, 11: 57-63.

[20] 蔡曉樺, 李月連. 基于I-TOPSIS 法的廣西北部灣經(jīng)濟區(qū)土地利用結(jié)構(gòu)研究[J]. 測繪與空間地理信息, 2012, 35(4): 42-45. Cai Xiaoye, Li Yuelian. Study of land use structure of Guangxi Beibu Gulf Economic Zone[J]. Geomatics & Spatial Information Technology, 2012, 35(4): 42-45.

[21] 邸向紅, 侯西勇, 吳莉. 中國海岸帶土地利用遙感分類系統(tǒng)研究[J]. 資源科學(xué), 2014, 36(3): 463-472. Di Xianghong, Hou Xiyong, Wu li. Land use classification system for China’s coastal zone based on remote sensing[J]. Resources Science, 2014, 36(3): 463-472.

[22] 吳莉, 侯西勇, 徐新良, 等. 山東沿海地區(qū)土地利用和景觀格局變化[J]. 農(nóng)業(yè)工程學(xué)報, 2013, 29(5): 207-216. Wu Li, Hou Xiyong, Xu Xinliang, et al. Land use and landscape pattern changes in coastal areas of Shandong province, China[J]. Transactions of the Chinese Society of Agricultural Engineering, 2013, 29(5): 207-216.

[23] 陳晴, 侯西勇. 集成土地利用數(shù)據(jù)和夜間燈光數(shù)據(jù)優(yōu)化人口空間化模型[J].地球信息科學(xué)學(xué)報, 2015, 17(11): 1370-1377. Chen Qing, Hou Xiyong. An improved population spatialization model by combining land use data and DMSP/OLS Data[J]. Journal of Geo-Information Science, 2015, 17(11): 1370-1377.

[24] 張華, 張勃. 國際土地利用|覆蓋變化模型研究綜述[J].自然資源學(xué)報, 2005, 20(3): 422-431. Zhang Hua, Zhang Bo. Review on land use and land cover change models [J]. Journal of Natural Resources, 2005, 20(3): 422-431.

[25] 陸汝成, 黃賢金, 左天惠, 等. 基于CLUE-S和Markov復(fù)合模型的土地利用情景模擬研究——以江蘇省環(huán)太湖地區(qū)為例[J]. 地理科學(xué), 2009, 29(4): 577- 581. Lu Rucheng, Huang Xianjin, Zuo Tianhui, et al. Land use scenarios simulation based on CLUE-S and Markov composite model: A case study of Taihu lake rim in Jiangsu province[J]. Scientia Geographica Sinica, 2009, 29(8): 577-581.

[26] 何丹, 金鳳君, 周璟. 基于Logistic-CA-Markov 的土地利用景觀格局變化—— 以京津冀都市圈為例[J]. 地理科學(xué), 2008, 31(8): 903-910. He Dan, Jin Fengjun, Zhou Jing. The changes of land use and landscape pattern based on Logistic-CA- Markov model: A case study of Beijing-Tianjin-Hebei metropolitan region[J]. Scientia Geographica Sinica, 2011, 31(8): 903-910.

[27] Wang N H, Daniel GBr, Li A, et al. Comparative performance of logistic regression and survival analysis for detecting spatial predictors of land-use[J]. International Journal of Geographical Information Science, 2013, 27(10): 1960-1982.

[28] 侯西勇, 常斌, 于信芳. 基于CA-Markov的河西走廊土地利用變化研究[J]. 農(nóng)業(yè)工程學(xué)報, 2004, 20(5): 286-291. Hou Xiyong, Chang Bin, Yu Xinfang. Land use change in Hexi corridor based on CA-Markov methods[J]. Transactions of the Chinese Society of Agricultural Engineering, 2004, 20(5): 286-291.

[29] 李強, 任志遠. 基于Logistic回歸分析的土地利用變化空間統(tǒng)計與模擬[J]. 統(tǒng)計與信息論壇, 2012, 27(3): 98-103. Li Qiang, Ren Zhiyuan. Spatial statistics and simulation of the land use change based on binary logistic regression[J]. Statistics & Information Forum, 2012, 27(3): 98-103.

[30] Pontius RG, Laura CS. Land-cover change model validation by an ROC method for the Ipswich watershed, Massachusetts, USA[J]. Agriculture, Ecosystems and Environment, 2001, 85 (1/2/3): 239-248.

[31] 黎夏, 葉嘉安. 知識發(fā)現(xiàn)及地理元胞自動機[J]. 中國科學(xué)D輯地球科學(xué), 2004, 34 (9): 865-872. Li Xia, Ye Jiaan. Mining cellular automata’s transition rules[J]. Science China (Series D, Earth Sciences), 2004, 34(9): 865-872.

[32] 楊俊, 解鵬, 席建超, 等. 基于元胞自動機模型的土地利用變化模擬——以大連經(jīng)濟技術(shù)開發(fā)區(qū)為例[J]. 地理學(xué)報, 2015, 70(3): 461-475. Yang Jun, Xie Peng , Xi Jianchao, et al. LUCC simulation based on the cellular automata simulation: A case study of Dalian Economic and Technological Development Zone[J]. Acta Geographica Sinica, 2015, 70(3): 461-475.

[33] Robinson, D. T., Brown, D. G.. Evaluating the effects of land-use development policies on ex-urban forest cover: An integrated agent-based GIS approach[J]. International Journal of Geographical Information Science, 2009, 23(9): 1211-1232.

[34] Hubacek K, Sun L X. A scenario analysis of China’s land use and land cover change incorporating biophysical information into input-output modeling[J]. Structural Change and Economic Dynamics, 2001, (12): 367-397.

[35] 戴聲佩, 張勃. 基于CLUE_S模型的黑河中游土地利用情景模擬研究——以張掖市甘州區(qū)為例[J]. 自然資源學(xué)報, 2013, 28(2): 336-348. Dai Shengpei, Zhang Bo. Land use change scenarios simulation in the middle reaches of the Heihe River Basin based on CLUE-S model—A case of Ganzhou district of Zhangye city[J]. Journal of Natural Resources, 2013, 28(2): 336-348.

(本文編輯: 康亦兼)

Spatio-temporal characteristics and scenario analysis of land-use change in the Beibu Gulf Economic Rim coastal area, China

SU Hong-fan1, HOU Xi-yong1, DI Xiang-hong2

(1. Yantai Institute of Coastal Zone Research, Chinese Academy of Science, Yantai 264003, China; 2. College of Resources, Environment and Planning, Dezhou University, Dezhou 253023, China)

Considering the Beibu Gulf Economic Rim coastal area as an example, land-use maps from 2000, 2005, and 2010 were created on the basis of remote sensing techniques. Moreover, the spatio-temporal characteristics of land-use change from 2000 to 2010 were revealed. Then, a Logistic-CA-Markov model was developed to simulate land-use change and forecast its future tendencies under different scenarios. Results show that from 2000 to 2010, farmland, grassland, coastal wetland, and unused land decreased, whereas forest, inland wetland, urban areas, rural settlements, and isolated industrial-mining areas increased. Farmland is the main source of land transformation and has been changed into forest, rural settlements, urban areas, and isolated industrial-mining areas. In addition, from 2005 to 2010, parts of the increased urban and isolated industrial-mining areas were reclaimed from coastal wetland and sea areas, which in turn, markedly changed the mainland shoreline. Under a business-as-usual scenario, land-use change will follow the trend as from 2005 to 2010; under an urbanization priority scenario, urban areas, rural settlements, and isolated industrial-mining areas will increase significantly and will occupy farmland, water bodies, or wetlands in both inland and coastal areas; under an ecological protection priority scenario, wetlands, natural conservation areas, vegetation cover, and farmlands will be effectively protected from excessive reclamation. Information from this study can be useful for strategies on urbanization optimization, coastal wetland conversation, and reasonable land-use planning in hotspot areas.

Land-use change; Beibu Gulf Economic Rim; Logistic Regression; CA-Markov model; Scenario Analysis

Jul. 8, 2015

S159

A

1000-3096(2016)09-0107-10

10.11759/hykx20150708001

2015-05-04;

2016-05-06

國家自然科學(xué)基金項目(31461143032), 中國科學(xué)院重點部署項目(KZZD-EW-TZ-15, KZZD-EW-14)

蘇紅帆(1986-), 女, 河北石家莊人, 碩士研究生, 主要研究方向: 海岸帶土地利用/覆蓋變化, E-mail: hf_su@sina.com; 侯西勇,通信作者, E-mail: xyhou@yic.ac.cn

[Foundation: National Natural Science Foundation of China, No.31461143032; Key Research Program of the Chinese Academy of Sciences, No. KZZD- EW-TZ-15 and No. KZZD-EW-14]