日本东北部潮下带海带Eisenia bicyclis的补充：多种环境因素的影响

《Oceans》：Recruitment of the Subtidal Kelp Eisenia bicyclis in Northeastern Japan: Effects of Multiple Environmental Factors

　　

为了确定影响海带Eisenia bicyclis幼体补充（juvenile recruitment）的因素，研究人员在受2011年地震引起的地震沉降（seismic subsidence）及随后防波堤修复影响的区域开展了一项为期七年的监测调查。幼体补充在2011年9月以及2013–2015年的5月至6月较高，但在2012年、2016年和2017年较低。对环境因子与幼体补充之间关系的分析表明，补充与光照强度（light intensity）、两个月前的较低水温（water temperature）以及四个月前营养盐（nutrients）的增加相关。冬季营养盐浓度的季节性增加可能受到季节性西北向沿岸流的影响。相比之下，尽管营养盐浓度升高，2016–2017年的补充相对较差。这可能归因于与防波堤建设相关的海床环境不稳定。研究人员对同一地点E. bicyclis幼体补充数量和环境因子的月度监测表明，营养盐浓度的时滞性增加和水温的降低与Eisenia bicyclis幼体补充相关。本研究提供了海带补充的基础信息，这将有助于预测干扰后海带群落的恢复以及环境变化下的补充动态。

### 论文解读：日本东北部潮下带海带Eisenia bicyclis补充机制的多环境因子研究

#### 1. 研究背景与问题
海带林（kelp forest）在浅海岩礁区提供高光合生产力和生物多样性，但近年来因全球变暖导致的海水温度升高、沿海城市化的栖息地丧失和水污染等多重压力而衰退。幼体补充（juvenile recruitment）是海带种群动态和干扰后恢复的关键过程，以往研究认为低温、高营养盐条件促进补充，高温和低营养盐抑制补充。然而，现有研究多基于实验室或短期观测，缺乏长期、原位同步监测环境因子与补充的野外证据。此外，2011年日本东北太平洋地震引发的地震沉降和后续防波堤修复工程改变了海底环境，为研究多因子交互作用对海带补充的影响提供了独特契机。本研究旨在同步测量与Eisenia bicyclis补充相关的环境因子（水温、营养盐、光照强度等），并明确与地震及修复工程影响下海带床恢复最密切的因素。

#### 2. 研究内容与结论
研究人员在宫城县牡鹿半岛西侧的狐崎（Kitsunezaki）设置永久样地（36 m2），从2011年7月至2017年7月每月监测E. bicyclis幼体（单叶、总叶状体长度>15 mm）和成体密度，同时连续记录海底水温、光照强度，每月采集水样测定营养盐（NH₄-N、NO₃-N、NO₂-N、PO₄-P），并利用CTD测量盐度、叶绿素a浓度、浊度，使用电子流速计监测潮汐流。采用零膨胀负二项模型（Zero-Inflated Negative Binomial Model, ZINB）分析环境因子（当月累积光照、时滞1–4个月的水温偏差和总营养盐浓度）对幼体补充密度的影响，通过赤池信息准则（AIC）选择最优模型。

主要结论：光照强度是幼体补充的必要条件；补充受到时滞环境因子的影响——当月光照强度越高、两个月前水温越低、四个月前营养盐浓度越高，幼体补充数量越多。冬季西北向沿岸流带来营养盐升高，但2016–2017年尽管营养盐充足，由于防波堤修复导致海床不稳定（物理扰动增加），补充仍较低。本研究首次在同一地点长期同步监测证实了这些时滞效应，为预测干扰后海带群落恢复和环境变化下的补充动态提供了基础数据。论文发表在《Oceans》。

#### 3. 主要技术方法
研究人员在日本宫城县牡鹿半岛西侧狐崎（38.355751° N, 141.423188° E）设立永久样地。每月人工计数并标记所有E. bicyclis个体（幼体与成体）。环境监测包括：使用HOBO数据记录器（Onset）每15分钟记录海底（水深0.8 m）水温和光照强度，光照值经现场标定转换为光子通量密度（mol m^-2 day^-1），并计算25天滑动平均值；每月采集海底水样，经0.2 μm滤膜过滤后用Autoanalyzer（QuAAtro 2HR）测定NH₄-N、NO₃-N、NO₂-N、PO₄-P浓度；自2013年4月起使用RINKO-Profiler（JFE Advantec）测量盐度、叶绿素a浓度、浊度和水温，并对水团进行分类（表层水、亲潮水等）；自2014年5月起使用INFINITY-EM电子流速计在样地外3.1 m水深监测海流速度和方向（每15分钟采集10个样本，取25天滑动平均）。统计上采用ZINB模型，以AIC选择最佳时滞模型。

#### 4. 研究结果

##### 3.1 幼体和成体密度的时序变化
通过每月计数发现：E. bicyclis幼体补充在2011年9月（1.56 ind. m^-2）及2013–2015年5–8月（峰值密度2.83–3.19 ind. m^-2）较高，而2012年、2016年和2017年无显著高峰（年补充密度0.47–0.69 ind. m^-2）。成体最大密度出现在2012年2月（0.92 ind. m^-2），之后逐渐下降至2016年8月的0.56 ind. m^-2，2016年10月降至0.36 ind. m^-2。

##### 3.2 水温和光照强度时序变化
月均水温在5.5–26.2 °C之间波动；2012年1–5月、2013年5–7月、2014年9月–2015年4月、2015年9–12月出现低于多年均值>2.0 °C的低温期；2012年9–11月、2015年5–8月、2017年7月出现高于均值>2.0 °C的高温期。25天滑动平均的日光子通量密度夏季高（最高21.8 mol m^-2 day^-1）、冬季低（最低2.0 mol m^-2 day^-1），且始终高于E. bicyclis的日补偿光量（0.3–1.5 mol m^-2 day^-1），表明光照条件总体充足。

##### 3.3 营养盐时序变化
营养盐浓度冬季升高，夏季常低于检出限（<0.01 μmol L^-1）。NH₄-N浓度0.22–3.76 μmol L^-1，在2011年8月、2013年12月和2017年5月出现峰值；NO₃-N（0.22–1.64 μmol L^-1）、NO₂-N（0.12–0.41 μmol L^-1）和PO₄-P（0.11–0.28 μmol L^-1）在每年10月至次年2月最高。溶解无机氮（DIN）总体冬季高、夏季低。

##### 3.4 水温、光照和营养盐对幼体补充密度的影响
通过ZINB模型比较17个候选模型（含1–4个月时滞），最佳模型（AIC=269.94）包含当月光照强度、两个月前水温和四个月前营养盐浓度。计数模型系数表明：当月光照强度具有显著正效应（p<0.001），两个月前水温有显著负效应（p<0.01），四个月前营养盐有显著正效应（p<0.01）。零膨胀模型显示光照强度负相关于零补充概率（p<0.05），即光照越高，零补充概率越低。

##### 3.5 水质时序变化
盐度28.2–33.7，冬季高夏季低；浊度和叶绿素a浓度夏季升高、冬季降低。基于温盐图的水团分类显示，研究点水流主要为表层水（Surface-layer water），但2014年2月出现亲潮（Oyashio）水团，2017年2–3月出现津轻暖流（Tsugaru Warm Current）水团。

##### 3.6 潮汐流时序变化
东西向流中，2013年11月–2014年2月和2014年7–9月以西向流为主；南北向流中，夏季南向流占主导，冬季北向流占主导。整体上5–9月为西南向流，10–2月为西北向流。

#### 5. 讨论与结论
讨论部分指出：光照强度是满足补充的充分必要条件，成体密度（<1 ind. m^-2）未造成遮阴限制，但可能影响孢子供应。水温降低单独不能完全解释补充高峰，因为2011年低温但补充不高，提示营养盐的时滞效应更为关键。营养盐高峰后四个月出现补充高峰，与实验室和野外研究中营养盐吸收和叶状体生长的时间滞后一致。冬季营养盐升高与西北向沿岸流相关，该流动由冬季西北风驱动的离岸流和湾内补偿流引起，将富含营养盐的沿岸水带入研究区。尽管营养盐在2016–2017年升高，但补充仍低，推测系因防波堤修复导致海床不稳定（底质扰动增加）所致，物理扰动是不稳定底质上群落结构的重要决定因素。

研究结论（翻译自摘要和讨论）：本研究通过在同一地点同步监测E. bicyclis幼体补充数量和环境因子，揭示了营养盐浓度的时滞性增加（四个月前）和水温的降低（两个月前）与幼体补充相关，光照强度是必要条件。这些发现为预测干扰后海带群落的恢复以及环境变化下的补充动态提供了基础信息。