模拟的卡廷加干旱森林采伐后土壤化学性质长期恢复轨迹的时间序列

《CATENA》：Modeled long-term recovery trajectories of soil chemical properties along a post-clearing Chronosequence in the caatinga dry Forest

　　

在卡廷加森林再生过程中，土壤化学性质以不同速率恢复。土壤有机碳（SOC）和碱饱和度（BS）沿时间序列逐渐增加。盐度相关属性和土壤酸度表现出较慢的长期稳定化。交换性Ca²⁺、Mg²⁺和K⁺增加，而Na⁺在再生过程中下降。短轮伐周期可能不足以在半干旱干旱森林中完全恢复土壤化学功能。

热带干旱森林是全球受威胁最严重的生态系统之一，皆伐（clear-cutting）破坏了养分循环、微气候稳定性和土壤功能。在卡廷加（Caatinga）干旱森林中，木质植被常作为“肥力岛”集中养分和有机质，清除植被会破坏这些积累核，延缓土壤恢复并增加空间异质性。水分限制、高蒸散和低风化速率进一步加剧退化，增加荒漠化风险。然而，以往研究多孤立评估土壤化学属性，缺乏跨再生时间序列（chronosequence）的综合分析，尤其在半干旱环境中，对采伐后不同化学属性稳定所需时间的信息有限。当前卡廷加地区约15年的轮伐周期仍被采用，但土壤化学功能是否以相似速率恢复尚不明确。因此，研究人员沿一条包含6个采伐后年龄（0.5、6、9、12、25、50年）和一片超过80年参考森林的时间序列，通过年龄特异性比较和非线性模型估计不同属性达到多个恢复阈值所需的时间，系统调查了土壤化学性质的恢复轨迹。该研究发表在《CATENA》。

研究人员在巴西伯南布哥州Floresta市的超旱生卡廷加森林区域选取7个样地，代表不同再生年龄。每个样地挖掘5个20×50×30 cm的沟槽（间隔约50 m），在干旱季节采集0–5、5–10、10–20 cm三个深度土壤样品，共104个样本。化学分析包括：饱和浸提液电导率（EC）、土壤pH（水提）、交换性Ca²⁺、Mg²⁺、Na⁺、K⁺（醋酸铵提取）、潜在酸度（H⁺+Al³⁺）、阳离子交换容量（CEC）、碱饱和度（BS）、总有机碳（TOC）、总氮（TN）和有效磷（P，Olsen法）。统计方法采用广义线性模型（GLM）评估再生年龄和深度效应，非线性逻辑模型和广义加性模型（GAM）拟合恢复轨迹，并估算50%、60%、70%、80%、90%、99%恢复阈值所需时间。

**3.1 再生时间序列中土壤化学性质的变化**
通过GLM分析，土壤pH随再生时间显著增加（χ²₍₆₎=42.64, p<0.001），从0.5年的5.5升至参考森林的6.6，深度效应不显著。总有机碳（TOC）显著增加（χ²₍₆₎=34.24, p<0.001），从9.77 g kg^-1升至16.15 g kg^-1，深度效应显著（χ²₍₂₎=268.95, p<0.001）。电导率（EC）呈下降趋势，深度特异性模型显示各层显著降低（p<0.001），从1.33 dS m^-1降至0.49 dS m^-1。交换性Ca²⁺、Mg²⁺、K⁺均随再生时间显著增加（p<0.001），分别从2.77、0.91、0.25 cmol_c kg^-1升至5.49、3.12、0.51 cmol_c kg^-1；而交换性Na⁺显著下降（p<0.001），从0.27降至0.08 cmol_c kg^-1。总和碱（SB）和碱饱和度（BS）显著增加（p<0.001），分别从4.19 cmol_c kg^-1和52.53%升至9.20 cmol_c kg^-1和85.69%。阳离子交换容量（CEC）缓慢增加（p<0.001），从8.02升至10.75 cmol_c kg^-1。总氮（TN）和有效磷（P）也显著增加（p<0.001），分别从1.61 g kg^-1和3.93 mg dm^-3升至2.78 g kg^-1和9.88 mg dm^-3。相关性分析表明，TOC、TN、Ca²⁺、K⁺与pH正相关，EC与这些属性负相关，Na⁺与EC正相关、与pH和TOC负相关。

**3.2 森林再生过程中土壤化学性质的恢复轨迹**
非线性模型显示，不同属性恢复速率差异近一个数量级。pH渐近估计值为7.07，90%恢复约需47年，99%恢复约需164年（外推）。TOC恢复较快，90%和99%恢复分别需22.5和54.8年，渐近值为15.71 g kg^-1。EC呈下降趋势，50%和90%减少分别需29.8和68年。交换性Ca²⁺渐近值为5.20 cmol_c kg^-1，90%和99%恢复需21.3和45.2年；Mg²⁺渐近值2.84 cmol_c kg^-1，需25.3和44.9年；K⁺渐近值0.50 cmol_c kg^-1，需29和62年。Na⁺下降，50%、90%、99%减少分别需30、67、78.7年。BS恢复最快，90%和99%恢复需16和36.1年，渐近值85.97%。潜在酸度下降，50%、90%、99%减少分别需9.6、19.8、23.5年。CEC恢复最慢，90%和99%恢复需117和>300年（外推）。TN和P分别需47.7和52.4年达到90%恢复，99%恢复需109.6和101.1年。

**3.3 长期恢复动态、研究局限及对干旱森林恢复的启示**
讨论指出，养分相关属性恢复快于调控属性（pH、EC、CEC），CEC最慢，受限于有机质稳定化。深度效应显著，EC和Na⁺在表层下降更快，Mg²⁺和K⁺在深层稳定更慢。当前15年轮伐周期远低于多数属性恢复所需时间，但应视为指示性而非确定性阈值。研究局限包括：每个年龄仅一个样地（空间替代时间设计），年龄与地点混淆；仅评估化学属性，未涉及物理、水文和生物过程；长期阈值属模型外推。

**结论**
土壤化学恢复在森林再生过程中以不同速率进行。与养分有效性和有机质积累相关的属性（Ca²⁺、Mg²⁺、SB、BS、TOC、TN）恢复较快，而pH、EC、Na⁺，尤其是CEC，恢复较慢。不同深度响应也存在差异。由于每个再生年龄仅由一个样地代表，这些模式描述了该时间序列内的趋势，较长的稳定化时间属于模型外推而非观测恢复时间。即使在此前提下，仅靠森林再生似乎不能保证半干旱环境中土壤化学功能的完全恢复，这支持将土壤指标纳入干旱热带森林管理指南。未来需结合样地重复以及物理、水文和生物指标，在变化气候条件下验证模拟的恢复轨迹。