王方怡，張嘉誠(chéng)，曹 彥,2，游巍斌，何東進(jìn)

(1.福建農(nóng)林大學(xué)林學(xué)院，福建 福州 350002；2.福建江夏學(xué)院金融學(xué)院，福建 福州 350108)

濕地被稱(chēng)為“地球之腎”，是寶貴的自然資源和重要的生態(tài)系統(tǒng)，具有不可替代的綜合功能[1-2]。濱海濕地是濕地的一種類(lèi)型，是陸地和海洋的交錯(cuò)過(guò)渡地帶[3-4]。濱海濕地具有提供水源、補(bǔ)充地下水、抗洪防澇、抵御風(fēng)暴等功能，對(duì)人類(lèi)社會(huì)的持續(xù)發(fā)展有重要貢獻(xiàn)[5-6]。1979年，我國(guó)從美國(guó)東海岸引進(jìn)互花米草，作為重要的生態(tài)工程材料使用[7]，取得了抗風(fēng)防浪、保灘護(hù)岸、促淤造陸、固碳等生態(tài)和經(jīng)濟(jì)效益，但另一方面互花米草破壞生物棲息環(huán)境、堵塞航道、影響海水交換能力，直接危及當(dāng)?shù)睾０渡鷳B(tài)系統(tǒng)，甚至對(duì)當(dāng)?shù)仞B(yǎng)殖業(yè)的經(jīng)濟(jì)發(fā)展產(chǎn)生威脅，2003年被我國(guó)國(guó)家環(huán)保總局列為首批16種外來(lái)入侵物種[8]。近年來(lái)，互花米草在我國(guó)東部沿海濕地正以驚人的擴(kuò)散速度入侵，已成為中國(guó)沿海灘涂地分布面積最廣的灘涂草本鹽沼植物[9]。黃河三角洲自然保護(hù)區(qū)互花米草的分布面積從2010年的62.6 hm2增至2015年的438.1 hm2，增幅近600%[10]。由此看來(lái)，亟需深入了解和評(píng)估互花米草入侵對(duì)濕地生態(tài)系統(tǒng)造成的危害。

互花米草入侵生態(tài)風(fēng)險(xiǎn)是指外來(lái)物種入侵發(fā)生、發(fā)展的可能性，以及給生態(tài)系統(tǒng)及其組分造成損失的可能性[11-12]。目前，關(guān)于互花米草入侵的研究多側(cè)重于互花米草發(fā)生、擴(kuò)散的自然過(guò)程以及入侵后造成的損失，而往往忽略入侵的驅(qū)動(dòng)力影響。研究區(qū)域主要集中在長(zhǎng)江口上海崇明東灘、九段沙以及三角洲區(qū)域，江蘇中南部海岸，浙江杭州灣，廣西北海海岸，福建閩江口、漳江口、羅源灣地區(qū)[13]。閩東濱海濕地為福建省3個(gè)重要濱海濕地之一，是我國(guó)紅樹(shù)林天然分布的最北緣，生境地位極其特殊[14]，并且是福建省互花米草入侵最嚴(yán)重的地區(qū)。但當(dāng)前對(duì)閩東濱海濕地的研究主要集中在生物學(xué)特性[15]、生態(tài)學(xué)特征[16]、對(duì)環(huán)境的影響[17]等方面，而針對(duì)互花米草生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)方面的宏觀(guān)研究還未見(jiàn)報(bào)道。

鑒于此，本研究以閩東濱海濕地為研究對(duì)象，引入驅(qū)動(dòng)力-壓力-狀態(tài)-影響-響應(yīng)(driving force-pressure-state-impact-response，DPSIR)概念模型，結(jié)合遙感影像、野外實(shí)測(cè)數(shù)據(jù)及氣象數(shù)據(jù)、社會(huì)經(jīng)濟(jì)統(tǒng)計(jì)資料，綜合評(píng)價(jià)互花米草入侵生態(tài)風(fēng)險(xiǎn)狀況，以期為閩東濱海濕地生態(tài)環(huán)境保護(hù)、經(jīng)濟(jì)建設(shè)及環(huán)境管理之間的協(xié)調(diào)有序發(fā)展提供參考價(jià)值。

1 研究區(qū)概況與數(shù)據(jù)收集

1.1 研究區(qū)概況

福建省寧德市別稱(chēng)閩東，位于福建省東北翼沿海，以寧德市霞浦縣(119°46′～120°26′E，26°25′～27°09′N(xiāo))作為研究區(qū)域，全縣陸地面積1708.4km2，海域面積29592.6km2，其中灘涂面積218.5km2，海岸線(xiàn)長(zhǎng)479.9km。該區(qū)氣候?qū)僦衼啛釒ШＱ笮约撅L(fēng)氣候，兼具山地氣候、盆谷地氣候等多種氣候特點(diǎn)，光照充足，春夏雨熱同期，秋冬光溫互利；氣候濕潤(rùn)，雨量充沛；年平均氣溫16～19℃。濱海灘涂由沙、粉砂質(zhì)泥組成，但因山麓迫岸，且有短小山溪注入，故在近岸潮灘上多見(jiàn)數(shù)10m的砂礫堆積帶，從而形成特殊的“礫泥灘”，灘面上生長(zhǎng)著稀疏的紅樹(shù)林和茂密的互花米草。

1.2 數(shù)據(jù)收集

統(tǒng)計(jì)數(shù)據(jù)主要來(lái)源于2006—2016年《霞浦縣統(tǒng)計(jì)年鑒》、霞浦縣林業(yè)局官網(wǎng)數(shù)據(jù)、福建省氣象站點(diǎn)監(jiān)測(cè)資料、地理空間數(shù)據(jù)云、外業(yè)調(diào)查數(shù)據(jù)等，鹽度數(shù)據(jù)源自室內(nèi)試驗(yàn)結(jié)果。

2 研究方法

2.1 遙感影像處理

收集覆蓋2006、2011和2016年3期的遙感數(shù)據(jù)，采用福建省地質(zhì)測(cè)繪院遙感中心提供的多源高分辨率衛(wèi)星影像數(shù)據(jù)，其中2016年遙感影像數(shù)據(jù)有高分一號(hào)(GF-1)、天繪一號(hào)(TH-1)，2011年的為RapidEye影像數(shù)據(jù)，2006年的為谷歌影像(Google Earth)共23景作為影像數(shù)據(jù)源。由于互花米草生長(zhǎng)在潮灘上，且生長(zhǎng)季在6—10月，為滿(mǎn)足覆蓋解譯的要求，故收集潮灘低潮位、無(wú)云覆蓋、6—10月時(shí)段的影像。影像數(shù)據(jù)均為全彩色影像，空間分辨為2.00～8.53m。另外收集用于校正的底圖1∶10萬(wàn)地形圖以及野外勘測(cè)GPS控制點(diǎn)。

首先對(duì)原始影像進(jìn)行幾何精校正、影像配準(zhǔn)、融合、鑲嵌及裁剪，然后通過(guò)現(xiàn)場(chǎng)野外考察9種類(lèi)型的實(shí)物現(xiàn)狀，獲取衛(wèi)星遙感影像解譯標(biāo)志，利用Arcgis10.0軟件平臺(tái)對(duì)研究區(qū)影像進(jìn)行目視解譯，數(shù)據(jù)矢量化后進(jìn)行融合，消除破碎的斑塊，確保保留準(zhǔn)確的有效信息，最后建立拓?fù)湟?guī)則，修正3期影像數(shù)據(jù)的邊界，使得研究范圍保持一致，最終得到閩東濱海濕地2006、2011、2016年土地利用類(lèi)型圖。計(jì)算得出分類(lèi)的精度分別為89.74%、91.20%、92.33%，綜合前人的研究可知,當(dāng)分類(lèi)精度大于85%時(shí)，分類(lèi)結(jié)果較好，因此本研究的遙感影像滿(mǎn)足要求。

由于存在著空間異質(zhì)性和不同程度的干擾強(qiáng)度，不同地區(qū)之間的風(fēng)險(xiǎn)水平具有一定的差異，為保證風(fēng)險(xiǎn)評(píng)價(jià)工作的準(zhǔn)確性，本研究把行政區(qū)的統(tǒng)計(jì)數(shù)分布在各單元格內(nèi)，以建立空間數(shù)據(jù)庫(kù)，一方面能夠簡(jiǎn)化不同來(lái)源數(shù)據(jù)的統(tǒng)計(jì)，另一方面可以提高計(jì)算結(jié)果和實(shí)際情況的符合程度。根據(jù)研究區(qū)面積大小及采樣工作量，采用3km×3km的正方形單元網(wǎng)格進(jìn)行等間距采樣，共劃分291個(gè)風(fēng)險(xiǎn)小區(qū)，空間疊加分析方法計(jì)算研究區(qū)生態(tài)風(fēng)險(xiǎn)值，并將這些風(fēng)險(xiǎn)值賦予到每個(gè)小區(qū)的中心點(diǎn)。

2.2 驅(qū)動(dòng)力-壓力-狀態(tài)-影響-響應(yīng)模型

1993年，OECD[18]在PSR模型基礎(chǔ)上正式提出了DPSIR模型，并被歐洲環(huán)境局(European Environment Agency)所采用，將表征自然系統(tǒng)的評(píng)價(jià)指標(biāo)從原來(lái)的3個(gè)擴(kuò)展為驅(qū)動(dòng)力、壓力、狀態(tài)、影響和響應(yīng)5個(gè)部分[19]。DPSIR模型的理論框架為：在經(jīng)濟(jì)發(fā)展、人口增多等潛在因子的驅(qū)動(dòng)下，產(chǎn)生破壞和干擾生態(tài)系統(tǒng)的直接因素，對(duì)生態(tài)系統(tǒng)施加一定的壓力，促使自然環(huán)境狀態(tài)的改變，進(jìn)而對(duì)環(huán)境系統(tǒng)產(chǎn)生各種影響，驅(qū)使人們做出響應(yīng)以阻止、補(bǔ)償、改善或適應(yīng)環(huán)境的變化[20]。該框架能在一定程度上揭示和深入剖析人與環(huán)境之間的因果關(guān)系，被廣泛應(yīng)用于分析和評(píng)估人類(lèi)和環(huán)境系統(tǒng)之間的相互作用的過(guò)程[21-25]。

模型中，驅(qū)動(dòng)力是指促使互花米草入侵的潛在原因。結(jié)合濱海濕地互花米草入侵化影響因子的作用機(jī)理可知，社會(huì)經(jīng)濟(jì)因素并不是直接作用于土地導(dǎo)致互花米草入侵的發(fā)生與擴(kuò)散，而是作為一種驅(qū)動(dòng)使得不合理的人為活動(dòng)的產(chǎn)生。另外，自然環(huán)境因素，如氣溫、降水量等驅(qū)動(dòng)互花米草在適宜的氣候定殖。壓力是指在驅(qū)動(dòng)力的作用之后，直接干擾生態(tài)系統(tǒng)的各種因素。自然環(huán)境因素會(huì)引發(fā)惡劣的生態(tài)問(wèn)題，如在全球變暖的大背景下，互花米草將逐漸向高緯度地區(qū)擴(kuò)散，極端氣候變化、海平面上升，地形高程降低，可能增加互花米草相對(duì)競(jìng)爭(zhēng)力，加速向高潮灘土著植物群落入侵。而人口、社會(huì)經(jīng)濟(jì)的驅(qū)動(dòng)會(huì)引起不合理的人為活動(dòng)，圍墾和過(guò)度養(yǎng)殖是引起濱海濕地互花米草入侵的人為因子。狀態(tài)是指在壓力作用下濕地環(huán)境的變化狀況。在互花米草入侵的地區(qū)，最明顯的環(huán)境特征就是濱海濕地互花米草入侵化程度，此外，自然和社會(huì)環(huán)境變化也可反映其狀態(tài)。生態(tài)環(huán)境的狀態(tài)越差，其入侵化生態(tài)風(fēng)險(xiǎn)越高。影響是指互花米草入侵對(duì)生態(tài)系統(tǒng)的自然環(huán)境和社會(huì)經(jīng)濟(jì)的影響。在入侵的過(guò)程中，必然會(huì)導(dǎo)致土地退化，對(duì)外界干擾的抵抗力及恢復(fù)力降低，進(jìn)一步加大了入侵生態(tài)風(fēng)險(xiǎn)。另一方面，土地退化使土地生產(chǎn)力下降，造成農(nóng)業(yè)減產(chǎn)，侵占養(yǎng)殖場(chǎng)的土地，造成漁業(yè)的損失。響應(yīng)是人類(lèi)為了預(yù)防、減輕或者消除不良的影響而采取的措施。如在互花米草防治的過(guò)程中提出化學(xué)、物理、生物方法控制互花米草入侵，治理互花米草的大力投入和產(chǎn)業(yè)結(jié)構(gòu)調(diào)整都可消除和替代不利的入侵化驅(qū)動(dòng)因素，從而防治互花米草入侵對(duì)生態(tài)環(huán)境造成損失的可能性。具體概念模型見(jiàn)圖1。

2.3 評(píng)價(jià)指標(biāo)體系構(gòu)建

建立的指標(biāo)體系框架將評(píng)價(jià)指標(biāo)分為3層。第1層是目標(biāo)層，即閩東濱海濕地互花米草入侵生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)；第2層為評(píng)價(jià)項(xiàng)目層，分別為驅(qū)動(dòng)力、壓力、狀態(tài)、影響、響應(yīng)；第3層為指標(biāo)層，即項(xiàng)目層每個(gè)指標(biāo)的具體評(píng)價(jià)指標(biāo)。通過(guò)大量文獻(xiàn)閱讀和實(shí)地調(diào)研，參考權(quán)威機(jī)構(gòu)發(fā)布的指標(biāo)體系和研究學(xué)者的研究成果[26]，篩選出16個(gè)具有代表性的指標(biāo)(表1)，構(gòu)建閩東濱海濕地互花米草入侵生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)的指標(biāo)體系。

圖1 互花米草入侵生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)概念模型

表1 互花米草入侵生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)指標(biāo)體系及權(quán)重

3 結(jié)果與分析

3.1 指標(biāo)權(quán)重的確定

評(píng)價(jià)指標(biāo)不同，其重要程度也不同。為將各個(gè)指標(biāo)的重要性做出定量分配，需對(duì)各指標(biāo)進(jìn)行權(quán)重的確定。本研究采用SAATY[27]提出的層次分析法，量化了客觀(guān)判斷的定性分析，將人們的思維過(guò)程數(shù)學(xué)化、系統(tǒng)化，利用數(shù)值化的方式來(lái)消除不同指標(biāo)帶來(lái)的差異，提高評(píng)價(jià)結(jié)果的科學(xué)性。層次分析法將復(fù)雜問(wèn)題中的各個(gè)指標(biāo)歸類(lèi)為互相關(guān)聯(lián)的有序?qū)?，通過(guò)對(duì)客觀(guān)和實(shí)際的大概性判斷，定量化每個(gè)層次中每個(gè)因素對(duì)結(jié)果的影響程度，是一種新型、簡(jiǎn)潔、實(shí)用的研究復(fù)雜系統(tǒng)的方法，已廣泛應(yīng)用于安全科學(xué)和環(huán)境科學(xué)領(lǐng)域[28-29]。結(jié)合專(zhuān)家咨詢(xún)和專(zhuān)家打分的方法，通過(guò)向35位霞浦鹽田林業(yè)站、福鼎林業(yè)局、福建農(nóng)林大學(xué)林學(xué)院等相關(guān)領(lǐng)域的專(zhuān)家、學(xué)者發(fā)放調(diào)查問(wèn)卷的方式，評(píng)價(jià)各指標(biāo)間的相對(duì)優(yōu)先級(jí)。即按一定的標(biāo)準(zhǔn)將每一指標(biāo)劃分為不同等級(jí)進(jìn)行賦分，并對(duì)所得分值標(biāo)準(zhǔn)化，經(jīng)過(guò)計(jì)算檢驗(yàn)，問(wèn)卷的矩陣一致性比率(consistencyratio,CR)均小于0.1，各矩陣均有較好的一致性，使用加權(quán)平均法[30-31]綜合計(jì)算各專(zhuān)家對(duì)各級(jí)指標(biāo)的權(quán)重值(表1)。

3.2 評(píng)價(jià)指標(biāo)等級(jí)劃分

本研究采用加權(quán)求和法獲得生態(tài)風(fēng)險(xiǎn)指數(shù)[32]。

式中:H為生態(tài)風(fēng)險(xiǎn)指數(shù)；Bi為第i參數(shù)的權(quán)重值；Wi為第i參數(shù)的標(biāo)準(zhǔn)化值。

評(píng)估數(shù)值僅是評(píng)估結(jié)果的體現(xiàn)，尚不能表達(dá)生態(tài)風(fēng)險(xiǎn)程度，因此，需要通過(guò)對(duì)數(shù)值限值進(jìn)行界定，才能表達(dá)出生態(tài)風(fēng)險(xiǎn)狀況[33]。為此，本研究依據(jù)歷史資料、實(shí)地考察、借鑒國(guó)內(nèi)外相關(guān)研究文獻(xiàn)中廣泛使用的標(biāo)準(zhǔn)和借助專(zhuān)家經(jīng)驗(yàn)，利用分位數(shù)對(duì)霞浦縣濱海濕地互花米草入侵生態(tài)風(fēng)險(xiǎn)進(jìn)行等級(jí)劃分，分為5等級(jí)，分別為極低生態(tài)風(fēng)險(xiǎn)(0.0≤生態(tài)風(fēng)險(xiǎn)指數(shù)<0.2)、低生態(tài)風(fēng)險(xiǎn)(0.2≤生態(tài)風(fēng)險(xiǎn)指數(shù)<0.4)、中生態(tài)風(fēng)險(xiǎn)(0.4≤生態(tài)風(fēng)險(xiǎn)指數(shù)<0.6)、高生態(tài)風(fēng)險(xiǎn)(0.6≤生態(tài)風(fēng)險(xiǎn)指數(shù)<0.8)和極高生態(tài)風(fēng)險(xiǎn)(0.8≤生態(tài)風(fēng)險(xiǎn)指數(shù)≤1.0)。各評(píng)價(jià)指標(biāo)對(duì)應(yīng)的等級(jí)標(biāo)準(zhǔn)見(jiàn)表2。

表2 互花米草入侵生態(tài)風(fēng)險(xiǎn)各評(píng)價(jià)指標(biāo)的各等級(jí)標(biāo)準(zhǔn)

3.3 生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)結(jié)果

根據(jù)DPSIR評(píng)價(jià)模型，采用Arc GIS 10.0軟件空間分析的反距離權(quán)重法(inverse distance weighted，IDW)插值模塊對(duì)劃分等級(jí)后的入侵生態(tài)風(fēng)險(xiǎn)離散點(diǎn)進(jìn)行插值，形成3期互花米草入侵生態(tài)風(fēng)險(xiǎn)綜合評(píng)價(jià)區(qū)劃圖(圖2)，分別對(duì)3期不同等級(jí)生態(tài)風(fēng)險(xiǎn)面積進(jìn)行統(tǒng)計(jì)(表3)。

圖2 閩東濱海濕地互花米草入侵生態(tài)風(fēng)險(xiǎn)綜合評(píng)價(jià)區(qū)劃圖

表3 閩東濱海濕地不同等級(jí)生態(tài)風(fēng)險(xiǎn)面積

結(jié)果表明，研究區(qū)域內(nèi)各風(fēng)險(xiǎn)等級(jí)空間分布格局呈西北低東南高，中間高，四周低，高、極高生態(tài)風(fēng)險(xiǎn)于研究區(qū)東南部最顯著，中生態(tài)風(fēng)險(xiǎn)集中在中部區(qū)域，低、極低生態(tài)風(fēng)險(xiǎn)以西北部為主，從西北到東南分布特點(diǎn)大致為：極低生態(tài)風(fēng)險(xiǎn)、低生態(tài)風(fēng)險(xiǎn)、中生態(tài)風(fēng)險(xiǎn)、高生態(tài)風(fēng)險(xiǎn)、極高生態(tài)風(fēng)險(xiǎn)、中生態(tài)風(fēng)險(xiǎn)、低生態(tài)風(fēng)險(xiǎn)梯度。2006年極高生態(tài)風(fēng)險(xiǎn)區(qū)主要位于長(zhǎng)春鎮(zhèn)西部、松港東部、下滸鎮(zhèn)北部、溪南鎮(zhèn)和沙江鎮(zhèn)的中部，在2011年蔓延至三沙鎮(zhèn)、鹽田畬族鄉(xiāng)、北壁鄉(xiāng)，2016年北壁鄉(xiāng)極高生態(tài)風(fēng)險(xiǎn)狀態(tài)消失；2006年高生態(tài)風(fēng)險(xiǎn)區(qū)主要位于三沙鎮(zhèn)、溪南鎮(zhèn)、下滸鎮(zhèn)、沙江鎮(zhèn)與長(zhǎng)春鎮(zhèn)交界處，零星分布于牙城鎮(zhèn)、北壁鄉(xiāng)、鹽田畬族鄉(xiāng)、松港，2011年已擴(kuò)散到鹽田畬族鄉(xiāng)與崇儒畬族鄉(xiāng)交界處，2016年則蔓延到牙城鎮(zhèn)，而鹽田畬族鄉(xiāng)與崇儒畬族鄉(xiāng)交界區(qū)域高生態(tài)風(fēng)險(xiǎn)狀態(tài)消失；中生態(tài)風(fēng)險(xiǎn)主要集中在牙城鎮(zhèn)、三沙鎮(zhèn)、松港東部、長(zhǎng)春鎮(zhèn)、下滸鎮(zhèn)、北壁鄉(xiāng)、沙江鎮(zhèn)西部、溪南鎮(zhèn)邊緣地區(qū)；2006年低生態(tài)風(fēng)險(xiǎn)區(qū)主要位于牙城鎮(zhèn)、水門(mén)畬族鄉(xiāng)、松城街道、松港街道北部、長(zhǎng)春鎮(zhèn)、北壁鄉(xiāng)、鹽田畬族鄉(xiāng)、崇儒畬族鄉(xiāng)南部，零星分布于下滸鎮(zhèn)、三沙鎮(zhèn)、柏洋鄉(xiāng)，2011年在鹽田畬族鄉(xiāng)和水門(mén)畬族鄉(xiāng)進(jìn)一步擴(kuò)散；2006年極低生態(tài)風(fēng)險(xiǎn)區(qū)主要位于柏洋鄉(xiāng)大部分區(qū)域、崇儒畬族鄉(xiāng)和鹽田畬族鄉(xiāng)的北部、水門(mén)畬族鄉(xiāng)西部，2011年鹽田畬族鄉(xiāng)極低生態(tài)風(fēng)險(xiǎn)狀態(tài)消失。極高、高風(fēng)險(xiǎn)區(qū)主要集聚在人口較多、社會(huì)經(jīng)濟(jì)相對(duì)發(fā)達(dá)的“兩洋三灣”(東吾洋、官井洋，三沙灣、福寧灣、牙城灣)沿岸地帶，極低、低風(fēng)險(xiǎn)區(qū)主要集聚在以山地、丘陵為主的西北部。

2006—2016年研究區(qū)互花米草的入侵增強(qiáng)，但后5 a呈現(xiàn)減緩的趨勢(shì)。主要表現(xiàn)為：極高生態(tài)風(fēng)險(xiǎn)面積不斷增加，2006—2011年間急劇增加了11 725.08 hm2，2011—2016年間共增加了857.92 hm2，10 a增加量達(dá)12 583.00 hm2；高生態(tài)風(fēng)險(xiǎn)區(qū)面積趨于減少，變化量相對(duì)較少，10 a共減少了2 509.92 hm2；中生態(tài)風(fēng)險(xiǎn)區(qū)范圍呈現(xiàn)擴(kuò)大的趨勢(shì)，面積共增加了4 433.04 hm2；低生態(tài)風(fēng)險(xiǎn)分布面積比例最大，10 a面積先大幅度減少5 593.50 hm2，后增加990.23 hm2，共減少4 603.28 hm2；極低生態(tài)風(fēng)險(xiǎn)面積不斷減少，2006—2011年面積減少5 944.22 hm2，2011—2016年面積減少3 967.98 hm2，10 a面積減少9 912.20 hm2?？梢?jiàn)，在整個(gè)研究期間，研究區(qū)生態(tài)風(fēng)險(xiǎn)格局無(wú)明顯變化，但不同等級(jí)的風(fēng)險(xiǎn)區(qū)面積變化存在明顯差異，中風(fēng)險(xiǎn)和極高風(fēng)險(xiǎn)區(qū)范圍不斷擴(kuò)張，極低風(fēng)險(xiǎn)和高風(fēng)險(xiǎn)區(qū)面積則逐年減少。

從不同生態(tài)風(fēng)險(xiǎn)等級(jí)之間面積相互轉(zhuǎn)化上看，各生態(tài)風(fēng)險(xiǎn)等級(jí)主要向鄰近的風(fēng)險(xiǎn)等級(jí)轉(zhuǎn)化。在整個(gè)研究期間，高風(fēng)險(xiǎn)和極高風(fēng)險(xiǎn)區(qū)面積比重由21.55%增加至26.50%又減至26.30%，而中等風(fēng)險(xiǎn)區(qū)面積上升，這表明研究區(qū)的入侵生態(tài)風(fēng)險(xiǎn)得到有效遏制，濱海生態(tài)環(huán)境有所改善。但極低風(fēng)險(xiǎn)向低風(fēng)險(xiǎn)轉(zhuǎn)化，低風(fēng)險(xiǎn)向中風(fēng)險(xiǎn)區(qū)轉(zhuǎn)化的趨勢(shì)也很明顯，反映出這10 a研究區(qū)互花米草潛在的入侵生態(tài)風(fēng)險(xiǎn)在加劇，沿海區(qū)域入侵生態(tài)風(fēng)險(xiǎn)形勢(shì)依然很?chē)?yán)峻。

3 結(jié)論與討論

2006—2016年間互花米草入侵生態(tài)風(fēng)險(xiǎn)時(shí)空格局變化不大，呈現(xiàn)一定的規(guī)律性。高風(fēng)險(xiǎn)區(qū)主要集中在溪南鎮(zhèn)、松港東部、下滸鎮(zhèn)和長(zhǎng)春鎮(zhèn)的西部、沙江鎮(zhèn)和三沙鎮(zhèn)的中部。低風(fēng)險(xiǎn)區(qū)主要分布在柏洋鄉(xiāng)、崇儒畬族鄉(xiāng)、水門(mén)畬族鄉(xiāng)、松城、牙城鎮(zhèn)西部，松港和鹽田畬族鄉(xiāng)的北部。

在整個(gè)研究期間，閩東濱海濕地互花米草入侵生態(tài)風(fēng)險(xiǎn)等級(jí)提高，但后5 a風(fēng)險(xiǎn)增加速度減緩，斑塊相對(duì)集中。極高生態(tài)風(fēng)險(xiǎn)面積不斷增加，10 a共增加12 583.00 hm2；高生態(tài)風(fēng)險(xiǎn)面積則不斷減少，共減少2 509.92 hm2；中生態(tài)風(fēng)險(xiǎn)面積10 a共增加4 433.04 hm2；低生態(tài)風(fēng)險(xiǎn)分布面積比例最大，面積先減少5 593.50 hm2后增加990.23 hm2，共減少4 603.28 hm2；極低生態(tài)風(fēng)險(xiǎn)面積不斷減少，10 a面積變化為9 912.20 hm2。

閩東濱海濕地互花米草入侵生態(tài)風(fēng)險(xiǎn)時(shí)空分布與產(chǎn)業(yè)結(jié)構(gòu)、人類(lèi)活動(dòng)和治理投入密切相關(guān)。因此，優(yōu)化區(qū)域產(chǎn)業(yè)結(jié)構(gòu)，避免不合理的沿海資源利用開(kāi)發(fā)方式和人為擾動(dòng)，積極落實(shí)互花米草的防控和治理，可有效防范研究區(qū)互花米草入侵生態(tài)風(fēng)險(xiǎn)，維持濱海濕地生態(tài)系統(tǒng)穩(wěn)定性。

至今，我國(guó)的生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)研究尚處于初始階段，統(tǒng)一、權(quán)威、公認(rèn)的生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)流程準(zhǔn)則還未出臺(tái)，學(xué)者多以有害生物風(fēng)險(xiǎn)分析作為理論研究基礎(chǔ)[34]。同時(shí)，對(duì)濕地互花米草入侵生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)指標(biāo)的研究還停留在概念性的設(shè)想和建議或簡(jiǎn)單的定量化階段，評(píng)價(jià)體系還有待進(jìn)一步完善。本研究中生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)指標(biāo)權(quán)重確定的方法為主觀(guān)性較強(qiáng)的層次分析法，在今后的研究中可選擇相對(duì)客觀(guān)的評(píng)價(jià)指標(biāo)方法，如灰色關(guān)聯(lián)、突變級(jí)數(shù)等。此外，在指標(biāo)的選擇上，不同因子層的指標(biāo)選擇存在部分的重合，如社會(huì)經(jīng)濟(jì)指標(biāo)、土地墾殖率、土地養(yǎng)殖率與農(nóng)漁牧業(yè)生產(chǎn)總值之間有明顯的相似性，未來(lái)研究中所選擇的指標(biāo)可盡量區(qū)分度明顯，以反映出不同的問(wèn)題。隨著人類(lèi)與生態(tài)環(huán)境之間相互聯(lián)系、相互依存、相互滲透的關(guān)系日益密切，評(píng)價(jià)指標(biāo)所涉及的范圍和尺度也亟需擴(kuò)大[35]，濕地互花米草入侵生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)指標(biāo)必將包括生態(tài)、社會(huì)、經(jīng)濟(jì)等綜合指標(biāo)，以準(zhǔn)確反映多壓力、多尺度、復(fù)雜的濕地生態(tài)系統(tǒng)。

互花米草入侵生態(tài)風(fēng)險(xiǎn)的空間分布格局受海拔、鹽度等自然因素的誘導(dǎo)及人為因素的影響較大。高、極高生態(tài)風(fēng)險(xiǎn)在研究區(qū)的東南部?jī)?nèi)海分布范圍最廣，該區(qū)土壤、水分、地形環(huán)境以及氣候環(huán)境相對(duì)適宜互花米草的生長(zhǎng)，且港灣眾多，在經(jīng)濟(jì)利益的驅(qū)使下，水域、潮間裸灘遭大規(guī)模圍墾和養(yǎng)殖，天然濕地面積大幅減少，生態(tài)系統(tǒng)的穩(wěn)定性被打破，致使互花米草入侵的機(jī)率增高。中生態(tài)風(fēng)險(xiǎn)主要集中在中部，該區(qū)域地勢(shì)復(fù)雜多樣，丘陵與低山、平原與盆谷錯(cuò)綜其間，適宜海灘潮間帶生長(zhǎng)的互花米草較難入侵該地，因此，中部地區(qū)可能是限制互花米草擴(kuò)張的緩沖區(qū)。低、極低生態(tài)風(fēng)險(xiǎn)集中于西北地區(qū)，西北地區(qū)峰巒聳峙，水分、土壤肥力等條件惡劣，其自然驅(qū)動(dòng)力因子較弱，且人口密度較低，人類(lèi)活動(dòng)不頻繁，該地區(qū)不易存在入侵生態(tài)風(fēng)險(xiǎn)。該結(jié)論與國(guó)內(nèi)眾多學(xué)者的相關(guān)研究結(jié)果基本一致，蘇少川等[1]、林立等[2]、劉明月[10]、歐文霞[36]認(rèn)為除自然環(huán)境對(duì)互花米草分布產(chǎn)生影響之外，人類(lèi)活動(dòng)對(duì)互花米草動(dòng)態(tài)變化影響亦較大，主要表現(xiàn)為人口因素、產(chǎn)業(yè)結(jié)構(gòu)調(diào)整、養(yǎng)殖池?cái)U(kuò)張、圍墾基建占用，只要人類(lèi)活動(dòng)壓力解除或減弱，生態(tài)系統(tǒng)尚能自我恢復(fù)。由此可見(jiàn)，需可持續(xù)開(kāi)發(fā)利用沿海資源，建立濱海濕地保護(hù)區(qū)，實(shí)現(xiàn)對(duì)耕地、裸灘、紅樹(shù)林群落景觀(guān)類(lèi)型的保護(hù)，維持濕地生態(tài)系統(tǒng)的穩(wěn)定。其次，應(yīng)加快閩東濱海濕地經(jīng)濟(jì)產(chǎn)業(yè)的轉(zhuǎn)移，加強(qiáng)發(fā)展第三產(chǎn)業(yè)，減少對(duì)第一產(chǎn)業(yè)的依賴(lài)。

治理互花米草是一個(gè)持久性、全面性的工程，霞浦縣通過(guò)采用物理、化學(xué)、生物防治的方法有效遏制了互花米草的蔓延，其中物理防治采用機(jī)械化處理，挖掘機(jī)作業(yè)除草；化學(xué)方法借鑒福建興華農(nóng)林高新技術(shù)研究所研制的“灘涂米草除控制”防治，效果達(dá)98%；生物防治采用種植秋茄抑制互花米草入侵。通過(guò)積極治理研究區(qū)互花米草入侵的高風(fēng)險(xiǎn)區(qū)，使減緩互花米草入侵速度，降低入侵生態(tài)風(fēng)險(xiǎn)達(dá)到一定成效。劉明月[10]、侯棟梁等[15]、張麗平等[37-38]的研究也證實(shí)了采取措施對(duì)互花米草進(jìn)行防控可達(dá)到較好的效果。因此，建議研究區(qū)充分考慮地區(qū)特色，因地制宜地開(kāi)展治理和防控，探索切實(shí)可行、科學(xué)的治理方式，降低互花米草的入侵危害。