基于栖息地适宜性指数模型及最大熵模型的弧边招潮蟹生境适宜性比较

周鑫,张丽源,何晨翔,胡成业,王晶,李艺,李雪丽,水柏年*

(浙江海洋大学 水产学院,浙江 舟山 316022)

摘要:为研究弧边招潮蟹(Uca arcuata)生境适宜性及最适分析模型,基于2020年7月(夏季)、10月(秋季)及12月(冬季)对浙江省温州市沿浦湾秋茄林大型底栖动物和若干环境因子的调查,采用栖息地适宜性指数(habitat suitability index,HSI)模型与最大熵(maximum entropy model,MaxEnt)模型对弧边招潮蟹生境适宜性进行比较。结果表明:采用HSI模型分析显示,各采样站位弧边招潮蟹的栖息地适宜性总体上呈现近岸段>中段>离岸段及东南部>西北部的空间分布特征,3个季节的栖息地适宜性两两之间存在显著性差异(P<0.05),总体上呈现夏季>秋季>冬季的季节分布特征;采用MaxEnt模型分析显示,各采样站位弧边招潮蟹的栖息地适宜性总体上呈现近岸段>中段>离岸段及东南部>西北部的空间分布特征,并呈现秋季>夏季>冬季的季节分布特征。研究表明,HSI模型较MaxEnt模型更适宜用于弧边招潮蟹适生区的调查分析,温度、底质、地质和潮汐等是影响弧边招潮蟹分布的关键环境因子,硝酸盐、亚硝酸盐和盐度等是影响弧边招潮蟹分布的关键理化因子,这些因子对弧边招潮蟹分布格局具有较强的驱动力,影响其适宜生境的选择。

关键词弧边招潮蟹;适宜性指数模型;最大熵模型;生境适应性

大型底栖动物作为红树林生态系统的重要组成部分,是物质循环和能量流动的重要消费者与转移者。同时,大型底栖动物活动能力有限,对逆境的逃避相对迟缓,但对红树林栖息地环境变化较为敏感[1]。其时空分布与群落结构不仅可以表征红树林栖息地环境质量状况,而且可以反映栖息地的适宜性特征[2],因此,大型底栖动物常作为评价红树林生态系统健康程度的指示物种[3-4]

关于生物对生境适应性的研究,国内外学者常采用物种分布模型(species distribution model,SDM),即以特定的算法将目标物种现有的分布信息及环境数据进行计算拟合的数学模型[5]。Roberto等[6]研究发现,物种分布模型中的最大熵(MaxEnt)模型在预测生物潜在分布区域时模拟精度较高,尤其是物种分布数据不全的情况下,相较于其他模型具有更好的性能表现。孙元超等[7]采用物种分布模型中的MaxEnt模型分析预估了纹沼螺(Parafossarulus striatulus)在中国的适生区,以及气候变化对纹沼螺分布的影响,指出海拔是影响纹沼螺分布的最主要环境因素,并提出了未来可建立针对螺类的系统监测体系等建议。总之,MaxEnt模型因精度较高、结果准确和操作方便等优点,已广泛应用于多种生物的适生区研究,且均证实其具有较高的预测性能,在保护物种多样性、防止外来物种入侵方面发挥着重要作用。

栖息地适宜性指数(habitat suitability index,HSI)模型最早由美国地理调查局国家湿地研究中心鱼类与野生生物署于20世纪80年代提出,主要用于描述野生动物栖息地质量,可预测和模拟生物的分布特征,评价栖息地环境变化对生物造成的影响。近年来,由于HSI模型在推断指示生物的最适环境条件、预测生物分布、估算生物量及评价栖息地适宜性等方面具有独特的优越性[8],其在物种管理、物种生境偏好评估[9]等研究中被广泛应用。

Vincenzi等[10]将GIS技术与栖息地适宜性模型相结合,构建菲律宾蛤仔(Ruditapes philippinarum)的HSI模型,对不同地区的菲律宾蛤仔栖息地适宜性开展了研究,有利于渔业生物资源的合理利用及有效分配。郑文浩等[11]和粟一帆等[12]分别对太子河及桥边河等内陆河流域的大型底栖动物栖息地适宜性开展了研究。综上,HSI模型在生物生境适宜性等方面具有独特的优越性,被广泛应用于内陆河流域生物栖息地适宜性的研究中,而针对海洋生物栖息地适宜性的研究主要围绕游泳动物展开,但关于红树林大型底栖动物的生境适宜性研究及报道较为罕见。

本文中,基于2020年7月(夏季)、10月(秋季)和12月(冬季)对浙江省温州市沿浦湾秋茄林大型底栖动物和环境因子的调查,开展3个季节大型底栖动物共有优势种弧边招潮蟹(Uca arcuata)生境适宜性研究,采用HSI模型及MaxEnt模型探究了影响其生境适宜性变化的关键环境因子,并比较了两种生境适宜性模型的预测性能、适用范围和局限性,以期为红树林生态系统物种分布及生物多样性保护等方面的应用研究提供理论支持,为生态环境治理与决策提供重要的科学依据。

1 材料与方法

1.1 研究区域

沿浦湾位于浙江省温州市苍南县沿浦镇(27.20°N、120.45°E),湾口东起霞关镇,西至虎头鼻,宽约为3 000 m,内延伸6 000 m,面积为19.235 km2,平均水深为2 m;属亚热带海洋性季风气候区,根据苍南县气象台多年统计资料,该区域多年平均气温为18.2 ℃,历史最冷月(1月)平均气温为7.9 ℃,累年极端最低气温为-2.2 ℃(1973年12月25日)。沿浦湾有下在河、沿浦河和岭尾河淡水注入,底质以软泥粉砂为主,2015年开始引种秋茄,至今已经引种达到100 km2以上,秋茄植株高度为2~3 m,长势良好。

1.2 方法

1.2.1 调查站位设置 依据《海洋调查规范》[13]和《海洋监测规范》[14],于2020年7月(夏季)、10月(秋季)和12月(冬季)在沿浦湾秋茄林区布设15个采样点,进行大型底栖动物调查(图1)。

图1 温州市沿浦湾采样点分布

Fig.1 Distribution of sampling sites in Yanpu Bay in Wenzhou City

1.2.2 大型底栖动物样品采集 采集的大型底栖动物样本用0.5 mm孔径的网筛淘洗后,装入样品瓶中,用体积分数为75%的乙醇溶液固定,带回实验室进行种类鉴定、计数及称重,并做好数据记录。

1.2.3 环境因子监测 根据《海洋监测规范》[14]的相关要求,采用便携式水质检测仪(YSI 2000)现场测定溶解氧、pH、温度、电导率和盐度,并做好数据记录。然后采集500 mL间隙水水样置于-20 ℃的车载冰箱中带回实验室,测定其氨氮、亚硝酸盐、硝酸盐及活性磷酸盐等水体生境指标。

1.2.4 HSI模型

1)栖息地指示因子选取。为避免各变量间存在多重共线性而导致模型过度拟合,在筛选环境变量前先对各变量进行相关性分析,剔除相关系数≥0.8的变量[15]

GAM模型是广义线性模型的非参数化扩展,其优点是能直接处理响应变量与多个解释变量间的非线性关系[16]。以溶解氧、pH、温度、电导率、盐度、氨氮、亚硝酸盐、硝酸盐及活性磷酸盐等9个环境因子作为自变量,以指示物种的丰度作为响应变量建立GAM模型,分析生物丰度与环境因子间的关系,进而评估沿浦湾秋茄林的栖息地适宜性。GAM模型的表达式为

(1)

式中:Y为丰度,为避免丰度为0的情况,本研究中采用丰度加1,再进行对数化处理的方法;α为截距项;n为解释变量的个数;fi(xi)为解释变量样条的平滑函数;ε为误差项。为避免过度拟合,所有的平滑被限制在4节[17]。赤池信息准则(akaike inforamation criterion,AIC)检验可明确添加不同环境因子后模型的拟合程度,AIC值越小,表示模型的拟合效果越好[18]

2)HSI模型建立。HSI模型的不同计算方法对其预测结果的影响不同,本研究中分别采用加权平均法(weighted average model,WAM)和几何平均法(geometric mean model,GMM)计算HSI模型,通过评价两种算法的性能,选取预测性效果最好的算法进行HSI模型的运算。计算公式如下:

(2)

(3)

式中:SHSI为HSI指数;wi为第i个栖息地环境因子的权重,权重总和为1,使用相对权重(relative weight)计算法确定,该方法是对所有子模型添加一个自变量引起的R2平均增加量的近似值;Si表示第i个栖息地因子的适宜度,表示大型底栖动物对该栖息地环境因子的适宜性程度,首先基于生物丰度计算适宜性指数,然后采用平滑回归模型拟合环境因子与丰度适宜性指数的关系,最后得出各环境因子的适宜度[19],取值为0~1,当其处于0.7~1.0时,对应的环境变量为最适环境范围。

栖息地的生境一般分为:极佳生境(0.75[20]。

3)HSI模型检验。采用交叉验证法评估HSI模型的预测性能。随机选取80%的调查数据建立模型,将剩余20%的数据代入模型进行交叉验证,该过程重复进行100次,从而评估HSI模型的预测准确性[21]。对采用WAM和GMM算法的HSI模型,使用GAM模型中的AIC值评估预测其综合性能,选取适用于本研究的最优模型,采用克里金插值绘制HSI时空分布图。分别采用SPSS 25、R 4.0.4和ArcGIS 10.7软件完成相关性分析、GAM模型与平滑回归模型的构建及克里金插值。

1.2.5 MaxEnt模型

1)MaxEnt模型。MaxEnt模型是一种基于生态位原理预测物种潜在分布的模型。该模型主要依据物种“出现点”的环境变量特征计算约束物种分布的条件,估算此约束条件下最大熵的可能分布,并据此推测目标物种的适宜栖息地分布[22]

对于目标地区XX其实为有限数量的网格(grid),记π(x)为π分布赋予各网格的概率,则∑π(x)=1,对π的估计分布熵为

(4)

熵是非负的,最多为X中元素数量的自然对数。最大熵原理可以被解释为上不应该存在毫无根据的约束[23]

2)MaxEnt模型数据处理。应用MaxEnt模型对大型底栖动物进行生境适宜性评价需要两方面的数据:一是大型底栖动物的地理分布数据;二是沿浦湾秋茄林区的环境变量数据。

大型底栖动物的物种分布数据来源于2020年7月、10月和12月的实地调查数据。将采集到生物样品的采样点数据均按照species、longitude和latitude格式依次输入Excel表格中,将其保存为MaxEnt模型支持的.csv格式文件备用。

沿浦湾秋茄林区的环境变量数据来源于实地调查及实验室测定数据,主要包括夏季、秋季及冬季的溶解氧、pH、温度、电导率、盐度、氨氮、亚硝酸盐、硝酸盐及活性磷酸盐含量数据。将所有环境因子数据采用ArcGIS 10.7软件转换为ASCII格式文件,并对结果进行评价[24]

3)MaxEnt模型建立。将大型底栖动物地理分布数据和环境因子数据导入MaxEnt模型中,随机分别选取75%的大型底栖动物分布数据作为训练集用于建立模型,剩余25%的分布数据作为测试集用于模型验证,模型重复运算10次,10次交叉检验的各项平均值作为模型最终的计算结果。选用Jackknife中的AUC(area under curve,AUC)评价指标对模型结果和环境变量的重要性及模型的准确度。

AUC为模型自带的受试者工作特征曲线(receiver operating characteristic,ROC)下的面积。其评价标准[22]:AUC=0.5~0.6为不及格;AUC=0.6~0.7为较差;AUC=0.7~0.8为一般;AUC=0.8~0.9为良好;AUC=0.9~1.0为优秀。

4)栖息地适宜性评价。将模型输出结果的.asc格式的图层导入ArcGIS并转换成栅格数据,栅格的数值代表大型底栖动物在此区域的生存概率。按照栅格数值的范围可将大型底栖动物的生境分布区划分为4个等级:0~0.25为不适宜生境;0.25~0.50为次适宜生境;0.50~0.75为适宜生境;0.75~1.0为极佳生境[2]

2 结果与分析

2.1 物种选择

根据2020年夏季、秋季、冬季3个季节大型底栖动物物种分布数据分析得出,弧边招潮蟹在3个季节中具有良好代表性,同时,弧边招潮蟹选择适应性生境的活动能力较强,故本试验中以弧边招潮蟹为例,采用HSI模型和MaxEnt模型进行融合分析,以揭示弧边招潮蟹的生境适应性特征。

2.2 基于HSI模型的弧边招潮蟹生境适应性

2.2.1 栖息地指示因子选取 以优势种弧边招潮蟹的丰度为因变量,不同环境因子组合为自变量,利用GAM模型筛选影响弧边招潮蟹适宜性的环境因子。分别将各环境因子组成不同的组不断添加至GAM模型中,观察不同环境变量组合形式下AIC值的变化情况,以具有较低AIC值的自变量组合作为解释变量。结果显示,硝酸盐+亚硝酸盐+盐度组合的AIC值最低(表1),因此,后续使用上述环境因子的组合构建HSI模型。

表1 所有环境因子逐级拟合后AIC最小组合的GAM模型检验

Tab.1 GAM model test of the combination with the minimum AIC value after fitting of all environmental factors step by step

注:本研究中GAM模型是对所有环境因子的所有组合形式的拟合,但因组合过多,本表中仅展示相同环境因子数量下AIC最小的组别,并将其进行比较。

Note:The GAM model in this study is a fitting of all combination forms of all environmental factors.However,only the group with the minimal AIC in the same number of environment factors are shown and compared due to too many combinations.

组别 groupAIC值 AIC value硝酸盐+亚硝酸盐 nitrate+nitrite105.219 90硝酸盐+亚硝酸盐+盐度 nitrate+nitrite+salinity98.271 68硝酸盐+亚硝酸盐+盐度+氨氮 nitrate+nitrite+salinity+ammonia nitrogen99.884 11

利用相对权重法,分析得到影响弧边招潮蟹生境适宜性的各主要环境变量权重值,其中,硝酸盐的权重最大(56.6%),盐度及亚硝酸盐等环境因子的权重分别为26.3%及17.1%。

2.2.2 栖息地适宜度曲线 根据单因子适宜度曲线,弧边招潮蟹最适硝酸盐质量浓度为1.8~2.5 mg/L,最适盐度≤16,最适亚硝酸盐质量浓度为0.08~0.021 mg/L(图2)。

图2 单环境因子适宜度曲线

Fig.2 Suitability curve of single environmental factor

2.2.3 HSI指数 从表2可见,采用交叉检验的方法评估WAM与GMM模型的预测性能发现,运用GMM算法构建的HSI模型的AIC值均较低,故选取该算法研究弧边招潮蟹的栖息地适宜性。

表2 WAM和GMM模型预测性能比较(弧边招潮蟹)

Tab.2 Comparison of predictive performance between WAM and GMM(Uca arcuata)

季节seasonWAMGMMAICPAICP夏季summer37.875 520.00028.666 970.000秋季autumn40.448 950.00039.476 940.000冬季winter38.427 520.00323.753 980.000

从图3可见,将各站位指示环境的单因子适宜度代入HSI模型,分析各采样站位的HSI值发现,各采样点弧边招潮蟹的栖息地适宜性总体上呈现近岸段>中段>离岸段,东南部>西北部,夏季>秋季>

图3 基于HSI模型的沿浦湾弧边招潮蟹栖息地适宜区分布

Fig.3 Distribution of suitable areas of Uca arcuata in Yanpu Bay based on HSI model

冬季的时空分布特征。

从图4可见:夏季、秋季和冬季的HSI指数之间均有显著性差异(P<0.05),夏季的HSI指数显著高于秋季及冬季(P<0.05),秋季的HSI指数极显著高于冬季(P<0.01);夏季、秋季和冬季的HSI指数分别为0.09~1.00、0.11~0.70和0~0.31。

*表示组间有显著性差异(P<0.05);**表示组间有极显著性差异(P<0.01)。

*means significant difference in the groups(P<0.05);**means very significant difference in the groups(P<0.01).

图4 夏季、秋季和冬季弧边招潮蟹HSI的差异检验

Fig.4 HSI difference test of Uca arcuata in summer,autumn and winter

2.3 基于MaxEnt模型的弧边招潮蟹生境适应性

2.3.1 MaxEnt模型适用性评价 由ROC曲线的评价结果可知:夏季弧边招潮蟹训练AUC平均值为0.750,表明MaxEnt模型的预测精度一般,其中,盐度、亚硝酸盐及氨氮对模型的贡献率较大;秋季弧边招潮蟹训练AUC平均值为0.774,表明MaxEnt模型的预测精度一般,其中,溶解氧、盐度及硝酸盐对模型的贡献率较大;冬季弧边招潮蟹训练AUC平均值为0.845,表明MaxEnt模型的预测精度良好,其中,氨氮、硝酸盐、活性磷酸盐、亚硝酸盐及盐度对模型的贡献率较大(图5)。

图5 夏季、秋季和冬季弧边招潮蟹生境分布预测结果的ROC曲线验证

Fig.5 ROC curve verification of habitat distribution prediction results of Uca arcuata in summer,autumn and winter

2.3.2 弧边招潮蟹生境适宜性分布 从图6可见,沿浦湾秋茄林弧边招潮蟹的栖息地适宜性总体上呈现近岸段>中段>离岸段,秋季>夏季>冬季,东南部>西北部的时空分布特征。

图6 基于MaxEnt模型的沿浦湾弧边招潮蟹栖息地适宜区分布

Fig.6 Distribution of suitable areas of Uca arcuata in Yanpu Bay based on MaxEnt model

利用ArcGIS计算弧边招潮蟹适宜性等级的面积,结果显示,在总面积1 km2以上秋茄林评价区域内,7月(夏季)适宜生境多分布于近岸段及东南部,约占秋茄林区总面积的17.92%,次适宜生境面积约占总面积的80.78%,不适宜生境面积约占总面积的1.30%;10月(秋季)极佳生境主要分布在近岸段及东南部,约占总面积的0.49%,适宜生境面积约占总面积的12.98%,次适宜生境面积约占总面积的86.53%;12月(冬季)适宜生境面积约占总面积的3.89%,次适宜生境面积约占总面积的80.89%,不适宜生境面积约占总面积的15.22%(图7)。

图7 弧边招潮蟹不同生境类型的面积预测

Fig.7 Prediction of the area of different habitat types of Uca arcuata

3 讨论

3.1 两种模型预测生物栖息地适宜性的性能比较

HSI和MaxEnt两种模型在分析生物对生境适宜性的研究上各有优势,MaxEnt模型在预测生物潜在分布区域时模拟精度较高,HSI模型更适合推断指示生物的最适环境条件、预测生物分布、估算生物量及评价栖息地适宜性等方面。

本研究中,将两种模型同时对弧边招潮蟹栖息地适宜性进行分析,结果显示,用HSI模型预测的弧边招潮蟹的栖息地适宜性总体上表现为近岸段>中段>离岸段,夏季>秋季>冬季,东南部>西北部的时空分布特征。用MaxEnt模型预测的弧边招潮蟹栖息地适宜性总体上表现为近岸段>中段>离岸段,秋季>夏季>冬季,东南部>西北部的时空分布特征。

实际调查得出,弧边招潮蟹丰度的时空分布总体呈现近岸段>中段>离岸段,夏季>秋季>冬季的特征,其空间分布特征与本研究中HSI模型及MaxEnt模型的预测结果均基本吻合,但两模型预测的弧边招潮蟹季节分布特征存在差异。其中,HSI模型的预测结果与实际分布规律基本相符,且生境分布的细节更符合实际。MaxEnt模型的预测结果与实际分布存在一定差异,而且其在夏季及秋季对弧边招潮蟹的预测结果中AUC值均为一般水平,说明MaxEnt模型对弧边招潮蟹的预测效果精度偏低。MaxEnt模型更适宜于大尺度范围适生区的调查研究,范围过小有可能出现过度拟合的问题,进而导致其预测精度降低。杨刚等[25]、刘振生等[26]的研究也得出了类似的结论。因此,本研究中采用HSI模型预测弧边招潮蟹生境适宜性比MaxEnt模型具有更大的可靠性。

3.2 弧边招潮蟹适宜栖息地分布

弧边招潮蟹所选取的适宜栖息地与其生活习性密不可分,在选取HSI模型作为最适合本研究的模型后,对弧边招潮蟹所适宜的栖息地进行分析,预测显示,弧边招潮蟹的时空分布特征主要表现为近岸段适应性高于中段,更高于离岸段,而夏季则高于秋季,更高于冬季。这可能是由弧边招潮蟹营穴居生活特性决定的。弧边招潮蟹挖筑洞穴的深度与被捕食风险密切相关,Klaassen等[27]研究发现,蟹类在基质中挖洞至少16 cm才可有效躲避鸟类捕食者。另外,涨潮期间,体型较大的雄性会向雌性求偶,雌蟹会选择拥有更合适洞穴的雄蟹,并在洞穴中交配繁殖。洞穴作为繁殖地其结构的完整性及稳定性很大程度上决定了繁殖的质量。但离岸段有氧层非常薄,不可能深入挖掘,近岸段和中段有氧层较厚,且受潮水浸没的时间较短,土质硬化程度较高,有助于弧边招潮蟹挖筑更深更稳固的洞穴,且可避免因潮汐淹没或由潮汐驱动的地下水上升而造成洞穴坍塌,故近岸段及中段更适宜弧边招潮蟹筑穴、栖息及繁衍。此外,弧边招潮蟹对水温变化较敏感,李明云[28]研究发现,冬季弧边招潮蟹几乎穴居洞内无外出活动,其活动能力会随着温度的升高而逐渐增强,温度高时比温度低时的出洞次数多、活动强,这可能是夏季弧边招潮蟹生境适宜性较高的原因。

3.3 弧边招潮蟹栖息地适宜性的影响因素

环境因子是影响弧边招潮蟹栖息地适宜性的重要因素。本研究中通过对溶解氧、pH、温度、电导率、盐度、氨氮、亚硝酸盐、硝酸盐及活性磷酸盐等9个环境因子的影响进行分析,HSI模型结果显示,夏、秋、冬季硝酸盐、亚硝酸盐及盐度是影响弧边招潮蟹生境适宜性的关键环境因素。张超文等[29]在洪泽湖的研究中也发现,硝酸盐是影响大型底栖动物分布的最关键环境因子。王航俊等[30]在乐清湾大型底栖动物的研究中发现,含氮营养盐是影响大型底栖动物分布的主要驱动因素。虽然亚硝酸盐与硝酸盐对弧边招潮蟹生境选择的影响机理尚不明确,但生长速率假说[31-32]指出,环境中营养盐的浓度变化会影响生物的生长速率,进而影响物种分布及栖息地选择。本研究区域沿浦湾附近围塘养殖和滩涂养殖等产业发达,致使间隙水中营养盐含量丰富,而且一些滩涂养殖户在蛏苗培育前违法使用农药杀灭所谓的“天敌”生物的事件常有发生。因此,水体中营养盐和有机污染物等不断变化影响了水中有机质含量及藻类的生长状态,进而通过食物链的传递效应对弧边招潮蟹的生境选择产生了影响。

本研究表明,盐度是影响弧边招潮蟹生存的关键因素,随着盐度的升高,弧边招潮蟹的分布概率呈现逐渐降低的态势。李庭古[33]研究表明,盐度主要通过影响生物的渗透压进而影响其活动能力,同时生物也可通过改变自身的代谢状态以适应不断变化的盐度环境。彭超[34]研究发现,弧边招潮蟹更偏爱低盐度,对高盐度的耐受性较低,盐度5~20为蟹类适宜生存范围,当盐度超过30时,弧边招潮蟹幼蟹便无法生存,本研究结果与之相吻合。

4 结论

1)弧边招潮蟹因其活动能力较强、活动范围较大和生境的选择能力较强,故采用HSI模型预测其生境适应性比MaxEnt模型更合适。

2)沿浦湾秋茄林区弧边招潮蟹的栖息地适宜性总体呈现近岸段>中段>离岸段的空间分布特征,以及夏季>秋季>冬季的季节分布特征。这与弧边招潮蟹对水温变化的敏感性及其营穴居生活特性有关。表明弧边招潮蟹因温度、底质、地质和潮汐等环境因子变化而驱动其选择适宜的生境。

3)硝酸盐、亚硝酸盐和盐度等是影响弧边招潮蟹物种分布概率的关键因素,而营养盐浓度过高会对弧边招潮蟹产生胁迫作用,影响其生长速率,进而影响其分布及栖息地选择。表明理化因子及其变化对弧边招潮蟹分布格局具有较强的驱动能力,影响其对适宜生境的选择。

参考文献:

[1] 徐兆东,石洪华,李乃成,等.庙岛群岛南部海域大型底栖动物群落结构及其与环境因子的关系[J].环境科学研究,2015,28(5):704-712.

XU Z D,SHI H H,LI N C,et al.Macrobenthic community structure and its relationship with environmental factors in the southern waters of the Miaodao archipelago[J].Research of Environmental Sciences,2015,28(5):704-712.(in Chinese)

[2] GORE J,LAYZER J,MEAD J.Macroinvertebrate in stream flow studies after 20 years:a role in stream management and restoration[J].Regulated Rivers Research &Management,2001,17(4/5):527-542.

[3] RÜGER N,SCHLÜTER M,MATTHIES M.A fuzzy habitat suitability index for Populus euphratica in the northern Amudarya delta(Uzbekistan)[J].Ecological Modelling,2005,184(2/3/4):313-328.

[4] SPRUNG M.Macrobenthic secondary production in the intertidal zone of the Ria Formosa—a lagoon in southern Portugal[J].Estuarine,Coastal and Shelf Science,1994,38(6):539-558.

[5] 郭彦龙,赵泽芳,乔慧捷,等.物种分布模型面临的挑战与发展趋势[J].地球科学进展,2020,35(12):1292-1305.

GUO Y L,ZHAO Z F,QIAO H J,et al.Challenges and development trend of species distribution model[J].Advances in Earth Science,2020,35(12):1292-1305.(in Chinese)

[6] ROBERTO M,ZAMORA R,RAMN J,et al.Predictive modeling of microhabitats for endemic birds in South Chilean temperate forests using Maximum Entropy (MaxEnt)[J].Ecological Informatics,2011,6(6):364-370.

[7] 孙元超,刘凯,姚晓燕,等.基于MaxEnt模型对我国纹沼螺潜在分布区的预估研究[J].中国媒介生物学及控制杂志,2022,33(4):555-561.

SUN Y C,LIU K,YAO X Y,et al.Maximum entropymodel-based projection of potential distribution areas of Parafossarulus striatulus in China[J].Chinese Journal of Vector Biology and Control,2022,33(4):555-561.(in Chinese)

[8] 龚彩霞,陈新军,高峰,等.栖息地适宜性指数在渔业科学中的应用进展[J].上海海洋大学学报,2011,20(2):260-269.

GONG C X,CHEN X J,GAO F,et al.Review on habitat suitability index in fishery science[J].Journal of Shanghai Ocean University,2011,20(2):260-269.(in Chinese)

[9] BRAMBILLA M,CASALE F,BERGERO V,et al.GIS-models work well,but are not enough:habitat preferences of Lanius collurio at multiple levels and conservation implications[J].Biological Conservation,2009,142(10):2033-2042.

[10] VINCENZI S,CARAMORI G,ROSSI R,et al.A GIS-based habitat suitability model for commercial yield estimation of Tapes philippinarum in a Mediterranean coastal lagoon (Sacca di Goro,Italy)[J].Ecological Modelling,2006,193(1/2):90-104.

[11] 郑文浩,渠晓东,张远,等.太子河流域大型底栖动物栖境适宜性[J].环境科学研究,2011,24(12):1355-1363.

ZHENG W H,QU X D,ZHANG Y,et al.Habitat suitability of macroinvertebrates in the Taizi River basin,northeast China[J].Research of Environmental Sciences,2011,24(12):1355-1363.(in Chinese)

[12] 粟一帆,李卫明,李金京,等.桥边河大型底栖动物生境适宜性[J].生态学报,2020,40(16):5844-5854.

SU Y F,LI W M,LI J J,et al.Habitat suitability of macroinvertebrates:a case study in Qiaobian River,a tributary of Yangtze River,China[J].Acta Ecologica Sinica,2020,40(16):5844-5854.(in Chinese)

[13] 中华人民共和国国家质量监督检验检疫总局,中国国家标准化管理委员会.《海洋调查规范》 (GB/T 12763—2007)[S].北京:中国标准出版社,2008.

General Administration of Quality Supervision,Inspection and Quarantine of the People’s Republic of China,Standardization Administration of the People’s Republic of China. Specifications for oceanographic survey (GB/T 12763—2007)[S]. Beijing:Standards Press of China,2008.(in Chinese)

[14] 中华人民共和国国家质量监督检验检疫总局,中国国家标准化管理委员会.《海洋监测规范》(GB 17378—2007) [S].北京:中国标准出版社,2008.

General Administration of Quality Supervision,Inspection and Quarantine of the People’s Republic of China,Standardization Administration of the People’s Republic of China. Marine monitoring specifications (GB 17378—2007)[S]. Beijing:Standards Press of China,2008.(in Chinese)

[15] 李国梁,汪文俊,李宝贤,等.基于MaxEnt模型和ArcGIS预测多肋藻在中国海域的适生分布特征[J].中国水产科学,2021,28(12):1588-1601.

LI G L,WANG W J,LI B X,et al.Potential geographic distribution of Costaria costata in China based on the MaxEnt model and ArcGIS[J].Journal of Fishery Sciences of China,2021,28(12):1588-1601.(in Chinese)

[16] 陈艺璇,张云雷,黄锘妍,等.海州湾鹰爪虾栖息地适宜性研究[J].海洋学报,2021,43(4):84-95.

CHEN Y X,ZHANG Y L,HUANG N Y,et al.Study on the habitat suitability of Trachypenaeus curvirostris in the Haizhou Bay[J].Haiyang Xuebao,2021,43(4):84-95.(in Chinese)

[17] ELITH J,LEATHWICK J.Species distribution models:ecological explanation and prediction across space and time[J].Annual Review of Ecology and Systematics,2009,40(1):677-697.

[18] ANDERSON R P.A framework for using niche models to estimate impacts of climate change on species distributions[J].Annals of the New York Academy of Sciences,2013,1297:8-28.

[19] CHANG Y J,SUN C L,CHEN Y,et al.Habitat suitability analysis and identification of potential fishing grounds for swordfish,Xiphias gladius,in the South Atlantic Ocean[J].International Journal of Remote Sensing,2012,33(23):7523-7541.

[20] 王昱,左一锋,冯起,等.黑河中上游大型底栖动物栖息地适宜度评估[J].生态学杂志,2021,40(4):1116-1127.

WANG Y,ZUO Y F,FENG Q,et al.Assessment of macrozoobenthos habitat suitability in the upper and middle reaches of the Heihe River[J].Chinese Journal of Ecology,2021,40(4):1116-1127.(in Chinese)

[21] TANAKA K R,CHEN Y.Spatiotemporal variability of suitable habitat for American lobster (Homarus americanus) in Long Island Sound[J].Journal of Shellfish Research,2015,34(2):531-543.

[22] 吴庆明,王磊,朱瑞萍,等.基于MAXENT模型的丹顶鹤营巢生境适宜性分析:以扎龙保护区为例[J].生态学报,2016,36(12):3758-3764.

WU Q M,WANG L,ZHU R P,et al.Nesting habitat suitability analysis of red-crowned crane in Zhalong Nature Reserve based on MAXENT modeling[J].Acta Ecologica Sinica,2016,36(12):3758-3764.(in Chinese)

[23] PHILLIPS S J,ANDERSON R P,SCHAPIRE R E.Maximum entropy modeling of species geographic distributions[J].Ecological Modelling,2006,190(3/4):231-259.

[24] 王庆莉,王茹琳,张利平,等.基于MaxEn模型的川西高原松茸气候生态适宜性与潜在分布[J]. 应用生态学报,2021,32(7):2525-2533.

WANG Q L,WANG R L,ZHANG L P,et al.Climatic ecological suitability and potential distribution of Tricholoma matsutake in wester Sichuan Plateau,China based on MaxEntmodel[J].Chinese Journal of Applied Ecology,2021,32(7):2525-2533.(in Chinese)

[25] 杨刚,张涛,庄平,等.长江口棘头梅童鱼幼鱼栖息地的初步评估[J].应用生态学报,2014,25(8):2418-2424.

YANG G,ZHANG T,ZHUANG P,et al.Preliminary assessment of habitat of juvenile Collichthys lucidus in the Yangtze Estuary[J].Chinese Journal of Applied Ecology,2014,25(8):2418-2424.(in Chinese)

[26] 刘振生,高惠,滕丽微,等.基于MAXENT模型的贺兰山岩羊生境适宜性评价[J].生态学报,2013,33(22):7243-7249.

LIU Z S,GAO H,TENG L W,et al.Habitat suitability assessment of blue sheep in Helan Mountain based on MAXENT modeling[J].Acta Ecologica Sinica,2013,33(22):7243-7249.(in Chinese)

[27] KLAASSEN M,ENS B J.Habitat selection and energetics of the fiddler crab (Uca tangeri)[J].Netherlands Journal of Sea Research,1993,31(4):495-502.

[28] 李明云.普陀及宁波滩涂弧边招潮蟹生物学特性的研究[J].动物学杂志,1996,31(6):32-35.

LI M Y.Study on biological characteristics of fiddler crabs along the arc in Putuo and Ningbo beaches[J].Chinese Journal of Zoology,1996,31(6):32-35.(in Chinese)

[29] 张超文,张堂林,朱挺兵,等.洪泽湖大型底栖动物群落结构及其与环境因子的关系[J].水生态学杂志,2012,33(3):27-33.

ZHANG C W,ZHANG T L,ZHU T B,et al.Community structure of macrozoobenthos and its relationship with environmental factors in Lake Hongze[J].Journal of Hydroecology,2012,33(3):27-33.(in Chinese)

[30] 王航俊,姚炜民,林义,等.乐清湾大型底栖动物群落及其与环境因子之间的关系[J].海洋学报,2020,42(2):75-86.

WANG H J,YAO W M,LIN Y,et al.Macrobenthic community and its relationship with environmental factors in the Yueqing Bay,Zhejiang Province,China[J].Haiyang Xuebao,2020,42(2):75-86.(in Chinese)

[31] ELSER J J,O’BRIEN W J,DOBBERFUHL D R,et al.The evolution of ecosystem processes:growth rate and elemental stoichiometry of a key herbivore in temperate and Arctic habitats[J].Journal of Evolutionary Biology,2000,13(5):845-853.

[32] 苏强.生长速率假说及其在浮游动物营养动力学中的研究进展[J].地球科学进展,2012,27(11):1204-1210.

SU Q.Growth rate hypothesis research progresses:implications for zooplankton[J].Advances in Earth Science,2012,27(11):1204-1210.(in Chinese)

[33] 李庭古.盐度对克氏原螯虾幼虾生长的影响[J].水产科学,2009,28(8):465-467.

LI T G.The effects of different salinity on survivals and growth of juvenile Procambarus clarkii[J].Fisheries Science,2009,28(8):465-467.(in Chinese)

[34] 彭超.氨氮对弧边招潮蟹的毒性效应及盐度对其存活、摄饵的影响[D].海口:海南师范大学,2013.

PENG C.Toxic effects of ammonia nitrogen on fiddler crabs along the arc and effects of salinity on their survival and feeding[D].Haikou:Hainan Normal University,2013.(in Chinese)

Comparison of habitat suitability of fiddler crab (Uca arcuata) based on habitat suitability index model and maximum entropy model

ZHOU Xin,ZHANG Liyuan,HE Chenxiang,HU Chengye,WANG Jing,LI Yi,LI Xueli,SHUI Bonian*

(College of Fisheries,Zhejiang Ocean University,Zhoushan 316022,China)

AbstractFor understanding of the habitat suitability and optimal analysis model of fiddler crab (Uca arcuata), key environment factors affecting changes in habitat suitability, and the predictive performance, scope and limitations of application of the two habitat suitability models including habitat suitability index (HSI) model and the maximum entropy (MaxEnt) model were compared with a view to species classification of mangrove ecosystems based on the survey on mangrove macrobenthos and several environmental factors in Yanpu Bay (27.20°N, 120.45°E) Wenzhou City, Zhejiang Province during July (summer), October (autumn) and December (winter) of 2020 to provide theoretical support for applied research on biodiversity conservation and scientific basis for ecological environment governance and decision-making. HSI model analysis showed that the descending order of the habitat suitability of the sampling stations was expressed as nearshore>middle>offshore section, with the southeast>the northwest in the spatial distribution, and significantly different habitat suitability in the three seasons (P<0.05), generally as summer>autumn>winter. The analysis of the maximum entropy (MaxEnt) model revealed that the descending order of habitat suitability of the fiddler crab in each sampling site as nearshore>middle section>offshore section, the southeast>the northwest, and autumn>summer>winter in the seasonal distribution. The HSI model was shown to be more suitable for the investigation of the suitable area of the fiddler crab than the MaxEnt model being; temperature, bottom, geology and tide were the environmental factors affecting the distribution; nitrate, nitrite and salinity were the key physical and chemical factors affecting the distribution of the fiddler crab, which have a strong driving force for the distribution pattern and influence the choice of suitable habitat.

Key wordsUca arcuata; suitability index model; maximum entropy model; habitat adaptability

中图分类号S 931

文献标志码:A

DOI10.16535/j.cnki.dlhyxb.2023-010

文章编号:2095-1388(2023)05-0819-09

收稿日期2023-01-17

基金项目浙江省自然科学基金(LQ21C030006);浙江省省属高校基本科研业务费项目(2021J003);国家重点研发计划重点专项(2017YFA0604902);“水产学”浙江省一流学科建设开放课题;浙江省软科学研究计划项目(2015C25018)

作者简介周鑫(1998—),男,硕士研究生。E-mail:1359489764@qq.com

通信作者水柏年(1965—),男,博士,教授。E-mail:shuibonian@163.com