大气二氧化硫浓度长期降低后桦树和针叶林步甲科甲虫群落的变化

《Forest Ecology and Management》：Changes in carabid beetle communities following long-term reductions in atmospheric sulphur dioxide concentrations in birch and conifer forests

　　

一项为期26年的步甲科甲虫研究追踪了极端二氧化硫(SO₂)污染后的森林恢复过程。随着排放量下降，林地群落取代了耐干扰的甲虫。森林演替和微气候，而非二氧化硫(SO₂)，主要驱动了功能变化。污染减轻后系统发育聚类加剧。恢复中的生态系统经历长期重组，而非简单的基线回归。

**论文解读：长期二氧化硫污染减少后步甲科甲虫群落的变化——基于26年监测数据的多维分析**

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

全球森林生态系统正受到空气污染、气候变化、土地利用历史等多重人为压力的共同影响，其中大气二氧化硫(SO₂)污染曾在中欧地区造成严重的土壤酸化和森林衰退。在捷克共和国与德国交界处的“黑三角”地区，20世纪晚期的高浓度SO₂（研究初期年均浓度持续超过40 μg/m3）导致挪威云杉（*Picea abies*）大面积死亡，超过40,000公顷生产林被毁，迫使大规模皆伐并改用先锋桦树（*Betula pendula*）等替代树种进行紧急造林。尽管自1980年代末期以来SO₂排放大幅下降，但生态系统恢复并非简单的线性回归，而是受到演替遗留效应和持续环境约束的复杂重组过程。步甲科甲虫（Carabidae）作为与土壤条件紧密相关的敏感指示生物，其群落结构的变化能反映这些动态。然而，目前缺乏长期、多维度的监测数据来检验污染减少如何重塑群落组成、功能及系统发育结构。为此，研究人员利用捷克共和国Dě?ín砂岩高地26年（1989–2015）的步甲科甲虫数据集，旨在检验三个假设：(i) 活动密度和物种丰富度沿SO₂梯度呈现显著非线性变化；(ii) 环境过滤减弱将增加整体功能多样性，但驱动孤立性状的方向性变化；(iii) 群落组成在政策区间（以2002年国家空气质量政策为界）发生离散的多变量移位，表现为系统发育聚类加剧，特征性林地物种取代耐干扰类群。该研究发表于《Forest Ecology and Management》。

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

研究人员在捷克共和国Dě?ín砂岩高地（50°47′ N，14°13′ E）的台地（海拔450–610 m）设置9个样地：6个桦树（*Betula pendula*）为主、1个欧洲落叶松（*Larix decidua*）、1个欧洲赤松（*Pinus sylvestris*）和1个挪威云杉（*Picea abies*）。自1989年至2015年连续26年，使用带顶盖的玻璃陷阱（0.7 L，含2–4%甲醛）采集步甲科甲虫，每年4月1日至10月30日每周收集一次，每个样地10个陷阱线性排列，样品合并处理。大气SO₂浓度数据来自捷克水文气象研究所（CHMI）监测站。统计方法包括：广义加性混合模型（GAMMs）分析环境驱动因素；主成分分析（PCA）评估性状共线性；非度量多维尺度分析（NMDS）和置换多元方差分析（PERMANOVA）检验群落组成变化；指示物种分析（IndVal）识别特征种；系统发育多样性采用Faith's系统发育多样性标准化效应值（SESpd）计算，基于两种零模型（丰富度约束和独立交换算法）。

**3. 研究结果**

**3.1 活动密度和物种丰富度对长期SO₂浓度的响应**
通过GAMMs分析（n=5455个采样事件），发现平均周SO₂浓度对活动密度（edf=3.94，F=2.17，p=0.044）和物种丰富度（edf=2.41，F=3.32，p=0.018）均产生显著非线性效应。此外，降水和风速也显著影响两者，且季节物候效应显著。优势树种作为生境过滤器，落叶松和松树显著提高物种丰富度和活动密度，而云杉显著降低物种丰富度。

**3.2 功能多样性（Rao's Q）和性状响应**
PCA显示性状空间由两个正交轴解释（PC1=73.8%，PC2=18.4%），分别代表体型与扩散能力的权衡以及营养策略。GAMMs表明，SO₂浓度显著降低整体功能多样性（Rao's Q，edf=1.011，F=5.079，p=0.025），但未对群落加权平均（CWM）性状（营养策略、扩散能力、体型）产生显著方向性过滤。相反，CWM性状和Rao's Q主要由当地生境（优势树种）、微气候（降水、风速）和物候驱动。

**3.3 群落组成和优势种结构在空气质量政策区间变化**
PERMANOVA显示，2002年前后两个时期群落组成显著分离（F=7.11，p<0.001，R2=0.025），年均SO₂浓度解释显著变异（F=5.26，p<0.001，R2=0.019）。NMDS排序（应力=0.116）支持这一移位，SO₂向量指向高污染期。指示物种分析表明，高污染期（2002年前）以耐干扰的广布种*Poecilus cupreus*（IndVal=0.46）、*Amara lunicollis*（0.56）和*Epaphius secalis*（0.31）为特征；低污染期（2002年后）以林地专化种*Abax parallelepipedus*（0.48）和旱生林地种*Amara makolskii*（0.36）为特征。优势种结构从高污染期以*Carabus hortensis*（31.8%）和*Carabus violaceus*（17.8%）等为主，转变为低污染期*C. hortensis*占42.2%，而*Poecilus cupreus*等耐干扰种完全退出前12位，森林专化种*Abax parallelepipedus*上升至第7位（3.1%）。

**3.4 系统发育结构**
尽管物种丰富度在政策区间稳定（均值11.4 vs 11.5），系统发育聚类显著加剧。丰富度约束的SESpd均值从-1.16降至-1.74（Wilcoxon检验，p<0.001），独立交换算法结果一致（p<0.001）。宏观时间尺度GAMMs（n=243）表明，年均SO₂浓度降低显著非线性驱动SESpd下降（丰富度：edf=2.29，F=6.36，p<0.001；交换：edf=2.39，F=7.40，p<0.001），表明污染减轻使群落从历史随机或过度分散状态转向显著的系统发育聚类。

**4. 讨论与结论**

讨论部分指出，SO₂浓度下降并未导致简单恢复，而是引发群落分类、功能和系统发育结构的复杂非线性重组。活动密度和丰富度对SO₂的非线性响应表明，步甲动态受森林演替和污染减轻的平行过程驱动。功能多样性受SO₂抑制，但性状过滤主要由当地生境和微气候主导，这体现了极端污染与功能维度之间的解耦。群落重组通过耐干扰种被林地专化种取代实现，而系统发育聚类加剧则反映了污染减轻后更强的谱系筛选。结论翻译如下：我们的26年分析表明，在严重大气SO₂污染后的生态系统恢复并非简单的、单调的回归到干扰前基线，而是一个复杂的、非线性的重组过程。尽管过去的高污染时期偏向耐干扰的广布种，但污染减轻时期以优势等级结构的显著重塑和特征性林地群落的重新建立为特征。这种转变是由广泛林地分类群（如*Abax parallelepipedus*和扩张的大型*Carabus*物种）以及稀有专化种（如*Amara makolskii*）的整合驱动的，后者利用了次生先锋林成熟过程中产生的复杂微气候镶嵌。此外，尽管物种丰富度无明显方向性波动，但系统发育聚类随着大气污染降低而显著加剧，表明在谱系水平上发生了强烈的环境过滤。这些发现强调了多维监测（整合分类、功能和系统发育指标）对于检测丰富度指标经常掩盖的进化历史和功能分化中隐藏变化的重要性。