地中海喀斯特地区黎巴嫩雪松（Cedrus libani）造林时间序列中土壤碳与微生物功能的多维格局

《Biology》：Soil Carbon and Microbial Functioning Reveal Multidimensional Patterns Across a Cedrus libani Afforestation Chronosequence in Mediterranean Karst

　　

：地中海喀斯特地区的造林评估往往侧重于植被和土壤碳，而对微生物活性的关注较少。研究人员考察了土耳其南部一个由四个样地组成的黎巴嫩雪松（Cedrus libani）时间序列的表层土壤性质。2024年10月，研究人员在一个退化参考样地和三个林龄分别约为9年、14年和24年的人工林样地采集了表土（0–10 cm深）样品（每个样地采集15个样品）。造林样地通常表现出更高的土壤有机碳（SOC）、全氮（TN）、团聚体稳定性和微生物生物量碳（MBC）。造林样地的平均MBC范围为288至536 mg C kg?1，而参考样地为180 mg C kg?1。基础呼吸（BR）较低（0.30–0.36对比0.96 μg CO2-C g?1 h?1），生物量特异性呼吸（qCO2：0.54–1.39对比5.41 mg CO2-C g?1 MBC h?1）也较低。微生物熵（qMic）无显著差异。合并的MBC–呼吸关联较弱，而样地校正后的关联为正。多变量格局和样地的非单调排序凸显了土壤异质性。由于造林林龄与样地身份完全混淆，这些单次观测结果无法确立因果性的时间效应。尽管如此，整合结构性、化学性和微生物指标可为表征这些喀斯特样地的土壤状况提供互补信息。

论文解读文章

一、研究背景与问题

在地中海地区，造林被广泛用作修复退化景观、增强土壤碳固存、减少侵蚀和提升生态系统韧性的基于自然的干预措施。然而，仅凭树木成功定植并不能推断恢复成效，因为地下功能的恢复可能滞后于植被发育。土壤结构、养分循环、碳稳定化和微生物功能共同决定了恢复植被能否转化为生态系统状况的持续改善。因此，评估恢复成效需要综合考量物理、生物地球化学和微生物过程的协同恢复，而非依赖单一土壤属性的变化。

地中海喀斯特地貌为功能恢复提供了极为复杂的背景。喀斯特地形约占全球陆地表面的12%，其特征是薄而空间不连续的土壤覆盖于碳酸盐基岩之上。土壤主要积聚在溶蚀坑、基岩裂隙和残余覆盖层中，但浅薄的剖面、高砾石含量、低有机质和有限持水能力严重制约了土壤发育和植被建立。由于碳酸盐风化产生的不可溶残余物极少，矿质土壤形成极为缓慢，这放大了有机质在维持土壤结构、养分保持和水分储存中的关键作用。因此，植被丧失及随之而来的有机质和团聚体稳定性损失，可能将这些脆弱生态系统推向严重且不可逆的退化。

在此背景下，本研究旨在评估土耳其南部退化地中海喀斯特地区四个样地组成的黎巴嫩雪松（Cedrus libani）造林时间序列中表层土壤功能格局。研究人员整合了结构性、生物地球化学和微生物指标，并明确评估了观测响应是否对碳酸盐含量和土壤pH的变化保持稳健。研究提出三个假设：（i）造林样地将表现出比退化参考样地更高的团聚体稳定性以及更高的SOC、TN和MBC；（ii）微生物代谢指标将在退化参考样地与造林样地之间存在差异，尤其表现为更低的生物量特异性呼吸；（iii）结构性和微生物属性对土壤异质性的敏感性将低于SOC和TN。该研究发表于《Biology》期刊。

二、主要技术方法

研究人员在土耳其南部梅尔辛省埃尔代姆利区Sorgun–Toros高地开展了观测性空间换时间时间序列研究。研究区包含四个样地：一个未造林退化参考样地（C0）和三个分别于2015年、2010年和2000年建立的人工林样地（C9、C14和C24），采样时林龄分别约为9年、14年和24年。每个样地设置15个空间分布采样点，共60个采样点，于2024年10月采集0–10 cm表层矿质土壤。物理化学分析包括：比重计法测定土壤颗粒组成、Walkley–Black湿氧化法测定SOC、凯氏消解法测定TN、Scheibler钙计法测定CaCO₃、湿筛法测定团聚体稳定性（AS）、环刀法测定容重（BD）。微生物分析包括：氯仿熏蒸–提取法（CFE）测定微生物生物量碳（MBC）、碱液吸收法测定基础呼吸（BR），并计算微生物熵（qMic）和代谢熵（qCO₂）。统计分析采用单因素方差分析（ANOVA）、Kruskal–Wallis检验、PERMANOVA、主成分分析（PCA）、Spearman相关分析和探索性协方差分析（ANCOVA）。

三、研究结果

3.1 样地相关的土壤物理化学性质变异

土壤物理化学性质在四个时间序列样地间存在显著差异。经Benjamini–Hochberg错误发现率（FDR）校正后，所有物理化学变量在样地间均保持显著差异。效应量从中等至极强，最强的响应出现在碳酸钙（CaCO₃，ω² = 0.640）、土壤pH（ε² = 0.620）、全氮（TN，ε² = 0.544）、电导率（EC，ε² = 0.532）和团聚体稳定性（AS，ε² = 0.510）。

3.1.1 土壤物理性质

颗粒组成差异主要源于C24样地的独特质地，其砂粒含量较高而粉粒和黏粒含量较低。C0、C9和C14在任一粒级上均无显著差异。容重（BD）在样地间存在差异（ω² = 0.102），C0与C24之间存在显著差异，C9和C14居中。团聚体稳定性（AS）差异更大（ε² = 0.510），C0显著低于各造林样地，而C9、C14和C24之间无显著差异。

3.1.2 土壤化学性质

土壤pH在样地间差异显著（ε² = 0.620），C9最高，C24低于其他所有样地。EC差异显著（ε² = 0.532），C14和C24低于C0。SOC差异显著（ε² = 0.313），C24超过C0和C14但未超过C9。TN差异显著（ε² = 0.544），C24超过C0和C14，C9超过C0。C/N比差异显著（ω² = 0.391），C14高于C9和C24。CaCO₃差异显著（ω² = 0.640），C24低于所有其他样地，C9低于C14。

3.2 样地相关的土壤微生物功能变异

四个微生物指标中有三个在FDR校正后差异显著。qCO₂效应量最大（ω² = 0.865），其次是BR（ω² = 0.566）和MBC（ω² = 0.337）。qMic不显著（p = 0.0801）。MBC在样地间差异显著（ω² = 0.337），C9和C24均超过C0，C24超过C14。BR差异显著（ω² = 0.566），C0超过所有三个造林样地，造林样地之间无显著差异。qCO₂差异极强（ω² = 0.865），C0最高，C24最低。总体而言，造林样地表现为更高的MBC、更低的BR和显著更低的qCO₂，而qMic变异较小。

3.3 整合多变量格局

基于11变量矩阵的PERMANOVA显示样地间多变量组成差异显著（Pseudo-F = 16.788，R² = 0.473，p = 0.0001），PERMDISP不显著，表明结果并非仅由组内离散度不等驱动。所有六组两两比较在Holm校正后均显著，C0–C24分离最大（R² = 0.507），C9–C14分离最小（R² = 0.114）。PCA中PC1解释41.2%、PC2解释20.0%的方差。C0和C24占据排序图相对区域，C9和C14居中且部分重叠。SOC、TN、团聚体稳定性和MBC向C24方向负载，pH、CaCO₃、EC和BR向相反方向负载。

3.4 土壤物理化学与微生物属性间的关系

Spearman相关分析显示SOC、TN和MBC之间存在强正相关（SOC–TN ρ = 0.936；SOC–MBC ρ = 0.923；TN–MBC ρ = 0.873）。pH与CaCO₃正相关（ρ = 0.680），与SOC、TN和MBC负相关。105对相关中有54对在FDR校正后显著。样地校正后的偏相关分析中，SOC–TN（偏ρ = 0.931）、SOC–MBC（0.881）、TN–MBC（0.798）和pH–CaCO₃（0.544）仍保持强相关，39对偏相关在FDR校正后显著。

3.5 林龄与基质

C24占据独特的碳酸盐–pH区域，CaCO₃与MBC在全部样本中呈负相关（Pearson r = ?0.494，p < 0.001），CaCO₃每增加一个百分点，MBC平均减少4.59 mg C kg^?1。C24显著的质地差异更可能被解释为预先存在的样地异质性而非24年造林发育的后果。探索性ANCOVA中，样地/类别×协变量交互项对log₁₀(MBC)、log₁₀(BR)和log₁₀(qCO₂)显著，表明共同斜率假设不成立，相应调整均值不可用于推断。

3.6 树木生长

24年生林分的树木平均胸径16.94 cm、树高7.00 m，约为两个较年轻林分的两倍。但生长并非随林龄单调递增：14年生林分的胸径（8.48 cm）和树高（4.08 m）均为最小，低于9年生林分，与该样地较高的碳酸盐含量和较低的有机碳一致。

四、讨论总结

造林样地与退化参考样地之间最清晰的物理差异是更高的团聚体稳定性，并伴随容重的适度降低。团聚体稳定性在控制CaCO₃和pH的敏感性分析中保持相对一致，表明这两个土壤变量本身不能轻易解释观测到的样地差异。颗粒组成差异主要反映固有的样地异质性，C24的砂质质地不太可能由24年森林发育造成。

SOC和TN在造林样地中较高，与凋落物输入、根系周转、根际沉积、干扰减少和团聚作用改善相关的有机质输入和保持增加一致。但这些差异不能独立于样地条件被解释为林龄效应。C24在质地、CaCO₃和pH上差异显著，而C14并未一致表现出高于C9的SOC、TN或MBC。探索性协变量调整降低了SOC和TN的样地/类别效应量，但这些模型无法区分造林林龄与样地身份。

MBC在造林样地中较高且与SOC强相关，与更大的底物可利用性一致。相反，退化参考样地的BR约为造林样地的2.7–3.3倍，尽管其SOC和MBC较低，表明底物可利用性本身无法解释呼吸格局。qCO₂在造林样地中较低，但其由BR/MBC计算得出，并非独立响应或碳利用效率的直接度量。qMic的总体检验不显著（p = 0.0801），这可能是由于MBC和SOC成比例差异导致其比值差异较小，而非证明微生物碳分配未发生变化。

多变量分析揭示了样地间的协调差异，但C9和C14部分重叠，多个变量未随名义林龄单调变化。14年生样地尤为关键：尽管名义林龄居中，但其CaCO₃浓度（47.56%）和C/N比最高，SOC、TN和MBC与退化参考样地无显著差异，树木生长最差。这一偏离单调序列的现象凸显了在解释时间序列格局时考虑样地异质性的重要性。

五、研究结论

黎巴嫩雪松（Cedrus libani）造林样地在表层土壤的结构性、生物地球化学和微生物属性上与退化参考样地存在显著差异。造林土壤通常具有更高的SOC、TN、团聚体稳定性和MBC，而基础呼吸和qCO₂较低。这些差异不随名义林龄单调变化，SOC、TN和MBC的变异部分与局部土壤条件相关。碳状况、土壤结构和微生物属性的协调差异表明，整合多种土壤指标比单独依赖SOC能更全面地表征地下状况。造林样地较低的qCO₂表明在标准化培养条件下生物量特异性呼吸较低，但这一格局背后的生理机制仍未阐明。

这些发现应在观测性时间序列设计的约束下进行解读。每个名义林龄等级仅由一个样地代表，测量仅限于0–10 cm表层矿质土壤，且仅于2024年10月采样一次。因此，样地身份无法与造林林龄分离，观测到的差异应被视为样地特异的时间序列相关格局，而非林分发育的普遍因果效应。需要在对比性喀斯特环境中开展重复时间序列和长期重复测量，以确定这些格局是否代表黎巴嫩雪松（Cedrus libani）造林的普遍时间响应。