張華玉, 秦年秀, 汪軍能, 李 敏, 楊愛霞, 盧芹芹

(1.北部灣環(huán)境演變與資源利用教育部重點(diǎn)實(shí)驗(yàn)室, 南寧 530001; 2.南寧師范大學(xué) 地理科學(xué)與規(guī)劃學(xué)院,南寧 530001; 3.南寧師范大學(xué) 自然資源與測繪學(xué)院, 南寧 530001; 4.北部灣大學(xué) 資源與環(huán)境學(xué)院, 廣西 欽州 535000)

海岸帶是陸海交匯、人類活動與自然系統(tǒng)交互作用最為頻繁、密集的地理空間之一[1]，資源豐富且生態(tài)系統(tǒng)結(jié)構(gòu)復(fù)雜、生態(tài)環(huán)境脆弱和敏感。在人類活動疊加自然變化的影響下，海岸帶地區(qū)環(huán)境和土地利用結(jié)構(gòu)發(fā)生了巨大變化，備受學(xué)術(shù)界關(guān)注，地理信息技術(shù)已廣泛運(yùn)用于海岸帶和區(qū)域土地利用變化研究中[2]。國內(nèi)外海岸帶土地利用研究主要集中在海岸帶岸線變化機(jī)理與用地變化效應(yīng)、土地利用/覆被與景觀格局演變及其驅(qū)動機(jī)制、景觀生態(tài)脆弱性與生態(tài)風(fēng)險(xiǎn)評價(jià)、生態(tài)系統(tǒng)服務(wù)價(jià)值等方面。Lam-Dao，Dada等[3-4]通過遙感手段分別對河口、海岸地區(qū)土地利用變化研究，發(fā)現(xiàn)海岸帶遭到不同程度侵蝕，養(yǎng)殖業(yè)發(fā)展和氣候變化對海岸帶土地利用影響較大。Bernard Ekumah等[5]通過強(qiáng)度分析方法評估了加納3個(gè)具有國際重要性的沿海城市濕地變化模式及其原因。Jayanta Saha等[6]分析了印度孟加拉沿海地區(qū)旅游業(yè)發(fā)展引致的土地利用變化及其環(huán)境效應(yīng)。吳琛璐[7]、高祥偉[8]等基于Landsat影像和3S技術(shù)分別對省、市尺度海岸帶展開研究，發(fā)現(xiàn)在社會經(jīng)濟(jì)因素驅(qū)動下，土地利用變化顯著。潘翔等[9]基于遙感和GIS技術(shù)及景觀格局指數(shù)法分析了廈門灣海岸帶景觀時(shí)空格局，發(fā)現(xiàn)研究期間農(nóng)用地景觀、城鎮(zhèn)用地與建設(shè)用地景觀變化顯著，農(nóng)用地和濱海濕地景觀轉(zhuǎn)出較多，受社會經(jīng)濟(jì)發(fā)展和政策調(diào)整影響，景觀格局向多樣化和破碎化演變。王奇等[10]運(yùn)用“壓力度—敏感度—恢復(fù)力”評價(jià)模型分析發(fā)現(xiàn)甬臺溫海岸帶景觀類型變化顯著，生態(tài)脆弱性增強(qiáng)。陳心怡等[11]對?？谑泻０稁?988—2017年景觀生態(tài)風(fēng)險(xiǎn)進(jìn)行評價(jià)，結(jié)果呈生態(tài)風(fēng)險(xiǎn)等級由低向高，高風(fēng)險(xiǎn)區(qū)向沿海蔓延的特征。寧立新等[12]對江蘇海岸帶地區(qū)的研究發(fā)現(xiàn)耕地減少、建設(shè)用地和人工塘增加，生態(tài)系統(tǒng)服務(wù)價(jià)值加速減少?？傮w上看，海岸帶土地利用研究體系不斷深化，方法也較為成熟。國外研究多從大尺度、國際重要海岸帶而開展[13-14]，中小尺度研究較少；國內(nèi)研究則多傾向于東部沿海地區(qū)，北部灣地區(qū)關(guān)注度有待提升。驅(qū)動機(jī)制分析中多采用社會經(jīng)濟(jì)指標(biāo)及少數(shù)易量化的自然因子指標(biāo)進(jìn)行統(tǒng)計(jì)分析[15-16]，而海拔、坡度、坡向、植被類型和土壤類型等對土地利用結(jié)構(gòu)及分布影響較大，且難以量化和標(biāo)準(zhǔn)化的自然因子，則未能進(jìn)行較好考量，以至于難以全面剖析土地利用變化的驅(qū)動機(jī)制。如何更全面系統(tǒng)探討海岸帶土地利用變化及其驅(qū)動機(jī)制，是當(dāng)前海岸帶土地利用變化研究遇到的新難題。

2008年廣西北部灣經(jīng)濟(jì)區(qū)發(fā)展規(guī)劃批復(fù)、2010年中國—東盟自由貿(mào)易區(qū)(CAFTA)全面啟動、2017年北部灣城市群規(guī)劃批復(fù)等一系列政策引導(dǎo)下，廣西海岸帶地區(qū)進(jìn)入快速發(fā)展階段，經(jīng)濟(jì)快速增長，人口不斷集聚，土地利用變化也進(jìn)入新階段。用于旅游、工業(yè)、港口和城鎮(zhèn)建設(shè)等填海面積不斷增加[17]，岸線縮短[18]，土地可持續(xù)利用與生態(tài)保護(hù)面臨新挑戰(zhàn)[19-20]。快速發(fā)展背景下迫切需要對廣西海岸帶土地利用變化與景觀格局演變特征及土地利用變化的驅(qū)動機(jī)制進(jìn)行全面、系統(tǒng)、綜合地分析研究。

研究采用2010年、2014年、2018年三期遙感影像解譯數(shù)據(jù)，運(yùn)用地理學(xué)和數(shù)理統(tǒng)計(jì)方法分析廣西海岸帶土地利用與景觀格局演化規(guī)律，并運(yùn)用地理探測器方法探究其背后的驅(qū)動機(jī)制。揭示人文與自然因子雙重作用下的海岸帶土地利用變化過程，為破解廣西沿海地區(qū)土地退化、生境功能減弱及人地矛盾尖銳等難題，合理開展國土資源利用與空間規(guī)劃服務(wù)，為優(yōu)化廣西海岸帶土地利用結(jié)構(gòu)、科學(xué)布局產(chǎn)業(yè)、科學(xué)利用和管理自然資源提供依據(jù)。

1 數(shù)據(jù)和方法

1.1 研究區(qū)概況

本研究包括沿廣西大陸海岸線的陸域與海域(圖1)。研究區(qū)位于我國西南部、廣西東南部(20°26′N—22°41′N，107°28′E—109°56′E)，自東向西由洗米河口至北侖河口，海岸線長1 595 km[21]；陸域包括廣西防城港、欽州、北海3個(gè)沿海城市及其下轄的各縣(市)，海域部分以海岸線向海延伸10 km作為緩沖區(qū)域[22]，面積約23 010 km2。廣西海岸帶地處亞熱帶海洋性季風(fēng)氣候區(qū)，雨熱同期，年平均氣溫在22℃以上，年降水量達(dá)2 100～2 760 mm，風(fēng)能資源豐富，霧少光足，夏季洪澇災(zāi)害頻發(fā)，是全區(qū)洪澇災(zāi)害高風(fēng)險(xiǎn)區(qū)[23]；土壤類型主要有磚紅壤、水稻土等；地勢北高南低，北部是以十萬大山為主的山地丘陵地貌，南部平原灘涂面積廣，是南流江、欽江、茅嶺江等河流入?？?；海岸曲折，紅樹林資源豐富，港灣優(yōu)良，養(yǎng)殖業(yè)發(fā)達(dá)。2019年廣西沿海三市GDP為3 358.30億元，常住人口約598.84萬人。廣西海岸帶自然資源豐富、生態(tài)系統(tǒng)復(fù)雜、人地互動高度活躍且區(qū)位獨(dú)特，是我國西南地區(qū)最便捷的出海通道，也是面向東盟的國際大通道、21世紀(jì)海上絲綢之路和絲綢之路經(jīng)濟(jì)帶有機(jī)銜接的重要門戶[24]。

圖1 研究區(qū)位示意圖

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

采用2010年、2014年、2018年三期共16景，空間分辨率為30 m的Landsat5 TM和Landsat8 OLI遙感影像。影像來源于美國國家航空航天局(https:∥www.nasa.gov/)；用ENVI 5.5軟件對影像進(jìn)行輻射定標(biāo)、大氣校正、波段運(yùn)算、影像鑲嵌、影像裁剪等預(yù)處理后，在432，764，652三組波段組合下，以Google Earth判讀和實(shí)地調(diào)查數(shù)據(jù)為參照，通過人機(jī)交互方法進(jìn)行解譯，三期遙感影像解譯總體精度和Kappa系數(shù)分別為98.08%，0.86，99.42%，0.95，98.51%，0.89?？傮w精度均在95%以上，Kappa系數(shù)均超過0.85，達(dá)到研究要求。

常住人口(X1)、人均GDP(X2)、第一產(chǎn)業(yè)增加值(X3)、第二產(chǎn)業(yè)增加值(X4)、第三產(chǎn)業(yè)增加值(X5)、城鎮(zhèn)化率(X6)、GDP(X7)和社會固定資產(chǎn)投資(X8)8個(gè)社會經(jīng)濟(jì)因子數(shù)據(jù)源自《廣西統(tǒng)計(jì)年鑒》(2011—2019年)、區(qū)內(nèi)縣(區(qū))年鑒及統(tǒng)計(jì)公報(bào)數(shù)據(jù)。

年均降水量(X9)、氣溫(X10)、濕度(X11)和蒸發(fā)量(X12)數(shù)據(jù)均源于中國氣象數(shù)據(jù)網(wǎng)；海拔(X13)資料源于91衛(wèi)圖助手，用ArcMap對海拔數(shù)據(jù)進(jìn)行表面分析得到坡度(X14)和坡向(X15)數(shù)據(jù)；土壤類型(X16)和植被類型(X17)數(shù)據(jù)源于中科院資源環(huán)境科學(xué)與數(shù)據(jù)中心(https:∥www.resdc.cn/)。文中對所有因子數(shù)據(jù)進(jìn)行了柵格化及離散化處理。

1.3 研究方法

隨機(jī)森林(random forest，RF)分類算法。該方法具有人工干預(yù)少，分類精度高等優(yōu)勢，在處理復(fù)雜數(shù)據(jù)源、數(shù)據(jù)噪聲及有限訓(xùn)練樣本方面有較好的表現(xiàn)。本文結(jié)合已有研究[25-26]和廣西海岸帶土地利用特征，構(gòu)建耕地、林草地、建設(shè)用地、水體、人工濕地、自然濕地及其他用地7個(gè)類型的廣西海岸帶土地利用分類體系，運(yùn)用ENVI 5.5軟件對3期遙感影像進(jìn)行隨機(jī)森林分類，并對分類結(jié)果進(jìn)行人工修正。

土地利用動態(tài)度和馬爾科夫模型[27]。土地利用動態(tài)度表征某一地區(qū)在特定時(shí)段內(nèi)的土地利用變化速率，文中采用單一土地利用動態(tài)度和綜合土地利用動態(tài)度來描述廣西海岸帶3個(gè)時(shí)期土地利用變化狀況。通過馬爾科夫模型構(gòu)建土地利用轉(zhuǎn)移矩陣，可直觀表達(dá)研究區(qū)內(nèi)土地利用類型之間的量化轉(zhuǎn)移情況，從而分析一定時(shí)期內(nèi)土地利用類型的轉(zhuǎn)變趨勢。

景觀格局指數(shù)法[28]。景觀格局演變是土地利用變化影響下的結(jié)果?；趶V西海岸帶3期土地利用類型數(shù)據(jù)，從類型水平(Class-level)和景觀水平(Landscape-level)兩個(gè)尺度，選取最大斑塊指數(shù)、斑塊數(shù)量、斑塊密度、破碎度指數(shù)、景觀分割度指數(shù)、分離度指數(shù)、香農(nóng)多樣性指數(shù)、香農(nóng)均勻度指數(shù)8個(gè)景觀格局指數(shù)，借助Fragstats 4.2軟件，分析廣西海岸帶景觀格局動態(tài)變化狀況。

地理探測器方法[29]，通過因子探測器模塊和交互作用探測器模塊識別單因子影響力及各驅(qū)動因子之間的交互作用。以土地利用程度綜合指數(shù)[30]作為因變量，公式如下：

(1)

式中:La為土地利用程度綜合指數(shù);Ai表示第i級土地利用程度分級指數(shù)(表1);Ci表示第i級土地利用程度分級面積百分比。

2 結(jié)果與分析

2.1 土地利用動態(tài)變化特征

2.1.1 時(shí)間變化特征 2010—2018年廣西海岸帶土地利用變化顯著。土地利用類型中優(yōu)勢地類是林草地和耕地(圖2)，人工濕地、自然濕地、其他用地占比較低。2010—2018年耕地減少了1 301.28 km2，年均減少162.66 km2，且減少趨勢在加快；建設(shè)用地和林草地呈增加態(tài)勢，建設(shè)用地由1 029.09 km2(4.46%)增加至1 844.24 km2(7.99%)，增加了約815.15 km2，年均增加101.89 km2，增速最快；林草地次之，由8 792.22 km2增加至9 391 km2，年均增加74.85 km2；人工濕地和自然濕地小幅增加，分別增加了7.69 km2，7.24 km2。水體縮減34.73 km2；其他用地減少92.85 km2。說明這一時(shí)期廣西海岸帶人類活動較為劇烈，城市化加速和經(jīng)濟(jì)快速發(fā)展，城鄉(xiāng)建設(shè)用地不斷擴(kuò)張，耕地遭侵蝕較為嚴(yán)重，保護(hù)力度有待加強(qiáng)；林草資源發(fā)展與保護(hù)較好，尤其是經(jīng)濟(jì)林發(fā)展不可小覷，紅樹林保護(hù)得到加強(qiáng)；除季節(jié)性降水影響外，填海造陸、水產(chǎn)養(yǎng)殖及過度用水等活動正侵蝕著水體。

表1 土地利用程度分級賦值

圖2 廣西海岸帶土地利用類型面積及變化

2.1.2 空間變化特征 廣西海岸帶土地利用變化空間差異較為明顯(圖3)。耕地集中分布在合浦縣、上思縣、欽北區(qū)和欽南區(qū)，2010—2018年靈山縣耕地縮減面積最大(373.77 km2)，耕地縮減處伴隨著林草地和建設(shè)用地?cái)U(kuò)張。建設(shè)用地主要分布在合浦縣、靈山縣、欽南區(qū)等行政區(qū)劃較大的區(qū)域及處于城市發(fā)展核心地帶的海城區(qū)、港口區(qū)、東興市等，2010—2018年靈山縣、浦北縣、銀海區(qū)及鐵山港區(qū)擴(kuò)張明顯。上思縣、靈山縣、浦北縣是林草地集中分布區(qū)，2010—2018年欽南區(qū)、合浦縣、鐵山港區(qū)林草地?cái)U(kuò)張面積大，其中欽南區(qū)林草地?cái)U(kuò)張面積最大(180.72 km2)，浦北縣林草地縮減最多(77.90 km2)。水體和自然濕地主要分布在緩沖區(qū)、合浦縣、欽南區(qū)和上思縣，2010—2018年緩沖區(qū)水體縮減面積最大(45.73 km2)，合浦縣水體擴(kuò)張面積最大(17.02 km2)。其他地類空間變化不明顯。

圖3 廣西海岸帶土地利用結(jié)構(gòu)分區(qū)

2.1.3 動態(tài)度變化特征 從土地利用單一動態(tài)度上看(圖4)，變化速度依次為建設(shè)用地>其他用地>耕地>自然濕地>林草地>人工濕地>水體，2010—2018年建設(shè)用地增幅最大，動態(tài)度為8.80%；2010—2014年的其他用地縮減幅度最大，動態(tài)度為-10.48%；2010—2018年水體減幅最小，動態(tài)度為-0.11%。表明該區(qū)域?qū)ζ渌玫剡M(jìn)行了大量開發(fā)和開荒造林活動，從而加速了其他用地的喪失，而建設(shè)項(xiàng)目劇增及城鄉(xiāng)發(fā)展需要，則導(dǎo)致建設(shè)用地急劇變化。土地利用綜合動態(tài)度隨時(shí)間變化而上升，2010—2014年、2014—2018年和2010—2018年3個(gè)時(shí)段的綜合動態(tài)度為分別為1.12%，1.37%和1.52%，總體屬土地利用快速變化型[31]，并有向土地利用急劇變化型演變趨勢。

圖4 廣西海岸帶土地利用動態(tài)度變化

2.1.4 土地利用轉(zhuǎn)移特征 通過構(gòu)建土地利用轉(zhuǎn)移矩陣分析，發(fā)現(xiàn)2010—2018年廣西海岸帶各土地利用類型間轉(zhuǎn)變較大，尤其是林草地、建設(shè)用地和耕地(表2)。耕地轉(zhuǎn)變最顯著，主要轉(zhuǎn)變?yōu)榱植莸?、建設(shè)用地和水體，其中轉(zhuǎn)變?yōu)榱植莸刈顬閯×?2 297.12 km2)；林草地主要轉(zhuǎn)變?yōu)楦?、建設(shè)用地、水體和其他用地，向耕地轉(zhuǎn)變的幅度最大(1 547.29 km2)；建設(shè)用地主要轉(zhuǎn)變?yōu)楦亍⒘植莸?、人工濕地和水體，向耕地轉(zhuǎn)變幅度最大(247.25 km2)；水體主要轉(zhuǎn)換為建設(shè)用地、耕地和人工濕地，向耕地轉(zhuǎn)化面積最大(92.99 km2)；其他用地主要轉(zhuǎn)變?yōu)榱植莸睾透?，轉(zhuǎn)變?yōu)榱植莸刈顬閯×?120.7 km2)；人工濕地主要轉(zhuǎn)為水體、建設(shè)用地和耕地；自然濕地轉(zhuǎn)變最不明顯，主要轉(zhuǎn)變?yōu)榱植莸睾退w。說明研究期內(nèi)耕地、林草地和建設(shè)用地3種地類間交替轉(zhuǎn)化較為頻繁和活躍，而耕地大面積轉(zhuǎn)為林草地和建設(shè)用地，存在耕地林化與建設(shè)占用現(xiàn)象，應(yīng)加強(qiáng)管控與監(jiān)測，減緩轉(zhuǎn)變趨勢；封山育林、開荒造林取得一定成效；水體漲落對人工濕地和自然濕地產(chǎn)生一定的影響。

表2 廣西海岸帶土地利用轉(zhuǎn)移矩陣 km2

2.2 景觀格局演變特征分析

2010—2018年，廣西海岸帶各景觀指數(shù)變化明顯。斑塊數(shù)量(198 504～240 731)、斑塊密度(8.602 7～10.432 7)、景觀分割度指數(shù)(0.951 8～0.956 0)、分離度指數(shù)(20.749 8～22.747 5)、香濃多樣性指數(shù)(1.285 9～1.332 9)和香濃均勻度指數(shù)(0.660 8～0.685 0)均呈波動上升的變化特征(表3)。表明廣西海岸帶地區(qū)景觀的異質(zhì)性和復(fù)雜性增加，愈加破碎和均勻。2010—2018年水體最大斑塊指數(shù)最高，且其斑塊數(shù)量、斑塊密度及破碎度指數(shù)波動下降，表明水體景觀隨時(shí)間變化而集聚性上升，斑塊片狀效果增強(qiáng)；耕地則相反，2010—2018年期間斑塊數(shù)量、斑塊密度和破碎度指數(shù)遞增，最大斑塊指數(shù)減小，表明其受干擾大，分布隨時(shí)間變化而散化；建設(shè)用地的最大斑塊指數(shù)、斑塊數(shù)量和斑塊密度增速明顯，破碎度指數(shù)略微下降，表明建設(shè)用地?cái)U(kuò)張迅速，建設(shè)用地的集聚特征更加突出(圖5)。

人工濕地、自然濕地兩種景觀的最大斑塊指數(shù)、斑塊數(shù)量及斑塊密度都增加，人工濕地破碎度指數(shù)上升，自然濕地破碎度指數(shù)下降，表明這兩種景觀出現(xiàn)一定程度擴(kuò)張，而人工濕地逐漸呈破碎化和分散化發(fā)展，自然濕地則更加均勻化和集中分布。

2.3 土地利用變化驅(qū)動因素分析

2.3.1 因子探測 對廣西海岸帶土地利用變化的各因子影響程度定量化，以3個(gè)年份的平均q值作為驅(qū)動因子貢獻(xiàn)力參考(表4)。GDP(0.202)、社會固定資產(chǎn)投資(0.195)和城鎮(zhèn)化率(0.181)等人文因子對其土地利用變化的影響較大，平均貢獻(xiàn)率為17.82%，說明廣西海岸帶土地利用變化主要受人類活動影響，人文因子在整個(gè)變化過程中起主導(dǎo)作用。自然因子的影響相對較小，平均貢獻(xiàn)率僅9.72%，其中植被類型(0.210)、坡度(0.156)、海拔(0.139)和土壤類型(0.132)四項(xiàng)因子在自然因子中貢獻(xiàn)率最高，主要體現(xiàn)在對土地類型空間分布和功能劃分的影響上；年平均降水量(0.033)、年平均氣溫(0.039)和蒸發(fā)量(0.050)等因子對其土地利用變化影響較小。

圖5 廣西海岸帶各地類景觀指數(shù)統(tǒng)計(jì)

表3 廣西海岸帶整體景觀指數(shù)統(tǒng)計(jì)

2.3.2 因子間交互作用探測 在交互探測結(jié)果中選取解釋力均在0.3以上的21組因子進(jìn)行分析。由表5可知，所選取的21組因子間交互作用類型以雙因子增強(qiáng)與非線性增強(qiáng)為主，人文因子與自然因子間交互作用較為強(qiáng)烈。人文因子與自然因子中的坡度、海拔兩個(gè)因子的交互最為強(qiáng)烈，其中GDP與坡度、第二產(chǎn)業(yè)增加值與海拔、GDP與海拔三組因子交互作用值最高，交互作用值均為0.36。說明廣西海岸帶地區(qū)土地利用變化是人文因子與自然因子協(xié)同作用的結(jié)果，人文因子占優(yōu)(平均單一因子力大)，且人文因子與自然因子交互作用強(qiáng)于人文因子間交互作用與自然因子間交互作用。自然因子的參與使得兩種因子共同影響力是各自單一影響力的將近兩倍，從而加速了廣西海岸帶土地利用變化。

3 討論與結(jié)論

3.1 討 論

(1) 人文因子在廣西海岸帶土地利用變化過程中起主導(dǎo)作用。廣西海岸帶地勢平坦，資源豐富，環(huán)境優(yōu)美，自古以來即適宜人類居住和開發(fā)建設(shè)。近年來，廣西北部灣經(jīng)濟(jì)區(qū)、中國—東盟自由貿(mào)易區(qū)、北部灣城市群的規(guī)劃、建設(shè)和新型城鎮(zhèn)化發(fā)展，使得廣西海岸帶在區(qū)域中地位快速提升，成為全區(qū)乃至全國重大項(xiàng)目建設(shè)的主要分布區(qū)。開發(fā)建設(shè)、人工填海等活動在促進(jìn)本地區(qū)經(jīng)濟(jì)發(fā)展的同時(shí)，也加大了土地開發(fā)、經(jīng)營和管理壓力，導(dǎo)致土地利用結(jié)構(gòu)和數(shù)量在一定程度上發(fā)生了變化。2010—2018年廣西海岸帶建設(shè)用地面積激增、耕地持續(xù)快速的縮減，而且還存在耕地撂荒和耕地種林的現(xiàn)象。相關(guān)部門應(yīng)該立足存量，積極開展城市更新和城鄉(xiāng)建設(shè)用地增減掛鉤等整治措施；加大宣傳力度，鼓勵(lì)居民拆舊建新；嚴(yán)格審批，有效控制建設(shè)用地有序增減。同時(shí)，應(yīng)加強(qiáng)對密集型農(nóng)業(yè)區(qū)域的監(jiān)管，加強(qiáng)耕地和永久基本農(nóng)田保護(hù)，防止違法占用及破壞耕作層導(dǎo)致質(zhì)量下降行為，嚴(yán)格落實(shí)耕地占補(bǔ)平衡機(jī)制，占優(yōu)補(bǔ)優(yōu)，確?；巨r(nóng)田數(shù)量不減少、質(zhì)量不降低。

(2) 近些年，在紅樹林生態(tài)修復(fù)與保護(hù)、退耕還林、營造生態(tài)防護(hù)林等生態(tài)文明建設(shè)措施驅(qū)動下，廣西海岸帶土地利用趨于健康發(fā)展，取得了良好的生態(tài)、經(jīng)濟(jì)和社會效益[32]。廣西山口紅樹林自然保護(hù)區(qū)、茅尾海紅樹林自然保護(hù)區(qū)及北海國家濕地公園的建設(shè)發(fā)展，提升了廣西海岸帶紅樹林資源修復(fù)與保護(hù)力度。廣西海岸帶林草地和自然濕地面積呈正增長趨勢，正成為區(qū)域可持續(xù)發(fā)展、海岸防護(hù)及綠色生態(tài)發(fā)展的強(qiáng)有力支撐。此外，大量增減掛鉤和舊房拆除等土地整治活動，退耕還林、還湖等生態(tài)建設(shè)措施，使得一些建設(shè)用地轉(zhuǎn)變?yōu)楦?、耕地變林草地或水體等。但是，廣西海岸帶的生態(tài)修復(fù)與保護(hù)任務(wù)依然艱巨，應(yīng)堅(jiān)持生態(tài)保育，加強(qiáng)監(jiān)測與審批，繼續(xù)落實(shí)林業(yè)體系建設(shè)與紅樹林保護(hù)政策，防止耕地林業(yè)化，促進(jìn)土地可持續(xù)利用。

表4 廣西海岸帶土地利用變化影響因子探測值

表5 因子交互作用q值及相應(yīng)交互關(guān)系

(3) 人類活動影響下，廣西海岸帶土地利用變化對氣候變化的影響不可忽視[33]。森林過度采伐、城市建設(shè)及工農(nóng)業(yè)生產(chǎn)等活動導(dǎo)致土地利用變化的同時(shí)，易造成溫室氣體增加、能量和水分輸送不平衡、形成局地微氣候等問題，影響區(qū)域乃至全球氣候。因此，海岸帶土地利用過程中應(yīng)注重各用地間的平衡協(xié)調(diào)和景觀格局的營建與優(yōu)化，從而一定程度緩解溫室氣體排放與氣候變化背景下海平面上升對海岸帶土地的侵蝕；同時(shí)應(yīng)加強(qiáng)對海岸帶地區(qū)極端天氣的監(jiān)控與預(yù)報(bào)，預(yù)防驟發(fā)性大規(guī)模氣象災(zāi)害等造成的不良影響。

3.2 結(jié) 論

(1) 2010—2018年廣西海岸帶土地利用變化顯著。耕地年均減少162.66 km2，建設(shè)用地年均增加101.89 km2，各地類單一動態(tài)度存在明顯差異，綜合動態(tài)度逐段上升。土地利用方式改變、經(jīng)濟(jì)社會發(fā)展等使得土地利用類型間發(fā)生不同程度轉(zhuǎn)變，林草地、耕地和建設(shè)用地間相互轉(zhuǎn)化最為劇烈，其他地類轉(zhuǎn)變幅度不大。

(2) 2010—2018年廣西海岸帶土地利用結(jié)構(gòu)不斷優(yōu)化，各地類配置趨于合理，景觀格局呈斑塊破碎化、多樣化和均勻化演變趨勢。景觀水平尺度上的各指數(shù)呈上升態(tài)勢。水體最大斑塊指數(shù)最高，斑塊數(shù)量、斑塊密度及破碎度指數(shù)波動下降，是廣西海岸帶優(yōu)勢景觀。耕地受干擾大，分散化明顯，建設(shè)用地?cái)U(kuò)張迅速且集聚特征突出。人工濕地、自然濕地呈一定程度擴(kuò)張，而自然濕地更加均勻和集中。

(3) 2010—2018年廣西海岸帶土地利用變化是人文因子與自然因子綜合作用的結(jié)果，人文因子是主導(dǎo)因子，自然因子的影響較小。因子間交互作用類型主要為雙因子增強(qiáng)和非線性增強(qiáng)兩種，人文因子與自然因子中的海拔、坡度兩個(gè)因子的交互最為強(qiáng)烈，人文因子占優(yōu)。土地利用變化抗干擾功能提升，驅(qū)動因子對土地利用變化的總體影響程度逐年下降。