DB1331T 090-2024 雄安新区农田土壤碳汇核算技术规程

ICS13.040.01CCSB04DB1331雄 安 新 区 地 方 标 准DB1331/T090-2024雄安新区农田土壤碳汇核算技术规程布 河北雄安新区管理委员会综合执法局发布河北雄安新区管理委员会公共服务局DB1331/T090—2024DB1331/T090—2024II目 录前 言 Ⅱ范围 1规范性引用文件 1术语与定义 1核算方法 3数据收集及质量管理 10附表A.1农田不同管理活动的相关库变化因子缺省值 11附表A.2化肥含氮量 11附表A.3农田施肥N2O排放因子 12附件B数据质量要求指南(规范性） 13DB1331/T090—2024DB1331/T090—2024ⅡⅡ前 言本标准由雄安新区农业农村标准化技术委员会归口管理并实施。本标准主要起草单位：雄安创新研究院、中国科学院生态环境研究中心。DB1331/T090—2024DB1331/T090—2024PAGEPAGE10雄安新区农田土壤碳汇核算技术规程范围本规程确定了河北雄安新区农田土壤碳汇核算方法，包括农田土壤碳储量和碳排放的测量、评估、监测和报告。壤碳汇核算可参照使用。下列文件对于本文件的应用是必不可少的。凡是注日期的引用文件，仅所注日期的版本适用于本文件。凡是不注日期的引用文件，其最新版本（包括所有的修改单）适用于本文件。GB/T21010土地利用现状分类GB/T33760基于项目的温室气体减排量评估技术规范通用要求GB/T36199土壤质量土壤采样程序设计指南GB/T36197GB/T32722DB42/T2204湖沼湿地温室气体通量监测技术规范NY/T1121.1土壤检测第1部分：土壤样品的采集、处理和贮存NY/T1121.4土壤检测第4部分：土壤容重的测定NY/T1121.66HJ658-滴定法HJ695土壤有机碳的测定燃烧氧化-非分散红外法HJ615土壤有机碳的测定重铬酸钾氧化-分光光度法术语与定义下列术语与定义适用于本文件。农田土壤碳汇 Farmlandsoilcarbonsink农田植物群落通过光合作用吸收大气中的二氧化碳，并将其固定在农田土壤中的碳储量扣除因管理措施而增加的非二氧化碳温室气体排放当量。农田土壤固碳 soilcarbonsequestrationinfarmland采用管理措施提高农田土壤的总碳含量，增加土壤有机碳库和无机碳库储量。农田土壤温室气体排放量 Greenhousegasemissionsfromfarmlandsoil由农田土壤边界内的人为活动或自然过程引起的、可测定的和可核查的温室气体排放当量。二氧化碳当量 carbondioxideequivalent（CO2e）在辐射强度上与某种温室气体质量相当的二氧化碳的量。土壤有机碳库 soilorganiccarbonpools30cm厚度耕层土壤中的有机碳储量。土壤有机质 soilorganicmatter土壤中形成的和外加入的所有动、植物残体不同阶段的各种分解产物和合成产物的总称。注：它包括高度腐解的腐殖物质、解剖结构尚可辨认的有机残体和各种微生物体。土壤有机碳 soilorganiccarbon土壤有机质中所含有的碳元素。土壤无机碳 soilinorganiccarbon3 一定深度内矿物质土和有机土（包括泥炭土）HCO-COCaCO3。3 土壤碳库变化因子 emissionfactorsofsoilcarbonpool量化单位面积管理、投入引起的农田土壤碳库增加的系数。土壤碳密度soilorganiccarbondensity单位面积固定深度（0-30cm）土壤有机碳碳含量，农田土壤耕种前的碳密度为土壤本底碳密度。土壤容重 soilbulkdensity田间状态下单位容积土体的干质量。秸秆还田 strawreturn把作物秸秆（小麦、玉米和水稻秸秆等）处理后归还土壤的一种方式。轮作 croprotation在同一块田地上，有顺序地轮换种植不同作物的种植方式。绿肥 greenmanure提供作物肥源或培肥土壤的作物。增施外源性碳 increasedexogenouscarbon增施外源的有机物质或有机物料，能够在达到稳定后提高土壤碳储量，包括增施符合标准的有机肥、生物有机肥、有机源土壤调理剂、有机源生物腐植酸肥料及外源秸秆等。基线年baselineyear核算期的起始年，一般为核算年20年前。核算年inventoryyear核算农田土壤碳汇的年份。非二氧化碳温室气体non-CO2greenhousegases除去二氧化碳外的其他温室气体，本标准中指农业生产过程中产生的甲烷和氧化亚氮两类温室气体。核算方法一般要求宜选取持续采用土壤有机质提升耕作管理措施3年及以上的农田进行核算。土壤碳汇核算方法分为实测法和估算法，前者依据检测的非二氧化碳温室气体排放量、土壤无机碳，以及土壤有机质或土壤有机碳含量计算，优先使用。后者依据已公开发表的不同土壤固碳及温室气体排放数据计算，如第三次全国土壤普查资料、中国土种数据库、《中国土种志》(1～6册)、《中国土壤》等，适用于非二氧化碳温室气体排放量、土壤无机碳和有机质含量数据无法获得的情况。土壤碳汇核算——实测法方法原理根据土壤有机碳、无机碳含量与土壤容重计算核算期内土壤耕层中总碳库的变化量，计算核算期内土壤碳库年度变化量。根据农田氧化亚氮、甲烷2种温室气体排放量，计算农田活动期间非二氧化碳温室气体排放当量，计算核算期内农田温室气体年度非二氧化碳排放当量。参数值测定土壤样品制备土壤样品采集时一般采用“S”形布点采样，采样点分布要保证均匀，从总体上控制整个采样区，避免耕作、施肥等农艺措施和田埂、沟边等特殊部位产生误差。样品采集应按GB/T36199-2018、GB/T36197-2018和NY/T1121.1―2006的规定执行，样品制备和保存应按GB/T32722-2016和NY/T1121.1―2006的规定执行。土壤容重的测定应按NY/T1121.4―2006的规定执行。土壤有机质含量的测定应按NY/T1121.6―2006的规定执行。土壤有机碳含量的测定应按HJ615―2011、HJ658―2013或HJ695―2014的规定执行。非二氧化碳温室气体排放量的测定监测、分析方法参照DB42/T2204-2024的第5、6节或其它相关方法的规定执行。农田土壤碳汇农田土壤碳汇由核算期内农田土壤固碳量减去土壤非二氧化碳温室气体排放当量计算得到。∆𝐶=−(1)式中：∆𝐶——农田土壤碳汇，单位为吨二氧化碳/年(tCO2e·a-1)；∆𝐶1——农田土壤总碳库年度变化量，单位为吨二氧化碳/年(tCO2e·a-1)；∆𝐶2——农田土壤非二氧化碳温室气体排放当量，单位为吨二氧化碳/年(tCO2e·a-1)；农田土壤固碳土壤真实碳密度土壤真实碳密度𝑆𝐶𝐷𝑎𝑐𝑡按式(2)或(3)计算，取平行分析结果的算术平均值为计算结果。𝑆𝐶𝐷𝑎𝑐𝑡，𝑖=×0.58+）𝑑××0.1 (2)或式中：

𝑆𝐶𝐷𝑎𝑐𝑡，𝑖=+）×𝑑××0.1 (3)𝑆𝐶𝐷𝑎𝑐𝑡，𝑖——第i年土壤真实碳密度（万tC·km-2）；（g·kg-1）；——第i年土壤有机碳含量（g·kg-1）；——第i年土壤无机碳含量（g·kg-1）；𝑑——土壤耕层深度，取值30，单位为厘米(cm)；𝐵𝐷𝑖——第i年土壤容重（g·cm-3）；0.58——土壤有机碳与土壤有机质的转化系数，无量纲；0.1——单位换算系数，无量纲。土壤固碳农田土壤总碳库年度变化量∆𝐶1按式(4)计算，取平行分析结果的算术平均值为计算结果。∆𝐶

=𝑆𝑇𝐶𝑇−𝑆𝑇𝐶0×44

(4)1 𝑇 12式中：

=𝑆𝐶𝐷𝑎𝑐𝑡，𝑖×𝐴×0.01 (5)∆𝐶1——农田土壤有机碳库年度变化量，单位为吨二氧化碳/年(tCO2e·a-1)；𝑆𝑇𝐶𝑇(tC)；𝑆𝑇𝐶0(tC)；𝑇——评估年与基线年间的时间跨度，一个单独核算期的年数，单位为年(a)；44——二氧化碳与碳的转换系数，无量纲；12𝐴——农田土地面积，单位为公顷(hm2)；0.01——单位换算系数，无量纲，1hm2=0.01km2；𝑆𝑇𝐶𝑖——第i年土壤总碳库，单位为吨碳(tC)；𝑆𝐶𝐷𝑎𝑐𝑡，𝑖——第i年土壤真实碳密度（万tC·km-2）。农田土壤非二氧化碳温室气体排放土壤非二氧化碳温室气体排放当量∆𝐶2按式(6)计算，取平行分析结果的算数平均值为计算结果。∆𝐶2

=𝑇𝑚

(6)=∑×𝑗

(7)式中：

𝑛=∑××(8)𝑖=1∆𝐶2——农田土壤非二氧化碳温室气体排放年度变化量，单位为吨二氧化碳/年(tCO2e·a-1)𝐴𝐷𝑃𝑙𝑎𝑛𝑡𝑖——第i种作物种植土地面积，单位为公顷(hm2)；𝐸𝐹𝑃𝑙𝑎𝑛𝑡——非二氧化碳温室气体排放二氧化碳当量，单位为(tCO2e)；𝑇——评估年与基线年间的时间跨度，一个单独核算期的年数，单位为年(a)；∆𝐸𝑗——第j种温室气体单位排放量，单位为(t)；𝐺𝑊𝑃𝑗——第j种温室气体全球增温潜势值，用来换算二氧化碳当量，单位为吨二氧化碳当量每吨温室气体(tCO2e/tEj)；𝐹𝑗,𝑖——第i种作物种植过程中第j种温室气体的平均排放通量，单位为(t·hm-2·a-1)；𝑇𝑖——评估年与基线年间的时间跨度内，第i种作物的种植时间，单位为(a)。土壤碳汇核算——估算法方法原理在核算期内，根据参考区域土壤碳库缺省值、农田不同管理活动相关库变化因子系数及每个时点(时间等于0和时间等于T)相应的面积，计算土壤总碳库𝑆𝑇𝐶𝑖，用核算期前后土壤总碳库变化量除以土壤总碳库达到稳定的时间20年，即为核算期内区域农田土壤固碳量。通过参考区域农田非二氧化碳温室气体排放量数据和农田不同管理活动变化因子，得到农田非二氧化碳温室气体排放碳当量，农田土壤固碳量和农田非二氧化碳温室气体排放碳当量之差，即为核算期内区域农田土壤碳汇的变化量。农田土壤碳汇农田土壤碳汇计算方法参见4.2.3，由核算期内农田土壤固碳量减去土壤非二氧化碳温室气体排放当量计算得到。土壤碳库变化因子农田土壤碳库变化因子计算方法数据获取土壤碳库变化因子数据可通过文献数据与长期定位监测数据结合的方式进行确定。数据选取的标准如下：1）监测持续时间长（>3年）；2）土样采自土层0-30cm；3）监测时段的起止年清楚；4）监测时段内土壤总碳的初始值和变化值明确；5）农田类型明确，为旱地、水浇地、果园中的一类。计算各数据源的土壤碳库年变化率𝑟=1×𝑆𝑇𝐶𝑡−𝑆𝑇𝐶𝑡0

(9)𝑡式中：𝑟——土壤碳库年变化率；𝑡——试验持续时间，单位为年；

𝑆𝑂𝐶𝑡0计算土壤碳库变化因子计算同一农田类型中，同一管理活动或退化程度的土壤碳库年变化率的均值和95%置信度区间，进而计算土壤碳库变化因子。𝐹=1+𝐷×(𝑟̅±𝜀) (10)式中：𝐹——土壤碳库变化因子；𝐷——评估年限20年；𝑟̅——土壤碳库的年变化率的均值；𝜀——土壤碳库的年变化率95%置信度区间上限与均值之差。农田类型土壤碳库变化因子按GB/T21010-2017将土地利用现状划分为旱地、水浇地和果园3种农田类型，依据4.3.3.1规定的步骤分别估算特定土地利用类型中农田系统或亚系统的碳库变化因子。农田管理土壤碳库变化因子按充分耕作、少耕和免耕3种耕作方式对农田管理进行划分，依据4.3.3.1规定的步骤分别估算每种农田类型下各管理措施的土壤碳库变化因子。有机质投入土壤碳库变化因子按低、中和高3种施肥方式对农田肥料投入进行划分，依据4.3.3.1规定的步骤分别估算每种施肥类型下的土壤碳库变化因子。农田土壤固碳土壤本底碳密度土壤本底碳密度通过土壤有机质或土壤有机碳含量、土壤无机碳含量，以及土壤容重计算得到。𝑆𝐶𝐷𝑟𝑒𝑓，𝑖=×0.58+）×𝑑××0.1 或𝑆𝐶𝐷𝑟𝑒𝑓，𝑖=+）×𝑑××0.1 式中：𝑆𝐶𝐷𝑟𝑒𝑓，𝑖——第i年被估算农田土壤本底碳密度（万tC·km-2）；（g·kg-1）；——第i年被估算农田土壤有机碳含量（g·kg-1）；——第i年被估算农田土壤无机碳含量（g·kg-1）；𝑑——土壤耕层深度，取值30，单位为厘米(cm)；𝐵𝐷𝑖——第i年被估算农田土壤容重（g·cm-3）；0.58——土壤有机碳与土壤有机质的转化系数，无量纲；0.1——单位换算系数，无量纲。土壤碳密度各农田类型不同管理和施肥模式下土壤碳密度由土壤本底碳密度和土壤碳库变化因子计算得到。式中：

=𝑆𝐶𝐷𝑟𝑒𝑓，𝑖×××(13)𝑆𝐶𝐷𝑖——第i年土壤碳密度（万tC·km-2）；𝑆𝐶𝐷𝑟𝑒𝑓，𝑖——第i年被估算农田土壤本底碳密度（万tC·km-2）；𝐹𝐿𝑈——特定土地利用中农田系统或亚系统的碳库变化因子，无量纲，亦可参见附表A.1相应碳库变化因子缺省值；——土地管理的碳库变化因子，无量纲，亦可参见附表A.1相应碳库变化因子缺省值；A.1相应碳库变化因子缺省值。土壤固碳根据土壤碳密度和被估算农田的土地面积核算土壤碳库，进而计算得到农田土壤固碳。∆𝐶

=𝑆𝑇𝐶𝑇−𝑆𝑇𝐶0×44

(14)1 𝑇 12式中：

=×𝐴×0.01 (15)∆𝐶1——农田土壤有机碳库年度变化量，单位为吨二氧化碳/年(tCO2e·a-1)；𝑆𝑇𝐶𝑇(tC)；𝑆𝑇𝐶0(tC)；𝑇——评估年与基线年间的时间跨度，土壤有机碳库达到稳定的时间，取20，单位为年(a)；44——二氧化碳与碳的转换系数，无量纲；12𝐴——被估算农田的土地面积，单位为公顷(hm2)；0.01——单位换算系数，无量纲，1hm2=0.01km2；𝑆𝑇𝐶𝑖——第i年土壤总碳库，单位为吨碳(tC)；𝑆𝐶𝐷𝑖——第i年土壤碳密度（万tC·km-2）。土壤固碳不确定性本规程的土壤固碳不确定性主要来自土壤碳库变化因子，采用不确定性传递法估算农田土壤固碳的不确定性范围。当某一估计值为n个估计值之和或差时，该估计值的不确定性采用下式计算。√∑(𝑈𝑠𝑛∙𝜇𝑠𝑛)2式中：

𝑈𝑐=

|∑

| (16)𝑈𝑐——n个估计值之和或差的不确定性（%）；𝑈𝑠𝑛——第n个相加减的估计值的不确定性（%）；𝜇𝑠𝑛——第n个相加减的估计值。当某一估计值为n个估计值之积时，该估计值的不确定性采用下式计算。=√∑𝑈𝑠𝑛2 农田非二氧化碳温室气体排放当量的计算方法数据获取农田土壤非二氧化碳温室气体排放变化因子数据可通过文献数据与长期定位监测数据结合的方式进行确定。数据选取的标准如下：1）监测持续时间长（>3年）；2）监测时段内农田土壤的非二氧化碳温室气体尽由作物种植和肥料使用产生；3）非二氧化碳温室气体类型明确，为甲烷、氧化亚氮的一类。农田土壤非二氧化碳温室气体排放=+式中：∆𝐶2——农田土壤非二氧化碳温室气体排放年度变化量，单位为吨二氧化碳/年(tCO2e·a-1)；∆𝐸𝑃𝑙𝑎𝑛𝑡——作物种植过程中非二氧化碳温室气体排放二氧化碳当量，单位为(tCO2e·a-1)；∆𝐸𝐹𝑒𝑟𝑡——作物种植过程中施肥产生的非二氧化碳温室气体排放二氧化碳当量，单位为(tCO2e·a-1)；作物生长非二氧化碳温室气体排放𝑚 𝑛=∑∑(𝐴𝐷𝑃𝑙𝑎𝑛𝑡𝑖×𝐹𝐸𝑃𝑖,𝑗)×

(19)式中：

𝑗𝑖=1∆𝐸𝑃𝑙𝑎𝑛𝑡——农田土壤作物种植非二氧化碳温室气体年度排放当量，单位为吨二氧化碳/年(tCO2e·a-1)；𝐴𝐷𝑃𝑙𝑎𝑛𝑡𝑖——第i种作物种植土地面积，单位为公顷(hm2)；𝐹𝐸𝑃𝑖,𝑗——第i种作作物种植过程中第j种温室气体排放因子，单位为(tCO2·hm2)；𝐺𝑊𝑃𝑗——第j种温室气体全球增温潜势值，用来换算二氧化碳当量，单位为吨二氧化碳当量每吨温室气体(tCO2e/tEj)。农田施肥N2O温室气体排放∆𝐸𝐹𝑒𝑟𝑡=++)×××(20)式中：∆𝐸𝐹𝑒𝑟𝑡——作物种植过程中施肥产生的N2O温室气体排放二氧化碳当量，单位为(tCO2e·a-1);𝐸𝐹1——农田施用化肥和有机肥料引起的N2O直接排放因子，单位为kgN2O-N/kgN，可参考附表A.3；𝐸𝐹2——基于氮挥发引起的N2O间接排放因子，单位为kgN2O-N/kgN-VOL，可参考附表A.3;𝐸𝐹3——基于氮淋溶发引起的N2O间接排放因子，单位为kgN2O-N/kgN-Lea，可参考附表A.3；𝑀𝑖——第i种肥料的使用量，单位为kg；𝑁𝑖——第i种肥料的含氮量，单位为%,可参考附表A.2；𝐺𝑊𝑃𝑁2𝑂——N2O的全球增温潜势值，取值298，单位为吨二氧化碳当量每吨氧化亚氮(tCO2e/tN2O)。数据收集及质量管理数据收集及质量要求数据收集及质量要求见表1。表 数据收集及质量要求方法数据数据收集数据质量要求实测法土壤有机质含量符合NY/T1121.6-2006等相关标准的规定土壤有机碳含量第三方检测数据、大样本普查数据、政府公布数据、科研数据符合HJ658-2013或HJ615-2011等相关标准的规定土壤容重符合NY/T1121.4-2006等相关标准的规定温室气体排放布数据、科研数据符合RB/T095-2022的规定被估算农田的面积有合法依据的统计数据合法合规估算法耕地类型现场勘查及调研、政府公布数据耕作方式现场勘查及调研、政府公布数据秸秆及有机肥投入现场勘查及调研、政府公布数据被估算农田的面积有合法依据的统计数据合法合规监测及数据质量管理监测计划及监测数据要求监测计划应按照GB/T33760--2017中5.10制定和执行。测量仪器/表精度应满足相关要求，每年定期校准，校准机构应具有测量仪器/表检定资质，校准相关要求应按照国家相关计量规程执行。技术应用方应确保监测计划有效实施，并记录、汇编和分析有关数据，对数据存档，保证测量管理体系符合质量和规范要求。数据质量管理数据质量管理要求按照GB/T33760--2017中5.11执行。数据质量应符合附件B的相关要求。附表A.1农田不同管理活动的相关库变化因子缺省值农田管理活动管理方式缺省值说明农田类型𝐹𝐿𝑈耕地旱地0.69指无灌溉设施，主要靠天然降水种植旱生农作物的耕地，包括没有灌溉设施，仅靠引洪淤灌的耕地。连续管理时间超过20年，主要种植一年生作物。水浇地0.69指有水源保证和灌溉设施，在一般年景能正常灌溉，种植旱生农作物(含蔬菜）的耕地。包括种植蔬菜的非工厂化的大棚用地。连续管理时间超过20年，主要种植一年生作物。园地果园1.00指种植果树的园地。长期生长多年生树种，如果树。耕作方式𝐹𝑀𝐺充分耕作1.00进行充分和/或一年中频繁耕作(如深翻等），对土壤产生大量干扰。在种植期，地表覆盖的残余物很少，通常低于30%。少耕1.08只进行一次和/或二次浅耕和不充分耕地，减少对土壤的干扰。在种植期，地表落叶残余物覆盖率通常高于30%。免耕1.15不经耕地直接进行播种，只在播种区最低限度干扰土壤。秸秆及有机肥投入𝐹𝐼低0.95(例如蔬菜），时不使用有机肥（秸秆还田、种植绿肥等），或不种植固氮作物。中1.17一年生作物秸秆还田、少耕和不施肥的管理模式。轮作中使用有机肥（秸秆还田、种植绿肥等）或种植固氮作物。高1.311.72增施外源性碳，包括有机肥、生物有机肥、有机源土壤调理剂、有机源生物腐植酸肥料、外源秸秆等。附表A.2化肥含氮量单位%化肥名称化学式含氮量碳酸氢铵NH4HCO330硝酸铵NH4NO335硫酸铵(NH)2SO421浓氨水NH3·H2O82硫硝铵硝酸铵和硫酸铵的混合物26磷酸二铵(NH4)2HPO418磷酸一铵(NH4)H2PO411尿素CH4N2O46.4硝酸钙Ca(NO3)215硝酸铵钙5Ca(N3)2NH4310H2O27附表A.3农田施肥N2O排放因子符号分类单位缺省值EF1农田施用化肥和有机肥料引起的N2O直接排放因子kgN2O-N/kgN0.01EF1FR稻田施用化肥和有机肥料引起的N2O直接排放因子kgN2O-N/kgN0.003EF2基于氮挥发引起的N2O间接排放因子kgN2O-N/kgN-VOL0.01EF3基于氮淋溶发引起的N2O间接排放因子kgN2O-N/kgN-Lea0.0075注：数据来源于IPCC《国家温室气体清单指南》附件B数据质量要求指南(规范性）活动水平数据监测数