保护区对阔叶落叶林林下植物多样性的对比效应

《Biodiversity and Conservation》：Contrasting effects of protected areas on understory plant diversity in broadleaved deciduous forests

　　

保护区对于保护森林生态系统和维持其生态完整性至关重要。然而，它们在意大利林下植物物种保护中的作用仍未得到充分研究。研究人员分析了1980年至2020年间在意大利阔叶落叶林中采样的5,423个植被样方。使用倾向性评分匹配（PSM），研究人员控制了潜在混杂因素，如海拔、气候、与道路的接近程度和样方大小，从而能够无偏评估保护状态对林下植物多样性的影响。研究人员比较了保护区和未保护区的alpha多样性和gamma多样性，同时考虑了总物种数以及分别针对森林栖息地或开放栖息地相关的林下植物物种。alpha多样性的差异使用广义线性模型（GLM）评估，而gamma多样性则通过量化保护区独有、未保护区独有以及两者共有的物种数量来评估。保护区通常倾向于支持更高的alpha和gamma多样性水平，尽管模式因栖息地类型而异。alpha多样性在保护区内普遍较高，但仅在部分栖息地中观察到统计学显著差异。同样，森林关联物种在所有栖息地类型中在保护区内显示出比外部更高的gamma多样性，而开放栖息地物种的模式变化较大，在保护区和未保护区之间并不一致。总体而言，研究结果表明，意大利保护区代表了林下植物多样性的重要部分，特别是对于森林关联物种。然而，栖息地特异性组合的存在表明，在保护区边界之外可能仍需要额外的保护努力。

### 论文解读：保护区对阔叶落叶林林下植物多样性的对比效应

#### 研究背景与问题
保护区是全球生物多样性保护的核心工具，旨在维持森林生态系统及其生态完整性。然而，在意大利，保护区对林下植物（understory plants）多样性的保护作用尚未得到充分研究。林下植物是温带森林维管植物多样性的主要组成部分（通常超过80%），对干扰响应迅速，并作为森林状况的可靠指示物种。现有研究常因数据非随机采样、缺乏对照设计而难以评估保护区的真实效应。此外，保护区选址常偏向偏远或低生产力区域，导致对生物多样性代表性的评估存在偏差。因此，需借助反事实分析（counterfactual analysis）和匹配技术（matching techniques），在控制环境混杂因素（如海拔、气候、道路距离）后，客观比较保护区内外的物种多样性。本研究利用意大利全国性植被样方数据集，结合倾向性评分匹配（Propensity Score Matching, PSM）和广义线性模型（Generalized Linear Models, GLMs），系统评估保护区对阔叶落叶林林下植物alpha多样性（局域尺度）和gamma多样性（区域尺度）的影响，并分别分析森林关联种（forest-associated species）和开放栖息地种（open-habitat species）的响应差异。研究假设：保护区局域尺度上森林关联种多样性更高，而区域尺度上开放栖息地种的物种库在保护区内外部仅部分共享。论文发表在《Biodiversity and Conservation》。

#### 关键技术方法
研究人员采用以下主要技术方法：
1. **数据来源与处理**：从意大利全国植被样方数据集（Alessi et al., 2023）中筛选出1980-2020年间采样的5,423个阔叶落叶林样方，样方面积100-400 m²。
2. **栖息地分类**：基于EUNIS专家系统（v2021-06-01）和硬噪声聚类模型（hard noise clustering）将样方归入8个EUNIS Level 3栖息地类型（如酸性栎林、山毛榉林等）。
3. **保护状态判定**：利用世界保护区数据库（WDPA）划定意大利保护区边界，将样方分为保护区内（3,638个）和保护区外（1,785个）。
4. **倾向性评分匹配（PSM）**：以保护状态为处理变量，匹配样方大小、海拔、年均温、年降水量、距森林道路最小距离和栖息地类型，采用最近邻法（ratio=1, caliper=0.25）消除选址偏倚。
5. **多样性分析**：alpha多样性采用Wilcoxon秩和检验与GLMs（泊松分布）分析；gamma多样性通过量化保护区内独有、保护区外独有及共有物种数评估，并绘制稀薄化曲线（rarefaction curves）检验采样完整性。

#### 研究结果
**Alpha多样性**
通过GLMs分析，保护区内森林关联种的alpha多样性普遍高于保护区外，仅在酸性栎林（acidophilous Quercus forests）中例外。开放栖息地种的alpha多样性在保护区内也呈较高趋势，但在河岸林（riparian forests）、温带硬木河岸林（temperate hardwood riparian forests）和栎类中生落叶林（Quercus mesic deciduous forests）中较低。Wilcoxon检验显示，仅峡谷林（ravine forests）、非酸性土山毛榉林（Fagus forests on non-acid soils）和暖温性落叶林（Thermophilus deciduous forests）中保护区内alpha多样性显著更高。

**Gamma多样性**
保护区内总gamma多样性普遍高于保护区外，但河岸林和栎类中生落叶林例外（如非酸性土山毛榉林：保护区内315种 vs. 保护区外277种）。物种组成上，21%的物种仅出现于保护区内，16%仅出现于保护区外，63%共有的。森林关联种在所有栖息地类型中均显示保护区内gamma多样性更高；开放栖息地种在多数栖息地中保护区内更高，但非酸性土山毛榉林和栎类中生落叶林中相反（保护区外独有种更多）。共记录32种意大利特有植物，其分布模式与总gamma多样性类似。

#### 讨论与结论
**讨论**：与先前研究（如Halpern & Spies, 1995）认为保护区或未管理森林中林下植物alpha多样性较低不同，本研究显示保护区内alpha多样性一般更高，主要由森林关联种驱动。这可能源于保护区对高价值地点的战略性选择，以及减少砍伐等人为干扰后森林结构完整性的维持。然而，高alpha多样性不代表高保护价值，因为历史土地利用、管理实践差异也可能贡献于观察到的模式。保护区对多样性的代表性因栖息地类型而异：酸性栎林和河岸林中保护区内alpha多样性较低，可能由于传统放牧停止后开放栖息地种扩张，或河岸林线性结构使外部受干扰地维持更多适应种。gamma多样性保护区内更高的趋势与北美洲及欧洲其他研究一致，表明保护区对区域物种库的维持作用。特有物种的保护区内出现（如Helleborus bocconei）进一步支持保护区的保护价值。开放栖息地种在保护区内更高独有性可能受自然干扰（如树倒、放牧）驱动，但需注意其存在也可能反映边缘效应或管理遗留，而非更高生态完整性。研究局限性包括：未评估稀有性、威胁状态或系统发育独特性，且采样强度可能不均衡，但PSM缓解了部分偏差。未来需整合时空数据以评估保护区的长期有效性。

**结论**（翻译自原文Conclusion部分）：
本研究表明，意大利保护区是生物多样性保护的基石，特别是对于森林林下植物物种。研究结果揭示，保护区通常支持比未保护区更高的alpha和gamma多样性。此外，森林关联物种在所有栖息地类型中在保护区内表现出更高的多样性，强化了保护区在保护森林生物多样性中的作用。然而，研究结果也突出了重要的例外情况。一些保护栖息地中开放栖息地相关物种的较低多样性表明，土地利用变化可能对在管理或半自然景观中繁盛的物种产生负面影响。总体而言，本研究强调了保护区在代表植物多样性方面的关键作用，同时指出需要更细致的保护方法，以确保不同栖息地类型中多样化植物群落的保护。