基于優(yōu)化模型的萊州灣三疣梭子蟹(Portunus trituberculatus)棲息地適宜性研究*

張孝民 王秀霞 李少文 楊艷艷 徐炳慶 李 凡

(山東省海洋資源與環(huán)境研究院 山東省海洋生態(tài)修復(fù)重點(diǎn)實(shí)驗(yàn)室 山東煙臺 264006)

棲息地適宜性是對棲息地特征與物種在該條件下生存質(zhì)量的定量描述(易雨君等, 2013), 是當(dāng)前海洋生態(tài)保護(hù)領(lǐng)域重要的研究內(nèi)容之一, 對于生態(tài)養(yǎng)護(hù)、資源評估、種群動態(tài)研究都具有重要的理論和現(xiàn)實(shí)意義(鄒易陽等, 2016)。在人類活動不斷加劇和氣候變化的大背景下, 全球生物多樣性正處于持續(xù)的快速喪失中(徐煒等, 2016), 其中棲息地的減少已經(jīng)成為物種和局部種群絕滅的重要原因(鄧文洪, 2009; Keinathet al, 2017)。三疣梭子蟹(Portunus trituberculatus)生命周期短、繁殖力強(qiáng)、經(jīng)濟(jì)價(jià)值高(林群等, 2015), 捕撈數(shù)量占全世界經(jīng)濟(jì)蟹類的四分之一左右(Liuet al, 2013)。萊州灣位于渤海南部, 是赫赫有名的“萊州梭子蟹”的產(chǎn)地, 更有黃河和小清河等眾多河流入海, 灣內(nèi)基礎(chǔ)餌料豐富, 水溫鹽度適宜, 是包括蟹類在內(nèi)的眾多經(jīng)濟(jì)性海洋生物的產(chǎn)卵和棲息場所。近幾十年來, 受環(huán)境污染、過度捕撈等人類活動的影響, 萊州灣海洋生境遭到破壞, 眾多生物的棲息地出現(xiàn)破碎化現(xiàn)象, 部分海區(qū)呈現(xiàn)荒漠化趨勢。與其他海洋生物一樣, 三疣梭子蟹的棲息地也在不斷減少, 資源量持續(xù)衰退(盧曉等, 2018)。棲息地適宜性評價(jià)在內(nèi)陸河流及湖泊生物資源中運(yùn)用較為廣泛, 在海洋中的運(yùn)用起步較晚(龔彩霞等, 2011), 近年來研究集中在大洋魚類生物資源, 對于近海生物資源有較大的應(yīng)用空間, 國內(nèi)有關(guān)三疣梭子蟹棲息地適宜性的研究尚未見報(bào)道, 亟待開展相關(guān)研究。

HSI 模型是進(jìn)行棲息地適宜性定量評價(jià)的經(jīng)典方法, 該方法最早由美國地理調(diào)查局國家濕地研究中心魚類與野生生物署于20 世紀(jì)80 年代初提出(U.S.Fish and Wildlife Service, 1981)。GAM 模型常用來分析生物豐度和環(huán)境變量的非線性關(guān)系(Woodet al,2002), 可以在HSI 模型中選擇環(huán)境變量。BRT 模型基于決策樹理論的學(xué)習(xí)方法, 可有效計(jì)算出變量因子對模型建立的貢獻(xiàn)大小(De’ath, 2007), 生成多重回歸樹, 獲取變量的權(quán)重分布。

本文構(gòu)建了四種HSI 模型, 包括普通HSI 模型;GAM 優(yōu)化HSI 模型; BRT 優(yōu)化HSI 模型; GAM 和BRT優(yōu)化HSI 模型, 對萊州灣三疣梭子蟹棲息地適宜性進(jìn)行研究, 以期摸清資源分布狀況, 明確適宜的環(huán)境特征, 全面了解其棲息地生存質(zhì)量, 為三疣梭子蟹棲息地保護(hù)和萊州灣生態(tài)修復(fù)提供幫助。

1 材料與方法

1.1 數(shù)據(jù)來源和預(yù)處理

三疣梭子蟹數(shù)據(jù)來源于2010～2020 年夏季底拖網(wǎng)調(diào)查, 調(diào)查區(qū)域?yàn)槿R州灣海域(119°00′～120°15′E,37°12′～37°40′N), 調(diào)查站位每年20 個(gè)(圖1)。調(diào)查船只功率260 kW, 調(diào)查網(wǎng)具為單船底拖網(wǎng), 網(wǎng)口周長30.6 m, 囊網(wǎng)網(wǎng)目20 mm, 拖曳時(shí)網(wǎng)口寬度約8 m。

圖1 萊州灣三疣梭子蟹調(diào)查站位Fig.1 Sampling stations of P. trituberculatus in the Laizhou Bay

環(huán)境數(shù)據(jù)利用美國 ORI-ON520M-01A 便攜式Y(jié)SI 水質(zhì)分析儀現(xiàn)場測定, 內(nèi)容為水深(depth,D)、底層水溫(bottom sea temperature, BST)、底層葉綠素(bottom sea chlorophyll, BSC)、底層鹽度(bottom sea salinity, BSS)。大型底棲生物(Macrozoobenthos, Mzb)取樣使用開口面積為0.05 m2箱式采泥器采集, 每個(gè)站位取2 個(gè)平行樣, 使用0.05 mm 的網(wǎng)篩分選大型底棲動物, 樣品的保存、處理、計(jì)數(shù)和稱量等均按《海洋調(diào)查規(guī)范》GB/T 12763.6—2007 (中華人民共和國國家質(zhì)量監(jiān)督檢驗(yàn)檢疫總局等, 2008)進(jìn)行。

以調(diào)查網(wǎng)具拖速3 kn、拖網(wǎng)時(shí)間1 h 為基準(zhǔn), 根據(jù)各站位實(shí)際拖網(wǎng)時(shí)間對調(diào)查數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化處理,將其換算為單位時(shí)間的生物量(biomass, Bio) (g/h)和個(gè)數(shù)(individual, Ind) (ind./h), 以表征三疣梭子蟹的資源豐度。

1.2 研究方法

1.2.1 變量因子的選擇 GAM 是一種非參數(shù)化的多元線性回歸模型, 可以用來表示響應(yīng)變量和解釋變量的非線性關(guān)系, 本文利用GAM 模型選擇HSI 模型中的變量。GAM 模型通常表示為:

式中,g為鏈接函數(shù),μi=E(Yi),Xi為第i個(gè)響應(yīng)變量的解釋變量,β為模型估計(jì)參數(shù),Yi為第i個(gè)響應(yīng)變量,fi為平滑函數(shù)。將變量因子逐漸加入GAM 模型中, 選取AIC 值最小的模型作為最佳模型(Damalaset al,2007)。AIC 值的計(jì)算(Reid, 2003)如下:

式中,m為模型中參數(shù)的個(gè)數(shù),n為觀察值(數(shù)據(jù)樣本)的個(gè)數(shù), RSS 為殘差平方和。

1.2.2 變量因子權(quán)重分析 BRT 模型基于變量因子的相對貢獻(xiàn)確定其在HSI 模型中的權(quán)重(Changet al, 2010; Yiet al, 2016), BRT 方法是在運(yùn)算過程中隨機(jī)抽取一定量的數(shù)據(jù), 分析自變量對因變量的影響程度, 而剩余數(shù)據(jù)用來對擬合結(jié)果進(jìn)行檢驗(yàn), 對生成的多重回歸樹取平均值輸出, 從而得出自變量對因變量的相對重要性(Keinathet al, 2017)。本文使用R 4.0.3 中“GBM”包的“GBM.STEP”功能實(shí)現(xiàn)BRT 模型,其中抽樣率(Train. Fraction)設(shè)置為 0.8, 重復(fù)循環(huán)計(jì)算1 000 次, 獲取每個(gè)環(huán)境變量的相對貢獻(xiàn), 即環(huán)境變量的權(quán)重分布。

1.2.3 HSI 模型 首先, 利用樣條平滑函數(shù)建立變量因子與適宜性指數(shù)(suitability index, SI)的關(guān)系模型,分別以生物量Bio 和個(gè)數(shù)Ind 建立SI, 調(diào)查數(shù)據(jù)中Bio和Ind 最小值為0, 三疣梭子蟹在調(diào)查海域資源密度不均, 最大值和最小值量級差距較大, 直接以Bio 和Ind 數(shù)據(jù)建模會降低SI 指數(shù)精度, 本文對環(huán)境變量進(jìn)行分組, 以各組距范圍內(nèi)Bio 和Ind 的總和建立每個(gè)組距下的SI。SI 的取值范圍為0～1, 一般認(rèn)為0.7～1之間是適宜的棲息環(huán)境范圍(Brooks, 1997; Xueet al,2017), 計(jì)算公式如下:

式中, SIi為i環(huán)境變量的適宜性指數(shù), SIi,Bio為以生物量Bio 為基礎(chǔ)建立的適宜性指數(shù), SIi,Ind為以個(gè)數(shù)Ind為基礎(chǔ)建立的適宜性指數(shù), Bioij為i環(huán)境變量的j組距的Bio, Bioi,max為i環(huán)境變量的最大Bio, Indij為i環(huán)境變量的j組距的Ind, Indi,max為i環(huán)境變量的最大Ind。

HSI 模型是由SI 模型構(gòu)建的綜合模型, 取值范圍為0～1, 代表?xiàng)⒌剡m宜性“差”到“好”, 本文采用算數(shù)平均法(arithmetic mean model, AMM)和幾何平均法(geometric mean model, GMM)來構(gòu)建HSI 模型, 計(jì)算公式分別為:

1.2.4 HSI 模型的比較和驗(yàn)證 采用2010～2018 年數(shù)據(jù)構(gòu)建HSI 模型, 包括四種類型, 分別為: 普通HSI 模型; GAM 優(yōu)化HSI 模型(GAM 對變量進(jìn)行篩選后建立的HSI 模型); BRT 優(yōu)化HSI 模型(BRT 對變量權(quán)重進(jìn)行優(yōu)化建立的HSI模型); GAM和BRT優(yōu)化HSI模型。通過Pearson 檢驗(yàn)2020 年實(shí)測值HSI 和模型預(yù)測值HSI 的相關(guān)性, 采用ArcGIS 10.2 普通Kriging 插值法繪制2020 年萊州灣三疣梭子蟹HSI 地圖。

2 結(jié)果

2.1 變量因子的選擇和權(quán)重分析

以生物量表征資源豐度的GAM 模型(表1), 依次選取了D、BST、Mzb 密度和BSS 組成AIC 值最小的模型, 累計(jì)解釋方差為49.3%,D和BST 與生物量極顯著相關(guān)(P＜0.01), Mzb 密度和BSS 與生物量顯著相關(guān)(P＜0.05)。以尾數(shù)表征資源豐度的GAM 模型, 依次選取了D、BST、BSS 和Mzb 密度組成AIC 值最小的模型, 累計(jì)解釋方差為43.5%,D、BST 和BSS與生物量極顯著相關(guān)(P＜0.01), Mzb 密度與生物量顯著相關(guān)(P＜0.05)。

表1 以生物量和尾數(shù)表征資源密度構(gòu)建的GAM 模型Tab.1 P. trituberculatus abundance shown in the GAM model constructed by using biomass and mantissa

由圖2 可知, 以生物量和尾數(shù)表征資源密度構(gòu)建的BRT 模型變量因子的相對貢獻(xiàn)順序相同, BSS 的相對貢獻(xiàn)最大, 分別為26.9% (生物量)和23.5% (尾數(shù)),BSC 的相對貢獻(xiàn)最小, 分別為 12.2% (生物量)和11.4% (尾數(shù)), 其他因子的相對貢獻(xiàn)大小依次為BST、D和Mzb 密度, 區(qū)間范圍為19.0%～23.3%。

圖2 以生物量和尾數(shù)表征資源密度構(gòu)建的提升回歸樹模型Fig.2 P. trituberculatus abundance shown in the BRT model constructed in biomass and mantissa

2.2 變量因子與生物量和尾數(shù)的適宜性指數(shù)關(guān)系模型

從變量因子的適宜性指數(shù)曲線(圖3)可以看出,以生物量和尾數(shù)表征資源豐度擬合的三疣梭子蟹適宜棲息環(huán)境范圍(SI＞0.7)基本相同。BST 的適宜范圍為27.2～28.5 °C, BSS 的適宜范圍為25.5～30.5,D的 適宜范圍為 9.0～10.5 m, Mzb 密度的適宜范圍為2 900～3 100 個(gè)/m2, BSC 濃度的適宜范圍為 13.5～14.5 mg/m3。

圖3 以生物量和尾數(shù)表征資源豐度與環(huán)境因子的適宜性指數(shù)曲線Fig.3 Fitted suitability index (SI) curves based on the relationship between P. trituberculatus abundance in terms of biomass and mantissa and environmental factors

2.3 HSI 模型的比較和驗(yàn)證

以生物量表征資源豐度的GAM 和BRT 優(yōu)化算術(shù)平均模型, 與 2019～2020 實(shí)際調(diào)查的生物量數(shù)據(jù)Pearson 相關(guān)系數(shù)最高, 為0.653, 表明預(yù)測值與實(shí)測值有較強(qiáng)的相關(guān)性, 最優(yōu)模型具有較強(qiáng)的預(yù)測能力?？傮w看來, 優(yōu)化后的HSI 模型均好于普通HSI 模型,BRT 優(yōu)化的模型略好于GAM 優(yōu)化的模型, GAM 和BRT 優(yōu)化的模型好于其他3 種模型(表2)。此外, 生物量表征資源豐度的模型好于尾數(shù)表征資源豐度的模型,算術(shù)平均模型的相關(guān)系數(shù)高于幾何平均模型。

2.4 萊州灣三疣梭子蟹棲息地空間分布

根據(jù)GAM 和BRT 優(yōu)化的算術(shù)平均HSI 模型預(yù)測結(jié)果可以看出, 伏季休漁結(jié)束前三疣梭子蟹主要分布在萊州灣南部、東南部和東北部海域, 中部和西南部海域分布較少, HSI 值較低(圖4)。模型預(yù)測結(jié)果與2020 年實(shí)際調(diào)查的三疣梭子蟹分布情況大體相符,以生物量表征資源豐度的HSI 模型預(yù)測的三疣梭子蟹最適海域(HSI＞0.7)范圍明顯大于以尾數(shù)表征資源豐度的HSI 模型, 以生物量表征資源豐度的HSI 模型對資源豐度較高區(qū)域的預(yù)測更加準(zhǔn)確。

圖4 2020 年萊州灣三疣梭子蟹棲息地優(yōu)化模型預(yù)測結(jié)果與實(shí)際調(diào)查情況分布圖Fig.4 Spatial distribution of observed abundance of P. trituberculatus in 2020 overlaid with the predicted HSI values based on optimized HSI models

3 討論

3.1 HSI 模型的優(yōu)化

生物的棲息地受多種因子的影響, 不同因子的影響大小不同, 合理選擇影響因子和設(shè)置各因子的權(quán)重是優(yōu)化HSI 模型的兩種主要方法。本研究分別采用GAM 模型對三疣梭子蟹影響因子進(jìn)行選擇和BRT模型設(shè)置因子權(quán)重, 結(jié)果(表2)顯示BRT 優(yōu)化模型好于GAM 優(yōu)化模型, 說明在對HSI 模型優(yōu)化的過程中,權(quán)重設(shè)置比因子選擇更加重要, 與Zhang 等(2020)研究結(jié)果一致。通過查閱文獻(xiàn)資料基本可以囊括研究對象棲息地的大部分影響因子, 加上專家知識的判斷篩選重要因子, 可以比較準(zhǔn)確的完成因子選擇的過程(龔彩霞等, 2011), 但權(quán)重設(shè)置受不同研究海域和不同研究對象的影響較大, 專家無法獲得足夠的信息對不同因子權(quán)重準(zhǔn)確設(shè)置(Tianet al, 2009), 本研究利用BRT 模型設(shè)置因子權(quán)重, 與非優(yōu)化模型相比更具科學(xué)性和針對性, 將因子影響程度大小定量表示, 這也是導(dǎo)致權(quán)重設(shè)置比因子選擇優(yōu)化效果更好的原因。但仍不能忽視GAM 模型對因子選擇的重要性, 綜合利用GAM 模型和BRT 模型才能構(gòu)建最佳的HSI 模型。

表2 不同HSI 模型預(yù)測值與2019～2020 年實(shí)際調(diào)查數(shù)據(jù)的Pearson 相關(guān)檢驗(yàn)結(jié)果Tab.2 Results of the Pearson correlation in different HSI model-predicted values and the actual survey data in 2019～2020

研究對象資源豐度指標(biāo)的選取也是影響HSI 模型建模效果的重要因素, 主要的資源豐度指標(biāo)包括漁獲量(范秀梅等, 2020)、漁獲尾數(shù)、作業(yè)次數(shù)、單位捕撈努力量漁獲量(catch per unit effort, CPUE) (劉瑜等, 2021)、生物量和漁獲密度(朱承之等, 2020)等,相關(guān)研究大都采用單一指標(biāo)描述資源豐度, 本文比較了生物量和尾數(shù)表征三疣梭子蟹資源豐度對構(gòu)建HSI 模型的影響, 發(fā)現(xiàn)以生物量表征資源豐度的模型好于尾數(shù), 在GAM 模型的累計(jì)解釋率也相對較高,可能的原因是三疣梭子蟹個(gè)體大小差距較大, 部分站位大型個(gè)體較多, 用尾數(shù)表征降低了該區(qū)域的資源豐度, 但在適宜性指數(shù)曲線擬合得到的適宜棲息環(huán)境范圍基本相同。范秀梅等(2020)分別基于漁獲量和作業(yè)次數(shù)表征日本鯖魚的資源豐度建立HSI 模型,發(fā)現(xiàn)兩種方法在預(yù)報(bào)漁場上精度相差不大, 原因是日本鯖漁船作業(yè)次數(shù)和漁獲量通常成正比關(guān)系, 作業(yè)次數(shù)反映了漁船的集中程度和漁民對產(chǎn)量的滿意度(Bordalo-Machado, 2006), 作業(yè)次數(shù)多表明該海域有很高的漁獲量。Tian 等(2009)比較了基于CPUE 和作業(yè)次數(shù)建立的西北太平洋柔魚HSI 模型, 結(jié)果表明,基于CPUE 的HSI 模型高估了柔魚的最佳棲息地范圍,基于作業(yè)次數(shù)的HSI 模型能更好地定義柔魚的棲息地范圍。因此, 在選取資源豐度指標(biāo)時(shí), 應(yīng)充分考慮研究對象的生物學(xué)以及作業(yè)方式等特性, 并采用多種指標(biāo)表征資源豐度, 以獲得更好的HSI 模型。此外,HSI 模型的計(jì)算方法中, 算術(shù)平均法AMM 普遍具有較高的效果(陳新軍等, 2009), 本文研究結(jié)果也證明了AMM 的優(yōu)勢, 但相關(guān)研究表明GMM 對SI 的極大值和極小值較敏感(Vinagreet al, 2006)。

3.2 三疣梭子蟹棲息地的影響因素

萊州灣三疣梭子蟹冬季在渤海10～30 m 水深較深海域越冬, 4 月開始生殖洄游, 到3～5 m 深的港灣或河口等淺海附近海域產(chǎn)卵, 8 月大多數(shù)三疣梭子蟹度過產(chǎn)卵期, 進(jìn)入索餌育肥期, 是三疣梭子蟹資源量最高的時(shí)期, 9 月伏季休漁結(jié)束后, 高強(qiáng)度的捕撈會使資源量迅速下降, 人為因素成為主導(dǎo)三疣梭子蟹棲息地分布的決定因素, 環(huán)境因素影響力下降(Tianet al, 2009), 因此8 月伏季休漁結(jié)束前是研究三疣梭子蟹棲息地適宜性的最佳時(shí)期。本文研究了影響萊州灣三疣梭子蟹的部分環(huán)境因子和生物因子, GAM 模型和BRT 模型結(jié)果均表明BST、BSS 和D對其棲息地影響較大, 與吳強(qiáng)等(2016a)研究結(jié)果相似。三疣梭子蟹是廣鹽性蟹類, 但其 BSS 的最適范圍為25.5～30.5 (圖3), 夏季雨水豐富, 萊州灣淡水入海量增加, 近海鹽度變化成為影響三疣梭子蟹棲息地的最主要因素(圖2)。三疣梭子蟹最適棲息地的BST 范圍為27.2～28.5 °C, 水溫不僅是決定三疣梭子蟹的季節(jié)性洄游的主要因素, 夏季是索餌育肥的重要時(shí)期,水溫過高會降低其體內(nèi)消化酶的活性, 影響攝食行為(廖永巖等, 2008), 進(jìn)而對棲息地分布產(chǎn)生影響。萊州灣大部分區(qū)域水深小于10 m, 最大水深18 m, 水深較淺的近岸水域是三疣梭子蟹產(chǎn)卵棲息的重要場所, 水深的變化會引起水溫、鹽度和營養(yǎng)鹽等的改變,對棲息地分布產(chǎn)生影響, 本文研究發(fā)現(xiàn)伏季休漁結(jié)束前三疣梭子蟹在萊州灣的最適水深為9.0～10.5 m。相關(guān)研究表明三疣梭子蟹分布與底棲生物的關(guān)系密切(吳強(qiáng)等, 2016b), 本文研究認(rèn)為 Mzb 密度為 2 900～3 100 個(gè)/m2時(shí)是三疣梭子蟹最適范圍, 底棲生物是大型甲殼類的重要食物來源, 三疣梭子蟹與日本蟳、口蝦蛄和中國明對蝦之間存在較強(qiáng)的食物競爭關(guān)系(李凡等, 2021), 導(dǎo)致底棲生物密度成為限制其棲息地分布的影響因素, 夏季萊州灣餌料食物較為豐富, 底棲生物與水溫和鹽度的影響相比較小, 在食物缺乏季節(jié)或海域, 底棲生物的影響力將增強(qiáng)。

3.3 夏季萊州灣三疣梭子蟹棲息地空間分布

結(jié)合實(shí)際調(diào)查數(shù)據(jù)以及HSI 模型的預(yù)測結(jié)果, 夏季萊州灣整體表現(xiàn)為近岸海域比離岸海域更適合三疣梭子蟹棲息, 這與徐炳慶等(2021)的研究結(jié)果一致,8～9 月份伏季休漁結(jié)束前水溫較高, 三疣梭子蟹尚未進(jìn)入越冬洄游階段, 集中分布于近岸海域。本研究調(diào)查船只為底拖網(wǎng)漁船, 無法在水深低于5 m 海區(qū)設(shè)置站位, 如采用小型船只或其他作業(yè)方式對淺海三疣梭子蟹資源進(jìn)行調(diào)查, 將使近岸資源好于離岸資源的趨勢更加明顯。2020 年萊州灣西南部近岸海域在HSI 模型預(yù)測中棲息地適宜性較差, 主要原因?yàn)?211站位和7214 站位底層平均水溫30 °C, 底層平均鹽度23.8, 均處于適宜性指數(shù)曲線SI 較低的位置, 導(dǎo)致HSI 值較小, 成為不適宜三疣梭子蟹棲息的海域。但在2020 年的實(shí)際調(diào)查中, 7214 站位生物量較高, 為2 196 g/h, 可能的原因是站位設(shè)置中環(huán)境數(shù)據(jù)是與三疣梭子蟹資源數(shù)據(jù)同步, 近岸海域溫鹽環(huán)境變化較大, 站位設(shè)置稀疏使得環(huán)境數(shù)據(jù)的獲取出現(xiàn)偶然性,在以后的調(diào)查研究中, 應(yīng)當(dāng)對近岸海域環(huán)境監(jiān)測站位設(shè)置更加緊密, 以獲取更加準(zhǔn)確地?cái)?shù)據(jù), 提高棲息地預(yù)測的精確性。

4 結(jié)論

本研究運(yùn)用GAM 模型和BRT 模型在建模過程中對HSI 模型進(jìn)行了優(yōu)化, 優(yōu)化模型對適宜棲息地預(yù)測效果均好于未優(yōu)化模型, 優(yōu)化模型對萊州灣三疣梭子蟹棲息地適宜性有較強(qiáng)的預(yù)測能力。優(yōu)化模型顯示BST、BSS 和D對三疣梭子蟹棲息地影響較大, 萊州灣南部、東南部和東北部是夏季三疣梭子蟹的適宜棲息海域。