無居民海島生態(tài)承載力動態(tài)評估
——以福建平潭大嶼島為例

田鵬，林麗萍，黃金良?，李湞

(1.中交上海航道勘察設(shè)計研究院有限公司，上海 200120；2.福建省海岸帶污染防控重點(diǎn)實(shí)驗(yàn)室/廈門大學(xué)，廈門 361102)

在我國海域范圍內(nèi)，面積為500m2及以上的海島約有6500 個，總面積達(dá)6600km2，其中有94%屬于無居民海島[1]。無居民海島大多面積狹小，地域相對獨(dú)立，生態(tài)系統(tǒng)結(jié)構(gòu)也比較單一，生物多樣性指數(shù)一般不高，穩(wěn)定性也較差，所以極易受到破壞[2]。同時，現(xiàn)行的海島開發(fā)工作管控力度不足，盲目性和破壞性開發(fā)情況突出，資源利用率低，生態(tài)環(huán)境問題嚴(yán)重，難以實(shí)現(xiàn)資源的可持續(xù)利用[3]。

在某個特定環(huán)境空間下，資源環(huán)境的數(shù)量與質(zhì)量以及可承受的人口與經(jīng)濟(jì)可持續(xù)發(fā)展需求的能力稱為生態(tài)承載力[4]。從生態(tài)承載力理論方面看，生態(tài)承載力的提出意義重大。目前，生態(tài)承載力的評估方法可歸納為“供需對比法”“指標(biāo)體系法”及“系統(tǒng)模型法”三大類；其中的“生態(tài)足跡法”[1，5－8]“狀態(tài)空間法”[9－11]以及“分級—綜合法”[12－15]是應(yīng)用最為廣泛的評估方法[16－17]。國外的海島生態(tài)承載力研究開展得較早，Jansson(1985)研究了瑞典Gotland 島的區(qū)域生態(tài)承載力并提出島群區(qū)域規(guī)劃[18]。池源等(2017)以廟島群島南部島群為研究區(qū)構(gòu)建了一套能夠反映陸海雙重特征及其空間分異性的海島生態(tài)系統(tǒng)承載力評估模型[19]。黃劍彬等(2017) 基于平潭島2005—2013 年土地利用數(shù)據(jù)，運(yùn)用生態(tài)足跡模型，對平潭島生態(tài)承載力進(jìn)行了分析[20]。李芬等(2016)從人—自然耦合的復(fù)合生態(tài)系統(tǒng)角度，基于經(jīng)濟(jì)社會發(fā)展與生態(tài)系統(tǒng)之間的交互作用，提出了包含海島基礎(chǔ)承載力、現(xiàn)實(shí)承載力兩個綜合指數(shù)的海島生態(tài)系統(tǒng)承載力評估方法，評估了相關(guān)海島生態(tài)承載力狀況[21]。目前，以生態(tài)環(huán)境脆弱的無居民海島為研究對象進(jìn)行系統(tǒng)的海島生態(tài)承載力研究還鮮見報道。現(xiàn)有海島生態(tài)承載力研究實(shí)踐以靜態(tài)評估為主[22－24]，且評價指標(biāo)及方法存在較大的局限性[25－28]，難以適應(yīng)生態(tài)系統(tǒng)資源開發(fā)與發(fā)展的要求。另外，也較為缺乏從生態(tài)系統(tǒng)服務(wù)和生態(tài)敏感性方面研究生態(tài)承載力的案例實(shí)踐。

本研究以平潭大嶼島這一無居民海島作為研究對象，基于以往學(xué)者們的研究基礎(chǔ)，分別從水土保持服務(wù)、海岸侵蝕、水土流失敏感性、植被初級凈生產(chǎn)力以及生物多樣性等單項指標(biāo)來分析無居民海島的生態(tài)條件，構(gòu)建生態(tài)承載力評估指標(biāo)體系，對大嶼島開發(fā)前后的生態(tài)承載力進(jìn)行綜合評估。期待著本研究結(jié)果能夠?yàn)榇髱Z島乃至我國無居民海島的保護(hù)與開發(fā)提供決策參考。

1 材料與方法

1.1 研究區(qū)概況

平潭大嶼島(以下簡稱大嶼島)屬于無居民海島，位于平潭海壇海峽中段東側(cè)，平潭海峽大橋南側(cè)，東距吉釣島約1.8km，東北距平潭娘宮碼頭約2.8km，西北距福清小山東碼頭約2.9km，具有優(yōu)越的地理位置、良好的資源稟賦、典型的生態(tài)系統(tǒng)和原始自然的海島風(fēng)貌[29]。2014 年大嶼島處于尚未開發(fā)階段，其總面積為30.76hm2。其中，未利用地面積13.52hm2，占總面積的43.94%，主要集中在該島嶼南部；林地面積7.60hm2，占總面積的24.71%；沙地面積9.29hm2，占總面積的30.20%；水域面積僅0.35hm2，占總面積的1.14%。2019 年大嶼島正在開發(fā)建設(shè)中，增加建設(shè)用地3.08hm2，占總面積的10.02%；林地面積增加0.7hm2，占總面積的26.99%；未利用地面積減少5.12hm2，占總面積的27.3%；沙灘面積減少0.37hm2，占總面積的29%；由于沙灘海岸侵蝕，導(dǎo)致水域面積增加1.71hm2，占總面積的6.68%，主要集中在沿海沙灘。圖1 為平潭大嶼島開發(fā)前后的土地利用示意圖。

圖1 平潭大嶼島土地利用示意圖

1.2 研究方法

1.2.1 評價指標(biāo)

本文根據(jù)科學(xué)性、系統(tǒng)性、實(shí)用性及數(shù)據(jù)可獲得性等原則[30]，結(jié)合自然資源部2019 年發(fā)布的資源環(huán)境承載能力和國土空間開發(fā)適宜性評價技術(shù)指南[31]和研究區(qū)實(shí)際情況，構(gòu)建生態(tài)承載力評價指標(biāo)體系(表1)。采用實(shí)用性較強(qiáng)的層次分析法(AHP，Analytic Hierarchy Process)確定評價指標(biāo)權(quán)重[32－33]。對比各個指標(biāo)的重要程度，根據(jù)遞階賦權(quán)定律可以得到最低層指標(biāo)相對于最髙目標(biāo)層的權(quán)重，具體見表1。

表1 生態(tài)承載力指標(biāo)的標(biāo)準(zhǔn)和權(quán)重

1.2.2 單因素評價方法

(1)水土保持功能評估:采用修正通用水土流失方程(RUSLE，Revised Universal Soil Loss Equation)的水土保持服務(wù)模型開展評價，公式[34]如下:

在式(1)中，AC為水土保持量(t/hm2a)；R為降雨侵蝕力因子(MJmm/hm2ha)，具體公式[35]如下；

在式(1)中，K為土壤可蝕性因子(thm2h/hm2MJmm)，可結(jié)合土地利用現(xiàn)狀，參考相關(guān)文獻(xiàn)[36－37]對相應(yīng)的土壤類型進(jìn)行賦值；LS為地形因子，S表示坡度因子，具體公式[38－41]如下:

L表示坡長因子，具體公式如[40，42]下:

在式(1)中，C為植被覆蓋因子；P為水土保持措施因子。根據(jù)相關(guān)文獻(xiàn)的研究方法[43]，可以得到各個土地利用類型的C值和P值。

(2)海岸侵蝕評估:基于ArcGIS 平臺，結(jié)合不同時期的岸線數(shù)據(jù)，根據(jù)數(shù)字岸線分析系統(tǒng)功能模塊來計算監(jiān)測岸段的侵蝕速率，并依據(jù)海岸侵蝕強(qiáng)度等級劃分(表2)，得到海岸侵蝕現(xiàn)狀等級[44]。

表2 海岸侵蝕強(qiáng)度等級劃分

(3)水土流失敏感性評估:根據(jù)通用水土流失方程的基本原理，公式如下:

式中，SSi為i空間單元水土流失敏感性指數(shù)。

(4)植被凈初級生產(chǎn)力:植被凈初級生產(chǎn)力(NPP，Net Primary Productivity)指綠色植物在單位面積、單位時間內(nèi)通過光合作用積累的有機(jī)物質(zhì)總量扣除自養(yǎng)呼吸后的剩余部分。有很多模型可估算NPP，其中朱文泉提出的CASA 模型運(yùn)用較為廣泛。本文將采用該方法求算區(qū)域的NPP，其中由于數(shù)據(jù)限制，大嶼島的太陽總輻射采用統(tǒng)一值。

(5)生物多樣性評估:生物多樣性維護(hù)功能是生態(tài)系統(tǒng)維持基因、物種、生態(tài)系統(tǒng)多樣性所發(fā)揮的作用，是生態(tài)系統(tǒng)提供的最主要功能之一?？捎蒙锒鄻有跃S護(hù)服務(wù)能力指數(shù)作為評估指標(biāo)。?生態(tài)保護(hù)紅線劃定指南?中其計算公式為:

式中，S為生物多樣性維護(hù)服務(wù)能力指數(shù)，NPPmean為年植被凈初級生產(chǎn)力平均值，F(xiàn)pre為年降水量，F(xiàn)tem為年氣溫，F(xiàn)alt為海拔因子。各因子均進(jìn)行歸一化。

將以上圖件均重采樣至10m×10m 柵格，并做歸一化處理。根據(jù)AHP 方法確定因子權(quán)重之后，將因子權(quán)重和因子得分進(jìn)行柵格計算，通過ArcGIS 10.2 中柵格分類的自然斷點(diǎn)法，生成2014 年與2019 年大嶼島生態(tài)承載力分布圖。

1.3 數(shù)據(jù)來源

結(jié)合研究區(qū)域具體情況以及數(shù)據(jù)可獲得性，選取相應(yīng)的指標(biāo)。氣象數(shù)據(jù)來自大嶼島附近的四個氣象站點(diǎn):平潭、福清、莆田和長樂；高程數(shù)據(jù)來源于福建大嶼生態(tài)示范規(guī)劃設(shè)計方案；土壤數(shù)據(jù)來源于全國第二次土壤普查結(jié)果；遙感影像圖來源于91 衛(wèi)圖助手的第20 級數(shù)據(jù)。2014 年遙感影像日期為2014 年1 月25 日，2019 年遙感影像圖日期為2019 年11 月17 日，采用ENVI 軟件，結(jié)合現(xiàn)場查勘情況進(jìn)行目視解譯。植被覆蓋圖是根據(jù)2014 年7 月14 日與2019 年8 月23 日Landsat OLI 衛(wèi)星(http://www.gscloud.cn/)生成的，分辨率為15m，具體信息見表3。

表3 研究區(qū)數(shù)據(jù)信息表

2 結(jié)果與討論

2.1 水土保持服務(wù)功能

圖2 為降雨侵蝕、土壤可蝕性、坡度坡長、植被覆蓋、水土保持措施等因子示意圖。

圖2 大嶼島水土保持服務(wù)功能相關(guān)因子示意圖

由圖2 可見，2014—2019 年，大嶼島降雨侵蝕力呈從西北地區(qū)向東南地區(qū)逐漸遞減的狀況(圖2a 和圖2b)，2019 年降雨侵蝕因子最高值為2014 年最高值的近2 倍，高達(dá)7.56，最低值也是同樣情況，高達(dá)7.54，這與降雨的年際變異性直接相關(guān)；土壤可蝕性因子K值為0 的主要分布在中部地區(qū)，K值為0.12 的主要分布在大嶼島的西北以及東南地區(qū)(圖2c)；坡長坡度因子整體呈現(xiàn)西北和東南地區(qū)較高，其他地區(qū)較低的情況，其中最高值高達(dá)72.88(圖2d)。由圖2e 和圖2f 可知，大嶼島的植被覆蓋因子和水土保持措施具有一定的空間差異性。

根據(jù)式(1)對圖2 各因子進(jìn)行運(yùn)算，獲得大嶼島水土保持服務(wù)功能空間分布圖(圖3a 和圖3b)。由圖3a 和圖3b 統(tǒng)計可得，2014—2019 年水土保持服務(wù)功能增強(qiáng)，其中分布面積最廣的是“一般”，2014—2019 年 從 62.47% 降低到53.76%，主要分布在中部以及西北地區(qū)，其次是“低度”，從2014 年的26.21%降低到2019 年的16.18%。2014—2019 年，水土保持服務(wù)功能中度、中高度以及高度所占比例增加，其中變化最大的是“中度”，從8.68%增加到17.84%，“中高度”從2.14%增加到8.22%，“高度”從0.5%增加到3.99%。

圖3 大嶼島水土保持服務(wù)功能空間分布示意圖(a 為2014 年、b 為2019 年)、大嶼島海岸侵蝕空間分布示意圖(c 為2014 年、d 為2019 年)、大嶼島水土流失敏感性空間分布示意圖(e 為2014 年、f 為2019 年)

隨著大嶼島的開發(fā)建設(shè)，水土保持服務(wù)功能也隨之變化。為此應(yīng)權(quán)衡兩者關(guān)系，探索兩者可以并舉的發(fā)展道路，可以種植某些草本植物來增強(qiáng)水土保持服務(wù)功能[45－46]，亦可從生態(tài)系統(tǒng)服務(wù)功能視角提出相應(yīng)生態(tài)管理措施。渠斌等(2019)利用最小累積阻力模型構(gòu)建了海壇島生態(tài)安全格局[47]。對水土保持服務(wù)等方面的研究有助于合理利用管控土地空間，合理分配自然資源，調(diào)節(jié)發(fā)展與生態(tài)環(huán)境保護(hù)間的相互平衡，為無居民海島的保護(hù)和開發(fā)提供理論依據(jù)和科學(xué)指導(dǎo)[48]。

2.2 海岸侵蝕

大嶼島海岸侵蝕評估結(jié)果見圖3c 和圖3d。從圖3c 和圖3d 可以看出，2014 年和2019 年海岸侵蝕穩(wěn)定級別比例最大，分別高達(dá)57.69%和63.54%，主要分布在西北地區(qū)和東南地區(qū)，具體見圖3。其中2014 年侵蝕級別所占比例僅次于穩(wěn)定級別，為26.22%，微侵蝕占比最小，僅為16.09%；2019 年侵蝕級別高于2014 年同一級別，高達(dá)36.46%。

文世勇等(2020)用該方法研究不同時段海南省的岸線侵蝕變化狀況并預(yù)測局部岸段岸線變化趨勢[49]，劉勇等(2016)基于多源數(shù)據(jù)利用GIS分析了東山島海岸線類型、位置與長度變化及灘面沖蝕演變[50]，Ahmed 等(2021)選擇13 項空間標(biāo)準(zhǔn)從物理和經(jīng)濟(jì)社會脆弱性兩方面，研究孟加拉國東部沿海地區(qū)海岸侵蝕脆弱性[51]，Sahin 等(2019)利用空間貝葉斯網(wǎng)絡(luò)預(yù)測太平洋小島的海平面上升引起的海岸侵蝕[52]。當(dāng)前短時間尺度的人類活動已成為海岸侵蝕的最主要因素[53－54]，研究及預(yù)測海岸侵蝕變化有助于在脆弱的沿海地區(qū)制定積極的緩解策略，以保護(hù)人類和資源免受破壞性沿海災(zāi)害的影響。

2.3 水土流失敏感性

根據(jù)式(6)，結(jié)合圖2 各子圖，開展大嶼島水土流失敏感性評估，結(jié)果如圖3e 和圖3f 所示。由圖3e 和圖3f 可見，2014 年和2019 年占比最大的就是極強(qiáng)水土流失敏感性，從2014 年的44.17%增加到2019 年的44.81%，主要分布在大嶼島的中部地區(qū)以及東南地區(qū)，主要是沙地和建設(shè)用地；其次是強(qiáng)度水土流失敏感性，從23.19%增加到25.07%。水土流失敏感性變化最大的為水土流失輕度敏感性，從 14.01% 降低到11.42%，主要分布在東南地區(qū)的未利用地。2014—2019 年大嶼島水土流失敏感性整體呈增加趨勢。趙舒騰(2016)基于RUSLE 模型運(yùn)用探索性空間分析法探究海壇島水土流失空間分布情況[55]，賴華燕等(2017)運(yùn)用GIS 研究平潭島建設(shè)發(fā)展過程中的水土流失狀況[43]，Lidia 等(2020)利用GIS 研究1960—2012 年西班牙大加那利島水土流失情況[56]，王維碧等(2020)使用帶有輔助數(shù)據(jù)的NSED－C 分析了2000—2015 年整個黑土地區(qū)的土壤侵蝕時空變化[57]。國內(nèi)外學(xué)者主要利用GIS 研究分析水土流失情況，可為該地區(qū)的土地合理利用以及海島開發(fā)利用采取的最優(yōu)防治水土流失措施提供參考。

2.4 植被凈初級生產(chǎn)力

圖4 為大嶼島月降水、植被覆蓋指數(shù)(NDVI，Normalized Difference Vegetation Index)、月氣溫等因子圖。

圖4 大嶼島植被凈初級生產(chǎn)力相關(guān)因子示意圖

由圖4 可見，2014—2019 年大嶼島月降雨量呈現(xiàn)從西北地區(qū)向東南地區(qū)逐漸遞減的狀況(圖4a和圖4b)，且月降雨量相差不大；NDVI 因子(圖4c和圖4d)整體呈現(xiàn)西北地區(qū)和東南地區(qū)較高、其他地區(qū)較低的狀況，2019 年最低值主要分布在沙灘；月氣溫從西北地區(qū)向東南地區(qū)遞減且相差不大(圖4e 和圖4f)。

基于ENVI 軟件，利用CASA 模型，可得到植被凈初級生產(chǎn)力空間分布圖(圖5)。由圖5 可見，2014—2019 年植被凈初級生產(chǎn)力持續(xù)增強(qiáng)，2014 年分布最廣的是中度植被凈初級生產(chǎn)力，主要分布在中部以及西北地區(qū)，其次是低度和高度植被凈初級生產(chǎn)力，所占面積比例分別為26%和26.48%。2019 年高度植被凈初級生產(chǎn)力分布最廣，從26.48%增加到39.64%，其中面積變化最大的是中度植被凈初級生產(chǎn)力，從31.12%降低到10.86%，主要與林地面積增加以及未利用地面積減少有關(guān)。

圖5 大嶼島植被凈初級生產(chǎn)力空間分布示意圖

植被凈初級生產(chǎn)力不僅是調(diào)節(jié)生態(tài)過程的重要指標(biāo)，也是碳源(匯)的判定因子，在碳平衡和應(yīng)對全球氣候變化中具有重要的作用。國外學(xué)者側(cè)重于相關(guān)理論及模型等方面研究[58－59]，國內(nèi)學(xué)者主要研究模型改進(jìn)及應(yīng)用等方面問題[60－61]。光能利用模型的相關(guān)參數(shù)可通過遙感手段獲取，并考慮了不同的植被覆蓋情況，精度較高且適用于小尺度地區(qū)。其中朱文泉(2005)改進(jìn)的CASA模型模擬了中國主要植被類型的最大光能利用率，根據(jù)區(qū)域蒸散模型來模擬水分脅迫因子，并以2002 年的內(nèi)蒙古自治區(qū)植被為例，研究了植被凈初級生產(chǎn)力及其時空分布[62－63]。

2.5 生物多樣性

圖6 為大嶼島年降雨、年氣溫、海拔等因子圖。

由圖6 可見，2014 年和2019 年大嶼島月降雨量呈現(xiàn)從西北地區(qū)向東南地區(qū)逐漸遞減的狀況(圖4a 和圖4b)；年氣溫因子整體呈現(xiàn)截然相反的情況，2014 年呈現(xiàn)從西北地區(qū)向東南地區(qū)逐漸遞減，2019 年呈現(xiàn)從西北向東南地區(qū)逐漸遞增的狀況(圖4c 和圖4d)；海拔因子整體呈現(xiàn)西北地區(qū)和東南地區(qū)較高、其他地區(qū)較低的狀況(圖4e)。

圖6 大嶼島生物多樣性相關(guān)因子示意圖

根據(jù)式(7)開展大嶼島生物多樣性評估，結(jié)果見圖7。由圖7 可知，2014—2019 年生物多樣性呈降低趨勢，分布最廣的是一般生物多樣性，主要分布在中部地區(qū)，2014 年分布面積其次的是高度生物多樣性，主要分布在西北地區(qū)；2019 年分布面積其次的是低度生物多樣性，從2014 年的21.68%增加到25.2%。其中面積變化最大的是高度生物多樣性，從2014 年的27.68%降低到3.24%，主要是人為因素造成的，與當(dāng)?shù)亟ㄔO(shè)開發(fā)導(dǎo)致建設(shè)用地增加密切相關(guān)。

圖7 大嶼島生物多樣性空間分布

生物多樣性受多種因素影響，比如土地、水以及兩者的共同作用。其中土地對生物多樣性的影響包括城市擴(kuò)張、景觀設(shè)計、農(nóng)業(yè)景觀/功能異質(zhì)性以及土地利用變化[64－67]。闡明生物多樣性—生態(tài)系統(tǒng)服務(wù)—人類社會福祉相互作用機(jī)理，辨析生物多樣性與其他生態(tài)系統(tǒng)服務(wù)之間的協(xié)同關(guān)系，采取適宜的土地利用政策手段豐富生物多樣性，實(shí)現(xiàn)生態(tài)系統(tǒng)服務(wù)的持續(xù)、穩(wěn)定供給，提升人類社會福祉是未來的研究重點(diǎn)之一[68]。

2.6 生態(tài)承載力

圖8 為基于上述五類子因子圖開展的大嶼島生態(tài)承載力計算結(jié)果圖。結(jié)合圖8 和表4 統(tǒng)計可得，2014 年較低生態(tài)承載力所占比例最大，高達(dá)32.19%，主要分布在中部地區(qū)；其次是中等承載力，所占比例為28.41%；再者是較高承載力，占比16.09%。經(jīng)過五年時間，2019 年占比最大的為較高承載力，高達(dá)25.13%，比2014 年多出9.04%，主要分布在東南地區(qū)；其次是低承載力和中等承載力，分別占23.99%和19.08%。其中變化最明顯是較低生態(tài)承載力，從2014 年的32.19%降低到2019 年的14.06%，主要分布在中部地區(qū)。2014—2019 年大嶼島生態(tài)承載力整體呈現(xiàn)持續(xù)增加的狀況，引致這一現(xiàn)象產(chǎn)生的原因是大嶼島的林地增加。該地區(qū)的凈初級生產(chǎn)力、水土保持服務(wù)呈遞增趨勢且差異較為明顯，水土流失敏感性雖呈遞增趨勢但差異不大，生物多樣性呈遞減趨勢。雖生態(tài)承載力整體呈持續(xù)增加但達(dá)到某一情況將會下降，因此需適度開發(fā)?，F(xiàn)有整體開發(fā)布局具有一定的合理性，但仍有超過一半的開發(fā)區(qū)域?qū)傧拗苹虿贿m宜開發(fā)區(qū)域，因此需將此部分開發(fā)區(qū)域內(nèi)的開發(fā)活動向承載力更高的地區(qū)進(jìn)行適當(dāng)轉(zhuǎn)移，減少島上限制或不適宜開發(fā)區(qū)域內(nèi)的建設(shè)面積，對于島上新的建設(shè)活動規(guī)劃也應(yīng)避開此類區(qū)域，使開發(fā)空間得到進(jìn)一步優(yōu)化。

圖8 大嶼島生態(tài)承載力空間分布

表4 生態(tài)承載力評估結(jié)果面積占比統(tǒng)計結(jié)果 單位/%

現(xiàn)有的有關(guān)生態(tài)承載力評估方法眾多，其中生態(tài)足跡、能值分析等方法[20，23，28]只能關(guān)注少數(shù)主要因素，通過忽略詳細(xì)過程簡化計算方法，對于生態(tài)過程機(jī)制的研究是不夠深刻的。如同渠斌等人(2019)從生態(tài)系統(tǒng)服務(wù)功能角度入手構(gòu)建海壇島生態(tài)安全格局[47]，本文從水土保持服務(wù)等生態(tài)系統(tǒng)服務(wù)、海岸侵蝕以及水土流失生態(tài)敏感性等方面構(gòu)建生態(tài)承載力評估指標(biāo)體系，可以更好地揭示生態(tài)結(jié)構(gòu)之間的關(guān)系，能夠產(chǎn)生更全面、更精確的生態(tài)承載力核算結(jié)果[69]。目前大多數(shù)生態(tài)承載力研究側(cè)重于靜態(tài)評估[22－24]，而本文研究2014—2019 年大嶼島生態(tài)承載力為動態(tài)評估。動態(tài)評估研究有助于深入動態(tài)預(yù)測研究，建立科學(xué)實(shí)用的預(yù)警機(jī)制。生態(tài)承載力作為可持續(xù)發(fā)展政策制定的重要參量，其未來的承載力預(yù)估則更為重要[70]，能夠?yàn)闊o居民海島合理開發(fā)與永續(xù)發(fā)展提供科技支撐。

3 結(jié)論

本研究從水土保持服務(wù)、海岸侵蝕、水土流失敏感性、植被初級凈生產(chǎn)力以及生物多樣性等方面分析無居民海島的生態(tài)條件，構(gòu)建生態(tài)承載力評估指標(biāo)體系，開展大嶼島開發(fā)前后的生態(tài)承載力評估。2014 年較低生態(tài)承載力所占比例最大，高達(dá)32.19%；其次是中等承載力，所占比例為28.41%；再者是較高承載力，占16.09%。經(jīng)過海島開發(fā)，2019 年所占比例最大為較高承載力，高達(dá)25.13%，比2014 年多出9.04%；其次是低承載力和中等承載力，分別占23.99%和19.08%。其中變化最明顯是較低生態(tài)承載力，從2014 年的32.19%降低到2019 年的14.06%，主要分布在中部地區(qū)。2014—2019 年大嶼島生態(tài)承載力整體呈一定的增加，擁有較好的國土空間自然生態(tài)本底條件，但未來需要加強(qiáng)對西北和東南兩區(qū)域的保護(hù)力度。本研究結(jié)果可為無居民海島的合理開發(fā)與永續(xù)發(fā)展提供科技支撐。