探索等效土壤质量替代方法以提升土壤碳监测、报告与核证的精度与稳健性

《Soil Use and Management》：Alternatives to Equivalent Soil Mass in Monitoring, Reporting and Verification of Changes in Soil Carbon

为提升土壤有机碳（SOC）储量估算的准确性与土壤碳信用体系的完整性，本文研究了土壤碳监测、报告与核证（MRV）中固定深度（FD）采样方法的局限性与等效土壤质量（ESM）方法的优化。研究人员对比了基于质量10%分位数（ESM10）、均值（ESMμ）以及新提出的等效容重（EBD10、EBDμ）四种SOC估算方法。基于对澳大利亚东部9个农场的密集采样数据，研究发现基于均值的EBDμ和ESMμ方法能提供更稳定、准确的SOC估算，相较于当前澳大利亚土壤碳方法中使用的ESM10，方差降低高达38%–86%。此研究为改进国际土壤碳监测方法提供了重要依据。

　　

土壤，这颗蓝色星球的皮肤，不仅是农作物生长的根基，更是应对气候变化的关键碳库。将大气中的二氧化碳“锁”进土壤——即土壤碳固存，被视为一项能够广泛应用且至关重要的二氧化碳移除技术。为此，全球各地建立了土壤碳信用体系，旨在激励农民采取有利于土壤碳积累的农地管理措施。然而，这些激励机制的公信力，完全依赖于一套能够精准监测、报告和核证土壤碳储量变化的科学方法。当前，最简单且广泛使用的估算方法是固定深度（FD）采样。但这种方法存在一个致命缺陷：它假设土壤的容重（单位体积土壤的干重）是恒定不变的。事实上，当耕作方式改变（例如从传统耕作转为免耕）时，土壤容重很可能会发生变化。这时，FD采样就会导致对土壤有机碳（SOC）储量的严重误判——容重增加会高估碳储量，容重减少则会低估。这无疑动摇了整个碳信用体系的科学基石。

为了纠正容重变化带来的误差，科学家们提出了等效土壤质量（ESM）的概念。与FD关注固定的采样体积不同，ESM的核心是确保每次采样都能获取“同等质量”的土壤颗粒。澳大利亚土壤碳方法（2021）以及国际上的Verra VM0042等方法学都已采纳ESM。然而，理想丰满，现实骨感。在真实的碳农业实践中，土壤采样常常面临障碍物，导致采集的土芯长度（即质量）差异巨大。这种采样长度的巨大变异性，给基于土芯总质量计算的ESM值，以及最终的SOC估算，带来了显著的不确定性。那么，有没有更稳定、更可靠的替代方案呢？

为了回答这个问题，并致力于改进澳大利亚土壤碳方法，研究人员开展了一项深入研究。他们的核心是比较四种不同的SOC估算策略：第一种是当前澳大利亚土壤碳方法中使用的、基于采样质量10%分位数的ESM₁₀；第二种是提议的、基于采样质量均值的ESM_μ；第三种和第四种则是本文新提出的、基于“等效容重”（EBD）概念的方法——分别基于土壤分层平均容重的10%分位数（EBD₁₀）和均值（EBD_μ）。研究人员假设，基于均值的方法（ESM_μ和EBD_μ）会比基于10%分位数的方法（ESM₁₀和EBD₁₀）更稳定，且EBD方法对“短样本”（采样长度不足）的敏感性低于ESM方法。

为验证这些假设，研究团队在2023年8月至2024年1月期间，对位于澳大利亚东部（昆士兰州、新南威尔士州、维多利亚州和南澳大利亚州）的九个农场进行了密集采样。这些农场覆盖了从低（年降雨量<500毫米）到高（>1000毫米）的不同降雨带，并包括种植和放牧两种管理方式。在每个农场，研究人员按照不规则均匀网格（包括一个密集采样的1公顷网格）采集了数十至上百个土芯，目标采样深度约为1米。每个土芯被分为0-30厘米和30-×厘米（×为实际采样深度）两层，随后在实验室中按照澳大利亚土壤碳方法的标准，测定了风干及烘干质量、砾石含量、含水量和SOC浓度等关键参数。

基于这些详实的实验数据，研究团队首先通过统计模拟，量化了不同样本量下四种估算方法（ESM₁₀, ESM_μ, EBD₁₀, EBD_μ）本身的变异性。结果表明，EBD估算值的变异性显著低于相应的ESM估算值，而基于均值的方法（ESM_μ和EBD_μ）又比基于10%分位数的方法（ESM₁₀和EBD₁₀）稳定得多。 例如，在样本量为9（澳大利亚土壤碳方法的最低要求）时，ESM₁₀的相对标准偏差在3.8%到21.0%之间波动，而EBD_μ的相对标准偏差仅为1.1%到5.4%。当样本量增加到50时，EBD_μ的变异性进一步降至0.5%到2.1%。研究还发现，小样本量会导致ESM₁₀和EBD₁₀被系统性高估，而这种偏差在基于均值的方法中并未出现。

接下来，研究人员分析了使用不同方法估算出的SOC储量。研究发现，基于均值的方法（ESM_μ和EBD_μ）估算出的SOC储量，普遍高于基于10%分位数的方法（ESM₁₀和EBD₁₀）。 与ESM₁₀相比，ESM_μ估算的SOC平均高出3.9%到23.4%，EBD_μ估算的则高出3.3%到25.9%。尽管EBD方法在估算“等效”参考值（EBD本身）时更稳定，但将其用于计算单个样本的SOC时，产生的SOC值方差与相应的ESM方法相似或略大。通过不确定性传播分析，研究揭示了ESM或EBD估算误差对最终SOC估算的影响是“衰减”的，即ESM/EBD变化x%，导致的SOC平均变化小于x%。同时，低估EBD对SOC估算的负面影响比低估ESM更严重。

最后，通过从拟合的土芯质量和长度分布中重复抽样，研究人员模拟了在真实抽样不确定性下，不同方法所得SOC平均值的分布。结果再次证实，基于均值的方法（尤其是EBD_μ）估算出的平均SOC值，其变异性最小。 在九个农场中的五个，EBD_μ的方差明显小于ESM_μ；在其余四个农场，两者差异很小。这强有力地支持了EBD_μ作为更优替代方案的潜力。

综上所述，本研究通过系统的理论分析和基于大量实地数据的模拟验证，得出明确结论：在土壤碳监测的背景下，与当前澳大利亚土壤碳方法所采用的ESM₁₀相比，基于采样质量均值的新方法ESM_μ，以及基于分层平均容重均值的新概念方法EBD_μ，能显著降低SOC储量估算的变异性，最高可减少38%至86%的方差，从而提供更稳定、更准确的估算结果。 特别是EBD_μ，其对短样本不敏感的特性使其在实际采样条件多变的场景下更具鲁棒性。

这项研究的意义深远。它不仅仅是对几个数学公式的优化，更是对土壤碳信用科学基础的一次重要加固。更稳定的估算方法意味着碳信用项目产生的结果更具可重复性和可信度，减少了因方法学不确定性导致的碳信用额波动风险，从而增强了投资者、政策制定者和公众对土壤碳市场整体的信心。本研究提出的EBD_μ和ESM_μ等替代方案，为澳大利亚乃至国际土壤碳监测、报告与核证（MRV）方法的演进提供了直接的、基于证据的改进路径，有助于推动土壤碳固存这一关键气候解决方案更加稳健、有效地实施。