《Journal of Hazardous Materials Advances》：Ecological restoration mitigates heavy metal pollution in rubber monocultures: Quantifying risk reduction and anthropogenic source control
编辑推荐：
橡胶（Hevea brasiliensis）单作扩张从根本上重塑了热带景观，然而生态修复在缓解相关土壤重金属（HMs）污染方面的功效仍未得到充分量化。本研究采用空间换时间替代法（space?for?time substitution approach），评估了
橡胶（Hevea brasiliensis）单作扩张从根本上重塑了热带景观，然而生态修复在缓解相关土壤重金属（HMs）污染方面的功效仍未得到充分量化。本研究采用空间换时间替代法（space?for?time substitution approach），评估了中国西南地区四种森林类型中八种HMs的污染特征、生态风险及来源，这四种森林类型构成了从橡胶单作（RM）向成熟雨林（MF）延伸的两个平行恢复演替序列：自然恢复（RM–自然恢复林，NF–MF）和人工恢复（RM–人工恢复林，AF–MF）。实证分析显示，RM土壤中汞（Hg）和镉（Cd）显著富集，主要由遗留的农用化学品投入驱动。生态修复导致潜在生态风险大幅缓解：相对于RM，综合潜在生态风险指数（Nemerow Integrated Risk Index, NIRI）在AF中下降了94.8%，在NF中下降了78.7%，其中AF达到了类似于顶极雨林参考基准的风险基线。通过正矩阵分解模型（Positive Matrix Factorization, PMF）进行的源解析确定，人为活动——特别是农业实践（31.5%）和交通排放（26.7%）——是表土HMs积累的主要驱动因素，超过了自然成土作用的贡献。此外，偏最小二乘结构方程模型（Partial Least Squares Structural Equation Modeling, PLS?SEM）阐明了风险缓解的双重路径机制，表明恢复主要依赖于切断外源性人为输入，其次才是增强土壤微环境的缓冲能力。因此，提出了一个空间差异化的“恢复?源控制”框架，将种植园内部的主动恢复与运输走廊沿线植被生物屏障的建立相结合。这些发现为生态修复作为热带农业生态系统强有力的修复策略提供了定量证据，前提是必须将其与有针对性的源头减排方案相结合。
该论文针对橡胶单作扩张导致热带土壤重金属（HMs）污染风险加剧且缺乏定量修复依据的现状，系统评估了自然与人工恢复途径对污染风险的缓解效能及驱动机制。研究人员在中国西双版纳热带植物园（XTBG）选取了橡胶单作（RM）、自然恢复林（NF）、人工恢复林（AF）和成熟热带雨林（MF）四种森林类型，构建了两个平行的恢复演替序列（RM–NF–MF和RM–AF–MF），采用空间换时间替代法，结合地统计指数、受体模型及结构方程模型，揭示了生态修复对土壤HMs的调控机理。研究发现，生态修复特别是人工恢复，能显著降低由历史农用化学品残留导致的Hg和Cd富集风险，其效果主要通过切断外源输入与改善土壤微环境双重路径实现。基于此提出的“恢复?源控制”空间差异化治理框架，为实现热带种植区生产力与环境安全的可持续平衡提供了科学依据。该研究成果发表于《Journal of Hazardous Materials Advances》。
为开展此项研究，研究人员主要采用以下关键技术方法：基于空间换时间替代法设置样地，采集0–20?cm表层土壤样品；利用微波辅助酸消解结合电感耦合等离子体质谱（ICP?MS）测定包括Hg、Cd在内的八种HMs总浓度；计算富集因子（EF）、地累积指数（Igeo）及综合潜在生态风险指数（NIRI）评估污染水平；应用正矩阵分解模型（PMF）进行源解析；并构建偏最小二乘结构方程模型（PLS?SEM）量化恢复过程中各变量间的因果路径系数。
3. 研究结果
3.1 土壤性质与HMs浓度
研究结果显示，不同森林类型的土壤理化性质存在显著差异，随着恢复进程，全氮（TN）和土壤有机碳（SOC）浓度增加，而土壤pH呈下降趋势。在HMs浓度方面，Ni、As和Pb普遍低于云南省背景值，而Cr、Cu、Zn、Cd虽有超标但幅度有限，Hg则表现出广泛的累积特征，62.5%的样本超过背景值。RM表现出最高的Hg和Cd平均浓度，其变异系数（CV）分析表明RM中Cd（56.3%）和Hg（18.8%）具有中到高度的变异性，反映了强烈的人为干扰特征。
3.2 污染水平与生态风险评估
富集因子（EF）和地累积指数（Igeo）揭示了明显的污染格局。对于Hg，超过半数的样本表现为“中度富集”至“显著富集”，这一现象在RM中最为突出，RM土壤的Hg污染等级达到“中度到重度”。相比之下，其他金属大多处于轻微富集或无污染状态。综合潜在生态风险指数（NIRI）分析进一步证实，RM属于高风险系统（均值?>?160），而生态修复显著降低了风险：与RM相比，NF的平均NIRI降低了78.7%，AF更是降低了94.8%，其风险水平已与顶极雨林（MF）相当。
3.3 多元统计分析及源解析
主成分分析（PCA）提取了三个主成分，解释了71.93%的总方差，分别对应Cr?Ni?Pb组合、Cd?Hg组合及Cu?Zn组合。PMF模型进一步将污染源解析为四个因子：因子1（生物地球化学循环）由Cr、Ni和少量Hg构成；因子2（交通排放）以Cu和Zn为主，与道路邻近性相关；因子3（自然成土作用）包含As、Cr、Ni和Pb，浓度均低于背景值；因子4（农用化学品遗留）对Hg和Cd的贡献最大，直接对应橡胶单作期间的农药和化肥施用历史。
3.4 HMs与土壤理化性质的关系
相关性分析揭示了土壤理化性质对HMs行为的控制作用。全磷（TP）作为化肥施用的指标，与Cd和Hg呈极显著正相关，验证了农业源的输入。土壤有机碳（SOC）与Cr、Cu呈正相关，但与Pb和Hg呈负相关，表明有机质的络合作用对不同金属的影响存在差异。土壤pH与Cd呈正相关，但在对比分析中发现，RM较低的pH并未导致其Cd含量最高，说明外部输入量是主导因素。
3.5 HMs风险降低的驱动机制
基于PLS?SEM的分析量化了风险降低的双重路径。结果表明，生态修复首先通过显著减少人为HMs输入（路径系数?p?
4. 讨论与结论
在讨论部分，研究人员指出橡胶种植留下的Hg和Cd富集主要源于长期施用含Cd磷肥及有机汞杀菌剂的历史。尽管土壤pH影响金属迁移，但外部输入量远超土壤性质的调控能力。关于生态风险，虽然NIRI指数显示AF和NF风险大幅下降，但需注意该指数基于总金属浓度，可能高估了有机碳含量高、金属被强烈吸附的森林土壤的实际生物有效性风险。
针对源解析的不确定性，PMF模型虽能有效区分化学特征，但仍属统计推断。研究明确了四种来源：自然成土作用（背景）、交通排放（路侧Cu、Zn）、农用化学品遗留（RM内部Hg、Cd）以及生物循环过程（植被?凋落物?土壤）。
基于上述结果，研究提出了“恢复?源控制”管理框架。鉴于AF比NF具有更高的风险削减效率，建议在种植园内部推广多样化农林复合系统，停止含Hg农药并合理化施肥，以切断外源输入；同时利用行间空间加速土壤微环境改良（SOC和氮积累），通过有机质官能团络合固定遗留污染物。对于种植园边界，由于内部恢复无法抵消扩散的路侧排放，必须建立5–10?m宽的植被缓冲区，利用耐金属乡土植物物理拦截车辆颗粒物。
研究局限性在于空间换时间替代法难以完全排除初始空间异质性，且仅关注旱季表土，未涵盖深层淋溶风险和季节性变化。未来需加强长期监测、深层剖面采样及化学形态分析。
结论
该研究定量证明，单作向多样化森林的转变显著缓解了计算的潜在生态风险，其中人工恢复实现的NIRI降低约95%，优于自然演替的79%。结构方程模型验证了风险缓解的双重路径机制：主要由切断人为输入驱动，其次为恢复的土壤微环境缓冲能力提升。然而，源解析表明人为输入仍是污染负荷的主导因素。因此，仅靠恢复不足以完全消除边界区的持续外部输入风险，须实施空间差异化治理：在内部实施合理化施肥，沿路侧建立植被缓冲区。最终，这些发现证实了将生态修复与针对性源头削减相结合的必要性，以实现农业生产力与环境安全的可持续平衡。
