深圳灣海岸帶景觀格局變化特征和趨勢(shì)預(yù)測(cè)

羅艷，婁全勝，謝健，李騫

(1.國(guó)家海洋局南海規(guī)劃與環(huán)境研究院 廣州 510300；2.中山大學(xué)地理科學(xué)與規(guī)劃學(xué)院 廣州 510275)

深圳灣海岸帶景觀格局變化特征和趨勢(shì)預(yù)測(cè)

羅艷1，婁全勝1，謝健1，李騫2

(1.國(guó)家海洋局南海規(guī)劃與環(huán)境研究院 廣州 510300；2.中山大學(xué)地理科學(xué)與規(guī)劃學(xué)院 廣州 510275)

文章以深圳灣海岸帶為研究對(duì)象，在遙感和GIS技術(shù)的支持下，應(yīng)用2001年和2009年遙感影像解譯信息，建立空間信息庫(kù)并疊加分析，獲取景觀類型動(dòng)態(tài)轉(zhuǎn)移數(shù)據(jù)矩陣，得出2001-2009年深圳灣海岸帶景觀空間分布規(guī)律及其變化特征，即紅樹林、建設(shè)用地和河流的面積有所擴(kuò)張，湖泊、灘涂沙洲、林地、淺海水域和綠地的面積有所減少。同時(shí)，應(yīng)用馬爾科夫模型對(duì)海岸帶景觀格局的變化趨勢(shì)進(jìn)行預(yù)測(cè)，得出淺海水域、灘涂沙洲、湖泊、林地和綠地的比重將繼續(xù)下降，建設(shè)用地和紅樹林的比重將持續(xù)上升，河流的比重基本不變的結(jié)論。由此提出深圳灣應(yīng)加強(qiáng)對(duì)填海造地工程的規(guī)劃，保護(hù)海洋生態(tài)環(huán)境，實(shí)現(xiàn)科學(xué)用海和生態(tài)用海。

海岸帶；景觀格局；趨勢(shì)預(yù)測(cè)；深圳灣；海洋生態(tài)環(huán)境

1 引言

海岸帶是潮間帶及其兩側(cè)一定范圍內(nèi)陸地和淺海的過渡地帶。作為第一海洋經(jīng)濟(jì)區(qū)，海岸帶也是地球環(huán)境動(dòng)力系統(tǒng)最復(fù)雜和地表生態(tài)系統(tǒng)最脆弱的地帶[1]。諸多內(nèi)外因素尤其是人類活動(dòng)易導(dǎo)致海岸帶景觀結(jié)構(gòu)和功能的復(fù)雜化[2-3]，基于生態(tài)系統(tǒng)的空間分布和隨時(shí)間變化的空間信息研究景觀的時(shí)空變化、分析其變化原因以及預(yù)測(cè)其變化趨勢(shì)，可為區(qū)域景觀的管理和規(guī)劃提供決策依據(jù)[4]。

隨著人口和產(chǎn)業(yè)向沿海城市不斷聚集，城市土地資源緊張、土地價(jià)格猛漲，大多數(shù)沿海城市通過填海造地增加土地資源。深圳市作為我國(guó)沿海首個(gè)經(jīng)濟(jì)特區(qū)，自1980年至今已開展過多次較大規(guī)模的填海造地工程，尤其20世紀(jì)90年代是深圳市填海造地高速增長(zhǎng)階段，其中以深圳灣較為典型[5]。大范圍的填海造地工程改變了深圳灣的自然岸線形態(tài)，同時(shí)可能給周圍的海陸生態(tài)系統(tǒng)造成不良影響，引起很多學(xué)者的關(guān)注。本研究采用深圳灣海岸帶2001年和2009年2個(gè)時(shí)期的衛(wèi)星影像數(shù)據(jù)進(jìn)行解譯，并開展空間疊加分析，反映該地區(qū)8年間海岸帶景觀格局的空間變化規(guī)律；同時(shí)，應(yīng)用馬爾科夫模型模擬動(dòng)態(tài)演變過程，預(yù)測(cè)景觀格局的變化趨勢(shì)，為深圳灣海岸帶的管理和保護(hù)提供科學(xué)依據(jù)。

2 研究區(qū)概況

深圳灣亦稱后海灣，是深圳3大海灣之一，整體呈SW—NE向延伸，口門外與伶仃洋東槽礬石水道-暗士頓水道相接，平均水深2.9 m，自然水深一般不超過5 m，是半封閉的小海灣[6]。

海岸帶是陸地系統(tǒng)和海洋系統(tǒng)的交界地帶，其范圍劃分尚無(wú)統(tǒng)一標(biāo)準(zhǔn)[7]。目前大部分采用的劃定原則是海岸線向陸延伸5 km或10 km，向海延伸5 km或10 km或至15 m等深線[8]。為保留海岸帶景觀的完整性，本研究擬向陸側(cè)以沿海城市干道為界，無(wú)干道處則以鄰近山體向陸側(cè)邊界為界；向海側(cè)以深圳市和香港特別行政區(qū)的分界線為界。由此界定研究區(qū)面積約為9 834.21 hm2。

3 研究方法

3.1 數(shù)據(jù)來(lái)源

研究數(shù)據(jù)來(lái)源于深圳市2001年Landsat-TM 影像和2009年SPOT影像。在ArcGIS10.0軟件下，將計(jì)算機(jī)解譯和目視解譯相結(jié)合提取景觀信息，將2001年最小斑塊(面積為0.27 hm2)作為2個(gè)時(shí)期影像識(shí)別斑塊的最小面積，以保證分辨率統(tǒng)一。

根據(jù)實(shí)地考察和影像分析確定解譯標(biāo)志，將海岸帶景觀分為淺海水域、紅樹林、灘涂沙洲、河流、湖泊、建設(shè)用地、林地和綠地8個(gè)類型。其中海水范圍內(nèi)均為淺海水域，紅樹林在海岸帶邊緣呈帶狀分布，深圳灣長(zhǎng)期填海造地產(chǎn)生的大片人造灘涂歸為灘涂沙洲，綠地主要是早期陸地未開發(fā)的綠地和城市綠化用地。

3.2 轉(zhuǎn)移矩陣

采用ArcGIS10.0軟件生成矢量數(shù)據(jù)和建立空間信息庫(kù)，并通過空間疊加功能生成景觀類型動(dòng)態(tài)變化圖，從而獲得動(dòng)態(tài)變化數(shù)據(jù)庫(kù)，由此得到2個(gè)時(shí)期深圳灣海岸帶景觀類型轉(zhuǎn)換矩陣，進(jìn)一步分析各景觀類型間的轉(zhuǎn)化特征。

3.3 馬爾科夫模型

馬爾科夫模型是關(guān)于時(shí)間發(fā)生的概率預(yù)測(cè)模型。馬爾科夫過程是指在一系列特定時(shí)間間隔下，某亞穩(wěn)態(tài)系統(tǒng)由t時(shí)刻狀態(tài)向t+1時(shí)刻狀態(tài)轉(zhuǎn)化的過程，這種轉(zhuǎn)化要求t+1時(shí)刻的狀態(tài)只與t時(shí)刻的狀態(tài)有關(guān)[9]。大量研究指出，景觀變化預(yù)測(cè)是由各景觀類型之間相互轉(zhuǎn)化的轉(zhuǎn)移概率矩陣決定的[10-12]，近年來(lái)學(xué)者們已經(jīng)開始采用馬爾科夫模型來(lái)預(yù)測(cè)景觀格局未來(lái)的變化趨勢(shì)[13-16]。

本研究以2001年和2009年深圳灣海岸帶景觀的解譯結(jié)果為基礎(chǔ)，通過馬爾科夫模型計(jì)算得到面積和概率轉(zhuǎn)移矩陣，從而預(yù)測(cè)2017年和2025年深圳灣海岸帶景觀變化情況。

4 結(jié)果和分析

4.1 海岸帶景觀動(dòng)態(tài)變化

2001—2009年研究區(qū)景觀類型轉(zhuǎn)移矩陣如表1所示。由表1可得出2001—2009年研究區(qū)各景觀類型面積變化情況(表2)。

表1 2001—2009年深圳灣景觀類型轉(zhuǎn)換矩陣 hm2

表2 2001—2009年深圳灣景觀類型面積變化情況 hm2

由表1和表2可以看出，2001—2009年深圳灣海岸帶的紅樹林、建設(shè)用地和河流的面積有所擴(kuò)張，而湖泊、灘涂沙洲、林地、淺海水域和綠地的面積有所減少；變化率最大的是湖泊，最小的是河流。

(1)紅樹林是增長(zhǎng)率最高的海岸帶景觀類型，增長(zhǎng)來(lái)源主要是淺海水域、灘涂沙洲和綠地。深圳市生態(tài)建設(shè)紅樹林修復(fù)示范工程自2007年正式啟動(dòng)，深圳灣海岸帶先后開展“深圳灣濱海休閑帶工程”“福田鳳塘河口紅樹林修復(fù)示范工程”和“深圳灣濱海紅樹林修復(fù)工程(一期)”等項(xiàng)目，紅樹林面積的擴(kuò)張?jiān)谝欢ǔ潭壬险f(shuō)明深圳灣紅樹林整治修復(fù)工程成效較顯著。

(2)建設(shè)用地的擴(kuò)張也非常明顯，新增面積最大，來(lái)源主要是綠地、淺海水域、林地和灘涂沙洲。根據(jù)《廣東省統(tǒng)計(jì)年鑒》歷年統(tǒng)計(jì)數(shù)據(jù)，2001—2009年深圳市GDP從2 482.49億元增長(zhǎng)到7 786.79億元，人口從724.57萬(wàn)增長(zhǎng)到954.28萬(wàn)。經(jīng)濟(jì)的迅速增長(zhǎng)和人口的快速增加需更多的城市建設(shè)用地支撐，在濱海范圍內(nèi)開展大規(guī)模填海造地工程是解決土地資源缺乏問題的有效途徑。廣東省填海造地活動(dòng)十分頻繁，填海面積較大，尤其是2004年和2007年[17]；曾輝等[18]在對(duì)深圳市1988-2007年濕地景觀動(dòng)態(tài)變化及其成因的分析中也指出，2000-2004年深圳市大量可開發(fā)的農(nóng)田和園地幾近消失，為增加建設(shè)用地而開展的填海造地工程導(dǎo)致濱海濕地迅速減少。

(3)湖泊減少態(tài)勢(shì)最為明顯，主要轉(zhuǎn)為建設(shè)用地和綠地；淺海水域面積減少最多，主要轉(zhuǎn)為建設(shè)用地和綠地，新增來(lái)源主要為灘涂沙洲；灘涂沙洲總體呈減少態(tài)勢(shì)，主要轉(zhuǎn)為淺海水域和建設(shè)用地，新增來(lái)源主要為淺海水域；河流面積有所增加，林地和綠地面積都有所減少。

4.2 海岸帶景觀變化的馬爾科夫過程模擬和預(yù)測(cè)

根據(jù)表1得到2001—2009年各類海岸帶景觀類型的轉(zhuǎn)移概率矩陣(表3)。

由2001年研究區(qū)海岸帶景觀類型的面積比，得到初始狀態(tài)概率向量P(n-1)，以8年為一步長(zhǎng)，在MATLAB7.1軟件下編程和運(yùn)算，當(dāng)n=1時(shí)預(yù)測(cè)年份為2017年，當(dāng)n=2時(shí)預(yù)測(cè)年份為2025年，以此類推，計(jì)算結(jié)果如表4所示。

表3 2001—2009年深圳灣景觀類型變化轉(zhuǎn)移概率矩陣

表4 深圳灣海岸帶景觀變化的馬爾科夫過程預(yù)測(cè)

通過預(yù)測(cè)結(jié)果可以看出，未來(lái)如按照現(xiàn)有海域開發(fā)規(guī)模和海域使用政策，淺海水域、灘涂沙洲、湖泊、林地和綠地的比重將繼續(xù)下降，但下降速度變緩；建設(shè)用地和紅樹林的比重將繼續(xù)上升，但上升趨勢(shì)減緩并趨于穩(wěn)定，其中建設(shè)用地在各景觀類型中的增長(zhǎng)速度最快；河流的比重將基本不變。

5 結(jié)語(yǔ)

2001—2009年深圳灣海岸帶的紅樹林、建設(shè)用地和河流的面積有所擴(kuò)張，湖泊、灘涂沙洲、林地、淺海水域和綠地的面積有所減少。從馬爾科夫模型預(yù)測(cè)結(jié)果看，如保持深圳市人口和經(jīng)濟(jì)的高增長(zhǎng)模式并以現(xiàn)有的海域使用和開發(fā)模式作為支持，須不斷占用海域海岸帶資源。深圳灣是深圳市最早發(fā)展的地區(qū)之一，作為深圳市和香港特別行政區(qū)重要的海洋資源，深圳灣還具有十分重要的海洋生態(tài)價(jià)值。深圳市經(jīng)濟(jì)的迅速發(fā)展離不開填海造地工程，這從很大程度上解決了深圳市發(fā)展用地不足的問題，但填海造地改變海域自然屬性，將給深圳灣海洋生態(tài)環(huán)境造成巨大壓力。因此，在填海造地工程實(shí)施前應(yīng)做好充分規(guī)劃，開展必要的生態(tài)環(huán)境影響和經(jīng)濟(jì)效益綜合評(píng)價(jià)，實(shí)現(xiàn)科學(xué)用海和生態(tài)用海，從而保障深圳灣在“粵港澳大灣區(qū)”戰(zhàn)略地位升級(jí)背景下的健康可持續(xù)發(fā)展。

利用馬爾科夫模型對(duì)于研究一定區(qū)域面積內(nèi)的海岸帶景觀動(dòng)態(tài)變化較適宜，這是由于一定區(qū)域內(nèi)不同海岸帶景觀類型之間具有相互轉(zhuǎn)化性。但海岸帶景觀類型之間的轉(zhuǎn)化過程包含較多尚難用函數(shù)關(guān)系準(zhǔn)確描述的事件(如空間轉(zhuǎn)化)，需結(jié)合國(guó)家和地方的相關(guān)政策變化進(jìn)一步綜合分析。

[1] 李振,黃海軍.膠州灣海岸帶土地利用/覆被變化研究[J].國(guó)土資源遙感，2010(4):71-76.

[2] 歐維新,楊桂山,李恒鵬,等.蘇北鹽城海岸帶景觀格局時(shí)空變化及驅(qū)動(dòng)力分析[J].地理科學(xué),2004,24(5)：610-615.

[3] 王穎,季小梅.中國(guó)海陸過渡帶：海岸海洋環(huán)境特征與變化研究[J].地理科學(xué),2011,31(2):129-135.

[4] 楊帆,趙冬至,索安寧.雙臺(tái)子河口濕地景觀時(shí)空變化研究[J].遙感技術(shù)與應(yīng)用,2008,23(1):51-56.

[5] 劉中民,張德民.海洋領(lǐng)域的非傳統(tǒng)安全威脅及其對(duì)當(dāng)代國(guó)際關(guān)系的影響[J].中國(guó)海洋大學(xué)學(xué)報(bào)(社會(huì)科學(xué)版),2004(4): 60-64.

[6] 宋紅,陳曉玲.基于遙感影像的深圳灣填海造地的初步研究[J].湖北大學(xué)學(xué)報(bào)(自然科學(xué)版),2004,26(3)：259-263.

[7] 黃寧,楊綿海,林志蘭,等.廈門市海岸帶景觀格局變化及其對(duì)生態(tài)安全的影響[J].生態(tài)學(xué)雜志,2012,31(12):3193-3202.

[8] 蘇勝金.七年全國(guó)海岸帶和海涂資源綜合調(diào)查綜述[J].海洋開發(fā)與管理,1988,5(2): 33-35.

[9] 施仁杰.馬爾科夫鏈基礎(chǔ)及其應(yīng)用[M].西安:西安電子科技大學(xué)出版社,1992:1-3.

[10] 羅艷,楊樹華,徐淑升,等.龍川江流域土地利用格局變化與趨勢(shì)預(yù)測(cè)研究[J].云南地理環(huán)境研究，2007,19(6)：29-33.

[11] 張海龍,蔣建軍,解修平.基于GIS和馬爾科夫模型的渭河盆地景觀動(dòng)態(tài)變化研究[J].干旱區(qū)資源與環(huán)境,2005,19(7)：119-124.

[12] 除多，張鐿鋰，鄭度.拉薩地區(qū)土地利用變化情景分析[J].地理研究,2005,24(6)：869-877.

[11] 王方雄,王博,杜研,等.基于GIS/RS的旅順濱海濕地景觀格局時(shí)空動(dòng)態(tài)變化研究[J].海洋開發(fā)與管理,2014,31(2):84-89.

[12] 汪文勇,向渝,董傳坤,等.用馬爾科夫模型優(yōu)化分布式最小連通支配集算法[J].電子學(xué)報(bào),2010,38(10):2441-2446.

[13] 賈勝韜.曹妃甸工業(yè)園區(qū)地表景觀格局變化預(yù)測(cè)分析[D].蘭州：蘭州交通大學(xué),2015：1-55.

[14] 蘇木蘭,戴文遠(yuǎn),黃華富,等.快速城鎮(zhèn)化背景下海島農(nóng)村居民點(diǎn)空間演變：以福建海壇島為例[J].熱帶地理,2016,36(6):1005-1018.

[15] 李利紅,張華國(guó),史愛琴,等.基于 RS/GIS的西門島海洋特別保護(hù)區(qū)灘涂濕景觀格局變化分析[J].遙感技術(shù)與應(yīng)用,2013,28(1):129-136.

[16] 何東進(jìn),林立,游巍斌,等.閩東濱海濕地景觀格局演化及其模擬[J].福建林學(xué)院學(xué)報(bào),2013.33(2):97-105.

[17] 張翠萍,謝健,婁全勝,等.廣東省填海造地的發(fā)展經(jīng)驗(yàn)及對(duì)策研究[J].海洋環(huán)境科學(xué),2013,32(2):311-315.

[18] 曾輝,啟輝,陳雪,等.深圳市1988-2007年間濕地景觀動(dòng)態(tài)變化及成因分析[J].生態(tài)學(xué)報(bào)，2010,30(10):2706-2714.

TheLandscapePatternCharacteristicsandForecastoftheCoastalZoneinShenzhenBay

LUO Yan1,LOU Quansheng1,XIE Jian1,LI Qian2

(1.South China Sea Institute of Planning and Environmental Research，SOA,Guangzhou 510300,China; 2.School of Geography and Planning,Sun Yat-Sen University,Guangzhou 510275,China)

The Shenzhen Bay coastal zone was regarded as the research object,supported by remote sensing and GIS technology,images of 2001’s and 2009’s were interpreted.The dynamic transfer data matrix of landscape types were obtained through the constructed space information base and overlaid the two kinds of vector maps.The spatial distribution and characteristics of landscape in Shenzhen Bay coastal zone were also obtained from 2001 to 2009.The results showed that the area of mangroves,construction lands and rivers expanded,the area of lakes,beaches,forests,shallow waters and green field decreased.At the same time,the developing tendency of landscape patterns were analyzed and forecasted with Markov Chain Model.The result also showed,the proportion of shallow waters,beaches,lakes,forests,green field would continue to fall,the proportion of construction land,mangroves would continue to rise,and the proportion of the rivers would stay unchanged essentially.Therefore,it is suggested that Shenzhen Bay should strengthen the planning of land reclamation,protect the marine ecological environment,and make use of the sea in a scientific and ecological way.

Coastal zone，Landscape pattern，Trend forecast，Shenzhen Bay，Marine ecological environment

2017-03-21；

2017-09-30

2012年度國(guó)家海洋局南海分局海洋科學(xué)技術(shù)局長(zhǎng)基金資助項(xiàng)目(1257).

羅艷，工程師，碩士，研究方向?yàn)樯鷳B(tài)環(huán)境和海洋環(huán)境管理

P748

A

1005-9857(2017)12-0046-05