《Canadian Journal of Fisheries and Aquatic Sciences》：The power of baseline data from fixed sampling sites: a simulation-based power analysis of marine harvest refugia
编辑推荐：
海洋捕捞庇护区（marine harvest refugia）有助于保护生物多样性和生态系统；然而，以足够的统计功效评估其有效性具有挑战性。量化具有足够功效的采样需求可改进数据收集方法、研究设计、精度并降低成本。传统功效方程并非为复杂生态系统模型而设计，而且与
海洋捕捞庇护区（marine harvest refugia）有助于保护生物多样性和生态系统；然而，以足够的统计功效评估其有效性具有挑战性。量化具有足够功效的采样需求可改进数据收集方法、研究设计、精度并降低成本。传统功效方程并非为复杂生态系统模型而设计，而且与稀疏的干预前数据相结合，导致庇护区有效性的证据有限。研究人员将空间显式高斯马尔可夫随机场（Gaussian Markov random field, GMRF）模型拟合至来自遥控潜水器（remotely operated vehicles, ROV）的鱼类密度数据，并利用基于模拟的功效分析比较两种采样方案。研究人员通过比较不同采样策略和努力水平下的早期-后期-对照-影响（Early-Later-Control-Impact, ELCI）设计与对照-影响（Control-Impact, CI）设计，量化检测丰度恢复的功效。结果表明：（1）早期样本具有很高价值，尤其是在庇护区内 vs. 外初始鱼类丰度与预期存在显著差异时；（2）小样本量会膨胀效应估计，提示应考虑类型M（量级）错误（type M (magnitude) errors）；（3）重复采样固定地点——而非在扩展采样域内部或跨域选择新的随机站点——可使功效提高至35%，并降低所需样本量。这些发现为设计能够可靠检测海洋捕捞庇护区恢复的成本有效监测计划提供了指导。
研究背景与问题：海洋捕捞庇护区（marine harvest refugia；限制人类利用以保护生物多样性和维持生态系统的海域）是保护生物多样性和维持海洋生态系统的重要工具，但其有效性评估受制于广阔、复杂或偏远海洋环境监测困难、边界执法困难、数据收集成本高和干预前基线数据稀疏等问题。传统统计功效方程基于简单模型，难以适配用于刻画鱼类分布的空间相关广义线性混合模型（generalized linear mixed model, GLMM；含随机效应的广义线性模型）等复杂生态模型；样本量不足、生态噪声和缓慢响应又会增加类型S错误（type S error；显著估计方向错误）和类型M错误（type M error；显著效应量被夸大/膨胀比错误）。因此，该研究以2018年Strait of Georgia的遥控潜水器（remotely operated vehicles, ROV）鱼类密度数据条件化空间高斯马尔可夫随机场（Gaussian Markov random field, GMRF；通过稀疏精度矩阵近似空间相关性的随机场）模型，采用基于模拟的功效分析，比较对照-影响（Control-Impact, CI）与早期-后期-对照-影响（Early-Later-Control-Impact, ELCI）设计及不同采样策略和努力水平，为平鲉保护区域（Rockfish Conservation Areas, RCAs；一类海洋捕捞庇护区）监测计划确定具有足够正确符号功效的采样需求，并评估固定站点重复采样的价值。
主要关键技术方法：基于2018年Strait of Georgia ROV调查数据（春季WDFW 49条、秋季DFO/WDFW/NOAA Fisheries 82条；共131条样带，54条在RCAs内、77条在外，4408个数据点）拟合空间GLMM-GMRF条件模型（负二项似然、log链接；纳入海底位置指数（BPI；描述局部地形相对位置的地形变量）、水深二次项、岩石、温度、RCA状态；随机偏微分方程（SPDE）/Matérn空间协方差）。再用该模型模拟鱼类密度，在CI与ELCI设计下变化样本量100–500和效应量1.1–1.4，重复400次，以Wald 95%置信区间判定显著性，计算正确符号功效、S错误和M错误，并比较固定站点、约束随机和未约束随机后期采样设计。
研究结果：3.1 条件模型表现（Conditioning model performance）：通过将空间GMRF-GLMM拟合至ROV数据，BPI与密度呈显著正相关（β=0.16，95% CI=0.06, 0.27），夏季平均温度呈显著负相关（β=-0.69，95% CI=-1.07, -0.32），岩石区域密度更高（β=1.10，95% CI=0.93, 1.27），RCAs内密度低于外（β=-0.50，95% CI=-0.74, -0.26）；深度关系在约50 m达峰但线性与二次项均不显著。空间随机场使预测值更极端，细尺度变异主要由固定效应决定。3.2 外推（Extrapolation）：模拟框架包含超出条件模型拟合范围的值；在BPI>466时预测鱼数增加，而在BPI<-376、夏季平均温度>14 °C、深度>166 m时预测值趋近零。该部分通过比较模拟协变量范围与拟合数据范围，量化了外推区预测行为。3.3 基于模拟的功效分析（Simulation-based power analysis）：3.3.1 对照-影响设计（Control-Impact design）：因条件模型中RCAs内密度低于外，CI设计在不同效应量和样带数下均产生100%符号错误，正确符号功效为0%；传统功效虽随样本量增加，但随效应量增大而下降，原因是估计越来越精确地指向错误效应。3.3.2 早期-后期-对照-影响设计（Early-Later-Control-Impact design）：ELCI设计的传统功效与正确符号功效相同（符号错误0%），并随样本量和效应量增加而增加；当RCAs内密度增加40%（效应量1.4）时，200–500条样带可达≥80%功效；增加30%（效应量1.3）时，400–500条样带可达；增加10%和20%时，即使500条样带也未达80%。M错误通常随样本量和功效增加而下降。3.3.3 采样设计分析（Sampling design analysis）：在后期样本量均为131条时，固定站点设计的功效高于约束随机设计，约束随机设计又略高于未约束随机设计；固定设计在密度增加30%（效应量≥1.3）时达到80%功效，而两种随机设计仅在增加40%时达到；固定设计的M错误率通常更低。
讨论部分指出，早期基线数据对于检测RCAs效应至关重要；如果没有基线数据，由于庇护区内初始密度较低，即使40%的增加也可能无法检出。固定地点重复采样可将功效提高至35%，并使M错误率降低至多8.2%，在早期数据质量高时可仅返回早期地点，以不足一半的采样努力达到80%功效。空间约束还可减少向未观测区域外推的不确定性。小但现实的效应量即使在大规模野外计划和固定站点下仍可能难以检测。M错误和S错误提醒，小样本和噪声系统中的显著结果可能被夸大或方向相反；CI设计在基线耗竭时可能失败。讨论也列出模拟假设与限制：仅一年数据，无法参数化潜在时空效应；未纳入溢出效应；恢复被假设为log空间线性效应，未包含种群动态。未来可检验不同监测方案，并在ROV样本不可行时考虑环境DNA（eDNA）、潜水员观察、声学监测或提取式采样，并通过方法特异检测概率整合多源数据。研究结论翻译：综合而言，这些结果突显了使用基于模拟的功效分析为复杂生态系统中海洋捕捞庇护区的成本有效且统计可靠的监测提供信息的价值。研究人员在此采用基于模拟的功效分析，在样本量和效应量范围内比较两种采样方案（即CI和ELCI）评估RCAs恢复的能力。研究表明，纳入早期数据对于检测RCA有效性至关重要；若起始密度与预期相反（即RCAs外鱼类密度高于RCAs内），S错误可普遍影响结果；高M错误影响低（但现实）样本量。该研究表明，周全的研究设计，尤其是纳入固定地点的早期数据，对于有效的海洋捕捞庇护区监测至关重要；在某些系统中，小效应量——即使庇护区内增加高达20%——在可行采样强度下可能无法被可靠检测。
