陆生蜗牛短期影响调查中的空壳纳入：方法学差异及其内在成因

《Ecology and Evolution》：Snail Inventories for Short-Term Impacts: Methodological Differences and Their Underlying Causes

　　

在陆生腹足类群落研究中，一个根本性的方法问题是：调查时是否应与活体个体一并纳入空壳。在这项研究中，研究人员检验了纳入新鲜空壳（具有完整角质层（periostracum））可能在分析中引入的偏差。研究人员首先重新分析了一个关于割草和碎草对草地蜗牛群落短期影响的数据集。原始分析仅使用活体标本数据，而这里研究人员使用了所谓的合并数据集（combined dataset），包括活体个体和新鲜空壳。活体数据揭示了割草对蜗牛群落的显著负面影响，而合并数据中这些影响无法被检测到。与此同时，研究人员观察到物种间过度代表值（overrepresentation value，即某物种合并丰度与其活体丰度的比值）存在大幅度的种间变异，这代表了总个体数在物种间分布的扭曲。为探究这一现象的内在原因，研究人员在实验室条件下测定了11个物种壳体的物理抗性和化学抗性，并使用多阶段种群模型（multistage population model）模拟了生活史特征和壳体降解速率对过度代表值的影响。与预期相反，研究结果表明，过度代表值的这种变异不由壳体抗性解释，而是由壳体大小解释：较小物种表现出显著更高的过度代表值，尽管它们的壳体更弱、可能分解更快。这一关系可能由较小物种更快的生命周期和更高的周转率驱动，这一发现得到了种群模型的印证。研究人员总结认为，纳入空壳可以显著扭曲对群落即时状态的感知，并掩盖短期时间变化。研究人员建议，在检测短期影响时优先使用活体数据，并在科学出版物中明确说明和论证所选的数据类型。

陆生腹足类群落调查中，动物遗骸常被用作物种出现和丰度的记录。陆生蜗牛（不含蛞蝓）的外壳在个体死亡后可长期存留，因此同一栖息地中活体与死亡空壳往往共存。当研究关心个体数量时，是否将空壳纳入分析是一个根本性方法学问题。空壳会膨胀物种丰度，且死个体不再参与种群过程；不同物种空壳降解速率可能不同，因而偏差可能具有物种特异性。尽管先前研究已意识到这些潜在偏差，但以比较实验方式系统检验其对短期影响研究影响的证据仍然缺乏。本研究正是为填补这一空白而开展。

研究人员利用Farkas等人开展的一项草地管理短期影响研究（割草和碎草处理）的原始数据，将原始仅含活体的数据集扩展为包含活体与新鲜空壳的合并数据集，重新运行影响评估分析。结果显示，活体数据能够检测出割草对蜗牛丰度和物种丰富度的显著负面影响，而合并数据却无法检测到这些效应。随后，研究人员计算各物种的过度代表值，发现物种间差异巨大，最大接近12倍。实验室抗性实验和种群模拟表明，造成这种差异的主要原因不是壳体抗性，而是壳体大小及其所反映的生活史策略：较小物种壳体更脆弱、抗性更低，但因生活史更快、周转率更高，产生空壳的速率更快，因此反而具有更高的过度代表值。该研究揭示，纳入空壳会像“扭曲的镜子”一样扭曲群落即时状态，掩盖短期时间变化。

主要技术方法（原始数据来自匈牙利北部Bükk西部两处湿草甸，处理为对照、割草和碎草，采样于2021年6月末和2022年8月中旬）：用广义线性混合模型（GLMM）分别分析活体与合并数据的总丰度、物种丰富度和Pielou均匀度指数；用合并丰度/活体丰度计算物种过度代表值，并用二项分布广义线性模型检验影响因素；在实验室测定11个物种壳体的物理抗性（压碎载荷）和化学抗性（醋酸溶解比例），测量壳体大小（长×宽）和形状（宽/长）；构建Leslie矩阵（Leslie matrix）多阶段种群模型，模拟生活史参数对成体过度代表的影响。

3.1 保护管理措施的影响：活体数据与合并数据的比较
基于活体数据的广义线性混合模型显示，割草处理与采样时间的交互项显著影响个体数（χ²₂=6.993，p=0.030）、物种数（χ²₂=11.318，p=0.004）和Pielou均匀度指数（χ²₂=6.275，p=0.043）；割草导致物种丰富度和丰度下降、均匀度上升。合并数据中这些交互项均不显著（p>0.05），表明纳入新鲜空壳会掩盖管理措施的短期影响。

3.2 物种丰度的差异
各物种过度代表值（合并丰度/活体丰度）在1至接近12之间，多数物种为2～6。该结果表明，合并数据中不同物种被空壳放大的程度不同，导致群落内物种相对丰度格局发生扭曲。

3.3 抗性与壳体特征的关系
实验室测定显示，壳体大小与物理抗性（t=6.74，df=9，p<0.001）和化学抗性（χ²=23.50，df=1，p<0.001）均呈显著正相关；壳体形状与物理抗性无显著线性关系，与化学抗性存在非线性关系，圆球形和长形壳体的化学抗性较低。

3.4 与过度代表值的关系
壳体大小和两种抗性均与过度代表值显著相关，但方向与预期相反：抗性越高的物种过度代表值越低，壳体越小的物种过度代表值越高。因此，物种间过度代表差异不能用壳体抗性解释。

3.5 生活史差异的影响
Leslie矩阵模型模拟表明，死壳存续概率（P_d）越高、成体存活率（P_a）越低、生活史越快（幼体到成体的转化率P_ya低、繁殖力F_y高），死壳在稳定阶段分布中的过度代表越高。模型结果支持“小物种因快速生活史和高周转率而具有更高空壳产生速率”这一解释。

讨论部分指出，本研究数据来自同一生境类型内两个相邻地点，推广性受局限；但核心发现——数据类型选择强烈影响群落变化评估——值得重视。研究人员认为，抗性差异不是过度代表种间差异的主要原因；Sólymos等人提出的微生境偏好假说在本研究相对均一的采样生境中不太可能解释所观察到的差异。研究结论翻译为：纳入空壳可显著扭曲对群落即时状态的感知，并掩盖短期时间变化；在检测短期影响时，应优先使用活体数据，并在科学出版物中明确说明和论证所选数据类型。