海洋大型藻类是海洋中重要的初级生产者，同时也是多种鱼类、贝类、底栖生物的重要栖息地，为其提供必不可少的遮蔽环境、食物来源、栖息场所、繁殖区域等环境条件。近年来，围绕山东半岛及辽东半岛藻场修复的研究有了很大进展，在受损海区生态环境功能修复上取得一定成效(涂忠等, 2007、2009)，并发挥出相应生态功能(Largo et al, 1994; Seaman, 2007; Jiao et al, 2014)。然而，海藻场的建设更多看重其生态功能，以吸引鱼虾生物及浮游生物作为主要目的(张磊等, 2012; 袁小楠等, 2017)，若能将富有工业价值的海洋大型藻类，如海带属(Laminaria)，作为主导藻种构建海藻床，在起到生态恢复功能的同时，也可以综合提高藻场的经济产出率。极北海带(Laminaria hyperborea)系多年生大型海藻，源于大西洋东北海域，是当地“海藻森林”的优势藻种(Steneck et al, 2002)。极北海带具有较高经济价值，其柄和叶片的褐藻胶含量分别为25%~30%和17%~33%(干重比)，属于褐藻胶含量较高的褐藻之一，且其褐藻胶凝胶强度大，因古罗糖醛酸(Guluronic acid, G)含量在其柄和叶片的褐藻胶中分别占比高达68%~70%和44%~55%，而在其他商用褐藻胶中的含量一般为10%~45% (Peter et al, 2015)，工业和商业价值显著。

在人工藻场的构建中，由于需要考虑生态恢复效果，提高底栖生物丰富度，海藻附着基通常设置在水深为15~25 m的硬质海底(许强等, 2011)。然而，我国常见的海带(Saccharina japonica)通常的垂直适生深度小于15 m (曾呈奎等, 1981)，在人工鱼礁深度的海水层中，光照强度往往无法达到海带正常生长所需的光强，使海带生长较慢，甚至死亡，导致生态功能恢复效果不理想。海带的适生深度主要受光照强度的影响，自然生长的极北海带对光照的需求较低，因此，在更深的海域也有分布。Joanna(1997)研究发现，在25 m以下深度的海底存在自然生长的极北海带。Luning(1979)对极北海带进行光照实验发现，极北海带生长所需光照显著低于其他海带属种类(L. saccharina, L. digitata)。因此，针对藻礁建设区域深度较深、光照不足的问题，极北海带可能是一个理想的适应品种。

Phillips等(2008)基于生态位理论，用最大熵原理作为统计推断工具，以物种在已知分布区的信息以及目标区的环境变量为基础，通过比较该物种在已知分布区的生态环境变量来预测其占有的生态位，并编写程序MaxEnt(Maxent entropy)。MaxEnt仅考虑目标物种在所在区域的“出现”情况(Presence-only)，综合考虑各种给定环境因子下物种出现的最大熵，从而对物种在预测区域的适生情况进行分析。该模型在预测物种适生性上有广泛应用，许多学者使用该模型对珍稀物种、入侵物种以及一些植物病虫害的潜在分布区域进行预测评价(曾辉等, 2008; 张飞等, 2017; 王运生等, 2007)。此外，该模型也被证实适合预测海洋藻类的适生范围(Martín-García et al, 2017)。

尽管是构建海藻场的优良藻种，极具开发潜力，但作为外来物种，极北海带在我国黄渤海海区的适应性有待研究，是否会成为入侵物种的风险也尚不明确。因此，本研究采用MaxEnt软件对极北海带在我国黄渤海海区进行适生范围预测，评价其生态入侵风险，为构建极北海带藻场的可行性提供参考依据。

本研究选用MaxEnt v3.3.3对我国近海海域模拟极北海带的适生性进行评估，选取logistic形式表达，其他算法选择模型自身默认值，并将所选用的环境因子做刀切法(JackKnife)检验，判断每个环境因子对预测结果的重要性。

本研究使用的极北海带空间数据从全球生物多样性信息网络GBIF (http://www.gbif.org/)中获取。收集1900~2018年极北海带出现位点共11894处。除去实验、标本存储及其他人为培养位点822处，得到极北海带自然扩散位点11072处的经纬度数据。考虑到模型运算量以及许多点位较为集中，使用Excel在11072个点位中随机选择300个位点作为极北海带的空间位点数据，将300个空间位点数据转换成*.csv格式导入MaxEnt模型。

于全球海洋生物扩散模型环境数据库Bio-ORICLE(http://bio-oracle.org/)中筛选出影响极北海带生长的21项关键环境因子，用于适生模型构建。其中包含12项水质参数，即：最高月均温度、溶解氧、氮浓度、磷浓度；最低月均温度、溶解氧、氮浓度、磷浓度，年平均温度、溶解氧、氮浓度、磷浓度，以及9项水文参数，即：最大海流流速、海区深度、光照强度；最小海流流速、海区深度、有效太阳辐射，年平均海流流速、海区深度、光照强度。

本研究采用MaxEnt模型中自带的刀切法，将上述21项因子进行筛选。将全球生物多样性信息网络中所收集到的极北海带空间位点经纬度数据转化成数据文本格式(.csv)，并于Bio-ORICLE中所收集的环境数据文件(.asc)同时导入MaxEnt模型，以logistic格式输出最终结果。模型共计算2次，第1次计算用于筛选对适生模型贡献度为0的环境因子，舍去这些因子，将剩余环境因子再次导入MaxEnt，进行第2次预测模型构建。

模型预测结果以*.asc格式输出。将模型输出结果导入ArcGIS 10.3，利用ArcGIS工具箱系统中的格式转换工具(conversion tools)，将预测结果数据转化为栅格数据(*.flt)格式。选定合适的阈值，将预测结果按适生概率等级划分为高适生区(> 0.5)、中适生区(0.3~0.5)、低适生区(0.1~0.3)、边缘适生区(0.05~0.1)、非适生区(< 0.05)(Martín-García et al, 2014; 齐国军等, 2012; 柯明剑等, 2017)。

本研究通过计算指受试者工作特征曲线(ROC)的曲线下面积(AUC)来评估模型准确性(Canran et al, 2005)。在进行模拟时，使用模型自身功能将所选取的300个位点分为2部分，随机选取75%的物种分布点用于建模，剩余25%用于模型的检验和校正(Phillips et al, 2006、2008)。在“presence-only”模型中，ROC曲线表示正确预测的出现位点概率与假信号出现概率的比值，ROC曲线下的面积(AUC)即表示预测结果的准确性。AUC数值范围为0~1，AUC数值越大，表示预测结果越准确，AUC值大于0.9时，表示模型预测结果有很高的可信度。

黄渤海海区统计共5213个栅格数据，最大值为0.472，最小值为0.000(精确到小数点后3位)。将适生概率按适生区域阈值划分，频数与频率见表 1。

表 1(Tab.1)

表 1 极北海带在我国黄渤海海域的适生区分布情况 Tab.1 The occurrence zones of L. hyperborean in Bohai and Huanghai Sea 适生区划分Occurrence zone 适生率范围Probability range 频数Frequency 频率Percentage (%) 非适生区Unbefitting zone 0~0.05 4614 88.51 边缘适生区Marginal zone 0.05~0.1 245 4.70 低适生区Low OPZ 0.1~0.3 267 5.12 中适生区Medium OPZ 0.3~0.5 87 1.67 高适生区High OPZ 0.5~1 0 0 合计Total 5213 注：适生频数及概率范围通过ArcGIS中的工具箱系统中的格式转换工具重新划分 Note: Occurrence probability zones (OPZs) and frequency are reclassified by the use of conversion tool from ArcGIS

表 1 极北海带在我国黄渤海海域的适生区分布情况 Tab.1 The occurrence zones of L. hyperborean in Bohai and Huanghai Sea

极北海带在中国黄渤海海域的适生性预测结果如图 1所示，适生概率从蓝色到红色依次上升。极北海带在山东半岛沿岸、渤海湾近岸海域以及辽东半岛沿岸均有一定的适生概率，低适生区多分布在渤海西北部及辽东半岛沿岸。山东半岛南岸有少量边缘适生区，再往南均无适生区。中、高适生区域主要分布在朝鲜及日本沿海地区，我国海域无分布。

图 1(Fig. 1)

图 1 MaxEnt模型预测极北海带在中国黄渤海海域的适生性概率范围 Fig.1 Occurrence probability zones of L. hyperborean in Bohai and Huanghai Sea

评估预测模型准确性的ROC曲线如图 2所示，其中，红色为用于模型检验校正，蓝色为模型建立。其AUC值分别为0.966与0.965。

图 1(Fig. 1)

图 2 MaxEnt预测模型的ROC曲线 Fig.2 ROC curve of MaxEnt

采用MaxEnt模型中自带的刀切法功能进行环境变量贡献度分析。第1次分析后，从所选21项环境变量中，将贡献度为0的12项变量删除，保留剩余9项环境变量。SPSS检测9项变量两两之间无显著相关性。图 3为9项有贡献度的环境变量刀切法分析结果。纳入考虑的9个环境变量分别为最大月光照强度Light(max)、年均光照强度Light(mean)、平均海流流速Current(mean)、年均溶解氧浓度DO(mean)、最低月均氮浓度C_N(min)、最低月均磷浓度C_P(min)、最大月均温度T (max)、年均温度T(mean)以及最小月均温度T(min)。光照强度和温度对极北海带适生概率值的贡献度最为显著。

图 1(Fig. 1)

图 3 环境因子对模型的贡献度 Fig.3 The contribution of each environmental variable

选取MaxEnt模型筛选出的对极北海带适生性影响较大的环境因子逐个分析(图 4)，仅考虑单一环境变量时的极北海带适生概率，用来反映模型在预测适生性时，对所选变量的依赖性以及所选变量和其他变量之间的相关性引起的依赖性。由图 4可见，极北海带适生概率随光照强度的增加而增加，在最大光照强度52 E/(m²∙d) (Photosynt. Avail. Radiation)、平均强度48 E/(m²∙d)时达到顶点。极北海带适生温度范围为7℃~18℃、在年均温度为13℃时达到最大适生概率。同时，极北海带的适生温度范围最小月均温度不低于–2℃，最大月均温度不超过25℃。氮、磷浓度的单因子适生概率分布图显示，极北海带在0~10 mg/L氮浓度范围、0~0.22 mg/L磷浓度范围内存在适生性。

图 1(Fig. 1)

图 4 单因子环境变量分析 Fig.4 Analysis of single environmental variable

本研究通过MaxEnt建立的极北海带适生性预测模型，训练集和测试集的AUC值均非常接近于1，说明结果预测准确(Phillips et al, 2006)。本研究通过刀切法筛选出对极北海带适生性贡献最大的9项环境因子也均符合极北海带的生长条件(Susan et al, 2011)。MaxEnt模型预测结果显示，除朝鲜西岸与日本主岛东部沿岸以外，在我国范围内仅辽东半岛东部大连长海县附近海域存在中适生区，其余地区均为低适生区、边缘适生区和非适生区。

温度是影响极北海带适生性的重要因子。本研究结果显示，极北海带适生温度年均值在6~17℃范围内，月均温度不超过25℃，最低不低于–2℃。大连长海县附近海域年均水温通常为–2~24℃，平均温度和最低温度在极北海带适生温度范围之内，但夏季月均最高温度接近于极北海带的最高耐受温度范围。因此，在夏季，温度有可能是在该海区内影响极北海带最重要的环境因子，在构建海带场时应重点关注。

光照强度是另一项影响极北海带分部主要的环境因子之一。模型自身的刀切法检测将最小光照强度对极北海带适生性贡献排除。该结果与预测结果相符，说明极北海带对低光照环境的适应性较好，相比与其他环境因子，自然的低光照环境对极北海带可能不构成影响。

氮、磷浓度对极北海带的适生范围也有一定的影响，单因子适生性分析结果显示，极北海带在氮、磷浓度极低的水平下均能达到较高的适生性，表明极北海带适生范围受氮、磷浓度影响相对较小。此外，最大及平均氮、磷浓度对该预测模型没有贡献值，也说明自然水体能达到的最高氮、磷浓度不足以成为极北海带生长的限制因子。

海洋底质对极北海带分布的影响没有纳入考虑范围。尽管海洋底质对大型海藻的自然分布影响很大，海带生长通常需要礁岩作为附着基。但藻场的建设通常包含人工投放礁石和海藻附着基，人为为海带生长提供适宜的海洋底质环境，或直接采用苗绳方法进行培养。因此，将底质类型排除在模型所考虑的环境因子之外。

引入外来物种需要严格把控其适生范围，防止其成为入侵物种。辽东半岛东部海域底质以沙质为主，不利于极北海带的自然生长。因此，在沙质底质建立引入的极北海带品种还可降低其成为入侵物种的风险，通过控制海洋牧场附着基质的投放区域，起到控制极北海带生长范围的作用。另一方面，适生概率超过0.5的高适生区有较高的入侵风险(Phillips et al, 2006)，而在我国近海范围内均无高适生区。因此，极北海带在我国海域内成为入侵物种的风险较低。需要注意的是，辽东湾西部和山东半岛北部海域底质以岩石为主，极北海带在此有一定的适生概率。然而，MaxEnt预测此海域属于低适生区，成为入侵物种的概率较低。其余低适生区例如渤海湾西岸，其底质以淤泥为主，不适合极北海带生长。

综上所述，引进极北海带作为海洋牧场的主导藻种在我国辽东半岛东岸的长海县附近海域有较高的可行性。同时，受到环境因子的限制，极北海带成为入侵物种的风险较低，但在辽东湾西部和山东半岛北部海域仍需要严格审查。