近30 年來萊州灣濱海濕地景觀格局變化特征研究

黃建濤 , 鄭 偉, 萬年新, 路景鈁 , 李 捷

(1. 青島理工大學環(huán)境與市政工程學院, 山東 青島 266033; 2. 自然資源部第一海洋研究所, 山東 青島266061; 3. 海洋科學與技術(shù)試點國家實驗室地質(zhì)過程與環(huán)境功能實驗室, 山東 青島 266200; 4. 山東省調(diào)水工程運行維護中心棘洪灘水庫管理站, 山東 青島 266111)

濱海濕地位于陸地與海兩大生境的連接處和過渡地帶, 因其具有較為復(fù)雜的生態(tài)環(huán)境, 形成生物多樣性和景觀格局的多樣化[1-2]。景觀格局是多種生態(tài)過程作用在不同尺度上所產(chǎn)生的結(jié)果[3-5], 景觀格局在不同時間下的變化可以直觀體現(xiàn)出土地利用的改變。研究區(qū)域生態(tài)系統(tǒng)時空演變規(guī)律, 揭示區(qū)域生態(tài)環(huán)境問題形成的驅(qū)動機制是當前的研究熱點[6-7]。萊州灣濱海濕地類型眾多且不同區(qū)位海岸帶開發(fā)利用方式差別較大, 景觀類型多樣、空間異質(zhì)性顯著。目前對萊州灣濱海濕地南岸部分的研究多集中在景觀格局變化、海岸侵蝕、海水入侵、濱海濕地退化和恢復(fù)等方面[8-13], 對萊州灣西岸和東岸研究相對較少且主要集中在重金屬以及污染負荷估算及風險評估等方面[14-17]。由于不同岸段景觀格局變化與驅(qū)動機制也不完全一致, 因此部分岸段的研究很難反映出濱海濕地整體的景觀類型演變及空間異質(zhì)性特征。本研究基于1990—2018 年轉(zhuǎn)移矩陣分析萊州灣濱海濕地景觀格局變化, 在此基礎(chǔ)上開展萊州灣景觀格局變化驅(qū)動機制分析, 探討不同岸段景觀格局變化的驅(qū)動機制; 同時基于 FRAGSTASTS 軟件,采用傳統(tǒng)景觀格局指數(shù)[18]計算方法對萊州灣濱海濕地的景觀格局整體特征和景觀異質(zhì)性空間特征進行研究。以期對萊州灣濱海濕地生態(tài)系統(tǒng)退化原因、合理利用和優(yōu)化土地空間布局、采取有效措施保護及恢復(fù)濱海濕地提供一定的技術(shù)支持[19-20]。

1 研究區(qū)及研究方法

1.1 研究區(qū)概況

萊州灣位于渤海南部, 山東半島西側(cè), 是山東省沿海面積最大的海灣[21-22], 地理位置處于36.96°—37.79°N, 118.73°—120.34°E。其中, 黃河口—羊角溝口為西岸區(qū), 羊角溝口—虎頭崖為南岸區(qū), 虎頭崖—屺姆角為東岸區(qū)。萊州灣年均氣溫11.9～12.6 ℃, 年均降水量為612.5～660.1 mm, 年平均風速4 m/s[23-24]。沿岸有小清河、維河、膠萊河等十余條河流注入, 行政區(qū)包括墾利、東營、壽光等9 個縣級市。參考濕地公約[25]以及《濱海濕地生態(tài)監(jiān)測技術(shù)規(guī)程》(HY/T080-2005)[26], 海水入侵危險區(qū)范圍為距離海岸20 km[10], 且人類活動和遙感影像信息提取難度隨著離岸距離的遠近會呈現(xiàn)一定程度的變化, 綜合萊州灣的入海河流上溯距離以及人類活動對濱海濕地的主要影響范圍, 確定以萊州灣大型海藻的生長區(qū)外緣向內(nèi)陸延伸20 km 作為本文的研究范圍(圖1)。

圖1 萊州灣研究區(qū)及野外觀測點位Fig. 1 Study area and field-observation points of the Laizhou Bay

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

從美國地質(zhì)調(diào)查局(USGS)和地理空間數(shù)據(jù)云(http://www.gscloud.cn/)中選取萊州灣1990 年7 月(Landsat5 TM)、2000 年8 月(Landsat7 ETM+)、2009 年8 月(Landsat5 TM)及2018 年8 月(Landsat8 OLI)四期遙感影像, 行列號120/34、121/34, 空間分辨率皆為30 m, 云層遮蓋度低于5%，滿足目視解譯的需求。利用ENVI 軟件對遙感圖像進行輻射定標、大氣校正和影像鑲嵌處理，并對遙感影像754、543、432 三種波段進行融合, 在ArcGIS10.0 中對萊州灣濱海濕地遙感影像進行目視解譯, 繪制1990—2018 年萊州灣濱海濕地景觀類型分布圖。

1.3 研究方法

1.3.1 分類體系

對濱海濕地類型進行合理劃分是分析其景觀格局變化的關(guān)鍵因素, 依據(jù)國際上對濱海濕地的分類方案, 結(jié)合現(xiàn)場踏勘掌握的情況以及目視解譯的可視性和可判讀性, 將濱海濕地分為兩級(表1)。

1.3.2 景觀類型轉(zhuǎn)移分析

通過轉(zhuǎn)移矩陣分析不同景觀類型的轉(zhuǎn)移方向和強度。景觀類型轉(zhuǎn)移矩陣公式為[27]:

式中,S代表面積;n代表景觀類型總數(shù);i、j(i,j=1, 2,…,n)分別代表轉(zhuǎn)移前與轉(zhuǎn)移后的景觀類型;表示景觀類型轉(zhuǎn)移前后的變化面積。

1.3.3 景觀格局指數(shù)分析

為反映萊州灣景觀格局特征, 在FRAGSTATS軟件的支持下, 選擇斑塊個數(shù)(NP)、景觀形狀指數(shù)(LSI)、最大斑塊指數(shù)(LPI)和斑塊所占景觀面積比例(PLAND)4 個指標, 選用傳統(tǒng)景觀指數(shù)計算法, 對萊州灣濱海濕地景觀斑塊類型水平進行了相關(guān)分析。景觀格局指數(shù)的計算方法見文獻[28-30]。

表1 萊州灣景觀類型及其解譯標志Tab. 1 Landscape types and interpretation signs of the Laizhou Bay

2 結(jié)果與分析

2.1 萊州灣濱海濕地類型及面積變化分析

如表 2 所示, 萊州灣濱海濕地面積在 1990—2018 年間呈上升趨勢, 28 年間面積約增加240.6 km2,平均每年增加約8.6 km2。其中天然濕地面積總體減少約80.79 km2, 人工濕地總體增加約321 km2。

表2 1990—2018 年萊州灣濱海濕地面積變化Tab. 2 Area changes in the coastal wetlands of the Laizhou Bay from 1990 to 2018

天然濕地中濕地植被呈現(xiàn)先減少后增加的趨勢,濕地植被主要分布在黃河三角洲保護區(qū)、各河流兩岸及入?？谔幩鶝_積出的灘涂上(圖2—圖5)。水域面積不斷增加, 2018 年約占天然濕地面積的36%, 水域和濕地植被為天然濕地中的主要類型(表3)。灘涂由于受到河口沖淤、年降水量變化和圍填?；顒拥挠绊? 不同年份的變化量較大, 至2018 年面積約占天然濕地的17%。

人工濕地中鹽田增加面積較大, 主要增加階段為2000—2009 年, 2009 年之后增加速度放緩, 至2018 年, 鹽田約占濱海濕地面積的43%。養(yǎng)殖池在2000 年之后面積基本趨于穩(wěn)定, 至2018 年, 養(yǎng)殖池面積約占濱海濕地面積的15%。

非濕地中, 農(nóng)田自1990 年起不斷向其余景觀類型轉(zhuǎn)化, 至2018 年, 農(nóng)田面積約為2 410 km2。建設(shè)用地在2009 年之后面積陡增。

2.2 遙感影像分類結(jié)果精度評價

利用2016—2017 年100 個現(xiàn)場踏勘站位數(shù)據(jù)信息, 與2018 年遙感影響解譯結(jié)果對比, 建立混淆矩陣對結(jié)果進行精度驗證。驗證結(jié)果顯示: 2018 年萊州灣遙感影像目視解譯總體精度為85%(表4), 滿足中分辨率遙感影像精度使用要求[31-32]。

圖2 1990 年萊州灣濱海濕地景觀類型分布Fig. 2 Landscape distribution of the coastal wetlands of the Laizhou Bay in 1990

2.3 景觀類型之間面積轉(zhuǎn)移情況

圖3 2000 年萊州灣濱海濕地景觀類型分布Fig. 3 Landscape distribution of the coastal wetlands of the Laizhou Bay in 2000

圖4 2009 年萊州灣濱海濕地景觀類型分布Fig. 4 Landscape distribution of the coastal wetlands of the Laizhou Bay in 2009

圖5 2018 年萊州灣濱海濕地景觀類型分布Fig. 5 Landscape distribution of the coastal wetlands of the Laizhou Bay in 2018

表3 1990—2018 年萊州灣濱海濕地景觀類型及面積變化Tab. 3 Landscape types and area changes in the coastal wetlands of the Laizhou Bay from 1990 to 2018

由圖2—圖5 可知, 萊州灣天然濱海濕地中, 10 年間灘涂向農(nóng)田和裸地轉(zhuǎn)化最多, 轉(zhuǎn)化面積分別達到約34 km2和29 km2(表5)。濕地植被主要分布在自然保護區(qū)和入海河流的兩岸, 由于鹵水開發(fā)和海水養(yǎng)殖的興起及大規(guī)模開發(fā), 濕地植被向養(yǎng)殖池和裸地轉(zhuǎn)化明顯。由于靠內(nèi)陸部分灌溉農(nóng)田的緣故, 約22 km2的農(nóng)田轉(zhuǎn)化為水域。在萊州灣南岸約30 km2的鹽田轉(zhuǎn)化為裸地, 用于建設(shè)用地的開發(fā)利用。此時間段內(nèi)人們的主要經(jīng)濟來源依然是農(nóng)業(yè)種植, 且養(yǎng)殖業(yè)和制鹽業(yè)發(fā)展迅速。

2000—2009 年較前一階段, 灘涂在萊州灣西岸向養(yǎng)殖池轉(zhuǎn)化和南岸部分向鹽田轉(zhuǎn)化明顯, 轉(zhuǎn)化面積分別達到約47 km2和37 km2(表6)。養(yǎng)殖池向鹽田轉(zhuǎn)化顯著, 主要是萊州灣南岸包括萊州市和昌邑市的北部地區(qū)等轉(zhuǎn)移較多, 總轉(zhuǎn)化面積約為122 km2。鹽田向建設(shè)用地轉(zhuǎn)化約76 km2。濕地植被除向裸地和建設(shè)用地轉(zhuǎn)化外, 向其余幾種景觀類型轉(zhuǎn)化并不明顯, 可以看出繼1990—2000 年階段后, 濕地植被仍被持續(xù)的開發(fā)利用為人工濕地或非濕地類型。農(nóng)田向建設(shè)用地轉(zhuǎn)化最多, 轉(zhuǎn)化面積約602 km2。

這一時期, 人類活動多種多樣, 工業(yè)化進程的加快促使人類生產(chǎn)生活重心由農(nóng)業(yè)向制鹽業(yè)以及海水養(yǎng)殖轉(zhuǎn)化, 濱海濕地生態(tài)環(huán)境也面臨較大壓力。

表4 2018 年萊州灣遙感影像目視解譯與野外驗證混淆矩陣Tab. 4 Visual interpretation and field-verification confusion matrix of the Laizhou Bay remote sensing image of 2018

表5 1990—2000 年萊州灣濱海濕地類型面積轉(zhuǎn)移矩陣Tab. 5 Land-use transfer matrix of coastal wetlands in the Laizhou Bay from 1990 to 2000

表6 2000—2009 年萊州灣濱海濕地土地利用轉(zhuǎn)移矩陣Tab. 6 Land-use transfer matrix of Laizhou Bay coastal wetlands from 2000 to 2009

2009—2018 年較上一階段而言, 主要轉(zhuǎn)移類型發(fā)生了改變。期間, 養(yǎng)殖池、鹽田轉(zhuǎn)出面積大于轉(zhuǎn)入面積,其中分別向建設(shè)用地轉(zhuǎn)化約12 km2和98 km2(表7),農(nóng)田向建設(shè)用地轉(zhuǎn)化最多, 轉(zhuǎn)出面積約226 km2。約26 km2的水域、34 km2的農(nóng)田和20 km2的灘涂向濕地植被轉(zhuǎn)化, 而濕地植被本身又向養(yǎng)殖池、農(nóng)田、水域和建設(shè)用地轉(zhuǎn)化, 總體轉(zhuǎn)出面積小于轉(zhuǎn)入面積。濕地植被向灘涂轉(zhuǎn)化約11 km2, 同時灘涂也向其余類型轉(zhuǎn)化, 其中向水域和濕地植被轉(zhuǎn)化最多, 轉(zhuǎn)化面積約31.35 km2和20.31 km2, 灘涂總體轉(zhuǎn)出面積大于轉(zhuǎn)入面積。

表7 2009—2018 年萊州灣濱海濕地土地利用轉(zhuǎn)移矩陣Tab. 7 Land-use transfer matrix of the Laizhou Bay coastal wetlands from 2009 to 2018

2.4 景觀破碎化分析

斑塊數(shù)量值得多少在反映景觀異質(zhì)性的同時,與破碎度也呈現(xiàn)了一定的相關(guān)性。各景觀類型中除山體植被斑塊數(shù)量基本不變之外, 其余景觀類型均呈現(xiàn)不同程度的增加(圖6), 且在2000—2009 年間達到最大增長速度。建設(shè)用地的斑塊數(shù)遠遠高于其他景觀類型(表8), 表明建設(shè)用地的景觀破碎度較高。鹽田和農(nóng)田所占研究區(qū)面積較大而斑塊數(shù)量并不是很多, 說明其中存在著大量的優(yōu)勢斑塊, 這和圖3 中的LPI 的分布情況以及濱海濕地的實際情況也比較吻合。

景觀形狀指數(shù)可以表示景觀破碎化的程度, 研究區(qū)內(nèi)各景觀類型LSI 均有不同程度的增長, 但養(yǎng)殖池、鹽田及山體植被變化幅度較小, 尤其是山體植被基本不變, 濕地植被、農(nóng)田和建設(shè)用地變化較為劇烈, 景觀形狀指數(shù)增長幅度在40%以上。表明研究區(qū)內(nèi)空間異質(zhì)性和時間異質(zhì)性加強，隨著人類活動強度增加的影響其破碎化程度加大。2009 年以后,景觀形狀指數(shù)均有下降可能與建設(shè)用地的集聚效應(yīng)有關(guān)。最大斑塊指數(shù)可以簡單地衡量景觀優(yōu)勢度,農(nóng)田最大斑塊指數(shù)由1990 年約49%下降到2018 年的23%, 表明農(nóng)田對萊州灣濱海濕地的組成影響程度降低, 建設(shè)用地LPI 增加較為明顯表明城鎮(zhèn)建設(shè)用地趨于集中化的變化趨勢, 這與景觀形狀指數(shù)的結(jié)果吻合, 建設(shè)用地的PLAND 也從不足10%增長到25%左右。

2.5 景觀格局變化原因分析

研究區(qū)景觀格局變化是自然環(huán)境和人類活動兩種因素共同干預(yù)的結(jié)果。自然環(huán)境對景觀類型轉(zhuǎn)變的影響是長期并且是較大范圍內(nèi)的，人類活動會使景觀類型的轉(zhuǎn)變更為徹底, 會使自然濕地向人工濕地和非濕地轉(zhuǎn)化, 然而研究區(qū)內(nèi)整體大范圍的濕地格局變化并不大, 主要是小范圍內(nèi)的景觀格局的改變, 因此, 人類活動對萊州灣濱海濕地的影響較大,為主要驅(qū)動力。下面分別對在萊州灣濱海濕地的西岸、南岸和東岸各取的部分典型變化區(qū)域進行分析。

圖6 萊州灣濱海濕地不同景觀類型的景觀指數(shù)對比Fig. 6 The comparison of landscape index of the different landscapes of coastal wetlands in the Laizhou Bay

表8 萊州灣濱海濕地1990—2018 年景觀類型指數(shù)表Tab. 8 Landscape-type index tables of the Laizhou Bay coastal wetlands from 1990 to 2018

黃河口及其他河流入?？谔幱捎诤涌跊_淤使灘涂面積不斷增長, 但灘涂總體呈現(xiàn)減少趨勢。1996 年黃河口改道是萊州灣西岸黃河口附近景觀格局改變的自然驅(qū)動力之一[33], 改道之后, 黃河故道處來水來沙條件改變, 來沙量減少, 海岸侵蝕加快, 另外受氣候變化影響, 黃河入海泥沙量呈階梯狀減少[34], 致使黃河故道處海岸線受侵蝕嚴重, 新黃河口處沖淤能力也逐步減弱。萊州灣西岸和南岸多為淤泥質(zhì)海岸, 風暴潮也會對灘涂和堤壩造成一定的侵蝕[35]。此外, 自然干旱、水庫引水供灌溉和城鎮(zhèn)供水等造成許多河流時空分布不均勻也對黃河口附近景觀格局改變產(chǎn)生了一定程度的驅(qū)動作用[36]。1990—2000 年灘涂和濕地植被主要轉(zhuǎn)化區(qū)域為研究區(qū)北部黃河沿岸, 向水田(本文中水田也劃為了農(nóng)田類型)和未利用地轉(zhuǎn)化, 轉(zhuǎn)化面積均超過60 km2。老黃河口南部向陸一側(cè)則主要是養(yǎng)殖池和鹽田等人工濕地增加迅速, 2000—2009 年間灘涂、農(nóng)田和裸地分別向人工濕地轉(zhuǎn)化約84、118、140 km2。從圖7 可以看出2000 年之前該區(qū)域以灘涂為主要景觀類型, 2009 年時已被規(guī)劃成養(yǎng)殖池以及少部分的鹽田,至2018 年, 已規(guī)劃但未投入使用的土地也投入使用, 且各類型之間界限分明。萊州灣西岸人工濕地中養(yǎng)殖池面積增加最多, 而養(yǎng)殖池總體面積增長不大而鹽田增加迅速的主要原因是萊州灣南岸的養(yǎng)殖池和其他土地利用類型向鹽田轉(zhuǎn)化較多。除此之外, 降水量[37]和人為利用的雙重影響以及近年來濕地保護政策的實施等原因都會引起萊州灣水域、濕地植被等面積的變化。

圖7 1990—2018 年萊州灣濱海濕地西岸部分岸段景觀變化Fig. 7 Landscape changes in the west bank of the Laizhou Bay coastal wetlands from 1990 to 2018

萊州灣南岸部分自1991 年山東省政府提出建設(shè)“海上山東”以來, 陸續(xù)實施“耕海牧漁”戰(zhàn)略, 灘涂養(yǎng)殖發(fā)展迅速。另外從鹵水中提取溴素所帶來的豐厚經(jīng)濟利益極大的刺激了對萊州灣南岸鹵水的開采,導(dǎo)致濱海濕地天然灘涂大面積消失轉(zhuǎn)為貯存提溴的鹵水的鹽池[38]。從圖8 中可以看出1990—2009 年在沿海灘涂建設(shè)了大量的鹽場和養(yǎng)殖池, 本文統(tǒng)計鹽田面積從1990年的476.2 km2到2009年增加至825.6 km2,增幅將近43%, 鹽田面積增加位置主要位于萊州灣南岸, 從斑塊數(shù)量上也可以看出該時間段內(nèi)養(yǎng)殖池和鹽田NP 增加明顯。根據(jù)中國海洋統(tǒng)計年鑒, 2010 年以來山東省海上鹽田和養(yǎng)殖池面積處于基本穩(wěn)定狀態(tài),增長速率較前兩個階段顯著下降, 表明萊州灣作為山東省的養(yǎng)殖和產(chǎn)鹽基地, 其總體已趨于飽和, 而大量的建設(shè)用地需求逐步取代已廢棄的鹽田和養(yǎng)殖池。從表8 可知, 鹽田面積從2009—2018 年增加了23.1 km2, 與前兩個階段相比增速大大下降, 而養(yǎng)殖池面積則減少了54.6 km2, 這與統(tǒng)計年鑒趨勢基本相符。“濰坊海上新城”集約用海區(qū)的建設(shè)是萊州灣南岸建設(shè)用地增加的主要原因, 從萊州灣景觀格局分布(圖2—圖5)和圖8 中2018 年部分景觀分布中可以看到, 白浪河兩岸的農(nóng)田大面積的轉(zhuǎn)化為建設(shè)用地, 白浪河入?？趦砂稙┩恳约敖逗Ｓ虮患谐善_發(fā)為建設(shè)用地和待利用地, 2009—2018 年鹽田和農(nóng)田分別向建設(shè)用地轉(zhuǎn)化98 km2和226 km2, 其中鹽田轉(zhuǎn)化部分基本全處于萊州灣南岸, 可見集約用海區(qū)建設(shè)對該地區(qū)的景觀格局產(chǎn)生了較大的影響。另外, 人口和社會經(jīng)濟發(fā)展作為持續(xù)不斷的外界壓力, 使得耕地及未利用土地不斷向工業(yè)區(qū)及城鎮(zhèn)用地轉(zhuǎn)化, 道路、水利設(shè)施(溝渠)、居住地以及工業(yè)區(qū)的建設(shè)使得萊州灣南岸破碎化加重。造成萊州灣南岸濱海濕地景觀格局變化的主要自然驅(qū)動力為海水入侵[39], 在其影響下, 南岸景觀類型由海向陸大致為灘涂-養(yǎng)殖池-鹽田-濕地植被-建設(shè)用地-農(nóng)田, 且海水入侵對工農(nóng)業(yè)生產(chǎn), 生態(tài)環(huán)境以及植被演替等方面有很大影響[13]。

圖8 1990—2018 年萊州灣濱海濕地南岸部分岸段景觀變化Fig. 8 Landscape changes in the south bank of the Laizhou Bay coastal wetlands from 1990 to 2018

圍海養(yǎng)殖以及填海造陸是萊州灣東岸沿岸景觀格局變化的主要人工驅(qū)動力, 從圖2—圖5 的景觀分布圖中可以看到圍海養(yǎng)殖主要集中在萊州市北部沿海, 從1993 年的海灣扇貝養(yǎng)殖業(yè)到2006 年之后“溫室大棚+深井海水”工廠化養(yǎng)殖和池塘養(yǎng)殖模式[40],養(yǎng)殖業(yè)利用了大面積的潮上帶土地; 填海造地主要發(fā)生在龍口灣附近, 尤其是2000 年之后, 填海造陸的規(guī)模不斷擴大, 2010 年龍口人工島群工程區(qū)域用海規(guī)劃獲得國家海洋局批復(fù), 人工島群的建設(shè)對于陸海統(tǒng)籌、集約用海等推進海洋文明建設(shè)有著極大的推動作用。集約用海區(qū)是在國家統(tǒng)一規(guī)劃下, 對同一區(qū)域內(nèi)多個圍填海項目集中開發(fā)的用海方式。集中集約用海對推進海陸統(tǒng)籌區(qū)域經(jīng)濟一體化起到了顯著作用, 是在“藍色”環(huán)保理念之下的一種用海新模式。2009—2018 年, “龍口灣海洋裝備制造業(yè)集聚區(qū)”和“萊州海洋新能源產(chǎn)業(yè)集聚區(qū)”等集約用海區(qū)的規(guī)劃開發(fā), 使生產(chǎn)活動的重心向制造業(yè)等偏移。萊州港和龍口港對萊州灣東岸的經(jīng)濟發(fā)展有著舉足輕重的作用, 港城經(jīng)濟帶動港口發(fā)展[41], 而近年來港口規(guī)模的擴大和集疏運體系的完善也反過來推動萊州市和龍口市經(jīng)濟快速發(fā)展。另外城鎮(zhèn)化規(guī)模的擴大是萊州灣東岸景觀格局變化的另一驅(qū)動力, 從不同年份之間景觀格局分布(圖2—圖5)可以看出2000—2009 年萊州市城區(qū)、龍口東、西城區(qū)城鎮(zhèn)規(guī)模增長較多, 該階段農(nóng)田向建設(shè)用地轉(zhuǎn)移602 km2, 除少部分轉(zhuǎn)移位于東營區(qū)、寒亭區(qū)和壽光市外, 其余大部分位于萊州灣東岸。從圖9 可以看出1990—2018 年龍口港及西城區(qū)規(guī)模不斷擴大, 原先零散分布的建設(shè)及居住用地逐漸趨于整合, 這與本文結(jié)果建設(shè)用地LPI 值增加的趨勢也吻合。2018 年建設(shè)用地的最大斑塊指數(shù)達到4.26%, 表明人類活動呈現(xiàn)出集聚效應(yīng)。

圖9 1990—2018 年萊州灣濱海濕地東岸龍口港景觀變化Fig. 9 Landscape changes in the east bank of the Laizhou Bay coastal wetlands from 1990 to 2018

1990—2009 年間景觀指數(shù)的增加表明萊州灣濱海濕地景觀破碎化加重, 景觀破碎化的增加主要是由于城市建設(shè)和鹽田的發(fā)展, 破壞了景觀斑塊的完整性。2009 年之后各項景觀指數(shù)不再增加且呈略微下降的趨勢, 景觀破碎趨勢變緩可能跟城市化進程加快、集約用海區(qū)建設(shè)工業(yè)區(qū)、港口以及城鎮(zhèn)居住地逐漸聯(lián)合為一體有關(guān)。隨著山東半島藍色經(jīng)濟區(qū)規(guī)劃的出臺, 4 個集中集約用海區(qū)位于萊州灣內(nèi)，待其建成后，可以預(yù)見其能長期影響萊州灣附近海岸及其濱海濕地景觀類型之間的相互轉(zhuǎn)化。

3 結(jié)論

本研究以遙感圖像為基礎(chǔ)，通過目視解譯等手段得到不同時期下萊州灣濱海濕地景觀格局分布,并通過景觀指數(shù)對研究區(qū)進行分析。有關(guān)結(jié)論如下:截至2018 年萊州灣濱海濕地面積總計約1 955 km2,其中天然濱海濕地面積約811 km2, 人工濱海濕地面積約1 143 km2。其近年來主要呈現(xiàn)天然濱海濕地減少，人工濱海濕地增加的總體趨勢。自2009 年起，在人工濱海濕地中養(yǎng)殖池和鹽田面積均接近飽和，且廢棄的養(yǎng)殖池和鹽田逐漸向建設(shè)用地轉(zhuǎn)化, 農(nóng)田向建設(shè)用地轉(zhuǎn)化顯著。萊州灣濱海濕地景觀格局變化的自然驅(qū)動力主要為河口沖淤、海水入侵、氣候變化等; 人工驅(qū)動力主要為圍海養(yǎng)殖、鹽田開發(fā)、城鎮(zhèn)化建設(shè)以及集約用海區(qū)的規(guī)劃開發(fā)等。1990—2018 年間萊州灣濱海濕地景觀破碎化加重, 景觀破碎化的增加主要是由于城市建設(shè)、港口碼頭、養(yǎng)殖業(yè)和鹽業(yè)的發(fā)展, 破壞了景觀斑塊的完整性,2009 年之后由于人口和建設(shè)用地的集聚效應(yīng), 破碎趨勢有所變緩。