西藏东南部样带海拔梯度下尼泊尔桤木(Alnus nepalensis)群落组成分化以物种周转为主导

《Diversity》：Compositional Differentiation of Alnus nepalensis Communities Along a Sampled Elevational Sequence in Southeastern Xizang, China, Is Dominated by Species Turnover

　　

摘要
山地研究常强调丰富度趋势，然而相似的丰富度可能掩盖物种更替。本研究旨在确定西藏东南部五个海拔带（2053–2463 m）中尼泊尔桤木(Alnus nepalensis)群落的组成分化究竟主要受丰富度丧失驱动还是受 assemblage 周转（assemblage turnover）驱动。研究人员分析了15个20 m × 20 m样地，采用存在–不存在Jaccard相异性（presence–absence Jaccard dissimilarity）分析，并以生活型标准化矩阵（life-form-standardized matrices）进行敏感性分析。群落组成在各海拔带间差异显著（PERMANOVA：F = 14.15，R2 = 0.850，p = 0.0001），而多元离散度无显著差异（PERMDISP：p = 0.617）；该模式在标准化生活型矩阵中依然存在（R2 = 0.910）。尽管样地丰富度随海拔升高而下降，周转（turnover）解释了相邻海拔带间平均样地对Jaccard相异性的71.0%–95.3%。Bootstrap置信区间支持所有四个过渡带中周转占主导地位，其中2053–2158 m处的不确定性最大（95% CI：54.6–87.3%）。2158 m与2247 m海拔带具有相同的合并丰富度，却显示出最高的周转比例（95.3%，95% CI：93.4–97.1%）。尼泊尔桤木在整个序列中均有分布，而与海拔相关的灌木和草本植物贡献了大部分分类学分化。这些结果表明，在该采样序列内，物种更替——而非单纯的丰富度丧失——是群落分化的主要驱动因素，尽管海拔与局地环境的混杂使得因果推断无法推广。

**西藏东南部样带海拔梯度下尼泊尔桤木群落组成分化的驱动机制——基于物种周转视角的解读**

**1. 研究背景与科学问题**

山地梯度能够在较短的地理距离内压缩气候、土壤和地形的变异，为生物多样性研究提供了理想的天然实验条件。然而，海拔并非单一的生态变量——温度、水分、土壤、坡度、干扰和林分历史往往协同变化，混淆了因果解释。海拔研究的核心贡献因此不在于简单呈现海拔与丰富度的相关性，而在于清晰描述组成如何沿样带变化以及推断的适用范围。

Alpha多样性（alpha diversity）无法单独捕捉组成变化的形式。两个地点可能具有相同的物种丰富度，但类群几乎完全不同；而物种贫乏的地点可能只是丰富地点的嵌套子集。将beta多样性（beta diversity）分解为周转（turnover）和丰富度差异（richness-difference）组分，能够区分这些替代情形。周转表示 assemblage 之间的物种更替；丰富度差异表示一个 assemblage 是另一个的子集。这一区分对保护具有实际意义：周转产生地点间的互补性，而丰富度差异则意味着冗余。

尼泊尔桤木(Alnus nepalensis D. Don)群落非常适合解决上述概念问题，因为该群落兼具三个关键特征：海拔与局地历史相混杂、冠层存在可能掩盖下层物种更替、以及该物种在森林恢复中的广泛使用意味着 assemblage 间究竟是相互更替还是仅形成嵌套子集对恢复目标的界定具有现实影响。西藏东南部位于东喜马拉雅–横断山生物地理过渡带附近，陡峭的地形和季风影响形成了精细尺度的生境异质性。然而，该区域尼泊尔桤木群落的组成组织方式尚不清楚，尤其是相邻海拔带间的差异主要反映渐进式物种丧失还是相互更替仍属未知。

基于此，本研究提出三个相互关联的假设：第一，在排除带内离散度不等的因素后，各海拔带的物种组成应存在差异；第二，周转应解释大部分相邻海拔带的beta多样性，包括丰富度相似的过渡带；第三，广泛分布的尼泊尔桤木冠层应与限定特定海拔带的窄域类群共存。

**2. 关键技术与方法**

研究样地位于中国西藏东南部林芝市巴宜区和波密县，沿雅鲁藏布江河谷设置五个海拔带（2053、2158、2247、2348和2463 m），每个海拔带设置3个20 m × 20 m样地，共15个样地，构成嵌套式采样设计（nested sampling design）。研究采用存在–不存在Jaccard相异性作为全群落分析的主要指标，以生活型标准化矩阵作为敏感性分析。核心分析方法包括：置换多元方差分析（PERMANOVA）用于检验组成差异，多元离散度置换分析（PERMDISP）用于检验组内离散度异质性，Baselga分解法将Jaccard相异性分解为周转与丰富度差异组分，Bootstrap重抽样（10,000次）量化周转比例的置信区间，以及Dufrêne–Legendre指示物种分析（IndVal）结合错误发现率（FDR）校正识别各海拔带关联类群。所有统计分析使用Python 3.12.13完成。

**3. 研究结果**

**3.1 采样完整性与丰富度背景**
样地丰富度总体随海拔升高而下降（r = ?0.544，R2 = 0.296，p = 0.036），但Shannon多样性无显著关联（r = ?0.310，p = 0.261）。最具信息量的是两个中间海拔带：2158 m和2247 m的合并丰富度完全相同（各16种），为检验“相似丰富度是否掩盖物种身份变化”提供了直接对照。Good's coverage范围为0.998–1.000，Chao1估计与观测丰富度接近，表明2463 m带仅6种的物种库并非采样不足所致。

**3.2 组成差异显著**
PERMANOVA结果表明，总Jaccard变异的84.98%归因于五个采样带间的差异（F = 14.15，R2 = 0.850，p = 0.0001）。PERMDISP不显著（p = 0.617），排除了离散度不等作为主要解释的可能性。主坐标分析（PCoA）显示样地按海拔带分离。

**3.3 相邻过渡带以物种周转为主导**
相邻海拔带间平均样地对Jaccard相异性为0.589–0.893，所有四个过渡带中周转均大于丰富度差异组分，占总相异性的71.0%–95.3%。2158–2247 m过渡带最具决定性意义：两带合并丰富度相同，但周转比例高达95.3%（95% CI：93.4–97.1%），表明相似的丰富度掩盖了assemblage 间的大量相互更替。

**3.4 分类学结构**
尼泊尔桤木在整个序列中持续存在，但伴生灌木和草本呈现较窄的海拔分布。九种类群通过FDR校正的指示物种阈值（q ≤ 0.05），其中Microstegium nudum和Berberis gyalaica的指示值达1.00。广泛分布的冠层物种与窄域伴生类群共同构成了“冠层连续、下层更替”的二层结构。

**3.5 敏感性分析与推断边界**
当定量记录按生活型标准化后，等权乔木–灌木–草本复合矩阵产生R2 = 0.910（p = 0.0001），各生活型单独分析R2为0.858–0.961，证实存在–不存在结果并非源于将不相容的测量单位强行合并。

**4. 讨论与结论**

研究表明，物种周转而非单纯丰富度丧失是该采样序列内群落分化的主要驱动力。2158 m与2247 m的对比提供了最有力的内部证据：相同的合并丰富度伴随着超过95%的周转比例，直接证明丰富度无法预测组成身份。该发现警示：将共享同一优势树种的林分视为组成上可互换的做法存在风险。

最高的2463 m带呈现特殊的模式——物种库收缩与更替并存。该带仅含6个合并物种，但记录个体数达1280，Good's coverage为1.000，Chao1与观测丰富度相等，表明其低丰富度并非采样稀缺所致。Jaccard相异性达到0.893，其中90.7%归因于周转。这种“收缩+更替”的组合意味着上带并非简单地是其相邻地带的贫化样本。

Beta多样性分解对保护具有重要意义。如果高海拔带仅仅是低海拔带的嵌套子集，保护丰富地点即可捕获大部分区域物种库。然而本研究中周转产生了互补的assemblage，特别在两个中间带和最高过渡带。因此，沿序列的空间代表性比依赖单一丰富度最大值更为合理。

研究的主要局限在于每个海拔带缺乏重复局地。虽然有15个植被样地捕获了局地内变异，但独立海拔级单位的有效数量仅为五个。海拔因此与局地特定条件不可分离。研究支持“采样带间存在系统性差异”的结论，但无法估计区域尺度的海拔效应或分离因果机制。Bootstrap重抽样量化的是观测样地的条件不确定性，而非可能局地间的设计不确定性。

结论：物种更替是西藏东南部五个采样尼泊尔桤木海拔带间组成分化的主导形式。基于存在–不存在数据的群落组成在各带间差异显著，且无可检测的全群落离散度异质性，该结果在乔木、灌木和草本数据分别标准化后依然成立。相邻带周转比例为71.0%–95.3%，包括丰富度相等海拔带间异常高的更替。广泛分布的尼泊尔桤木维持了冠层连续性，而窄域伴生类群区分了各采样群落；最高带还表现出物种库收缩。这些发现证明，在采样序列中，仅凭丰富度低估了空间生物多样性变异。由于海拔与局地相混杂，研究结论是组成性的而非因果性的：数据确立了五个局地间的分化和物种更替，但不支持推广为普遍的海拔机制。