导 读
土壤有机碳作为土壤肥力形成的基础，不但影响土壤质量、功能和粮食产量，而且在全球气候变化中扮演重要角色。在我国土壤资源同时面临保障粮食安全、发挥生态系统服务功能和应对气候变化等多重挑战的背景下，准确把握中国农田土壤固碳潜力及速率，对于实现土壤资源合理利用和农业可持续发展具有重要意义。文章首先介绍了对中国农田土壤有机碳变化速率和土壤固碳潜力的基本认识以及研究中面临的挑战，而后从基础研究、土壤信息平台、方法体系及研究成果与国家农业管理决策支撑方面提出了研究建议。
文/赵永存 徐胜祥 王美艳 史学正（中国科学院南京土壤研究所 土壤与农业可持续发展国家重点实验室）
来源：中国科学院院刊（2018年第2期）
土壤是陆地生态系统的核心，是人类赖以生存的重要自然资源。土壤有机碳（soil organic carbon，SOC）作为土壤肥力形成的基础，不但影响土壤质量和功能，而且在全球气候变化中扮演重要角色。SOC 是土壤肥力的决定性因素，其含量高低、质量好坏直接影响土壤肥力属性，即土壤有效持水量、保肥能力、养分利用效率、土壤微生物数量和活性，进而显著影响作物产量。同时，作为土壤碳库的重要组成部分，SOC 通过土壤微生物分解释放二氧化碳（CO2），而大气中的 CO2 则通过光合作用被固定到植物体，植物根系、凋落物及人为归还使得植物体中的部分碳再次归还到土壤中。因此，SOC 具有一定的大气 CO2 浓度调节功能。地球上 SOC 储量巨大且较为活跃，因而其微小变化就可能对大气 CO2 浓度产生重大影响，进而影响全球气候变化。
我国人多地少，耕地土壤质量总体不高，随着工业化和城市化进程的高速发展，人地、人粮矛盾日益突出，土壤资源正同时面临着保障粮食安全、发挥生态系统服务功能和应对气候变化等多重挑战。而农田作为受人为管理措施影响最为强烈的土壤利用方式，其 SOC 库最为活跃。同时，农田 SOC 库也是唯一可在较短时间尺度上通过合理利用而进行适度调节的碳库。因此，准确把握农田 SOC 变化速率及固碳潜力对于实现我国土壤资源高效利用及农业可持续发展战略，意义十分重大。
1 我国农田土壤固碳潜力及速率的基本认识
国家尺度农田SOC变化速率估算主要采用 Meta 分析、土壤调查数据差减和过程模型模拟3类方法。Meta分析采用已发表文献中的 SOC 数据，计算 SOC 变化速率；调查数据差减法通过两期土壤调查采样的SOC实测数据直接差减计算变化速率；过程模型模拟则采用SOC周转机理模型，在气候、土壤、农业管理措施等因子驱动下，实现SOC变化速率估算。Meta分析和两期调查数据差减，在SOC变化速率估算中没有整合SOC周转过程，主要用于估算“过去”的SOC变化速率，并不能获得SOC动态的逐年“演变”规律，也难以预测“未来”的固碳潜力。过程模型则整合了SOC 周转的机理过程，广泛用于“未来”不同假定情景下的固碳潜力预估。此外，基于长期试验数据建立SOC速率与管理措施的统计模型或获取SOC饱和水平并外推至全国尺度也被用于估算农田土壤固碳潜力。
我国农田SOC变化速率及固碳潜力估算中所采用的方法、农田面积、土层深度、土壤及管理措施数据来源及分辨率等均存在一定差异。但总体来看，目前对我国农田土壤固碳潜力及速率的认识基本形成了以下4点共识。
1.1 近30年来中国农田表层SOC库总体增加，发挥了碳汇功能
早期基于特定年份农业管理措施数据及 DNDC 过程模型②的估算认为，中国农田土壤丢碳，随着不同尺度SOC观测及农田管理措施数据收集的日趋完善，Meta分析、两期调查数据差减及过程模型等方法估算的结果均证明（表 1），近30年来中国农田表层 SOC 库总体增加，发挥了碳汇功能。其中，农田20cm 深度的土壤年固碳量在9.6—25.5Tg之间，30cm 深度在11—36.5Tg之间。单位耕地面积的固碳速率，20cm 深度每年74—184kgC/ha，30 cm深度则为每年85—281kg C/ha。
1.2 相对较低的初始SOC含量是中国农田土壤固碳的先决条件
初始SOC含量是影响SOC变化的重要因素之一。初始SOC含量对 SOC变化的影响通常具有负基线效应，即初始SOC含量高的地区土壤容易丢碳，而初始SOC含量低的地区则容易实现固碳。与西方发达国家相比，中国农田土壤除东北地区外，20世纪80年代的初始SOC库总体偏低，这不仅与我国农田土壤上千年的耕种和持续农业利用有关，也与盐碱土开垦耕地（如华北平原）、土壤侵蚀（如黄土高原）等因素有关。
基于第二次全国土壤普查数据估算的中国农田20世纪80年代初始 SOC库（20cm深度）基本处于26.6—32.5tC/ha之间，远低于美国农田的平均值43.7tC/ha和欧洲农田的平均值40.2tC/ha。中国农田相对较低的初始SOC含量，为利用化肥增加作物干物质产量从而提高土壤碳输入、促进农田土壤固碳提供了先决条件。
1.3 作物增产及秸秆还田大面积推广在中国农田土壤碳汇形成中发挥了关键作用
影响农田SOC变化的因素主要包括气候、土壤理化性质等自然因子，以及农业管理措施等人为因子。气候影响SOC的积累和分解过程，而土壤质地、pH值、通气性和土壤结构等理化性质则是影响SOC库容量的重要因素。施肥、耕作、秸秆管理等农业措施影响土壤碳输入和输出；同时，相对于温度、降水、土壤黏粒含量等因子而言，施肥、耕作、秸秆管理等农业措施促进农田土壤固碳是现实途径。因此，从这一点来看，中国近30年来农田土壤总体固碳必然与农业管理措施的变化历程密切相关。
从我国农田施肥和秸秆利用的阶段性特点来看，20世纪80年代早期化学肥料投入非常有限，而农家肥使用较为普遍；作物秸秆则主要用于燃烧和饲料。随后，农田化肥投入和粮食产量均快速增长，而农家肥投入逐年降低，但作物秸秆还田依然非常有限。90 年代后期，我国出现秸秆田间焚烧现象。1999年，原国家环保总局联合农业部、财政部等部门下发《秸秆禁烧和综合利用管理办法》的通知，禁止田间焚烧秸秆，引导秸秆综合利用。随后，秸秆还田在全国大规模推广。基于这些阶段性特点，从农田土壤碳输入变化角度来看，早期粮食产量的快速增长导致作物根系生物量增加从而显著增加了土壤碳输入；而随后秸秆还田的大规模推广则促进了作物地上生物量归还到土壤中。
中科院战略性先导科技专项项目课题“中国农田土壤固碳潜力与速率研究”覆盖全国农田的11470个点位的作物秸秆还田比例调查结果也表明（图 1），2000年之前我国水稻、小麦和玉米三大作物的秸秆还田比例均不足25%；但截至2010年，还田比例则分别达到了36.4%、39.7%和23.5%，从而在中国农田土壤碳汇形成中发挥了巨大作用。
1.4 中国农田SOC库当前水平依然总体偏低，未来还有进一步增加的潜力
“中国农田土壤固碳潜力与速率”课题基于2011年全国58个典型县4060个点位的SOC及容重等实测数据，估算的我国农田表层20 cm深度SOC库为32.34—33.47tC/ha，过程模型模拟的2010年全国农田平均SOC上限为36.0tC/ha，均还远低于欧盟同期的平均值 46.8tC/ha。
此外，全国95个水稻土长期施肥和秸秆还田试验、102个旱地长期施肥试验以及157个农田保护性耕作试验的Meta分析也表明，氮、磷、钾配施有机肥的固碳年速率可达0.37±0.04g/kg，而传统耕作+秸秆还田、免耕以及免耕+秸秆还田也均能显著提高我国农田的SOC含量，固碳年速率分别可达0.22g/kg、0.35g/kg和 0.52g/kg。
因此，在当前中国农田SOC依然总体相对较低的条件下，通过优化农业管理措施，我国农田土壤固碳潜力仍然有进一步提升的空间。目前的长期试验外推、过程模型预估等结果也表明，尽管中国农田土壤固碳潜力估算结果差异还较大，但可以确定的是SOC未来还有进一步增加的潜力。
2 面临的挑战
2.1 现有机理过程模型有待进一步改进
相对于两期调查数据差减及长期试验数据直接外推等方法的诸多局限性，过程模型整合了影响SOC动态的更多机理，能够反映气候、土壤、耕作管理等的空间变异性对 SOC 动态的影响，并且易于与土壤数据库、地理信息系统（GIS）、遥感等手段结合，因此已经成为大尺度 SOC 变化速率及固碳潜力研究的主流方法。
多库土壤碳周转模型（如 RothC、Century 等）在揭示不同尺度 SOC 库的变异性以及由于管理措施和土地利用变化而导致的土壤碳库变化方面均被证明是有效的，但在 SOC 对气候变化响应的适用性方面还有争议，主要问题是离散分库划分中的难分解/惰性碳库对温度变化的响应是否像假定的那样敏感性很小。同时，多库土壤碳周转模型中，微生物作物分解者的作用也被隐含在分解速率常数中，没有明确表达。
而目前新出现的微生物模型则直接把土壤碳周转和微生物生物量及生理机能耦合，对于探讨全球变暖和 SOC 间的反馈很有用。但是，该方法还存在震荡效应及对碳输入响应的不敏感性等问题。
因此，对农田土壤来说，在模型方面如何能同时有效的反映管理措施及气候变化对SOC的影响，还面临巨大挑战。
2.2 土壤基础数据分辨率低、更新慢、缺乏生物学指标
中国土壤资源丰富多样、空间分异明显。同时，在农业高强度利用、工业化及城镇化进程快速发展的背景下，土壤资源受人为活动影响强烈，土壤属性的变异性极强。土壤基本属性信息（SOC、容重、机械组成、pH值等）是固碳潜力估算模型的最基本驱动数据，同时SOC等直接观测数据也是模型参数校正及估算结果验证的关键。
然而，在初始土壤条件信息方面，目前可用的最精细的国家尺度土壤数据也仅为基于20世纪80年代第二次全国土壤普查成果建立的1:1000000土壤数据库。相较我国土壤的变异性而言，该数据的分辨率依然偏低。在模型参数校正及估算结果验证方面，目前可用的主要是不同地区零散的土壤调查及SOC数字制图数据。这些数据对我国典型农田区的空间覆盖不够完整，覆盖的时段则更少。此外，土壤基础信息的更新非常慢，或者说在国家尺度上基本没有更新。同时，表征土壤生物属性的微生物量、微生物多样性等空间分布基础信息极度缺乏，严重制约了国家尺度农田土壤固碳潜力研究。
2.3 缺乏不确定性的定量评估，制约了应用效果
由于土壤有机质形成机理和周转过程认识、观测数据误差、时空分辨率及可获取性等因素限制，农田 SOC变化速率及固碳潜力估算结果中必然包含一定的不确定性。这些不确定性会通过随后的分析及应用过程加以传递，从而对土壤资源可持续管理决策产生重大影响。
然而，目前土壤评估结果所提供的不确定性定量信息非常有限，综合考虑模型结构、内部参数及外部驱动数据不确定性的评估更少。中国国家尺度农田土壤固碳潜力及速率估算的不确定性空间分布信息则基本没有，在很大程度上制约了相关研究成果在国家农田土壤碳分区管理及指导土壤资源可持续利用中的实际应用效果。
3 行动建议
3.1 加强土壤有机质分解对温度变化响应的敏感性等基础研究
目前，关于土壤有机质分解对温度变化响应的敏感表性尚存很大争议，这也是未来气候变化条件下，土壤对大气CO2正负反馈的综合效应估算存在很大不确定性的主要原因之一。温度直接影响有机质分解速率，同时，环境因子也对土壤有机质分解产生影响，而环境因子本身也受气候变化影响。例如，气候和管理措施影响团聚体的形成，而团聚体对土壤有机质具有物理保护作用；温度影响有机质在矿物表面的吸附和解吸的化学过程；气候驱动的水平衡决定了土壤水膜的厚度，而可溶性有机碳和胞外酶通过土壤水膜进行扩散等。
因此，加强土壤有机质分解的温度敏感性等基础研究工作，对于改进现有SOC周转过程模型、准确估算SOC变化速率和固碳潜力具有重要意义。
3.2 强化数字土壤制图研究，丰富土壤基础数据，建立国家土壤信息服务平台
基于土壤空间变异先验信息建立土壤采样优化设计方案，整合遥感、近感、高光谱分析等技术，建立土壤基本属性快速获取方法，有效地降低土壤数据获取成本。加强土壤空间变异驱动因子定量识别研究，建立高强度人为作用影响下平原区农田土壤属性三维预测制图方法，为固碳潜力估算提供精细的土壤驱动数据。
建立土壤生物学指标的空间分布图以提供更全面的验证数据；建立历史土壤图的降尺度及空间和属性数据融合新方法，实现土壤遗产数据的高效利用，为固碳潜力研究提供高分辨率的初始土壤条件信息；建立包含多元土壤信息库、多服务模式的国家土壤信息服务平台，有效整合农田长期试验、土壤及农业管理措施区域调查等观测数据，实现土壤基础信息的有效共享，推进中国农田土壤固碳潜力研究。
3.3 建立系统化的研究方法体系，强化研究成果对土壤资源管理决策的支持作用
以土壤有机质形成机理和周转过程研究为基础，改进现有SOC过程模型，以数字土壤制图获取的 SOC、微生物量碳等多时段空间分布数据为依托，结合数据同化方法，实现模型内部参数校正以及固碳潜力估算结果的不确定性定量评估。整合模型结构改进、参数敏感性分析、参数校正、模型结构、参数及外部驱动数据的不确定性，建立系统化的中国农田土壤固碳速率及潜力研究方法体系。加强以SOC变化与相关政策关联性为切入点的研究工作，强化中国农田土壤固碳与相关政策之间的反馈研究，为我国土壤资源可持续利用管理提供强大的战略决策支撑。
