温室气体自愿减排项目方法学造林碳汇(CCER—14—001—V01)1引言造林碳汇项目可通过增加森林面积和森林生态系统碳储量实现二氧化碳清除,是减缓气候变化的重要途径。本方法学属于林业和其他碳汇类型领域方法学。符合条件的造林碳汇项目可按照本文件要求,设计和审定温室气体自愿减排项目,以及核算和核查温室气体自愿减排项目的减排量。2适用条件本文件适用于乔木、竹子和灌木造林,包括防护林、特种用途林、用材林等造林,不包括经济林造林、非林地上的通道绿化、城镇村及工矿用地绿化,使用本文件的造林碳汇项目必须满足以下条件:a)项目土地在项目开始前至少三年为不符合森林定义的规划造林地;b)项目土地权属清晰,具有不动产权属证书、土地承包或流转合同;或具有经有批准权的人民政府或主管部门批准核发的土地证、林权证;c)项目单个地块土地连续面积不小于400m2。对于2019年(含)之前开始的项目,土地连续面积不小于667m2;d)项目土地不属于湿地;e)项目不移除原有散生乔木和竹子,原有灌木和胸径小于2cm的竹子的移除比例总计不超过项目边界内地表面积的20%;f)除项目开始时的整地和造林外,在计入期内不对土壤进行重复扰动;g)除对病(虫)原疫木进行必要的火烧外,项目不允许其它人为火烧活动;h)项目不会引起项目边界内农业活动(如种植、放牧等)的转移,即不会发生泄漏;i)项目应符合法律、法规要求,符合行业发展政策。3规范性引用文件本文件引用了下列文件或其中的条款。凡是注日期的引用文件,仅注日期的版本适用于本文件。凡是未注日期的引用文件,其有效版本(包括所有的修改单)适用于本文件。GB/T26424森林资源规划设计调查技术规程GB/T38590森林资源连续清查技术规程CH/T8024机载激光雷达数据获取技术规范LY/T1812林地分类BD420073全球卫星导航系统(GNSS)定向设备性能要求及测试方法14术语和定义GB/T15776、GB/T18337.1、GB/T26423、LY/T1812、LY/T2647、LY/T2736和TD/T1055界定的以及下列术语和定义适用于本文件。4.1造林forestation在不符合森林定义的规划造林地上,通过人工措施营建或恢复森林的过程。[来源:GB/T15776-2023,3.1,有修改]4.2规划造林地planningafforestationland依据全国国土调查及其最新年度国土变更调查成果数据,综合考虑降水、积温、地貌、海拔、坡度、坡向、地表基质、土壤类型等自然条件,在各级国土空间规划中明确的,可用于造林绿化的用地空间。[来源:《造林绿化落地上图技术规范(试行)》(办生字〔2021〕87号),3.1,有修改]4.3森林forest包括乔木林、竹林和国家特别规定的灌木林。其中,国家特别规定的灌木林按照国家相关行业主管部门的规定。[来源:中华人民共和国森林法,第八十三条]4.4乔木林arborealforests由乔木(含因人工栽培而矮化的)树种组成,郁闭度≥0.2的片林或林带。其中,乔木林带行数应在2行以上且行距≤4m或林冠冠幅水平投影宽度在10m以上。[来源:GB/T26423-2010,6.68,有修改]4.5竹林forestsofbamboos附着有胸径2cm以上的竹类植物,郁闭度≥0.2的片林或林带。[来源:GB/T26423-2010,6.69;LY/T1812-2021,表1,有修改]4.6通道绿化passagewaygreening在铁路、公路沿线及两侧,以及河道两岸、渠道两侧,以保护铁路、公路、河岸、渠岸等生态安全、改善生态景观为主要目标的造林绿化。[来源:LY/T2647-2016,3.1]4.72城镇村及工矿用地towns,villages,industrialandminingarea城乡居民点、独立居民点以及居民点以外的工矿、国防、名胜古迹等企事业单位用地,包括城市、建制镇、村庄、盐田及采矿用地、特殊用地。[来源:TD/T1055-2019,表A.3,有修改]4.8生态公益林non-commercialforest为维护和改善生态环境,保持生态平衡,保护生物多样性等满足人类社会的生态、社会需求和可持续发展为主体功能,主要提供公益性、社会性产品或服务的森林、林木、林地。生态公益林按事权等级划分为国家级公益林和地方级公益林。国家级公益林区划界定执行《国家级公益林区划界定办法》的相关规定。[来源:GB/T18337.1-2001,2.1;《国家级公益林区划界定办法》(林资发〔2017〕34号),有修改]4.9经济林non-woodproductforest以生产果品,食用油料、调料、饮料,工业原料,药材和生物质能源为主要目的的林种。包括:以生产各种干、鲜果品为主要目的的果品林(如香榧、枣、苹果、梨、桃等),以生产食用油料、饮料、调料、香料等为主要目的的食用原料林(如咖啡、茶树、椰子等),以生产工业油料、树脂、木栓、单宁等非木质林产化工原料为主要目的的林产工业原料林(如油茶、小桐子、棕榈等),以生产药材、药用原料为主要目的的药用林(如杜仲、厚朴、肉桂等),以及以生产其他林副(特)产品为主要目的的其他经济林。[来源:LY/T2736-2016,2.1,有修改]4.10碳库carbonpools生态系统中碳储存的形式或场所,包括地上生物质、地下生物质、枯落物、枯死木、土壤有机碳和木(竹)产品。4.11地上生物质abovegroundbiomass土壤层以上所有活体植物的生物质,包括茎干、桩、枝、皮、叶、花、果和繁殖体等。4.12地下生物质belowgroundbiomass土壤层以下所有植物活根的生物质,通常不包括难以从土壤有机成分或枯落物中区分出来的直径≤2mm的细根。4.13生物量biomass地上生物质和地下生物质总的干物质质量。4.143枯落物litter枯落物是土壤层以上,直径≤5cm、处于不同分解状态的所有死有机质,包括凋落物、腐殖质,以及难以从地下生物质区分出来的细根。4.15枯死木deadwood枯落物以外的所有死有机质,包括枯立木、枯倒木以及直径>5cm的枯枝、死根和树桩。4.16土壤有机碳soilorganiccarbon一定深度内(通常为30cm)矿质土和有机土(包括泥炭土)中的有机碳,包括难以从地下生物质中区分出来的直径≤2mm的细根。4.17木(竹)产品harvestedwoodorbambooproducts由项目产生的、从项目边界内移出的木材(或竹材)加工而成,在项目计入期结束后仍然在用或进入到垃圾填埋的木制(或竹制)产品。4.18湿地wetlands全年(或一年中大部分时间,如泥炭土)被水淹没或土壤水分处于饱和状态的土地,且不属于森林、农田、草地和居住用地的范畴。5项目边界、计入期、碳库和温室气体排放源5.1项目边界造林碳汇项目区域可包括若干个不连续的地块,每个地块应有特定的地理边界。项目边界内不包括宽度大于3m的道路、沟渠、坑塘、河流等不符合适用条件的土地。项目边界可采用下述方法之一确定:a)利用北斗卫星导航系统(BDS)、全球定位系统(GPS)等卫星定位系统,直接测定项目地块边界的拐点坐标,单点定位误差不超过±5m;b)利用空间分辨率不低于5m的地理空间数据(如卫星遥感影像、航拍影像等)、林草资源“一张图”、造林作业设计等,在地理信息系统(GIS)辅助下直接读取项目地块的边界坐标。5.2项目计入期5.2.1项目计入期为可申请项目减排量登记的时间期限,从项目业主申请登记的项目减排量的产生时间开始,最短时间不低于20年,最长不超过40年。项目计入期须在项目寿命期限范围之内。5.2.2项目寿命期限应在项目业主对项目边界内土地的所有权(或使用权)或项目边界内林木的所有权(或经营权)的有效期限之内。项目寿命期限的开始时间即项目边界内首次实施整地、播种或植苗的项目开工日期。45.3碳库和温室气体排放源的选择项目边界内选择或不选择的碳库如表1所示。表1碳库的选择情景碳库是否选择理由基准线情地上生物质否地下生物质在计算项目清除量时扣除景否枯死木否根据适用条件,该碳库的清除量所占比例小,忽略不计项目枯落物否情景土壤有机碳是根据适用条件,土地处于稳定或退化状态,忽略不计木(竹)产品根据适用条件,该碳库的清除量所占比例小,忽略不计地上生物质是或否地下生物质主要碳库是或否枯死木是相比基准线情景,造林项目通常会增加枯死木碳储量;如果否项目存在移除枯死木的情形,基于保守性原则不选择该碳库枯落物相比基准线情景,造林项目通常会增加枯落物碳储量;如果项目存在移除枯落物的情形,基于保守性原则不选择该碳库土壤有机碳木(竹)产品造林项目会引起土壤有机碳储量发生变化按照保守性原则,忽略不计项目边界内选择或不选择的温室气体排放源与种类如表2所示。表2温室气体排放源的选择情景温室气体温室气体是否选择理由基准线排放源种类否否按照保守性原则,忽略不计情景火灾或人为火烧CO2、CH4是和N2O否生物质燃烧导致的CO2排放已在碳储量项目使用车辆、机械设备变化中考虑情景等过程中化石燃料燃在项目设计阶段计为0;如果项目计入期烧产生的排放内发生森林火灾或人为的火烧活动,则使用石灰、污泥施肥过程中产生的排放必须选择该排放源火灾或人为火烧CO2相对于基准线情景,排放量的变化量不显著,忽略不计使用车辆、机械设备CH4和等过程中化石燃料燃N2O烧产生的排放CO2、CH4使用石灰、污泥施肥和N2O过程中产生的排放6项目减排量核算方法6.1基准线情景识别本文件规定的造林碳汇项目基准线情景为:维持造林项目开始前的土地利用与管理方式。56.2额外性论证6.2.1免予论证以保护和改善人类生存环境、维持生态平衡等为主要目的的公益性造林项目,在计入期内除减排量收益外难以获得其他经济收入,造林和后期管护等活动成本高,不具备财务吸引力。符合下列条件之一的造林项目,其额外性免予论证:a)在年均降水量≤400mm的地区开展的造林项目;注:年均降水量≤400mm的地区可参考《国家林业局关于颁发<“国家特别规定的灌木林地”的规定>(试行)的通知》(林资发〔2004〕14号)。b)在国家重点生态功能区开展的造林项目;注:国家重点生态功能区可参考《国务院关于印发全国主体功能区规划的通知》(国发〔2010〕46号)、《国务院关于同意新增部分县(市、区、旗)纳入国家重点生态功能区的批复》(国函〔2016〕161号)。c)属于生态公益林的造林项目。6.2.2一般论证其他造林项目按照《温室气体自愿减排项目设计与实施指南》中“温室气体自愿减排项目额外性论证工具”对项目额外性进行一般论证。6.3项目碳层划分6.3.1为提高碳储量变化量计算的精度,并在一定精度要求下精简监测样地数量,应按照不同的分层因子将项目边界内的地块划分为不同的层次,包括项目设计阶段的碳层划分和项目实施阶段的碳层划分。6.3.2项目设计阶段划分的碳层用于预估碳储量变化量,综合考虑项目边界内土地在造林前的立地条件(如土壤类型、坡度坡向、海拔等),以及拟实施的项目造林时间、造林树种、造林密度等因素划分项目碳层,将无显著差别的造林地块划分为同一碳层。6.3.3项目实施阶段划分的碳层用于计算碳储量变化量,主要基于项目设计阶段碳层的划分,结合造林活动的实际情况进行调整确定。若存在自然因素(如立地条件、火灾、病虫害等)或人为干扰(如火烧、采伐等)导致原有碳层的异质性增加,或土地利用发生变化,须对项目碳层进行调整。6.4基准线清除量计算基准线情景下原有植被的生物质碳储量变化量在项目清除量的计算中给予考虑。根据表1和表2,项目开始后第t年的基准线清除量计为0,即:Δ𝐶BSL,𝑡=0(1)式中:Δ𝐶BSL,𝑡——项目第t年的基准线清除量,单位为吨二氧化碳当量每年(tt——CO2e·a-1);自项目开始以来的年数,t=1,2,3……,无量纲。6.5项目清除量计算项目开始后第t年的项目清除量按照公式(2)计算:6Δ𝐶PROJ,𝑡=Δ𝐶BiomassPROJ,𝑡+Δ𝐶DOMPROJ,𝑡+∆𝑆𝑂𝐶PROJ,𝑡−𝐺𝐻𝐺PROJ,𝑡(2)−∆𝐶BiomassPE,𝑡式中:——项目第t年的项目清除量,单位为吨二氧化碳当量每年(t∆𝐶PROJ,𝑡——CO2e·a-1);∆𝐶BiomassPROJ,𝑡——项目第t年的生物质碳储量变化量,单位为吨二氧化碳当量∆𝐶DOMPROJ,𝑡——每年(tCO2e·a-1),采用本文件附录A计算;∆𝑆𝑂𝐶PROJ,𝑡——项目第t年的死有机质碳储量变化量,单位为吨二氧化碳当𝐺𝐻𝐺PROJ,𝑡——量每年(tCO2e·a-1),采用本文件附录B计算;∆𝐶BiomassPE,𝑡项目第t年的土壤有机碳储量变化量,单位为吨二氧化碳当——量每年(tCO2e·a-1),采用本文件附录C计算;t项目第t年因火烧引起的温室气体排放量,单位为吨二氧化碳当量每年(tCO2e·a-1),采用本文件附录D计算;项目第t年原有植被(乔木、竹子和灌木)的生物质碳储量变化量,单位为吨二氧化碳当量每年(tCO2e·a-1),采用本文件附录A计算;自项目开始以来的年数,t=1,2,3……,无量纲。在项目设计阶段,火烧引起的温室气体排放通常无法预料,因此项目情景下火烧引起的温室气体排放量计为0。在项目实施阶段,通过监测项目边界内实际火烧发生情况,计算项目温室气体排放量。6.6项目泄漏计算根据本文件适用条件,项目不考虑泄漏。6.7项目减排量核算项目开始后第t年的项目减排量按照公式(3)核算:𝐶𝐷𝑅𝑡=(Δ𝐶PROJ,𝑡−Δ𝐶BSL,𝑡−𝐿𝐾𝑡)×(1−𝐾RISK)(3)式中:𝐶𝐷𝑅𝑡——项目第t年的项目减排量,单位为吨二氧化碳当量每年(tCO2e·a-1);Δ𝐶PROJ,𝑡——项目第t年的项目清除量,单位为吨二氧化碳当量每年(tCO2e·a-1);Δ𝐶BSL,𝑡——项目第t年的基准线清除量,单位为吨二氧化碳当量每年(tCO2e·a-1);𝐿𝐾𝑡——项目第t年的泄漏量,单位为吨二氧化碳当量每年(tCO2e·a-1);根据适用条件,𝐿𝐾𝑡=0;𝐾RISK——项目的非持久性风险扣减率,单位为百分比(%);t——自项目开始以来的年数,t=1,2,3……,无量纲。7监测方法7.1项目设计阶段确定的参数和数据项目设计阶段需确定的参数和数据的技术内容和确定方法见表3—表26。表3𝑲RISK的技术内容和确定方法7数据/参数名称𝐾RISK应用的公式编号公式(3)造林碳汇项目可能会由于自然因素(如火灾、病虫害、雨雪冰冻、风灾等)或数据描述人为干扰(如非法采伐和破坏等)原因导致项目清除的温室气体重新释放到大气中,即非持久性风险。在核算减排量时须按照项目非持久性风险扣减率,扣数据单位除一定比例的项目减排量。非持久性风险扣减率采用历史火灾、病虫害等灾害导致的森林蓄积量或森林面积的损失比例计算确定数据来源%本表默认值,根据《中国林业统计年鉴》统计全国及各省(区、市)因火灾引数值起的蓄积损失量占当年森林蓄积增长量的比例,以及病虫害重度危害面积占森数据用途林面积的比例,1999-2018年全国年均因灾损失率约4.98%,出于保守性原则取值10%数据/参数名称10%应用的公式编号用于计算项目减排量的非持久性风险数据描述表4𝑨𝒊,𝒋,𝒕的技术内容和确定方法数据单位数据来源𝐴𝑖,𝑗,𝑡公式(A.2)、公式(A.3)、公式(A.4)、公式(B.2)数值第t年时,第i项目碳层树种j的森林面积数据用途hm2项目设计文件及审定确认的项目碳层各树种的森林面积数据/参数名称应用的公式编号/用于预估项目清除量数据描述数据单位表5𝑪𝑭的技术内容和确定方法数据来源𝐶𝐹公式(A.2)、公式(A.3)、公式(A.4)、公式(B.2)、公式(D.3)数值生物量含碳率数据用途tC·(td.m.)-1数据/参数名称项目业主须按照如下优先顺序选择:应用的公式编号a)地方标准;b)国家或行业标准中适用于项目区的数据;数据描述c)本文件及附录中推荐的缺省值;数据单位d)项目区当地、相邻地区或相似生态条件下的调查统计数据。须基于5篇以数据来源上国内外核心期刊发表的研究结果或总数不少于30个样本的调查数据的整合分析,且经过同行专家评议表A.10用于将生物量转换为生物质碳储量表6𝑹𝑺𝑹的技术内容和确定方法𝑅𝑆𝑅公式(A.8)、公式(A.13)、公式(A.14)、公式(A.16)、公式(A.22)、公式(A.24)乔木、竹子或灌木的地下生物量占地上生物量的比例无量纲项目业主须按照如下优先顺序选择:a)地方标准;b)国家或行业标准中适用于项目区的数据;c)本文件及附录中推荐的缺省值;d)项目区当地、相邻地区或相似生态条件下的调查统计数据。须基于5篇以上国内外核心期刊发表的研究结果或总数不少于30个样本的调查数据的8数值整合分析,且经过同行专家评议数据用途表A.8、表A.9、表A.12、表A.15用于利用地上生物量计算地下生物量数据/参数名称应用的公式编号表7𝒇(𝑫𝑩𝑯,𝑯)的技术内容和确定方法数据描述𝑓(𝐷𝐵𝐻,𝐻)数据单位公式(A.10)单株乔木全株(或地上)生物量与胸径和(或)树高的相关方程数据来源kgd.m.·stem-1数值项目业主须按照如下优先顺序选择:数据用途a)地方标准;b)国家或行业标准中适用于项目区的数据;数据/参数名称c)本文件及附录中推荐的缺省值;应用的公式编号d)项目区当地、相邻地区或相似生态条件下的调查统计数据。须基于5篇以数据描述上国内外核心期刊发表的研究结果或总数不少于30个样本的调查数据的数据单位整合分析,且经过同行专家评议表A.1、表A.2、表A.3数据来源用于利用胸径和树高计算乔木林单株生物量数值表8𝒇(𝑽AF,𝒕)的技术内容和确定方法数据用途𝑓(𝑉AF,𝑡)数据/参数名称公式(A.6)应用的公式编号乔木林单位面积全林(或地上)生物量与蓄积量的相关方程数据描述td.m.·hm-2数据单位项目业主须按照如下优先顺序选择:a)地方标准;数据来源b)国家或行业标准中适用于项目区的数据;c)本文件及附录中推荐的缺省值;数值d)项目区当地、相邻地区或相似生态条件下的调查统计数据。须基于5篇以数据用途上国内外核心期刊发表的研究结果或总数不少于30个样本的调查数据的整合分析,且经过同行专家评议表A.5用于利用单位面积蓄积量计算乔木林全林(或地上)生物量表9𝒇(𝑨𝒈𝒆AF,𝒕)的技术内容和确定方法𝑓(𝐴𝑔𝑒AF,𝑡)公式(A.7)乔木林单位面积蓄积量与林龄的相关方程m3·hm-2项目业主须按照如下优先顺序选择:a)地方标准;b)国家或行业标准中适用于项目区的数据;c)本文件及附录中推荐的缺省值;d)项目区当地、相邻地区或相似生态条件下的调查统计数据。须基于5篇以上国内外核心期刊发表的研究结果或总数不少于30个样本的调查数据的整合分析,且经过同行专家评议表A.11用于利用单位面积蓄积量计算乔木林全林生物量表10𝑩𝑪𝑬𝑭的技术内容和确定方法9数据/参数名称𝐵𝐶𝐸𝐹应用的公式编号公式(A.13)、公式(A.14)乔木林生物量转换与扩展因子,即全林(或地上)生物量与蓄积量的比值数据描述数据单位td.m.·m-3项目业主须按照如下优先顺序选择:数据来源a)地方标准;b)国家或行业标准中适用于项目区的数据;数值c)本文件及附录中推荐的缺省值;数据用途d)项目区当地、相邻地区或相似生态条件下的调查统计数据。须基于5篇以数据/参数名称上国内外核心期刊发表的研究结果或总数不少于30个样本的调查数据的应用的公式编号整合分析,且经过同行专家评议表A.5、表A.6数据描述用于将乔木林蓄积量转化为全林生物量或地上生物量数据单位表11𝑩𝑬𝑭的技术内容和确定方法数据来源𝐵𝐸𝐹数值公式(A.17)数据用途乔木林生物量扩展因子,即地上生物量与树干生物量的比值无量纲数据/参数名称项目业主须按照如下优先顺序选择:应用的公式编号a)地方标准;b)国家或行业标准中适用于项目区的数据;数据描述c)本文件及附录中推荐的缺省值;数据单位d)项目区当地、相邻地区或相似生态条件下的调查统计数据。须基于5篇以数据来源上国内外核心期刊发表的研究结果或总数不少于30个样本的调查数据的整合分析,且经过同行专家评议数值表A.7数据用途用于将乔木林树干生物量转化为地上生物量数据/参数名称表12𝑺𝑽𝑫的技术内容和确定方法应用的公式编号𝑆𝑉𝐷数据描述公式(A.17)数据单位乔木树种的基本木材密度,即单位体积木材的干物质重量td.m.·m-3项目业主须按照如下优先顺序选择:a)地方标准;b)国家或行业标准中适用于项目区的数据;c)本文件及附录中推荐的缺省值;d)项目区当地、相邻地区或相似生态条件下的调查统计数据。须基于5篇以上国内外核心期刊发表的研究结果或总数不少于30个样本的调查数据的整合分析,且经过同行专家评议表A.10用于将林木蓄积量转换为树干生物量表13𝒇𝒗(𝑫𝑩𝑯,𝑯)的技术内容和确定方法𝑓𝑣(𝐷𝐵𝐻,𝐻)公式(A.15)乔木单株材积与树高和(或)胸径的相关方程m3·stem-110数据来源项目业主须按照如下优先顺序选择:a)地方标准;数值b)国家或行业标准中适用于项目区的数据;数据用途c)本文件及附录中推荐的缺省值;d)项目区当地、相邻地区或相似生态条件下的调查统计数据。须基于5篇以数据/参数名称应用的公式编号上国内外核心期刊发表的研究结果或总数不少于30个样本的调查数据的整合分析,且经过同行专家评议数据描述表A.4数据单位用于计算乔木林单位面积全林生物量数据来源表14𝑨𝑮𝑩BF,𝑻𝒃的技术内容和确定方法数值𝐴𝐺𝐵BF,𝑇𝑏数据用途公式(A.20)、公式(A.22)竹林成熟稳定后的单位面积地上生物量数据/参数名称应用的公式编号td.m.·hm-2项目业主须按照如下优先顺序选择:数据描述a)地方标准;数据单位b)国家或行业标准中适用于项目区的数据;c)本文件及附录中推荐的缺省值;数据来源d)项目区当地、相邻地区或相似生态条件下的调查统计数据。须基于5篇以数值上国内外核心期刊发表的研究结果或总数不少于30个样本的调查数据的数据用途整合分析,且经过同行专家评议表A.12用于计算竹林单位面积地上生物量表15𝒇AGB,BF(𝑫𝑩𝑯,𝑯)的技术内容和确定方法𝑓AGB,BF(𝐷𝐵𝐻,𝐻)公式(A.21)竹子单株生物量与胸径和(或)竹高的相关方程kgd.m.·stem-1项目业主须按照如下优先顺序选择:a)地方标准;b)国家或行业标准中适用于项目区的数据;c)本文件及附录中推荐的缺省值;d)项目区当地、相邻地区或相似生态条件下的调查统计数据。须基于5篇以上国内外核心期刊发表的研究结果或总数不少于30个样本的调查数据的整合分析,且经过同行专家评议表A.13、表A.14用于计算竹林单位面积地上生物量数据/参数名称表16𝑨𝑮𝑩𝐒𝐅的技术内容和确定方法应用的公式编号𝐴𝐺𝐵SF数据描述公式(A.24)数据单位灌木林成熟稳定时的单位面积地上生物量数据来源td.m.·hm-2项目业主须按照如下优先顺序选择:a)地方标准;b)国家或行业标准中适用于项目区的数据;c)本文件及附录中推荐的缺省值;11数值d)项目区当地、相邻地区或相似生态条件下的调查统计数据。须基于5篇以数据用途上国内外核心期刊发表的研究结果或总数不少于30个样本的调查数据的整合分析,且经过同行专家评议数据/参数名称应用的公式编号表A.15用于计算灌木林地上生物量数据描述数据单位表17𝒇SF(𝒙1,𝒙2,𝒙3……)的技术内容和确定方法数据来源𝑓SF(𝑥1,𝑥2,𝑥3……)公式(A.23)数值灌木单株生物量与灌木测树因子(如基径、灌高、冠幅、灌径等)的相关方程数据用途kgd.m.·stem-1项目业主须按照如下优先顺序选择:数据/参数名称a)地方标准;应用的公式编号b)国家或行业标准中适用于项目区的数据;c)本文件及附录中推荐的缺省值;数据描述d)项目区当地、相邻地区或相似生态条件下的调查统计数据。须基于5篇以数据单位数据来源上国内外核心期刊发表的研究结果或总数不少于30个样本的调查数据的整合分析,且经过同行专家评议数值表A.16、表A.17数据用途用于计算灌木林地上生物量、地下生物量和全林生物量数据/参数名称表18𝑪𝑫PE,𝒊,𝒕𝟎的技术内容和确定方法应用的公式编号𝐶𝐷PE,𝑖,𝑡0数据描述公式(A.25)数据单位项目开始前,各项目碳层样地内原有散生木(竹)的平均冠层盖度,即树冠投影面积与林地面积的比值数据来源无量纲,以小数计数值数据用途根据项目开始前的森林资源调查数据确定;或利用分辨率不超过2m的高清卫星影像数据,借助GIS等工具进行分析确定数据/参数名称实际测量用于扣减项目碳层原有散生木(竹)继续生长产生的清除量表19𝑫𝑭𝐋𝐈的技术内容和确定方法𝐷𝐹LI公式(B.3)枯落物生物量与森林地上生物量的比例无量纲项目业主须按照如下优先顺序选择:a)地方标准;b)国家或行业标准中适用于项目区的数据;c)本文件及附录中推荐的缺省值;d)项目区当地、相邻地区或相似生态条件下的调查统计数据。须基于5篇以上国内外核心期刊发表的研究结果或总数不少于30个样本的调查数据的整合分析,且经过同行专家评议表B.1用于计算枯落物生物量表20𝑫𝑭DW的技术内容和确定方法𝐷𝐹DW12应用的公式编号公式(B.4)数据描述枯死木生物量与森林地上生物量的比例数据单位无量纲项目业主须按照如下优先顺序选择:数据来源a)地方标准;b)国家或行业标准中适用于项目区的数据;数值c)本文件及附录中推荐的缺省值;数据用途d)项目区当地、相邻地区或相似生态条件下的调查统计数据。须基于5篇以数据/参数名称上国内外核心期刊发表的研究结果或总数不少于30个样本的调查数据的应用的公式编号整合分析,且经过同行专家评议表B.2数据描述用于计算枯死木生物量数据单位表21𝜹𝑺𝑶𝑪的技术内容和确定方法数据来源𝛿𝑆𝑂𝐶数值公式(C.2)数据用途造林后土壤有机碳密度平均年变化率数据/参数名称tC·hm-2·a-1应用的公式编号项目业主须按照如下优先顺序选择:a)地方标准;数据描述b)国家或行业标准中适用于项目区的数据;数据单位c)本文件及附录中推荐的缺省值;d)项目区当地、相邻地区或相似生态条件下的调查统计数据。须基于5篇以数据来源上国内外核心期刊发表的研究结果或总数不少于30个样本的调查数据的数值整合分析,且经过同行专家评议数据用途表C.1用于计算土壤有机碳储量变化数据/参数名称应用的公式编号表22𝑪𝑶𝑴𝑭的技术内容和确定方法数据描述𝐶𝑂𝑀𝐹数据单位公式(D.2)、公式(D.4)数据来源燃烧因子(针对不同的植被类型)无量纲项目业主须按照如下优先顺序选择:a)地方标准;b)国家或行业标准中适用于项目区的数据;c)本文件及附录中推荐的缺省值;d)项目区当地、相邻地区或相似生态条件下的调查统计数据。须基于5篇以上国内外核心期刊发表的研究结果或总数不少于30个样本的调查数据的整合分析,且经过同行专家评议表D.1用于计算森林火烧引起的温室气体排放量表23𝑬𝑭CH4的技术内容和确定方法𝐸𝐹CH4公式(D.2)、公式(D.4)CH4排放因子gCH4·(kgd.m.)-1项目业主须按照如下优先顺序选择:a)地方标准;b)国家或行业标准中适用于项目区的数据;13数值c)本表缺省值;数据用途d)项目区当地、相邻地区或相似生态条件下的调查统计数据。须基于5篇以上国内外核心期刊发表的研究结果或总数不少于30个样本的调查数据的整合分析,且经过同行专家评议热带森林:6.8;其它森林:4.7用于计算森林火烧引起的CH4排放量数据/参数名称表24𝑬𝑭N2O的技术内容和确定方法应用的公式编号𝐸𝐹N2O数据描述公式(D.2)、公式(D.4)数据单位N2O排放因子数据来源gN2O·(kgd.m.)-1数值项目业主须按照如下优先顺序选择:数据用途a)地方标准;b)国家或行业标准中适用于项目区的数据;c)本表缺省值;d)项目区当地、相邻地区或相似生态条件下的调查统计数据。须基于5篇以上国内外核心期刊发表的研究结果或总数不少于30个样本的调查数据的整合分析,且经过同行专家评议热带森林:0.20;其它森林:0.26用于计算森林火烧引起的N2O排放数据/参数名称表25𝑮𝑾𝑷CH4的技术内容和确定方法应用的公式编号𝐺𝑊𝑃CH4数据描述公式(D.2)、公式(D.4)数据单位100年时间尺度下CH4的全球增温潜势数据来源无量纲本表默认值,参考IPCC第五次评估报告数值28数据用途将CH4排放量转化为CO2e数据/参数名称表26𝑮𝑾𝑷N2O的技术内容和确定方法应用的公式编号𝐺𝑊𝑃N2O数据描述公式(D.2)、公式(D.4)数据单位100年时间尺度下N2O的全球增温潜势数据来源无量纲本表默认值,参考IPCC第五次评估报告数值265数据用途将N2O排放量转化为CO2e7.2项目实施阶段需监测的参数和数据项目实施阶段需确定的参数和数据的技术内容和确定方法见表27—表34。表27𝑨𝒊,𝒋,𝒕的技术内容和确定方法数据/参数名称𝐴𝑖,𝑗,𝑡应用的公式编号公式(A.2)、公式(A.3)、公式(A.4)、公式(B.2)14数据描述第t年时,第i项目碳层树种j的森林面积数据单位数据来源hm2监测点要求监测仪表要求野外测定监测程序与方法要求所有实际实施造林活动的项目地块及其拐点坐标实时动态差分技术(RTK)、GPS、BDS等导航设备、高分辨率卫星影像和监测频次与记录要求大比例尺地形图质量保证/质量控制程核对实际实施的项目地块及其拐点坐标与项目设计是否一致,针对不一致的序要求地方:数据用途a)位于项目设计边界之外的部分,不得纳入项目边界内;b)在监测时,项目设计边界内尚未实际实施造林的部分地块,如果面积≥400m2,须单独纳入新的碳层或移出项目边界外(如改变土地用途),并重新测定相关部分的项目边界坐标自首次核查后,一般每5年至少监测一次。须有项目及碳层边界坐标的.shp或.kml文件采用森林资源规划设计调查技术规程(GB/T26424)和森林资源连续清查技术规程(GB/T38590)使用的质量保证和质量控制(QA/QC)程序用于计算项目清除量表28𝑫𝑩𝑯的技术内容和确定方法数据/参数名称𝐷𝐵𝐻应用的公式编号公式(A.10)、公式(A.15)、公式(A.21)数据描述乔木或竹子的胸径数据单位数据来源cm监测点要求野外测定监测仪表要求所有野外监测样地胸径测量仪、皮尺;罗盘、RTK、GPS、BDS等定位和导航设备监测程序与方法要求采用森林资源规划设计调查技术规程(GB/T26424)和森林资源连续清查技术规程(GB/T38590)使用的标准操作程序(SOP)监测频次与记录要求自首次核查后,一般每5年至少监测一次。精确到小数点后一位质量保证/质量控制程采用森林资源规划设计调查技术规程(GB/T26424)和森林资源连续清查技序要求术规程(GB/T38590)使用的质量保证和质量控制(QA/QC)程序数据用途用于计算监测样地的单位面积生物量表29𝑯的技术内容和确定方法数据/参数名称𝐻应用的公式编号公式(A.10)、公式(A.15)、公式(A.21)、公式(A.23)数据描述乔木(或竹子、灌木)的高度数据单位数据来源m监测点要求野外测定监测仪表要求所有野外监测样地测高仪、皮尺;罗盘、RTK、GPS、BDS等定位和导航设备监测程序与方法要求采用森林资源规划设计调查技术规程(GB/T26424)和森林资源连续清查技术规程(GB/T38590)使用的标准操作程序(SOP)监测频次与记录要求自首次核查后,一般每5年至少监测一次。精确到小数点后一位质量保证/质量控制程采用森林资源规划设计调查技术规程(GB/T26424)和森林资源连续清查技序要求术规程(GB/T38590)使用的质量保证和质量控制(QA/QC)程序数据用途用于计算监测样地的单位面积生物量15数据/参数名称表30𝑺𝑪BF,𝒕的技术内容和确定方法应用的公式编号数据描述𝑆𝐶BF,𝑡数据单位公式(A.22)数据来源竹林t年时,累计择伐地上生物量占其单位面积地上生物量的比例(如株数监测点要求比例)监测仪表要求无量纲野外测定监测程序与方法要求所有野外监测样地测高仪、皮尺;罗盘、RTK、GPS、BDS等定位和导航设备监测频次与记录要求采用森林资源规划设计调查技术规程(GB/T26424)和森林资源连续清查技质量保证/质量控制程术规程(GB/T38590)使用的标准操作程序(SOP)序要求每次择伐时记录数据用途采用森林资源规划设计调查技术规程(GB/T26424)和森林资源连续清查技术规程(GB/T38590)使用的质量保证和质量控制(QA/QC)程序数据/参数名称用于计算监测样地的单位面积生物量应用的公式编号数据描述表31𝑪𝑪SF,𝒕的技术内容和确定方法数据单位数据来源𝐶𝐶SF,𝑡监测点要求公式(A.24)监测仪表要求第t年时,灌木林覆盖度,用小数表示(例如覆盖度10%记为0.10)无量纲监测程序与方法要求野外测定灌木冠幅和株数,换算成单位面积灌木覆盖度所有野外监测样地监测频次与记录要求皮尺;罗盘、RTK、GPS、BDS等定位和导航设备质量保证/质量控制程采用森林资源规划设计调查技术规程(GB/T26424)和森林资源连续清查技序要求术规程(GB/T38590)使用的标准操作程序(SOP)数据用途自首次核查后,一般每5年至少监测一次。精确到小数点后一位采用森林资源规划设计调查技术规程(GB/T26424)和森林资源连续清查技数据/参数名称术规程(GB/T38590)使用的质量保证和质量控制(QA/QC)程序应用的公式编号用于计算监测样地的单位面积生物量数据描述数据单位表32𝑨BURN,𝒊,𝒕的技术内容和确定方法数据来源监测点要求𝐴BURN,𝑖,𝑡监测仪表要求公式(D.2)、公式(D.3)第t年时,第i项目碳层发生燃烧的面积监测程序与方法要求hm2野外测定监测频次与记录要求所有野外监测样地皮尺;罗盘、RTK、GPS、BDS等定位和导航设备采用森林资源规划设计调查技术规程(GB/T26424)和森林资源连续清查技术规程(GB/T38590)使用的标准操作程序(SOP)。对于在坡地上的样地,须进行坡度校正火灾发生后当年监测。火烧地块位置应用GPS或BDS定位并记录经纬度坐标(以度表示的坐标至少保留6位小数)、位置(县、乡、村和小地名)以及地块的形状和大小16质量保证/质量控制程采用森林资源规划设计调查技术规程(GB/T26424)和森林资源连续清查技序要求术规程(GB/T38590)使用的质量保证和质量控制(QA/QC)程序数据用途用于计算火烧引起的温室气体排放量数据/参数名称表33𝑹BURN,𝒊,𝒕的技术内容和确定方法应用的公式编号数据描述𝑅BURN,𝑖,𝑡数据单位公式(D.4)数据来源第t年时,第i项目碳层烧除的病原疫木的数量占比监测点要求无量纲监测仪表要求野外测定所有野外监测样地监测程序与方法要求胸径测量仪、皮尺、测高仪、罗盘、RTK、GPS、BDS等定位和导航设备采用森林资源规划设计调查技术规程(GB/T26424)和森林资源连续清查技监测频次与记录要求术规程(GB/T38590)使用的标准操作程序(SOP)质量保证/质量控制程伐除前监测。记录伐除木树种、胸径、树高、单位面积株数等序要求采用森林资源规划设计调查技术规程(GB/T26424)和森林资源连续清查技数据用途术规程(GB/T38590)使用的质量保证和质量控制(QA/QC)程序用于计算火烧引起的温室气体排放量数据/参数名称表34𝒙𝒌𝐥的技术内容和确定方法应用的公式编号数据描述𝑥𝑘l数据单位公式(G.4)数据来源第k株单木的激光雷达特征参数监测点要求m或无量纲监测仪表要求野外测定监测程序与方法要求所有野外监测样地监测频次与记录要求无人机或有人机、机载激光雷达、RTK、GPS、BDS等定位和导航设备采用机载激光雷达数据获取技术规范(CH/T8024)使用的标准操作程序质量保证/质量控制程(SOP)序要求与样地监测频次同时进行,首次核查后,一般每5年至少监测一次。无人机数据用途点云密度不得少于100个/m2采用机载激光雷达数据获取技术规范(CH/T8024)使用的标准操作程序(SOP)用于计算乔木或竹林生物量7.3项目实施及监测的数据管理要求7.3.1一般要求7.3.1.1项目业主应采取以下措施,确保监测参数和数据的质量:a)遵循项目设计阶段确定的数据监测程序与方法要求,制定详细的监测方案;b)建立可信且透明的内部管理制度和质量保障体系,包括但不限于可靠的外业测定、外业测定的互检互核、内业数据的输入、计算和核实等;c)明确负责部门及其职责、具体工作要求、数据管理程序、工作时间节点等;d)指定专职人员负责项目边界、项目实施情况、测树因子、火烧等数据的监测、收集、记录和交叉核对。7.3.1.2鼓励项目减排量收益至少不低于90%返给具体实施了造林活动,并拥有林木的所17有权或经营权的项目实施主体。7.3.2项目边界监测要求7.3.2.1在项目设计阶段,项目业主须明确项目计划造林的地块边界,并提供所有项目地块边界的矢量数据文件。在项目实施阶段,项目业主须测量项目实际造林的地块边界。7.3.2.2在计入期内,项目业主须根据监测方案对项目边界进行监测,检查项目实际边界是否与项目设计文件一致。如果实际边界位于项目设计文件描述的边界之外,则边界以项目设计文件为准;如果实际边界位于项目设计文件描述的边界之内,则以实际边界为准,并提供新的项目边界矢量数据文件。7.3.2.3如果项目边界发生任何变化,例如土地利用类型发生变化,应测定被征占地块的地理坐标和面积,将这部分地块调出项目边界,并在后续减排量核算报告中予以说明,之后不再纳入项目边界。7.3.3项目实施情况监测要求项目实施阶段,主要监测和记录项目边界内所发生的造林、管护以及与温室气体排放有关项目活动的实施情况,并判断是否与项目设计文件及监测方案一致。主要内容包括:a)造林活动:造林时间、造林地块、造林树种、造林密度、苗木成活率和保存率、整地清林方式等;b)管护活动:巡护、补植、采伐、有害生物防治和森林火灾预防措施等;c)项目边界内自然灾害(如火灾、病虫害、台风、干旱等)和人为干扰(如土地利用变化等)的发生情况(如时间、地点、面积、边界、损害强度等)。7.3.4项目碳层划分要求项目实施阶段,如果项目边界内出现下述情形之一,项目业主须在每次监测前对上一次划分的碳层进行调整:a)项目实际活动与项目设计文件不一致,并影响了项目碳层内的均一性,如造林时间、树种选择、造林面积以及边界等发生变化;b)因自然因素(如立地条件、火灾、病虫害等)或人为干扰(如火烧、采伐等)导致碳层内的变异性增加;c)因土地利用类型变化等造成碳层边界发生变化。若上一次监测发现,两个或多个碳层具有相近的碳储量及变化,则可将这些不同的碳层合并成一个碳层,以降低监测工作量。7.3.5抽样设计要求本文件要求对项目生物质碳储量进行抽样监测,监测应达到90%可靠性水平下90%的精度要求。项目业主须按照附录E步骤计算获得抽样监测所需的样地数量及在各碳层中的分布。7.3.6样地设置要求项目生物质碳储量的变化可采用固定样地连续监测。项目业主须按照附录E步骤,采用随机起点、系统布点的方法设置样地。187.3.7监测频率与时间要求项目业主应在项目设计阶段确定固定样地监测频率,一般每5年至少监测一次。首次监测时间不早于项目申请登记时间。7.3.8项目生物质碳储量监测与计算要求项目业主须按照附录F步骤,通过在项目样地监测得到的平均单位面积年生物质碳储量,计算项目边界内生物质碳储量的年变化量。7.3.9数据精度控制与校正要求基于样地的生物质碳储量的抽样调查(见附录F),要求达到90%可靠性水平下90%的精度要求。如果测定精度低于该值,项目业主可通过增加样地数量进行补测,从而使测定结果达到精度要求;或选择扣减一定比例清除量的方式进行校正。Δ𝐶BiomassPROJ,𝑡=Δ𝐶Biomass,𝑡×(1−𝐷𝑅)(4)式中:Δ𝐶BiomassPROJ,𝑡——校正后第t年的项目生物质碳储量变化量,单位为吨二氧Δ𝐶Biomass,𝑡化碳当量每年(tCO2e·a-1);——监测的第t年的项目生物质碳储量变化量,单位为吨二氧DR化碳当量每年(tCO2e·a-1);t——扣减率,单位为百分比(%);——自项目开始以来的年数,t=1,2,3……,无量纲。扣减率可从表35获得,不确定性的计算见附录F。表35样地监测生物质碳储量变化量的扣减率不确定性扣减率(DR)小于或等于10%0%大于10%但小于或等于20%6%大于20%但小于或等于30%11%大于30%须增加样地数量,直至测定结果达到精度要求为保证监测数据的可靠性,在现有监测方法的基础上,本文件鼓励使用机载激光雷达技术作为林木生物质碳储量监测的辅助校核工具。在拟开发项目区域首次应用机载激光雷达技术开展监测时,须与样地监测数据比照,充分检验该技术在当前场景下的适用性(见附录G)。在确保能够满足最低精度要求后,方可采用其计算项目范围内林木生物质碳储量。若基于样地监测方法与基于激光雷达方法的计算结果相比误差小于10%,可基于保守性原则,在首次监测时选取较小值作为项目当年的林木生物质碳储量。在满足上述精度与误差要求的基础上,后续可选用机载激光雷达方法进行监测,并复用首次监测时建立的激光雷达生物量模型进行计算。若首次监测时,基于激光雷达方法计算结果与基于样地监测方法相比较大,则须采用首次监测时基于样地监测方法与基于激光雷达方法计算结果的比值,向下修正基于激光雷达方法的计算结果。l𝐶Bmiomass,𝑇1𝐶Biomass,𝑡=𝐶Biomass,𝑡×l(5)𝐶Biomass,𝑇1式中:19𝐶Biomass,𝑡——项目第t年时的森林生物质碳储量,单位为吨碳(t𝐶Bliomass,𝑡——C);𝐶Bliomass,𝑇——第t年时,使用激光雷达方法监测的森林生物质碳储——量,单位为吨碳(tC);1——首次监测时(𝑇1),使用激光雷达方法监测的森林生物质碳储量,单位为吨碳(tC);𝐶Bmiomass,𝑇首次监测时(𝑇1),使用样地监测方法监测的森林生1物质碳储量,单位为吨碳(tC);l通过激光雷达获取或计算的值,无量纲;m——样地实测直接获取或计算的值,无量纲。在计算最终生物质碳储量变化量时,根据激光雷达生物量模型精度检验结果,采用扣减的方法进行校正,扣减率见表36。表36激光雷达监测的生物质碳储量变化量的扣减率精度检验R2值扣减率(DR)大于或等于0.90%大于或等于0.8但小于0.96%大于或等于0.7但小于0.811%小于0.7须优化模型或额外增加样地数量,直至测定结果达到精度要求7.3.10项目火烧排放监测和计算要求7.3.10.1项目业主须详细记录项目边界内每一次森林火灾发生的时间、面积、地理边界等信息,根据火烧前最近一次核查确定的单位面积生物量,计算项目边界内由于森林火灾燃烧地上生物质和死有机质所引起的温室气体排放量。7.3.10.2如果发生有计划的人为清除病原疫木引起的火烧,项目业主须在伐前监测各样地内病原疫木的树种、胸径、树高、数量等,计算清除的病原疫木的生物质碳储量,进而计算人为火烧引起的温室气体排放量。7.3.11数据管理与归档要求7.3.11.1对于收集到的监测数据,项目业主应建立数据、信息等原始记录和台账管理制度,妥善保管监测数据、原始记录(植被调查、生态综合调查数据)、证明材料(权属证明文件、土地合格性证明)相关的书面文件等。原始记录和台账应明确数据来源、数据获取时间及填报台账的相关责任人等信息。7.3.11.2项目监测的所有数据均应进行电子存档,在该温室气体自愿减排项目最后一期减排量登记后至少保存10年,确保相关数据可被追溯。8项目审定与核查要点8.1项目适用条件的审定与核查要点审定与核查机构应基于项目设计文件,对本文件适用条件进行逐条分析,重点确定以下内容:a)核实项目是否符合法律、法规要求,符合行业发展政策。可查阅《中华人民共和国森林法》《中华人民共和国森林法实施条例》《中华人民共和国土地管理法》等法律及20造林相关的法规和政策,确认项目不违反有关法律法规和政策;b)核实项目地块的合格性:——可通过项目开始前项目边界内的土地利用现状图、现场照片或高分辨率的地理空间数据(如卫星影像、航片)、林地“一张图”或林草资源“一张图”、森林资源规划调查数据、造林作业设计、同行评议的研究报告等确认项目边界内地块在项目开始前符合本文件的适用条件;——对于无法通过上述方式确定造林地块合格性的,应采用现场走访的方式确认,并以影像等可证明的方式记录赴现场开展审核活动的点位和线路。c)核实项目边界内的土地权属。在项目审定和减排量核查时,须核对项目全部的土地所有权(或使用权)或林木所有权(或经营权)的证据,如不动产权属证书、土地承包或流转合同;或经有批准权的人民政府或主管部门批准核发的土地证、林权证等;如果项目业主不是项目边界内土地(或林木)权属所有人,项目业主应取得权属人授予的相关权利,并提供相关证明文件;d)核实造林活动对项目地块的影响。可基于造林作业设计、卫星影像和实地走访等方式核查项目区域原有植被的移除比例、土壤的扰动程度、清林方式、造林模式等。8.2项目开始时间的审定与核查要点审定与核查机构须通过对比项目开始前后的卫星遥感影像或现地走访(如查数针叶树轮枝数或年轮数确定树木年龄),并结合证据文件核实等方法,验证项目开始时间的真实性。项目业主可选择提供下列材料之一,说明项目的开始时间:a)经县级(含)以上行业主管部门批复的作业设计和(或)出具的验收报告;b)项目业主与施工方签署的施工合同和相关付款证明;c)其他具有法律效力的、注明项目开始日期的文件(如项目监理报告)。8.3项目边界的审定与核查要点审定与核查机构须根据项目业主提供项目边界的矢量数据文件(如.shp文件或.kml文件,并细化到地块),重点开展以下工作:a)从每个碳层中以随机方式选取至少1个地块(或总共不少于5个地块),利用BDS或GPS系统,直接测定项目地块的全部多边形边界及其拐点坐标,核实单点定位误差是否超过±5m。根据重要拐点坐标定位,计算选取项目地块的面积,与项目业主的测定结果进行对比,核实项目边界面积误差是否超过±5%;b)通过遥感影像或实地走访,确认项目边界内是否包含宽度大于3m的道路、沟渠、坑塘、河流等不符合适用条件的土地;c)通过项目所在地遥感影像、造林作业设计、造林验收报告等资料,核对实际造林地块的边界与项目设计文件中计划实施的边界是否一致,识别项目实施与项目设计是否出现偏移,并确认出现的偏移是否按照6.3和7.3.4要求调整碳层划分;d)核实项目边界内土地利用类型是否发生变化。对土地利用方式已经发生变化的地块,需要从项目边界内调出。8.4项目减排量核算的审定与核查要点8.4.1审定与核查机构须核实项目减排量核算过程符合本文件的要求,项目实施阶段每次监测和计算方法一致,参数选择合理,计算结果准确且符合保守性原则。8.4.2在项目设计阶段,基准线清除量经过审定后,在整个计入期内都是有效的。在项目实施阶段,项目业主可选择不对基准线清除量进行监测,核查机构也不需要对基准线清除21量进行核查。8.4.3如果项目业主选择对项目开始时的原有散生木(竹)进行标记,在项目实施阶段,项目核查机构须随机抽取不少于10个样地,核对原有散生木(竹)的标记情况,株数误差不超过1株,并确认没有人为移除原有散生木(竹)的情况发生。如果项目业主未选择对项目开始时的原有散生木(竹)进行标记,在项目设计阶段或项目实施阶段,审定与核查机构须根据项目业主提供的项目开始前的森林资源调查数据、或利用高清卫星影像(分辨率不超过2m)来确定项目边界内原有散生木(竹)的冠层盖度,与项目业主测定值误差不超过10%。8.4.4审定与核查机构须对项目业主自选参数的真实性和保守性进行核实。自选参数必须是基于5篇以上国内外核心期刊发表的或总数不少于30个样本的调查数据的整合分析结果,且经过同行专家评议。整合分析须采用剔除异常值之后的样本数据,计算参数平均值和95%置信区间,或构建模拟模型且决定系数(R2值)不低于0.7,并将参数或模型的不确定性整合到碳储量变化计算结果的不确定性。不确定性的整合方法可参考《IPCC2000不确定性管理优良做法指南》。8.5样地监测的审定与核查要点确认项目是否按照本文件的要求制定了监测计划并实施,重点审定与核查以下要点:a)确认监测计划是否包含了监测实施的组织形式和职责分工,监测方法、程序和频次,数据记录与收集程序,抽样方案等;b)确认项目碳层划分、抽样设计和样地设计是否满足90%可靠性水平下90%的精度要求,是否满足7.3.9节的要求;c)确认项目监测阶段项目碳层调整与地块生物质碳储量异质性变化的符合性,可使用项目开始时和发生干扰时的卫星影像进行对比,确定项目实施阶段项目碳层调整的合理性;d)确认固定样地的布设是否根据7.3.6节执行;e)现场核查须在监测完成后6个月之内完成。核查机构须从项目所有监测样地中随机选择至少10个样地,且每个碳层至少抽1个监测样地(以数量多的为准),进行现场测定核查。首先须核实监测样地与所属碳层样地外的造林与经营措施的一致性,确定无误后开展样地测定核查。测定内容包括:样地位置、面积以及每木检尺的株数、胸(地)径、树(株)高或灌木覆盖度,并与项目业主的测定结果进行对比。在误差允许范围内,使用业主的测量值;在误差允许范围之外,项目业主须重新监测和核算。样地监测的平均允许误差如下:——样地位置:样地中心点复位误差不超过±5m;——样地面积:样地面积与核算报告描述面积一致;——株数:胸(地)径≥2cm的检尺株数测量误差不超过±5%,最多不超过±3株;——胸(地)径:样地平均胸(地)径测量误差不超过±5%;——灌木覆盖度:平均覆盖度误差不超过±10%。8.6参数的审定与核查要点及方法参数的审定与核查要点及方法见表37。表37参数的审定与核查要点及方法序号内容审定要点及方法核查要点及方法1𝐶𝐹确认参数的来源并评估其适用a)查阅项目减排量核算报告中的参222𝑅𝑆𝑅性:数取值是否与项目设计文件一3𝐵𝐶𝐸𝐹a)若参数来源为地方标准、致、准确;4𝐵𝐸𝐹行业标准、国家标准或本b)须对参数选取的适用性、准确文件及附录中推荐的缺省性、保守性进行核查,核查是否5𝑆𝑉𝐷值,则须审定参数是否适符合7.1对参数的技术内容和确6𝐴𝐺𝐵BF,𝑇𝑏用于项目所在地和所选树定方法的要求。7𝐴𝐺𝐵SFb)种若;参数为项目业主自选参8𝐷𝐹LI数,则须判断参数选择是9𝐷𝐹DW否符合以下要求:10𝛿𝑆𝑂𝐶——参数为项目区当地、相邻地区或相似生态条件11𝐶𝑂𝑀𝐹下的调查统计数据;——参数为基于5篇以上12𝐸𝐹CH413𝐸𝐹N2O国内外核心期刊发表的研14𝑓(𝐷𝐵𝐻,𝐻)究结果的整合;或经过同行专家评议的不少于3015𝑓(𝑉AF,𝑡)个样本的调查数据的整合须对参数选取的适用性、准确性、保结果。守性进行核查,核查是否符合7.1对参16𝑓AGB,BF(𝐷𝐵𝐻,𝐻)数技术内容和确定方法的要求。17𝑓SF(𝑥1,𝑥2,𝑥3……)a)确认项目碳层边界的准确性:a)确认项目碳层边界的准确——若利用BDS或GPS系统,性:直接测定项目地块边界的拐点坐——通过查阅项目设计文标,应现场核查拐点情况,并确件、造林作业设计等,确认拐点坐标采集设备精度不低于认项目边界的实际位置;5m;——确认空间数据分辨率——若利用地理空间数据(如卫不低于5m,保障项目边星遥感影像、无人机航拍影界的准确性;像),在GIS辅助下直接读取项18𝐴𝑖,𝑗,𝑡b)确认项目碳层面积的准确目地块的边界坐标,应确认空间性:数据分辨率不低于5m,同时现——通过查阅项目设计文场核查边界坐标的实际情况;件、造林作业设计等,确——比较项目边界和碳层核查面认碳层划分是否符合6.3积与项目业主的监测结果误差。节要求;b)确认项目碳层调整的准确性:——确认空间数据分辨率——查阅项目监测记录,确认项不低于5m,保障地块边目碳层调整是否符合7.3.4节要求;界的准确性。——确认碳层调整后,项目减排量核算是否符合7.3.9节要求。a)查阅项目减排量核算报告中的参根据项目开始前的森林资源调数取值是否与项目设计文件一查数据,或分辨率不超过2m的致、准确;19𝐶𝐷PE,𝑖,𝑡0高清卫星影像数据,审定该数b)须对参数选取的适用性、准确据是否准确。性、保守性进行核查,核查是否符合7.1对参数的技术内容和确定方法的要求。20𝐷𝐵𝐻确认测树因子监测的准确性:a)从监测样地中,随机选择至少1021𝐻个样地,且每个碳层至少抽1个监测样地,以数量多的为准进行现场测定,包括样地位置、面积22𝑆𝐶BF,𝑡项目审定阶段不包括该参数。以及株数、胸(基)径和株高;b)比较核查测树因子与项目业主的23𝐶𝐶SF,𝑡监测结果误差是否符合8.5节的要求。24𝐴BURN,𝑖,𝑡核查管理日志,确定项目边界内火灾发生的时间、面积、地理边界等信息2325𝑅BURN,𝑖,𝑡与项目业主的监测结果相比是否在误差范围内。26𝑥𝑘l核查病原疫木伐除记录,确定各样地内疫木的树种、胸径、树高、数量等与项目业主的监测结果相比是否在误差范围内。核查机载激光雷达数据的获取是否符合《CH/T8024-2011机载激光雷达数据获取技术规范》。9方法学编制单位在本文件编制工作中,中国林业科学研究院森林生态环境与自然保护研究所,以及大自然保护协会北京代表处、北京林业大学生态与自然保护学院、清华大学能源环境经济研究所、中国绿色碳汇基金会、扬州大学植物保护学院生态学系、中国林业科学研究院资源信息研究所、北京天德泰科技股份有限公司、中国林业集团有限公司、大兴安岭林业集团公司、腾讯科技(北京)有限公司、国家林业和草原局林草调查规划院、中国林业科学研究院林业研究所、生态环境部环境与经济政策研究中心、中国内蒙古森林工业集团有限责任公司、北京汇智绿色资源研究院、中国计量科学研究院、中国林业与环境促进会、北京碳宝科技有限公司、中建生态环境集团有限公司、中国林场协会、江南大学国家安全与绿色发展研究院、北京和碳环境技术有限公司、中国民贸一乡一品产业促进中心、海南斯兰低碳研究中心有限公司、广西森工集团股份有限公司、广西壮族自治区国有高峰林场、广西力源宝科技有限公司、广西壮族自治区国有三门江林场、广西壮族自治区国有黄冕林场、广西壮族自治区国有维都林场等单位作出了积极贡献。24附录A(资料性附录)森林生物质碳储量变化计算方法A.1森林生物质碳储量变化的计算本文件采用“储量变化法(StockDifferenceMethod)”计算项目边界内的森林生物质碳储量在一段时期内的年均变化量:Δ𝐶Biomass,𝑡=𝐶Biomass,𝑡2−𝐶Biomass,𝑡1×44(A.1)𝑡2−𝑡112式中:Δ𝐶Biomass,𝑡——项目开始第t年的森林生物质碳储量的年变化量,单位为吨二氧化碳当量每年(tCO2e·a-1);𝐶Biomass,𝑡1——第t1年时,森林生物质碳储量,单位为吨碳(tC);𝐶Biomass,𝑡2——第t2年时,森林生物质碳储量,单位为吨碳(tC);44/12——二氧化碳与碳的相对分子质量之比,无量纲;t——自项目开始以来的年数,t=1,2,3……,无量纲;t1,t2——项目开始后的第t1年和第t2年,单位为年(a),且t1≤t≤t2。A.2森林生物质碳储量的计算首先分别计算各碳层内各树种(含竹子和灌木)的全林生物质碳储量。当无法直接计算全林生物质碳储量时,可选择分别计算各碳层内各树种的地上和地下生物质碳储量,然后加总得到全林生物质碳储量:𝐶Biomass,𝑡=∑∑(𝐴𝑖,𝑗,𝑡×𝐵Total,𝑖,𝑗,𝑡×𝐶𝐹Total,𝑖,𝑗)(A.2)(A.3)𝑖𝑗(A.4)(A.5)𝐶AGB,𝑡=∑∑(𝐴𝑖,𝑗,𝑡×𝐴𝐺𝐵𝑖,𝑗,𝑡×𝐶𝐹AGB,𝑖,𝑗)𝑖𝑗𝐶BGB,𝑡=∑∑(𝐴𝑖,𝑗,𝑡×𝐵𝐺𝐵𝑖,𝑗,𝑡×𝐶𝐹BGB,𝑖,𝑗)𝑖𝑗𝐶Biomass,𝑡=𝐶AGB,𝑡+𝐶BGB,𝑡式中:𝐶Biomass,𝑡——第t年时,森林生物质碳储量,单位为吨碳(tC);𝐶AGB,𝑡——第t年时,森林地上生物质碳储量,单位为吨碳(tC);𝐶BGB,𝑡——第t年时,森林地下生物质碳储量,单位为吨碳(tC);𝐴𝑖,𝑗,𝑡——第t年时,第i项目碳层树种j的森林面积,单位为公顷(hm2);𝐴𝐺𝐵𝑖,𝑗,𝑡——第t年时,第i项目碳层树种j的单位面积地上生物量,单位为吨每公顷(td.m.·hm-2);𝐵𝐺𝐵𝑖,𝑗,𝑡——第t年时,第i项目碳层树种j的单位面积地下生物量,单位为吨每公顷(td.m.·hm-2);𝐵Total,𝑖,𝑗,𝑡——第t年时,第i项目碳层树种j的单位面积全林生物量,单位为吨每公顷(td.m.·hm-2);25𝐶𝐹AGB,𝑖,𝑗——第i项目碳层树种j的地上生物量含碳率,单位为吨碳每吨(tC·(td.m.)-1);𝐶𝐹BGB,𝑖,𝑗——第i项目碳层树种j的地下生物量含碳率,单位为吨碳每吨(t𝐶𝐹Total,𝑖,𝑗——C·(td.m.)-1);i——第i项目碳层树种j的全林生物量含碳率,单位为吨碳每吨(tC·(td.m.)-1);j——项目碳层,i=1,2,3……,无量纲;t——树种,j=1,2,3……,无量纲;自项目开始以来的年数,t=1,2,3……,无量纲。A.3乔木林生物量的计算推荐优先使用生物量方程法,其次为生物量转换因子法,并可采用激光雷达估测法进行校核。当选择生物量方程法或生物量转换因子法时,优先选择全林生物量的计算方法,其次选择分别计算地上生物量和地下生物量的方法。在项目实施阶段,每次监测和核算均须使用同一类方法,以保证项目减排量核算结果的可比性。A.3.1生物量方程法在项目设计阶段,采用乔木林全林(或地上)生物量与单位面积蓄积量的相关方程(表A.5)计算乔木林全林(或地上)生物量。其中,乔木林单位面积蓄积量采用单位面积蓄积量与林龄的相关方程(表A.11)计算:𝐵Total,AF,𝑡=𝑓(𝑉AF,𝑡)(A.6)𝑉AF,𝑡=𝑓(𝐴𝑔𝑒AF,𝑡)(A.7)式中:𝐵Total,AF,𝑡——第t年时,乔木林单位面积全林生物量,单位为吨每公顷(td.m.·hm-2);𝑓(𝑉AF)——乔木林单位面积全林(或地上)生物量与单位面积蓄积量的相关方程,单位为吨每公顷(td.m.·hm-2);𝑓(𝐴𝑔𝑒AF)——乔木林单位面积蓄积量与林龄的相关方程,单位为立方米每公顷(m3·hm-2);𝑉AF,𝑡——第t年时,乔木林单位面积蓄积量,单位为立方米每公顷(m3·hm-2);𝐴𝑔𝑒AF,𝑡——第t年时,乔木林的林龄,无量纲。如果采用通过地上生物量计算地下生物量的方法,在得到乔木林地上生物量后,利用乔木林地下生物量与地上生物量的比值计算乔木林地下生物量,并加和得到乔木林全林生物量。𝐵𝐺𝐵AF=𝐴𝐺𝐵AF×𝑅𝑆𝑅AF(A.8)𝐵Total,AF=𝐴𝐺𝐵AF+𝐵𝐺𝐵AF(A.9)式中:——乔木林单位面积全林生物量,单位为吨每公顷(td.m.·hm-2);𝐵Total,AF26𝐴𝐺𝐵AF——乔木林单位面积地上生物量,单位为吨每公顷(td.m.·hm-2);𝐵𝐺𝐵AF——乔木林单位面积地下生物量,单位为吨每公顷(td.m.·hm-2);𝑅𝑆𝑅AF——乔木林地下生物量与地上生物量的比值,无量纲。在项目实施阶段,须使用单株林木的全株(地上或地下)生物量(表A.1、表A.2和表A.3)与林木胸径和(或)树高的相关方程,再根据单位面积林木数量,计算乔木林全林生物量。𝐵Total,AF,𝑖,𝑡=∑[𝑓(𝐷𝐵𝐻AF,𝑖,𝑗,𝑡,𝐻AF,𝑖,𝑗,𝑡)×𝑁𝑖,𝑗,𝑡]×10−3(A.10)(A.11)𝑗𝐵Total,AF,𝑡=∑𝐵Total,AF,𝑖,𝑡𝑖式中:𝐵Total,AF,𝑡——第t年时,乔木林单位面积全林生物量,单位为吨每公顷(td.m.·hm-2);𝐵Total,AF,𝑖,𝑡——第t年时,第i项目碳层乔木林单位面积全林生物量,单位为吨每公顷(td.m.·hm-2);𝑓(𝐷𝐵𝐻AF,𝐻AF)——乔木全株(或地上)生物量与胸径(单位为厘米,cm)和(或)树高(单位为米,m)的相关方程(单位为千克每株,kgd.m.·stem-1);𝐷𝐵𝐻AF,𝑖,𝑗,𝑡——第t年时,第i项目碳层乔木林第j树种的平均胸径,单位为厘米(cm);𝐻AF,𝑖,𝑗,𝑡——第t年时,第i项目碳层乔木林第j树种的平均树高,单位为米(m);j——树种,j=1,2,3……,无量纲;i——项目碳层,i=1,2,3……,无量纲;𝑁𝑖,𝑗,𝑡——第t年时,第i项目碳层乔木林树种j的单位面积株数,单位为株每公顷(stem·hm-2);t——自项目开始以来的年数,t=1,2,3……,无量纲;10-3——将千克转换为吨的常数。在项目实施阶段,如果项目首次监测时间(T1)是在项目计入期开始(t0)之后,则可根据经核查的首次监测结果,采用如下办法计算计入期开始时(第t0年)的单位面积生物量:——根据项目设计阶段预估的单位面积生物量(公式A.6和公式A.7),计算项目开始至首次监测间每一年的单位面积生物量与上一年的比值(𝑏Total,AF,𝑡⁄𝑏Total,AF,𝑡−1);——通过首次监测的项目实际单位面积生物量,计算前一年的单位面积生物量,直至计入期开始时。𝐵Total,AF,𝑇1−1=𝐵Total,AF,𝑇1⁄(𝑏Total,AF,𝑇1⁄𝑏Total,AF,𝑇1−1)(A.12)式中:𝐵Total,AF,𝑇1——经核查的首次监测时的乔木林单位面积生物量,单位为吨每公顷(td.m.·hm-2);𝐵Total,AF,𝑇1−1——首次监测前一年时,乔木林单位面积生物量,单位为吨每公顷(td.m.·hm-2);27𝑏Total,AF,𝑇1——项目设计阶段预估的第T1年时,乔木林单位面积生物量,单位为吨每公顷(td.m.·hm-2);𝑏Total,AF,𝑇1−1——项目设计阶段预估的第T1-1年时,乔木林单位面积生物量,单位为吨每公顷(td.m.·hm-2)。A.3.2生物量转换因子法首先选择利用生物量转换与扩展因子法,将乔木蓄积量(或单株材积)转换为乔木林(或单木)的全林(或地上)生物量。若计算的是全林生物量,则利用表A.5提供的全林生物量方程计算。若计算的是地上生物量,则利用乔木林地下生物量与地上生物量的比值转化为乔木林全林生物量,见公式(A.13)、公式(A.14)和公式(A.15):(A.13)𝐵Total,AF,𝑡=𝑉AF,𝑡×𝐵𝐶𝐸𝐹×(1+𝑅𝑆𝑅AF)𝐵Total,AF,𝑡=∑∑[𝑣𝑗,𝑠,𝑡×𝐵𝐶𝐸𝐹𝑗×(1+𝑅𝑆𝑅AF,𝑗)]⁄𝐴AF,𝑡(A.14)(A.15)𝑗𝑠𝑣𝑗,𝑠,𝑡=𝑓𝑣(𝐷𝐵𝐻𝑗,𝑠,𝑡,𝐻𝑗,𝑠,𝑡)式中:𝐵Total,AF,𝑡——第t年时,乔木林单位面积全林生物量,单位为吨每公顷(td.m.·hm-2);𝑉𝐴𝐹,𝑡——第t年时,乔木林的单位面积蓄积量,单位为立方米每公顷(m3·hm-2);𝑣𝑗,𝑠,𝑡——第t年时,乔木林内树种j第s株树的立木材积,单位为立方米(m3·stem-1);𝐵𝐶𝐸𝐹——基于林分的乔木林地上生物量转换与扩展因子,单位为吨每立方米(td.m.·m-3);𝐵𝐶𝐸𝐹𝑗——基于单木的乔木林树种j的地上生物量转换与扩展因子,单位为吨每立方米(td.m.·m-3);𝑅𝑆𝑅AF——基于林分的乔木林地下生物量与地上生物量的比值,无量纲;𝑅𝑆𝑅AF,𝑗——基于单木的乔木林树种j的地下生物量与地上生物量的比值,无量纲;𝐴AF,𝑡——第t年时,乔木林面积,单位为公顷(hm2);𝑓𝑣(𝐷𝐵𝐻𝑗,𝐻𝑗)——乔木林内树种j的单株立木材积方程,单位为立方米(m3),见表A.4;j——树种,j=1,2,3……,无量纲;s——乔木林内树种j的单木,s=1,2,3……,单位为株(stem);t——自项目开始以来的年数,t=1,2,3……,无量纲。其次,可利用基本木材密度(表A.10)和生物量扩展因子(表A.7)将乔木林蓄积量转化为乔木林地上生物量,再利用乔木林地下生物量与地上生物量的比值计算全林生物量:𝐵Total,AF,𝑡=𝐴𝐺𝐵AF,𝑡×(1+𝑅𝑆𝑅AF)(A.16)𝐴𝐺𝐵AF,𝑡=𝑉AF,𝑡×𝑆𝑉𝐷×𝐵𝐸𝐹(A.17)式中:𝐵Total,AF,𝑡——第t年时,乔木林单位面积全林生物量,单位为吨每公顷(td.m.·hm-2);28𝑉AF,𝑡——第t年时,乔木林的单位面积蓄积量,单位为立方米每公顷(m3·hm-2);𝑅𝑆𝑅AF——基于林分的乔木林地下生物量与地上生物量的比值,无量纲;𝐴𝐺𝐵AF,𝑡——第t年时,乔木林单位面积地上生物量,单位为吨每公顷(td.m.hm-2);𝑆𝑉𝐷——乔木林的基本木材密度,单位为吨每立方米(td.m.m-3);𝐵𝐸𝐹——乔木林的生物量扩展因子,是地上生物量与树干生物量的比值,无量纲;t——自项目开始以来的年数,t=1,2,3……,无量纲。在项目设计阶段可采用乔木林单位面积蓄积量随林龄变化的模型(表A.11)计算从项目开始以来逐年的单位面积蓄积量(公式A.7)。在项目实施阶段则须通过监测乔木林密度、平均树高和(或)胸径等(见附录F)来计算单位面积蓄积量。A.3.3机载激光雷达估测法利用机载激光雷达数据进行单木分割,获取单木位置、树高和冠幅,并与实际中单木位置进行精准匹配。通过实测不同径级的单木样株生物量,构建单木生物量模型。基于激光雷达单木分割获取的每株单木的激光点云数据,计算每株单木的地上生物量,以及单位面积地上生物量:𝐴𝐺𝐵=∑𝑖∑𝑗∑𝑠𝑓AGB,AF(𝑥1,𝑖,𝑗,𝑠,𝑡,𝑥2,𝑖,𝑗,𝑠,𝑡,𝑥3,𝑖,𝑗,𝑠,𝑡……)×10−3(A.18)AF,𝑡𝐴𝐴𝐹,𝑡式中:𝐴𝐺𝐵AF,𝑡——第t年时,乔木林单位面积地上生物量,单位为吨每公顷(td.m.·hm-2);𝑓AGB,AF(𝑥1,𝑥2,𝑥3……)——乔木单株地上生物量与激光点云数据的相关方程,单位为千克每株(kgd.m.·stem-1)。其中,x1,x2,x3……为树高、冠幅、枝下高等;𝐴AF,𝑡——第t年时,乔木林面积,单位为公顷(hm2);i——项目碳层,i=1,2,3……,无量纲;j——树种,j=1,2,3……,无量纲;s——第i项目碳层内第j类树种的单木,s=1,2,3……,单位为株(stem);t——自项目开始以来的年数,t=1,2,3……,无量纲。基于上述方法获得乔木林单位面积地上生物量,再利用乔木林地下生物量与地上生物量的比值转换成单位面积全林生物量。A.4竹林生物量的计算单位面积竹林生物量包括地上生物量和地下生物量:𝐵Total,BF=𝐴𝐺𝐵BF+𝐵𝐺𝐵BF(A.19)式中:——竹林单位面积全林生物量,单位为吨每公顷(td.m.·hm-2);𝐵Total,BF29𝐴𝐺𝐵BF——竹林单位面积地上生物量,单位为吨每公顷(td.m.·hm-2);𝐵𝐺𝐵BF——竹林单位面积地下生物量,单位为吨每公顷(td.m.·hm-2)。A.4.1竹林地上生物量的计算方法新造竹林生长发育可分为发育阶段(一般大径散生竹林1-9年、小径散生竹林1-5年、丛生竹1-5年、混生竹1-6年)和成林稳定阶段(一般大径散生竹林从第10年开始、小径散生竹林从第6年开始、丛生竹从第6年开始、混生竹从第7年开始)。在发育阶段,竹林的生物量、株数、平均胸径、平均竹高等都会发生明显的变化。而达到成林稳定阶段后,由于择伐或自然枯损以及新竹的生长,竹林地上生物量基本上处于动态平衡状态。假定竹林达到成熟稳定的竹龄为Tb,在竹龄(tb)达到成熟稳定年龄之前(tb≤Tb),可采用下列方法计算竹林地上生物量。竹林达到成林稳定后(tb>Tb),则等于第Tb年时的竹林地上生物量。在项目设计阶段,采用平均生长量的方法计算竹林地上生物量:𝐴𝐺𝐵BF,𝑇𝑏×𝑡𝑡≤𝑇𝐴𝐺𝐵BF,𝑡𝑏={𝑇𝑏𝑏𝑏𝑏(A.20)𝐴𝐺𝐵BF,𝑇𝑏𝑡𝑏>𝑇𝑏式中:𝐴𝐺𝐵BF,𝑡𝑏——第𝑡𝑏年时,竹林单位面积地上生物量,单位为吨每公顷(td.m.·hm-2);𝐴𝐺𝐵BF,𝑇𝑏——竹林达到成熟稳定后的单位面积地上生物量,单位为吨每公顷(td.m.·hm-2),见表A.12;𝑇𝑏——竹林达到成熟稳定的年龄,无量纲;𝑡𝑏——项目竹龄,无量纲。在项目实施阶段,监测竹林单位面积株数、平均胸径和(或)竹高等(见附录F),采用单株地上生物量与胸径和(或)竹高的相关方程计算单位面积地上生物量。𝐴𝐺𝐵BF,𝑡𝑏=𝑓AGB,BF(𝐷𝐵𝐻𝑡𝑏,𝐻𝑡𝑏)×𝑁𝑡𝑏×10−3(A.21)式中:𝐴𝐺𝐵BF,𝑡𝑏——第tb年时,竹林单位面积地上生物量,单位为吨每公顷(td.m.·hm-2);𝑓AGB,BF(𝐷𝐵𝐻,𝐻)——竹子单株地上生物量与胸径和(或)竹高的相关方程,单位为千克每株(kgd.m.·stem-1),见表A.13和表A.14;𝐷𝐵𝐻𝑡𝑏——第tb年时,竹林平均单株胸径,单位为厘米(cm);𝐻𝑡𝑏——第tb年时,竹林平均单株高度,单位为米(m);𝑁𝑡𝑏——第tb年时,竹林的单位面积株数,单位为株每公顷(stem·hm-2);𝑡𝑏——项目竹龄,无量纲;10-3——将千克转换为吨的常数。30A.4.2竹林地下生物量的计算方法由于竹林经营通常只移除地上部分,而地下部分(竹蔸、竹根和竹鞭)仍会在较长时间内留存于林地中,竹林地下生物质碳储量通常还会继续增加,即竹林地下生物量与地上生物量的比值会随着竹林年龄的增加而增加,呈现动态变化关系。在竹林达到成熟稳定的年龄前,通过竹林地下生物量与地上生物量的比值,结合竹林地上生物量的变化,计算竹林地下生物量。当竹林成熟稳定后,经过一段时间的经营,考虑到地下生物量生长也存在上限,本文件保守地假定tb>2Tb时,地下生物量不再增长。𝐴𝐺𝐵BF,𝑡𝑏×𝑅𝑆𝑅BF,𝑡𝑏≤𝑇𝑏𝐵𝐺𝐵BF,𝑡𝑏=𝐴𝐺𝐵BF,𝑇𝑏×𝑅𝑆𝑅BF+𝑆𝐶BF,𝑡𝑏×𝐴𝐺𝐵BF,𝑇𝑏×𝑅𝑆𝑅BF,𝑇𝑏<𝑡𝑏≤2𝑇𝑏(A.22)𝑡𝑏>2𝑇𝑏{𝐴𝐺𝐵BF,𝑇𝑏×𝑅𝑆𝑅BF+𝑆𝐶BF,2𝑇𝑏×𝐴𝐺𝐵BF,𝑇𝑏×𝑅𝑆𝑅BF,式中:𝐵𝐺𝐵BF,𝑡𝑏——第tb年时,竹林单位面积地下生物量,单位为吨每公顷(td.m.·hm-2);𝐴𝐺𝐵BF,𝑡𝑏——第tb年时,竹林单位面积地上生物量,单位为吨每公顷(td.m.·hm-2);𝑅𝑆𝑅BF——基于林分的竹林地下生物量与地上生物量的比值,无量纲,见表A.12;𝐴𝐺𝐵BF,𝑇𝑏——竹林达到稳定成熟年龄时的单位面积地上生物量,单位为吨每公顷(td.m.·hm-2);𝑆𝐶BF,𝑡𝑏——第tb年时,竹林累计择伐地上生物量占𝐴𝐺𝐵BF,𝑇𝑏的比例(如株数比例),无量纲。A.5灌木林单位面积生物量的计算灌木林单位面积生物量计算方法,可按照下列方法的优先顺序进行选择。在项目实施阶段,每次监测和核算均须使用同一类方法,以保证项目减排量核算结果的可比性。A.5.1生物量方程法通过构建灌木林单株全株生物量、地上(或地上部各器官)或地下生物量与灌木测树因子(如基径、灌高、冠幅等)的相关方程,再结合单位面积灌木株数进行计算。在选择生物量方程时,优先选择全株生物量方程,其次选择地上和地下生物量方程(表A.16和表A.17)。𝐵Total,SF=∑∑𝑓SF,𝑗,𝑝(𝑥1,𝑥2,𝑥3……)×𝑁𝑗×10−3(A.23)𝑗𝑝式中:𝐵Total,SF——灌木林单位面积全林生物量,单位为吨每公顷(td.m.·hm-2);𝑓SF,𝑗,𝑝(𝑥1,𝑥2,𝑥3……)——灌木林中第j类灌木器官p的生物量与测树因子(x1,x2,x3……,如基径、灌高、冠幅等)的相关方程,单位为千克每株(kg·stem-1);𝑁𝑗——灌木林中第j类灌木的单位面积株数,单位为株每公顷(stem·hm-2);31j——灌木种类,j=1,2,3……,无量纲;p——灌木的器官,可分为叶、枝、茎、根等;也可分为地上与地下部分;也可是全株,无量纲;10-3——将千克转换为吨的常数。A.5.2缺省值法对于灌木造林,当灌木盖度<0.05时,灌木生物量可忽略不计,计为0。当灌木盖度≥0.05时,按照下列方式计算:𝐵Total,SF,𝑡=𝐴𝐺𝐵SF×𝐶𝐶SF,𝑡×(1+𝑅𝑆𝑅SF)(A.24)式中:𝐵Total,SF,𝑡——第t年时,灌木林单位面积全林生物量,单位为吨每公顷(td.m.·hm-2);𝐴𝐺𝐵SF——灌木林成熟稳定时的平均单位面积地上生物量,单位为吨每公顷(td.m.·hm-2),见表A.15;𝐶𝐶SF,𝑡——第t年时,灌木林盖度,用小数表示(例如盖度10%记为0.10),无量纲;𝑅𝑆𝑅SF——基于林分的灌木林地下生物量与地上生物量的比值,无量纲,见表A.15。A.6原有植被生物质碳储量变化量的计算项目业主须选择以下方式之一,确定项目原有植被生物质碳储量变化量:a)在项目设计阶段和项目实施阶段,项目业主可选择对每个项目碳层内原有的散生木(竹)进行逐一标记,记录编号和数量,并确保项目实施阶段不会人为移除原有散生木(竹)。项目业主仅对造林项目新增的林木(竹)进行监测,同时忽略原有植被的生物质碳储量变化量,即∆𝐶BiomassPE,𝑡=0;b)如果项目开始时项目业主未对原有散生木(竹)进行标记,则在项目首次监测时须补充标记原有散生木(竹)。项目业主须对新增林木进行监测,并计算项目各碳层的生物质碳储量变化量,再对原有植被的生物质碳储量变化量进行扣除。原有植被的生物质碳储量变化量按照公式(A.25)进行计算。∆𝐶BiomassPE,𝑡=∑(∆𝐶BiomassPROJ,𝑖,𝑡×𝐶𝐷PE,𝑖,𝑡0)(A.25)𝑖式中:∆𝐶BiomassPE,𝑡——项目原有散生木(竹)第t年的生物质碳储量变化量,单位为吨二氧化碳当量每年(tCO2e·a-1);∆𝐶BiomassPROJ,𝑖,𝑡——第i个项目碳层第t年的生物质碳储量变化量(包括原有植被和新增林木),单位为吨二氧化碳当量每年(tCO2e·a-1),𝐶𝐷PE,𝑖,𝑡0——采用本文件附录A计算;项目开始时第i个项目碳层原有散生木(竹)的平均冠层盖度,无量纲。须利用项目开始前的调查数据或高清遥感影像资料进行确定。以小数表示,例如盖度10%计为0.10,若无散生木(竹)则计为0;i——项目碳层,i=1,2,3……,无量纲。c)如果基准线和项目情景下会发生灌木的移除,可根据移除比例(如灌木盖度变化),采用公式(A.24)计算项目边界内原有灌木生物量的减少量。32A.7生物量模型与参数表A.1中国主要乔木林树种单株生物量与胸径或树高的相关方程(行业标准汇总)适用范生物一元模型二元模型树种围量应用条件模型参数R2模型参数R2油松湿地ababc松晋冀鲁AGBDBH≥5cm0.086112.461570.950.067772.180500.43610.95云南京津蒙AGBDBH<5cm0.429371.463290.398351.079940.4361松马尾辽豫甘BGBDBH≥5cm0.010932.664780.870.010902.661840.00460.87松青宁陕BGBDBH<5cm0.109311.233820.139011.079940.0046杉木AGBDBH≥5cm0.083892.440910.940.047442.103590.63110.97落叶闽赣湘AGBDBH<5cm0.309761.629280.290020.978680.6311松粤桂BGBDBH≥5cm0.043572.228770.910.035262.103590.23420.90云杉BGBDBH<5cm0.071441.921510.071991.660110.2342川滇藏AGBDBH≥5cm0.094922.356670.920.070232.103920.41120.95BGBDBH≥5cm0.016542.344900.780.014362.226370.19280.81AGBDBH≥5cm0.099492.408590.950.066622.093170.49760.96苏浙徽AGBDBH<5cm0.147692.16312闽粤赣桂湘贵BGBDBH≥5cm0.008112.695050.890.117271.741790.4976BGBDBH<5cm0.033841.807540.008832.73828-0.08030.880.043671.74485-0.0803AGBDBH≥5cm0.137922.343590.910.092352.028170.49760.92鄂川AGBDBH<5cm0.189192.147210.181671.607780.4976BGBDBH≥5cm0.011252.630050.580.012242.67327-0.80260.75BGBDBH<5cm0.043351.791630.067661.61084-0.8026AGBDBH≥5cm0.076372.403930.960.065392.017350.49430.97湘鄂粤AGBDBH<5cm0.212771.76730桂苏浙川徽贵BGBDBH≥5cm0.015502.444210.920.190711.352260.4943BGBDBH<5cm0.032991.975050.016392.52941-0.11740.920.346552.06398-0.1174AGBDBH≥5cm0.043632.545890.980.032722.110930.60210.98赣闽AGBDBH<5cm0.208481.574050.173821.073220.6021BGBDBH≥5cm0.008862.586170.760.008202.62298-0.00960.84BGBDBH<5cm0.032321.781790.315861.78495-0.0096AGBDBH≥5cm0.112702.395820.960.068482.015490.59150.97黑吉AGBDBH<5cm0.182542.09620辽东、部蒙BGBDBH≥5cm0.042582.370530.940.145831.545810.5915BGBDBH<5cm0.016712.951760.044412.40255-0.04980.940.016323.02446-0.0498AGBDBH≥5cm0.073022.472980.920.062332.015490.59150.92晋冀、AGBDBH<5cm0.142142.05910蒙中部西BGBDBH≥5cm0.028292.364030.900.170511.390240.5915BGBDBH<5cm0.022752.499380.028672.40255-0.04980.900.039192.20824-0.0498AGBDBH≥5cm0.117952.336120.920.058512.015490.59150.96新AGBDBH<5cm0.142362.219230.122911.554320.5915BGBDBH≥5cm0.020552.375560.700.021802.40255-0.04980.70BGBDBH<5cm0.078521.542590.053251.84756-0.0498AGBDBH≥5cm0.111612.328030.900.055772.015490.59150.95川滇藏AGBDBH<5cm0.163412.091180.156781.373320.5915BGBDBH≥5cm0.021362.376230.710.022652.40255-0.04980.71BGBDBH<5cm0.036632.041270.031452.19867-0.4980AGBDBH≥5cm0.093072.432150.970.080702.259570.25660.98黑吉AGBDBH<5cm0.136952.192110.129761.964480.2566BGBDBH≥5cm0.020722.474310.810.025132.70802-0.34750.80BGBDBH<5cm0.775341.654340.087511.93277-0.3475甘青、AGBDBH≥5cm0.148652.289060.910.128902.098280.25660.91新疆天AGBDBH<5cm0.317961.816640.307981.557100.2566山BGBDBH≥5cm0.046322.288360.810.056182.54672-0.34750.8533适用范生物一元模型二元模型树种围量应用条件模型参数R2模型参数R2冷杉ababc柳杉栎树BGBDBH<5cm0.097031.828930.094152.22593-0.3475AGBDBH≥5cm0.155592.258770.830.134922.073310.25660.89新疆阿AGBDBH<5cm0.170812.200790.174921.912010.2566尔泰山BGBDBH≥5cm0.034082.270610.840.041342.52175-0.34750.85BGBDBH<5cm0.066221.857950.066902.22700-0.3475AGBDBH≥5cm0.129612.309610.890.112392.122570.25660.89川AGBDBH<5cm0.308801.770200.310041.492080.2566BGBDBH≥5cm0.026792.317710.970.032492.57101-0.34750.97BGBDBH<5cm0.692861.727320.059142.19889-0.3475AGBDBH≥5cm0.105542.389260.960.091522.210600.25660.98滇AGBDBH<5cm0.187462.032300.169231.828660.2566BGBDBH≥5cm0.024132.417100.940.029272.65904-0.34750.93BGBDBH<5cm0.030622.269070.038122.49485-0.3475AGBDBH≥5cm0.181432.218040.910.157322.036230.25660.92藏AGBDBH<5cm0.136922.392930.126912.169730.2566BGBDBH≥5cm0.031442.238460.850.038132.48467-0.34750.80BGBDBH<5cm0.026422.346460.032782.57872-0.3475AGBDBH≥5cm0.090882.417620.950.069452.057530.50840.96黑吉AGBDBH<5cm0.140352.147610.125451.690120.5084BGBDBH≥5cm0.026082.301980.810.031182.54141-0.33800.82BGBDBH<5cm0.022302.399040.025162.67482-0.3380AGBDBH≥5cm0.097512.428780.890.074512.057530.50840.95甘青、AGBDBH<5cm0.131282.24404新疆山天BGBDBH≥5cm0.038632.294570.890.102791.857620.5084BGBDBH<5cm0.013602.943040.046202.54141-0.33800.900.015703.21188-0.3380AGBDBH≥5cm0.084862.409850.880.064842.057530.50840.93川AGBDBH<5cm0.168961.981950.157651.505520.5084BGBDBH≥5cm0.024072.307150.850.028792.54141-0.33800.83BGBDBH<5cm0.025192.278930.027222.57625-0.3380AGBDBH≥5cm0.080182.410490.990.061272.057530.50840.99滇AGBDBH<5cm0.181751.902050.194061.341220.5084BGBDBH≥5cm0.029422.306730.870.035182.54141-0.33800.90BGBDBH<5cm0.050961.965300.052572.29188-0.3380AGBDBH≥5cm0.081162.424110.940.062022.057530.50840.95藏AGBDBH<5cm0.103662.272090.091591.815260.5084BGBDBH≥5cm0.030352.297670.940.036302.54141-0.33800.94BGBDBH<5cm0.010412.962710.012633.19747-0.3380AGBDBH≥5cm0.154832.171000.940.093111.811740.60680.97浙鄂湘AGBDBH<5cm0.227541.931800.197801.343550.6068川渝BGBDBH≥5cm0.014062.571930.870.016792.69756-0.21220.86BGBDBH<5cm0.062791.931800.058841.91838-0.2122AGBDBH≥5cm0.091352.489540.920.061492.143800.58390.95黑吉AGBDBH<5cm0.119632.32194辽东、部蒙BGBDBH≥5cm0.045882.300790.810.109181.787050.5839BGBDBH<5cm0.086461.907050.051832.40723-0.17970.800.293521.32981-0.1797AGBDBH≥5cm0.093932.546080.940.075092.326370.33020.95豫晋冀AGBDBH<5cm0.204842.06167陕青宁甘BGBDBH≥5cm0.052092.301300.890.154191.879370.3302BGBDBH<5cm0.127301.746120.069892.58980-0.43350.880.177022.01240-0.4335AGBDBH≥5cm0.115202.424240.910.078062.063210.57390.92川滇AGBDBH<5cm0.298131.833420.229991.391830.5739BGBDBH≥5cm0.048902.207300.770.055612.32664-0.18970.78BGBDBH<5cm0.140671.550770.156211.68493-0.189734适用范生物一元模型二元模型树种围量应用条件模型参数R2模型参数R2桦树ababc木荷AGBDBH≥5cm0.213602.304160.910.131881.828920.71120.97枫香鄂湘赣AGBDBH<5cm0.225862.26960徽浙贵渝BGBDBH≥5cm0.110602.057300.790.154441.730820.7112柏木BGBDBH<5cm0.062712.409820.114602.09235-0.05250.75高山松0.053442.56640-0.0525思茅黑吉、AGBDBH≥5cm0.102982.440220.940.068072.108500.52020.96松蒙东部AGBDBH<5cm0.143052.236030.104081.844700.5202樟子松(白BGBDBH≥5cm0.055112.254640.810.044092.075910.28030.79桦)BGBDBH<5cm0.056122.243340.047292.032420.2803黑吉、AGBDBH≥5cm0.095882.425640.920.063382.108500.52020.94蒙东部AGBDBH<5cm0.138632.196530.110291.764260.5202(其他BGBDBH≥5cm0.049162.246780.710.039332.075910.28030.71桦)BGBDBH<5cm0.059982.123250.054091.878010.2803蒙中西AGBDBH≥5cm0.111462.429830.940.073672.108500.52020.95部、豫AGBDBH<5cm0.306461.801360.202541.480130.5202晋陕青BGBDBH≥5cm0.050842.249040.860.040682.075910.28030.85甘宁新BGBDBH<5cm0.088751.902900.075661.690310.2803AGBDBH≥5cm0.096152.418610.880.063562.108500.52020.90川滇AGBDBH<5cm0.130972.226600.089071.898780.5202BGBDBH≥5cm0.041592.243000.820.033282.075910.2803BGBDBH<5cm0.020042.696730.017232.485040.28030.82滇闽赣AGBDBH≥5cm0.176852.263140.960.120452.064460.38270.96贵浙粤湘BGBDBH≥5cm0.064082.197840.810.081772.32395-0.24290.82桂闽赣AGBDBH≥5cm0.106152.466500.950.089092.255640.30410.96湘徽贵渝浙BGBDBH≥5cm0.095522.141900.920.120522.42178-0.40370.90AGBDBH≥5cm0.153412.316960.920.124831.914890.61520.95京津冀AGBDBH<5cm0.362741.782240.270861.433570.6152晋蒙豫鲁BGBDBH≥5cm0.059612.166570.880.043811.767500.64280.900.095071.286180.6428BGBDBH<5cm0.140961.631850.089471.914890.61520.90AGBDBH≥5cm0.133132.253590.850.338431.088250.6152陕甘青AGBDBH<5cm0.483311.452480.033701.767500.64280.86BGBDBH≥5cm0.054872.103200.810.127490.940850.6428BGBDBH<5cm0.199211.302090.094331.914890.61520.98AGBDBH≥5cm0.141792.329280.890.160901.583140.6152川滇藏AGBDBH<5cm0.239072.004720.024531.767500.64280.93BGBDBH≥5cm0.039172.178900.770.041841.435740.6428BGBDBH<5cm0.066031.854340.100141.914890.61520.95AGBDBH≥5cm0.147342.345120.930.120881.797910.6152浙鄂湘AGBDBH<5cm0.191172.183320.035181.767500.64280.91粤渝贵BGBDBH≥5cm0.056172.194730.870.042471.650520.6428BGBDBH<5cm0.072882.032930.089092.255640.30410.94AGBDBH≥5cm0.103872.371220.900.223061.623700.2747藏川滇AGBDBH<5cm0.244171.840150.015862.191830.32570.79BGBDBH≥5cm0.021742.371220.770.039581.623700.32570.79BGBDBH<5cm0.050241.840150.770.021062.523980.41630.97普洱、AGBDBH≥5cm0.027422.803630.960.103971.531800.41630.97西双版AGBDBH<5cm0.111441.932470.960.004022.634450.22520.93纳、德BGBDBH≥5cm0.006362.687310.94宏、临0.016241.767610.22520.93沧BGBDBH<5cm0.013792.206060.94蒙东AGBDBH≥5cm0.075992.425390.940.054602.234120.34780.94部、黑AGBDBH<5cm0.075992.425390.940.175771.507700.34780.94吉辽BGBDBH≥5cm0.012812.506590.780.011312.57232-0.03980.7635适用范生物一元模型二元模型树种围量应用条件模型参数R2模型参数R2椴树榆树ababc黄山BGBDBH<5cm0.081051.360510.780.033431.89875-0.03980.76松黑吉辽AGBDBH≥5cm0.071112.458530.930.037982.128250.61120.96杨树BGBDBH≥5cm0.021582.515090.860.016712.70249-0.13120.86注黑吉辽AGBDBH≥5cm0.185272.175220.920.102661.985200.49070.95注冀蒙BGBDBH≥5cm0.071872.140110.870.043232.060560.33720.91浙徽闽AGBDBH≥5cm0.094212.416670.930.059522.121830.52380.95赣AGBDBH<5cm0.182972.004250.165911.484910.2524黑吉辽AGBDBH≥5cm0.073572.496460.940.045022.122230.59160.96蒙AGBDBH<5cm0.133072.128260.084281.732680.5916(天然BGBDBH≥5cm0.024072.362070.920.022432.308200.08520.92山杨)BGBDBH<5cm0.041522.023440.057761.720500.0852陕甘青AGBDBH≥5cm0.075292.511650.950.046072.148920.59160.95宁新AGBDBH<5cm0.159112.046690.091071.725500.5916(天然BGBDBH≥5cm0.025122.387110.900.023412.334890.08520.90山杨)BGBDBH<5cm0.095391.558130.075081.610700.0852AGBDBH≥5cm0.095492.428420.930.058432.051860.59160.96川滇藏AGBDBH<5cm0.099772.401180.065301.982760.5916(山天杨然)BGBDBH≥5cm0.026682.292040.870.024862.237830.08520.87BGBDBH<5cm0.022172.407040.022082.311390.0852黑吉辽AGBDBH≥5cm0.090842.372090.920.055592.008610.59160.94蒙AGBDBH<5cm0.126382.166930.095851.670050.5916(人工BGBDBH≥5cm0.033892.246900.670.031572.194580.08520.66杨树)BGBDBH<5cm0.032012.282210.039612.053640.0852鲁豫蒙AGBDBH≥5cm0.076852.507310.920.047032.124870.59160.94冀晋京AGBDBH<5cm0.101822.332510.067411.901150.5916(人工BGBDBH≥5cm0.026062.365900.910.024282.310850.08520.91杨树)BGBDBH<5cm0.035392.175730.036692.054210.0852苏徽鄂AGBDBH≥5cm0.047922.673460.970.029332.276340.59160.98湘渝贵AGBDBH<5cm0.138162.015590.098371.524350.5916(人工BGBDBH≥5cm0.012322.519480.940.011482.462320.08520.94杨树)BGBDBH<5cm0.058511.551300.059661.438200.08521:一元模型方程的表达式为𝑀=𝑎⋅𝐷𝐵𝐻𝑏,二元模型方程的表达式为𝑀=𝑎⋅𝐷𝐵𝐻𝑏⋅𝐻𝑐。其中,M为单株生物量,单位为千克(kgd.m.);AGB为地上生物量,单位为千克(kgd.m.);BGB为地下生物量,单位为千克(kgd.m.);DBH为胸径,单位为厘米(cm);H为树高,单位为米(m);a、b、c为模型参数;R2为决定系数。2:在造林密度和立地条件等相对一致的情况下,应优先选择基于胸径的一元模型,否则应选择基于胸径和树高的二元模型。本表中未包含的树种,可参考表A.2和表A.3中的方程。本表未包含的树种,可采用文献值。表A.2中国主要乔木林树种单株生物量与胸径的相关方程(文献混合建模)树种(组)器官abR2胸径范围0.11122.36890.9261.0~95.0地上0.15332.33770.9171.0~95.0针叶树0.06222.52890.9331.0~150.00.02772.75180.9441.0~150.0整株0.28482.13440.9621.0~95.00.34922.11030.9631.0~95.0地上0.13802.28850.9801.0~80.0阔叶树0.03802.73300.9961.0~50.00.00763.28950.9642.0~40.0整株0.01303.19950.9772.0~40.036地上云冷杉整株地上红松整株地上栲类整株树种(组)器官abR2胸径范围桉树地上0.18822.19160.9532.0~19.7整株0.18982.24070.9372.0~19.7橡胶地上0.17842.31970.9842.0~38.0整株0.21132.31170.9912.0~38.0青冈地上0.22352.23110.9773.2~37.5整株0.19302.35900.9983.2~37.5注:方程表达式为𝑌=𝑎⋅𝐷𝐵𝐻𝑏。其中,Y为单株生物量,单位为千克(kgd.m.);DBH为胸径,单位为厘米(cm);a、b为模型参数;R2为决定系数。表A.3中国主要乔木林树种单株生物量与胸径和树高的相关方程(文献混合建模)树种(组)器官abR2胸径范围树高范围针叶树地上32.63350.94720.9331.0~80.00.7~36.00.92120.9231.0~80.00.7~36.0整株51.5232阔叶树地上119.16320.85420.9211.0~150.00.7~66.80.87930.9311.0~150.00.7~66.8整株102.1363云冷杉地上227.36270.76600.8531.0~46.02.3~33.40.67030.7931.0~46.02.3~23.3整株707.3147红松地上164.16890.79980.9711.0~50.01.5~25.5整株141.86790.84270.9851.0~50.01.5~25.5水杉地上131.98140.75890.9672.0~32.05.8~22.7整株191.42210.73850.9642.0~32.05.8~22.7侧柏地上270.77370.73260.8612.0~20.02.4~11.5整株72.50150.93340.9852.0~20.02.4~9.8桉树地上10.68701.12190.8842.0~23.04.1~29.0整株11.60781.13560.8582.0~23.04.1~29.0刺槐地上56.55560.92820.9762.0~34.02.1~18.0整株103.80420.87640.9692.0~34.02.1~18.0金合欢地上321.25320.71950.8752.0~26.03.4~24.00.86340.9932.0~26.03.4~24.0整株114.6927樟树地上5.76801.18460.9852.0~26.02.8~19.1整株7.43091.18270.9852.0~26.02.8~19.1注:方程表达式为𝑌=𝑎⋅(𝐷𝐵𝐻2⋅𝐻)𝑏⋅10−3。其中,Y为单株生物量,单位为千克(kgd.m.);DBH为胸径,单位为厘米(cm);H为树高,单位为米(m);a、b为模型参数;R2为决定系数。表A.4中国主要乔木林树种单株材积与胸径或树高的相关方程(行业标准汇总)树适用应用条件一元模型二元模型模型参数R2模型参数R2种范围ababc晋冀鲁DBH≥5cm0.128672.450510.96590.085471.970570.744670.9876油京津蒙松辽豫甘DBH<5cm青宁陕0.319781.884860.96590.281351.230260.744670.9876湿闽赣湘DBH≥5cm0.163552.378330.94570.077881.939300.821360.9791地粤桂DBH<5cm松0.287092.028720.94570.263511.181950.821360.9791川滇藏DBH≥5cm0.172512.375940.88000.096531.888860.792420.973337树适用应用条件一元模型二元模型模型参数R2模型参数R2种范围ababc云南松DBH<5cm0.172512.375940.88000.096531.888860.792420.9733苏浙徽DBH≥5cm0.138812.484920.95520.066941.911400.904850.9854马闽赣湘尾粤桂贵DBH<5cm0.214332.214990.95520.140911.448880.904850.9854松鄂川DBH≥5cm0.146442.484920.94480.070631.911400.904850.9831DBH<5cm0.181422.351840.94480.168511.371000.904850.9831湘鄂粤DBH≥5cm0.093012.618380.94130.067991.837980.997760.9934杉桂浙徽DBH<5cm木苏川贵0.231492.051850.94130.185601.213980.997760.9934赣闽DBH≥5cm0.092262.603140.95540.055611.837981.059220.9934DBH<5cm0.209962.092130.95540.152491.211181.059920.9934黑吉DBH≥5cm0.176942.435830.96170.076691.800350.992460.9878辽、蒙DBH<5cm0.166542.476160.96170.114261.552620.994260.9878东部落晋冀、DBH≥5cm0.094772.568010.95740.072671.800350.992460.9868叶蒙中西DBH<5cm松部0.108442.484280.95740.147151.361940.994260.9868新DBH≥5cm0.247732.338350.88570.076401.800350.992460.9669DBH<5cm0.208642.445050.88570.163051.329360.994260.9669川滇藏DBH≥5cm0.240672.324780.88080.075141.800350.992460.9760DBH<5cm0.255512.287620.88080.238361.083070.994260.9760黑吉DBH≥5cm0.126772.489370.93940.070891.785621.046510.9876DBH<5cm0.164502.327500.93940.132021.399241.046510.9876甘青、DBH≥5cm0.118012.563620.86160.065991.785621.046510.9903新疆天DBH<5cm0.263852.063670.86160.231691.005261.046510.9903山新疆阿DBH≥5cm0.117472.541900.79190.065691.785621.046510.9680云尔泰山DBH<5cm0.204712.196800.79190.225531.109201.046510.9680杉DBH≥5cm0.121232.548380.83130.067791.785621.046510.9576川DBH<5cm0.292522.001090.83130.297380.866951.046510.9576滇DBH≥5cm0.126032.514190.85040.070471.785621.046510.9847DBH<5cm0.266082.049850.85040.175301.219411.046510.9847藏DBH≥5cm0.117732.527040.80660.065831.785621.046510.9807DBH<5cm0.222632.131160.80660.163351.220961.046510.9807黑吉DBH≥5cm0.108712.551150.96500.062651.813411.041580.9895DBH<5cm0.209662.143040.96500.166591.205741.041580.9895甘青、DBH≥5cm0.113902.540100.86230.065641.813411.041580.9818新疆天DBH<5cm山0.173542.312370.86230.105141.520681.041580.9818冷DBH≥5cm0.109802.535240.93390.063281.813411.041580.9842杉川DBH<5cm0.225582.087850.93390.195741.111751.041580.9842滇DBH≥5cm0.119562.536540.90810.06891.813411.041580.9647DBH<5cm0.200972.213880.90810.229851.064861.041580.9647藏DBH≥5cm0.114212.564440.93840.065821.813411.041580.9794DBH<5cm0.175142.298770.93840.135921.362811.041580.9794柳浙鄂湘DBH≥5cm0.155982.431250.94940.066501.829091.017030.9902杉川渝DBH<5cm0.209442.248160.94940.165611.262171.017030.9902黑吉DBH≥5cm0.119832.448280.91000.064051.901180.923940.9780栎辽、蒙DBH<5cm树东部0.120702.443780.91000.104431.597390.923940.9780DBH≥5cm0.120072.516080.92890.064181.901180.923940.974138树适用应用条件一元模型二元模型模型参数R2模型参数R2种范围ababc豫晋冀陕青甘DBH<5cm宁0.421401.736010.92890.190291.225840.923940.9741川滇DBH≥5cm0.124662.482380.91720.066631.901180.923940.9671DBH<5cm0.196492.199670.91720.129391.488790.923940.9172鄂湘赣DBH≥5cm0.128812.518600.92590.068841.901180.923940.9837徽浙渝DBH<5cm0.180932.307460.92590.110421.607640.923940.9837贵黑吉、DBH≥5cm0.164152.406880.94820.080771.838640.891080.9797蒙东部(白DBH<5cm桦)0.195002.299870.94820.113091.629520.891080.9797黑吉、DBH≥5cm0.164152.381900.88750.080771.838640.891080.9444蒙东部桦(其他DBH<5cm树桦)0.315781.975370.88750.213431.234900.891080.9444蒙中西DBH≥5cm0.164152.389090.86170.080771.838640.891080.9283部、豫晋陕青DBH<5cm甘宁新0.338941.938600.86170.166731.388330.891080.9283川滇DBH≥5cm0.164152.369860.94620.080771.838640.891080.9676DBH<5cm0.161592.379630.94620.083491.818090.891080.9676木滇闽赣荷贵浙湘DBH≥5cm粤0.185092.341760.95420.092501.982900.691150.9814枫桂闽赣香湘贵浙DBH≥5cm徽渝0.128992.490440.96440.077521.877440.884180.9922京津冀DBH≥5cm0.105522.431350.95130.080770.072121.688930.9838晋蒙鲁DBH<5cm0.324641.733100.95130.113090.189311.089290.9838豫柏陕甘青DBH≥5cm0.151312.314330.86470.080770.072641.688930.9772DBH<5cm0.261011.975560.86470.213430.135181.303000.9772木川滇藏DBH≥5cm0.139762.454110.79160.080770.065851.688930.9699DBH<5cm0.187662.270990.79160.166730.090331.492540.9699浙鄂湘DBH≥5cm0.136612.483350.93910.080770.066951.688930.9731粤渝贵DBH<5cm0.187582.286320.93910.083490.080471.574670.9731高DBH≥5cm0.127972.454100.88140.079891.890950.862270.9848山藏川滇DBH<5cm松0.282321.962470.88140.212551.282970.862270.9848普洱、DBH≥5cm0.099162.577130.96560.056481.981280.887090.9890思西双版茅松纳宏、、德临DBH<5cm沧0.222192.075850.96560.191671.222130.887090.9890樟蒙东DBH≥5cm0.175232.357850.88350.075711.872230.882880.9759子松部吉、辽黑DBH<5cm0.206292.256450.88350.191511.295590.882880.9759椴黑吉辽DBH≥5cm树0.178602.372430.94030.077251.931030.816810.9840榆黑吉辽DBH≥5cm树冀蒙0.259772.131010.87190.098381.818490.807120.9494DBH≥5cm0.146872.455810.92210.067991.961270.878510.980439树适用应用条件一元模型二元模型模型参数R2模型参数R2种范围ababc黄浙徽闽山赣DBH<5cm松0.222782.196930.92210.189051.325800.878510.9804黑吉辽DBH≥5cm0.134472.494610.061241.895210.94759蒙天然DBH<5cm0.176982.323950.94620.085151.690370.947590.9830山杨陕甘青DBH≥5cm0.134272.476180.061141.895210.94759宁新(天然DBH<5cm山杨)0.218502.173630.94370.089401.659190.947590.9848川滇藏DBH≥5cm0.135602.498340.061751.895210.94759(天然DBH<5cm0.150052.435460.94470.076111.765300.947590.9848山杨)杨黑吉辽DBH≥5cm0.127632.477390.058121.895210.94759树蒙(人工杨DBH<5cm树)0.199932.198510.93260.128401.402680.947590.9700鲁豫蒙DBH≥5cm0.122962.507740.055991.895210.94759冀晋京(人工DBH<5cm杨树)0.209462.176770.96240.108201.485890.947590.9889苏徽鄂DBH≥5cm0.119912.531250.054601.895210.94759湘渝贵(人工DBH<5cm杨树)0.179002.282290.94890.103901.495470.947590.9849注1:一元模型方程的表达式为𝑉=𝑎⋅𝐷𝐵𝐻𝑏⋅10−3,二元模型方程的表达式为𝑉=𝑎⋅𝐷𝐵𝐻𝑏⋅𝐻𝑐⋅10−3。其中,V为单株材积,单位为立方米(m3);DBH为胸径,单位为厘米(cm);H为树高,单位为米(m);a、b、c为模型参数;R2为决定系数。注2:在造林密度和立地条件等相对一致的情况下,应优先选择基于胸径的一元模型,否则应选择基于胸径和树高的二元模型。表A.5中国主要乔木树种林分水平全林(或地上)生物量方程与生物量转换与扩展因子森林类型abcR2冷杉林4.10950.59760.81810.871云杉林1.38790.71680.79910.9130.69860.82620.77940.936落叶松林0.57430.71200.81090.936杉木林0.45451.53410.80070.947柏木林1.42820.97680.82970.9530.21121.12350.80980.966马尾松林0.91580.65010.84530.973油松林1.07210.83030.80410.9460.71961.29480.79300.917云南松林0.75071.01180.78210.898其他针叶林0.29450.89780.82460.9340.38641.54990.78330.933栎类林0.33301.17400.77930.962桦树林0.24010.86790.79810.989杨树林刺槐林桉树林橡胶林40其他软阔类0.34541.21300.78800.933其他硬阔类0.65340.99200.79540.897针叶混交林3.36120.83480.80880.919针阔混交林2.25750.95600.79610.892阔叶混交林1.39441.10860.79380.916注1:全林生物量的方程表达式为𝐵Total=𝑎+𝑏∙𝑉,地上生物量的方程表达式为𝐴𝐺𝐵=𝑐⋅(𝑎+𝑏∙𝑉)。其中,𝐵Total为单位面积全林生物量,单位为吨每公顷(td.m.·hm-2);𝐴𝐺𝐵为单位面积地上生物量,单位为吨每公顷(td.m.·hm-2)。注2:全林生物量转换与扩展因子的方程表达式为𝐵𝐶𝐸𝐹Total=𝑎⁄𝑉+𝑏,地上生物量转换与扩展因子的方程表达式为𝐵𝐶𝐸𝐹AGB=𝑐⋅(𝑎⁄𝑉+𝑏)。其中,𝐵𝐶𝐸𝐹Total为全林生物量转换与扩展因子,单位为吨每立方米(td.m.·m-3);𝐵𝐶𝐸𝐹AGB为地上生物量转换与扩展因子,单位为吨每立方米(td.m.·m-3)。𝑉为单位面积蓄积量,单位为立方米每公顷(m3·hm-2);a、b、c为模型参数;R2为决定系数。表A.6中国主要乔木林树种的地上生物量转换与扩展因子(行业标准汇总)一元模型二元模型模型参数树种适用范围应用条件模型参数baba0.20994c0.79289-0.15032-0.30857油松晋冀鲁京津蒙DBH≥5cm0.669250.011061.415870.16429-0.30857湿地松辽豫甘青宁陕DBH<5cm0.60913-0.20327-0.19028云南松DBH≥5cm1.34270-0.421571.100580.21506-0.19028马尾松闽赣湘粤桂DBH<5cm0.727540.18177-0.38122DBH≥5cm0.512940.062580.995190.29290-0.40722杉木川滇藏DBH≥5cm0.832200.11677-0.40722苏浙徽闽赣湘DBH<5cm1.07895-0.399441.307740.23678-0.40722落叶松DBH≥5cm1.078050.17937-0.40722粤桂贵DBH<5cm0.55024-0.019270.961760.13827-0.50351DBH≥5cm1.027530.27295-0.50351鄂川DBH<5cm0.71673-0.076330.58839-0.13795-0.45710DBH≥5cm1.139910.21510-0.45710湘鄂粤桂浙徽DBH<5cm0.68906-0.051870.89292-0.00682-0.40100苏川贵DBH≥5cm1.276300.21514-0.40100DBH<5cm0.94183-0.141340.857790.02830-0.40100赣闽DBH≥5cm1.158710.21514-0.40100DBH<5cm1.04284-0.204640.765840.22496-0.40100黑吉辽、蒙东DBH≥5cm0.753830.21514-0.40100部DBH<5cm0.82107-0.214450.742260.29025-0.40100DBH≥5cm0.657740.47395-0.40100晋冀、蒙中西DBH<5cm0.91913-0.284551.138430.56524-0.78988部DBH≥5cm0.982870.31266-0.78988DBH<5cm0.47291-0.057251.953420.55184-0.78988新DBH≥5cm1.329270.28769-0.78988DBH<5cm0.99295-0.518152.053970.89281-0.78988DBH≥5cm0.775590.33695-0.78988DBH<5cm0.63696-0.042711.657890.62514-0.78988DBH≥5cm1.04260-0.78988DBH<5cm1.09608-0.379970.77052-0.095031.31083-0.425700.47612-0.002240.68235-0.22583川滇藏0.463760.003250.63954-0.19644黑吉0.73413-0.05722甘青、新疆天0.83254-0.13538山1.25968-0.27456新疆阿尔泰山云杉1.20508-0.247031.32452-0.283130.834400.00399川1.06911-0.238771.05564-0.2308941一元模型二元模型模型参数树种适用范围应用条件模型参数baba0.42498c1.298660.60925-0.78988滇DBH≥5cm0.83745-0.124930.965360.25061-0.78988DBH<5cm2.389780.94877-0.78988藏DBH≥5cm0.70453-0.017540.776900.24412-0.78988DBH<5cm1.108530.48438-0.53319黑吉DBH≥5cm1.54107-0.309000.753040.24412-0.53319DBH<5cm1.135180.33694-0.53319甘青、新疆天DBH≥5cm0.614990.261770.977660.24412-0.53319山DBH<5cm1.024760.39377-0.53319DBH≥5cm0.83602-0.133530.805430.24412-0.53319川DBH<5cm0.889260.37636-0.53319DBH≥5cm0.669400.004570.844300.24412-0.53319滇DBH<5cm0.942280.45244-0.53319DBH≥5cm0.85612-0.145230.67386-0.01734-0.53319藏DBH<5cm1.400020.08138-0.41026DBH≥5cm0.75645-0.068331.194350.24261-0.41026浙鄂湘川渝DBH<5cm0.960000.18966-0.34004冷杉DBH≥5cm0.77285-0.125391.045400.42518-0.34004柳杉黑吉辽、蒙东DBH<5cm1.170120.65353-0.59379栎树部DBH≥5cm0.74899-0.105900.810260.16203-0.59379DBH<5cm1.17159-0.09695-0.35001桦树豫晋冀陕青甘DBH≥5cm0.67066-0.126051.77745-0.07227-0.35001木荷宁DBH<5cm1.915680.12318-0.21275枫香DBH≥5cm0.90439-0.311831.398660.26986-0.21275柏木川滇DBH<5cm0.843790.21518-0.37089DBH≥5cm0.71064-0.140340.920320.26986-0.37089鄂湘赣徽浙渝DBH<5cm0.784620.52936-0.37089贵DBH≥5cm0.59184-0.026670.516740.26986-0.37089DBH<5cm0.912110.09180-0.37089黑吉、蒙东部DBH≥5cm0.99262-0.260251.214800.26986-0.37089(白桦)DBH<5cm0.786880.08069-0.37089DBH≥5cm1.08643-0.316361.06992-0.37089黑吉、蒙东部DBH<5cm(其他桦)0.762330.041261.30222DBH≥5cm蒙中西部、豫0.99116-0.12184晋陕青甘宁新DBH≥5cm0.782320.03000川滇DBH≥5cm滇闽赣贵浙湘DBH<5cm0.486100.32566DBH≥5cm粤DBH<5cm0.92409-0.05814桂闽赣湘贵浙DBH≥5cmDBH<5cm1.51730-0.36625徽渝DBH≥5cm京津冀晋蒙鲁DBH<5cm1.65835-0.21443豫1.24832-0.03796陕甘青0.627380.03334川滇藏0.73359-0.06384浙鄂湘粤渝贵0.584070.043740.438990.221160.678980.040750.90420-0.137240.587670.048750.81051-0.153040.95545-0.078620.08156-0.308510.82290-0.023941.149260.37820-0.580041.45379-0.114391.730940.22596-0.520741.117380.049141.430790.34429-0.520740.87984-0.060751.231750.22596-0.520741.85167-0.523082.50363-0.21475-0.520741.01457-0.124821.432570.22596-0.520741.27394-0.266271.781270.09060-0.520741.07857-0.138231.495670.22596-0.520741.01913-0.103011.502230.22324-0.5207442一元模型二元模型树种适用范围应用条件模型参数模型参数ababc高山松藏川滇DBH≥5cm0.81169-0.082891.118970.30088-0.58760DBH<5cm0.86488-0.122321.049450.34073-0.58760思茅松普洱、西双版DBH≥5cm0.276560.226500.372810.54270-0.47075纳、德宏、临DBH<5cm0.50156-0.143380.542470.30967-0.47075沧樟子松蒙东部、黑吉DBH≥5cm0.433670.067550.721230.36189-0.53513辽DBH<5cm0.87738-0.370290.917830.21211-0.53513椴树黑吉辽DBH≥5cm0.398150.086100.491680.19723-0.20564榆树黑吉辽冀蒙DBH≥5cm0.713220.044211.043560.16671-0.31638黄山松浙徽闽赣DBH≥5cm0.64147-0.039140.875510.16057-0.35477DBH<5cm0.82127-0.192670.877560.15911-0.35477黑吉辽蒙DBH≥5cm0.547130.001850.735230.22702-0.35596(天然山杨)DBH<5cm0.75192-0.195690.989760.04231-0.35596陕甘青宁新DBH≥5cm0.560700.035470.753470.25371-0.35596(天然山杨)DBH<5cm0.72820-0.126941.018710.06631-0.35596杨树川滇藏DBH≥5cm0.70416-0.069910.946250.15665-0.35596(天然山杨)DBH<5cm0.66491-0.034280.858020.21747-0.35596DBH≥5cm0.71175-0.105300.956450.11340-0.35596黑吉辽蒙DBH<5cm0.63212-0.031580.746510.26737-0.35596(人工杨树)鲁豫蒙冀晋京DBH≥5cm0.62504-0.000430.839920.22966-0.35596(人工杨树)DBH<5cm0.486130.155740.623030.41527-0.35596苏徽鄂湘渝贵DBH≥5cm0.399680.142210.537090.38113-0.35596(人工杨树)DBH<5cm0.77182-0.266690.946810.02887-0.35596注1:一元模型方程的表达式为𝐵𝐶𝐸𝐹=𝑎⋅𝐷𝐵𝐻𝑏,二元模型方程的表达式为𝐵𝐶𝐸𝐹=𝑎⋅𝐷𝐵𝐻𝑏⋅𝐻𝑐。其中,BCEF为地上生物量转换与扩展因子,是地上生物量与树干蓄积量的比值,单位为吨每立方米(td.m.·m-3);DBH为胸径,单位为厘米(cm);H为树高,单位为米(m);a、b、c为模型参数。注2:在造林密度和立地条件等相对一致的情况下,应优先选择基于胸径的一元模型,否则应选择基于胸径和树高的二元模型。本表未包含的树种,可采用文献值。表A.7中国主要乔木林树种的生物量扩展因子(文献整合分析)森林类型应用条件:公顷蓄积≤100m3·hm-2应用条件:公顷蓄积>100m3·hm-2BEFBEF云冷杉林1.82751.4048落叶松林1.45111.2224温性针叶林1.81001.3405油松林1.83591.4081马尾松林1.55651.2063暖性针叶林1.71191.39711.90851.2875杉类1.70291.3593柏木林1.36941.26931.38891.2416栎类1.56701.3104桦木林其它硬阔类43杨树林1.55581.4184桉树林1.24131.1266其它软阔类1.47191.3335针叶混1.61661.3033阔叶混1.40421.3587针阔混1.67131.3725注:BEF为生物量扩展因子,是地上生物量与树干生物量的比值,无量纲。表A.8中国主要乔木林树种的地下生物量与地上生物量的比值(行业标准汇总)树种适用范围应用条件一元模型二元模型模型参数模型参数DBH≥5cmDBH<5cmababcDBH≥5cm0.160820.48134-0.43154油松晋冀鲁京津蒙DBH<5cm0.126880.203210.348970.00000-0.43154湿地松辽豫甘青宁陕DBH≥5cm0.743310.00000-0.39689云南松DBH≥5cm0.25457-0.229470.248240.68143-0.39689马尾松闽赣湘粤桂DBH<5cm0.204510.12245-0.21837DBH≥5cm0.51937-0.212140.132520.64510-0.57789杉木川滇藏DBH<5cm0.372430.00307-0.57789苏浙徽闽赣湘DBH≥5cm0.230640.292230.132520.64510-0.57789落叶松DBH<5cm0.372430.00307-0.57789粤桂贵DBH≥5cm0.17427-0.011770.250590.52050-0.61168DBH<5cm0.181710.71173-0.61168鄂川DBH≥5cm0.081540.286460.250590.52050-0.61168DBH<5cm0.181710.71173-0.61168湘鄂粤桂浙徽DBH≥5cm0.22916-0.355580.648490.38706-0.64126苏川贵DBH<5cm0.111921.47865-0.64126DBH≥5cm0.081540.286460.459900.38706-0.64126赣闽DBH<5cm0.229850.81800-0.64126DBH≥5cm0.22916-0.355580.372530.38706-0.64126黑吉辽、蒙东DBH<5cm0.433250.29325-0.64126部DBH≥5cm0.203010.040280.406080.38706-0.64126DBH<5cm0.200570.82535-0.64126晋冀、蒙中西DBH≥5cm0.155050.207750.311390.44844-0.60415部DBH<5cm0.6744-0.03171-0.60415DBH≥5cm0.203010.040280.435860.44844-0.60415新DBH<5cm0.30570.66884-0.60415DBH≥5cm0.155050.207750.30640.44844-0.60415DBH<5cm0.382450.31069-0.60415DBH≥5cm0.37783-0.025290.289120.44844-0.60415DBH<5cm0.190760.70681-0.60415DBH≥5cm0.091540.855560.319790.44844-0.60415DBH<5cm0.225250.66619-0.60415DBH≥5cm0.38739-0.108950.242380.44844-0.604150.258270.40899-0.604150.160050.440270.449020.48388-0.846420.174210.039440.55155-0.67664川滇藏0.191410.048210.22415-0.04991黑吉0.222620.04217甘青、新疆天0.56614-0.53777山新疆阿尔泰山0.31161-0.00070川0.305160.012290.219060.01184云杉0.38768-0.342840.206700.008100.22437-0.04288滇0.228630.027840.163340.23677藏0.173280.020420.19298-0.04647冷杉黑吉0.28691-0.1156344树种适用范围应用条件一元模型二元模型模型参数模型参数DBH<5cmDBH≥5cmababcDBH<5cm0.200540.9847-0.84642DBH≥5cm0.158920.251430.619970.48388-0.84642DBH<5cm0.152761.35426-0.84642甘青、新疆天DBH≥5cm0.39615-0.134210.443950.48388-0.84642山DBH<5cm0.172641.07073-0.84642DBH≥5cm0.103630.699010.574160.48388-0.84642川DBH<5cm0.270870.95066-0.84642DBH≥5cm0.38367-0.102700.585250.48388-0.84642滇DBH<5cm0.137861.38221-0.84642DBH≥5cm0.149090.296980.180340.88582-0.81895藏DBH<5cm0.297480.57483-0.81895DBH≥5cm0.36688-0.103760.84290.26343-0.76364浙鄂湘川渝DBH<5cm2.68849-0.45725-0.76364DBH≥5cm0.280410.063250.930740.26343-0.76364黑吉辽、蒙东DBH<5cm1.148090.13303-0.76364部DBH≥5cm0.37396-0.126440.71240.26343-0.76364DBH<5cm0.679190.29309-0.76364豫晋冀陕青甘DBH≥5cm0.100400.690620.868940.26343-0.76364宁DBH<5cm0.83558-0.76364柳杉DBH≥5cm0.090770.400930.346-0.03259-0.23991栎树川滇DBH<5cm0.647670.18772-0.23991DBH≥5cm0.27596-0.289930.45432-0.03259-0.23991桦树鄂湘赣徽浙渝DBH<5cm0.620630.11375-0.23991木荷贵DBH≥5cm0.50224-0.188750.4904-0.03259-0.23991枫香DBH<5cm0.552110.31018-0.23991黑吉、蒙东部0.72273-0.414890.37354-0.03259-0.23991(白桦)DBH≥5cm0.523560.58626-0.239910.55458-0.244790.19338黑吉、蒙东部DBH≥5cm(其他桦)0.62147-0.315540.67887DBH≥5cm蒙中西部、豫0.42448-0.21694晋陕青甘宁新DBH<5cm0.47183-0.28265川滇DBH≥5cm0.51776-0.24687滇闽赣贵浙湘DBH<5cm粤0.277650.14032DBH≥5cm桂闽赣湘贵浙0.53510-0.18558徽渝DBH<5cm0.392290.00731京津冀晋蒙鲁DBH≥5cm豫0.51276-0.17885DBH<5cm0.43263-0.07328DBH≥5cm0.45615-0.180790.289580.101530.43256-0.175610.153000.470130.36234-0.0652990.25949-0.62553柏木陕甘青0.89984-0.32461.35270.16614-0.70784川滇藏0.38859-0.150390.35097-0.147400.38859-0.150390.35097-0.14740.02765浙鄂湘粤渝贵0.41217-0.150390.37672-0.1474010.41217-0.150390.37672-0.14740高山松藏川滇0.27622-0.150390.26005-0.147400.027650.27622-0.150390.26005-0.147410.38121-0.150390.35131-0.147400.38121-0.150390.35131-0.14740.027650.205760.177451000.0276510.0276510.0276510.0276510.0276510.05105445一元模型二元模型树种适用范围应用条件模型参数模型参数思茅松普洱、西双版DBH≥5cmababc纳、德宏、临DBH<5cm0.23177-0.116320.191050.11047-0.19119樟子松DBH≥5cm0.123740.273590.156150.23581-0.19119椴树沧DBH≥5cm蒙东部、黑吉0.171860.0750650.21810.36695-0.440640.303490.0565630.440030.57424-0.74241辽黑吉辽榆树黑吉辽冀蒙DBH≥5cm0.3879-0.035110.421120.075352-0.15356黄山松浙徽闽赣DBH≥5cm0.158250.066210.201020.36734-0.46549DBH<5cm0.34655-0.420850.39221-0.04794-0.46549黑吉辽蒙DBH≥5cm0.32719-0.134390.498190.18597-0.50647(天然山杨)DBH<5cm0.31198-0.104820.68533-0.01218-0.50647陕甘青宁新DBH≥5cm0.33369-0.124540.508080.18597-0.50647(天然山杨)DBH<5cm0.59949-0.488570.82444-0.11480-0.50647杨树川滇藏DBH≥5cm0.27936-0.136390.425370.18597-0.50647(天然山杨)DBH<5cm0.22220.005860.338110.32863-0.50647DBH≥5cm0.37302-0.125190.567980.18597-0.50647黑吉辽蒙DBH<5cm0.253310.115280.413240.38359-0.50647(人工杨树)鲁豫蒙冀晋京DBH≥5cm0.33904-0.141410.516230.18597-0.50647(人工杨树)DBH<5cm0.34753-0.156780.544310.15306-0.50647苏徽鄂湘渝贵DBH≥5cm0.25702-0.153980.391350.18597-0.50647(人工杨树)DBH<5cm0.42352-0.464290.60641-0.08615-0.50647注1:一元模型方程的表达式为𝑅𝑆𝑅AF=𝑎⋅𝐷𝐵𝐻𝑏,二元模型方程的表达式为𝑅𝑆𝑅AF=𝑎⋅𝐷𝐵𝐻𝑏⋅𝐻𝑐。其中,RSRAF为乔木林地下生物量与地上生物量的比值,无量纲;DBH为胸径,单位为厘米(cm);H为树高,单位为米(m);a、b、c为模型参数。注2:在造林密度和立地条件等相对一致的情况下,应优先选择基于胸径的一元模型,否则应选择基于胸径和树高的二元模型。本表未包含的树种,可采用文献值。表A.9中国主要乔木林树种的地下生物量与地上生物量的比值(文献整合分析)森林类型RSRAF冷杉林云杉林0.2223落叶松林0.2514杉木林0.2830柏木林0.2332马尾松林0.2489油松林0.2053云南松林0.2349其他针叶林0.1830栎树林0.2436桦树林0.2610杨树林0.2786刺槐林0.2127桉树林0.2767橡树林0.2832其他软阔类0.25300.269046其他硬阔类0.2572针叶混交林0.2364针阔混交林0.2561阔叶混交林0.2598注:RSRAF为乔木林地下生物量与地上生物量的比值,无量纲。表A.10中国主要乔木林树种的基本木材密度和生物量含碳率(文献整合分析)森林类型SVDCFTotalCFAGBCFBGB云冷杉林0.35970.49310.49310.4933落叶松林0.40590.48930.48950.4884温性针叶林0.38970.49610.49670.4955油松林0.42430.51650.51840.5093马尾松林0.44760.52520.52540.5082暖性针叶林0.39420.50340.50450.4912杉类0.30980.49900.50030.4880柏木林0.50100.48470.48460.4851栎类0.57620.48020.48270.4678桦木林0.48480.48720.48970.4779其它硬阔类0.52570.47110.47340.4637杨树林0.41770.47050.47280.4644桉树林0.38480.47300.47500.4685其它软阔类0.38480.47300.47500.4685针叶混0.38280.50050.50140.4943阔叶混0.49670.47180.47410.4652针阔混0.43970.48610.48770.4797注:SVD为基本木材密度,单位为吨每立方米(td.m.·m-3);CFTotal为全树生物量含碳率,单位为吨碳每吨(tC·(td.m.)-1);CFAGB表示地上生物量含碳率,单位为吨碳每吨(tC·(td.m.)-1);CFBGB表示地下生物量含碳率,单位为吨碳每吨(tC·(td.m.)-1)。表A.11中国主要乔木林树种(组)单位面积蓄积量随林龄的Richards生长方程区域树种(组)abcR2东北落叶松100.9212.1760.0720.739284.5501.6220.0140.819华北其他针叶树171.9602.0290.0320.776栎类156.9172.3020.0260.746白桦204.5351.3260.0190.701224.0472.3110.0310.814其他阔叶树236.0511.6290.0190.738针叶混交类117.8902.1890.0440.612针阔混交类207.7301.7450.0140.743102.0440.2880.0010.591落叶松193.5861.2080.0080.649油松173.6041.5870.0200.661其他针叶树55.55302.3190.1770.628栎类205.4191.1750.0090.665桦木杨树其他阔叶树47区域树种(组)abcR2针阔混交类200.4470.9310.0060.68落叶松293.6101.4980.0120.767云杉331.6771.3950.0100.621其他针叶树105.3191.8880.0350.488西北栎类135.7390.9390.1790.520杨树107.4181.1200.0330.500其他阔叶树133.5161.3170.0160.702阔叶混交类120.3231.9470.0400.675针阔混交类146.2211.3640.0230.659马尾松151.5181.5790.0290.616杉木149.5441.7860.0560.693其他针叶树73.3520.9810.0360.558栎类154.6101.5750.0290.721中南杨树92.6942.4660.2100.673其他阔叶树115.4121.2050.0310.603针叶混交类126.8781.5780.0460.624阔叶混交类475.2541.1340.0070.670针阔混交类216.5711.1010.0150.615马尾松142.3381.7520.0390.641杉木183.2071.3550.0370.58其他针叶树161.8921.2520.0310.602栎类265.6962.0280.0310.704东南沿海桉树70.3552.7130.4570.658其他阔叶树269.4051.2180.0160.655针叶混交类266.0601.6560.0250.657阔叶混交类320.2611.6440.0220.742针阔混交类286.9271.5020.0210.735冷杉395.0172.0920.020.533云南松172.7131.7310.030.578马尾松176.9272.0450.0450.740杉木123.3373.3520.140.759西南其他针叶树465.2040.8570.0040.647栎类267.1221.1850.0110.712其他阔叶树492.4650.8370.0340.653针叶混交类513.4581.1230.0080.678阔叶混交类259.7591.3150.0190.724针阔混交类375.5391.2460.0120.740注:方程表达式为𝑉=𝑎∙(1−𝑒−𝑐∙𝐴𝑔𝑒𝑡)𝑏。其中,V为单位面积蓄积量,单位为立方米每公顷(m3·hm-2);𝐴𝑔𝑒𝑡为林龄,无量纲;a、b、c为模型参数;R2为决定系数。表A.12中国主要竹林类型的生物量参数(文献整合分析)指标毛竹其他竹子林分地上生物量(td.m.·hm-2)63.423734.1104林分地下生物量(td.m.·hm-2)30.514416.521148全林生物量(td.m.·hm-2)99.070740.2236林分地下/地上生物量的比值0.51100.7224表A.13中国主要竹林类型单株生物量与胸径的相关方程(文献混合建模)生长型器官abR2胸径范围散生竹地上0.16972.08120.9121.9~17.0整株2.10030.9541.9~17.00.1782丛生竹地上0.47231.79280.9480.5~7.0整株1.89210.9210.5~7.00.4117混生竹地上0.33821.91560.9891.0~5.5整株1.28320.8981.0~5.51.0491注:方程表达式为𝑀=𝑎∙𝐷𝐵𝐻𝑏。其中,M为单株生物量,单位为千克每株(kg·stem-1);DBH为胸径,单位为厘米(cm);a、b为模型参数;R2为决定系数。表A.14中国主要竹林类型单株生物量与胸径和竹高的相关方程(文献混合建模)生长型器官abcR2胸径范围竹高范围地上0.00190.12393.18700.9371.9~17.04.9~20.1散生竹0.03900.51121.89560.9451.9~17.04.9~20.1整株地上0.99932.0499-0.56150.8810.5~7.01.5~16.0丛生竹0.32590.61721.18980.9880.5~7.01.5~16.0整株地上0.27563.4175-1.03150.5191.0~5.51.5~13.8混生竹1.86625.9533-3.29900.6611.0~5.51.5~13.8整株注1:方程表达式为𝑀=𝑎∙𝐷𝐵𝐻𝑏∙𝐻𝑐。其中,M为单株生物量,单位为千克每株(kg·stem-1);DBH为胸径,单位为厘米(cm);H为竹高,单位为米(m);a、b、c为模型参数;R2为决定系数。注2:在造林密度和立地条件等相对一致的情况下,应优先选择基于胸径的一元模型,否则应选择基于胸径和竹高的二元模型。表A.15中国主要灌木林的生物量参数(文献整合分析)指标人工灌木林天然灌木林均值均值林分地上生物量(td.m.·hm-2)8.7383林分地下生物量(td.m.·hm-2)13.47049.3568全林生物量(td.m.·hm-2)6.6772林分地下/地上生物量的比值24.635817.27000.65901.3838表A.16中国主要灌木林单株生物量模型(基于文献汇总)优势灌木物种研究区域或器官自变量模型形式模型系数R2地点abλ沙针、坡柳、石株、虾地上𝐷2𝐻𝑀=𝑎⋅(𝐷2𝐻)𝑏0.0916430.67110.976𝑀=𝑎⋅𝐷𝑏0.0378792.3410.614子花、朝天罐、水锦云南省临沧地下D树、悬钩子、铁刀木、市,膏桐种算盘子、酸藤子、余甘植区整株𝐷2𝐻𝑀=𝑎⋅(𝐷2𝐻)𝑏0.0135940.7370.875子腾格里沙地上𝑉𝑐𝑀=𝑎⋅𝑉𝑐𝑏0.0350.6940.693驼绒藜、盐爪爪、珍珠漠,荒漠生地下𝑉𝑐𝑀=𝑎⋅𝑉𝑐𝑏0.0230.6560.653猪毛菜、红砂𝑀=𝑎⋅𝑉𝑐𝑏0.0690.6720.718态系统整株𝑉𝑐49优势灌木物种研究区域或器官自变量模型形式模型系数R2地点abλ四川红淡、格药柃、栀子、满树星、美丽胡枝江西千烟洲子木、、白山株莓、、杜长鹃托、裁盐英肤、区红,壤林丘下陵灌整株𝑉𝑐𝑀=𝑎⋅𝑉𝑐𝑏⋅𝜆0.1655150.8171.36山矶、白檀、轮叶蒲木层桃、南烛、牡荆大白杜鹃、白背杜鹃、叶𝐷2𝐻𝑀=𝑎+𝑏⋅𝐷2𝐻0.0430.0040.828-0.2110.0640.885金顶杜鹃、小鞍叶羊蹄四川、云南茎𝐷2𝐻𝑀=𝑎+𝑏⋅𝐷2𝐻-0.2160.0820.968甲地下𝐷2𝐻𝑀=𝑎+𝑏⋅𝐷2𝐻千里香杜鹃、密枝杜青海、四叶𝐷2𝐻𝑀=𝑎+𝑏⋅𝐷2𝐻0.0950.020.979鹃、裂毛雪山杜鹃、草川、云南茎𝐷2𝐻𝑀=𝑎+𝑏⋅𝐷2𝐻0.1990.0370.933原杜鹃、头花杜鹃、金地下𝐷2𝐻𝑀=𝑎+𝑏⋅𝐷2𝐻0.3930.0610.969背杜鹃千果榄仁、香合欢、岗叶𝐷2𝐻𝑀=𝑎+𝑏⋅𝐷2𝐻0.0890.0110.9000.260.0960.901柃、胡颓子、潺槁木姜四川、云南茎𝐷2𝐻𝑀=𝑎+𝑏⋅𝐷2𝐻-0.2260.070.966子、盐肤木地下𝐷2𝐻𝑀=𝑎+𝑏⋅𝐷2𝐻小果蔷薇、绢毛蔷薇贵州、四叶𝐷2𝐻𝑀=𝑎+𝑏⋅𝐷2𝐻-0.0010.0270.934川、云南茎𝐷2𝐻𝑀=𝑎+𝑏⋅𝐷2𝐻0.170.060.924地下𝐷2𝐻𝑀=𝑎+𝑏⋅𝐷2𝐻0.0830.0280.955一担柴、银柴、银叶叶𝐷2𝐻𝑀=𝑎+𝑏⋅𝐷2𝐻0.0540.0030.908锥、锐齿斛株、唐梨四川、云南茎𝐷2𝐻𝑀=𝑎+𝑏⋅𝐷2𝐻0.0640.0180.9890.2370.0150.867地下𝐷2𝐻𝑀=𝑎+𝑏⋅𝐷2𝐻叶𝑉𝑐𝑀=𝑎+𝑏⋅𝑉𝑐0.2831.4660.909高山柏、香柏四川、西藏茎𝑉𝑐𝑀=𝑎+𝑏⋅𝑉𝑐-0.4663.860.905地下𝑉𝑐𝑀=𝑎+𝑏⋅𝑉𝑐0.0656.2650.897叶𝑉𝑐𝑀=𝑎+𝑏⋅𝑉𝑐0.0120.1810.858红棕杜鹃、单色杜鹃西藏茎𝑉𝑐𝑀=𝑎+𝑏⋅𝑉𝑐-0.0273.150.932地下𝑉𝑐𝑀=𝑎+𝑏⋅𝑉𝑐-0.051.970.976注:M为单株生物量,单位为千克每株(kg·stem-1);D为基径,为地面高度5cm处的树干直径,单位为厘米(cm);H为植株高度,单位为米(m);D2H为基径平方与株高的乘积;Vc为冠幅投影体积,单位为立方米(m3);a、b为模型系数,λ为幂函数模型中由对数单位(lnM)转换为测量单位(M)时的标准误修正因子;R2表示调整后的决定系数(adj-R2)。表A.17中国主要灌木林单株生物量模型(基于调查的混合物种建模)植被分区枝干形态器官自变量模型形式模型系数R2地上𝐷2𝐻𝑀=𝑎⋅(𝐷2𝐻)𝑏⋅𝜆abλ0.7751地下𝐷2𝐻𝑀=𝑎+𝑏⋅𝐷2𝐻0.03230.8642温带针叶、落叶整株𝐷2𝐻𝑀=𝑎+𝑏⋅𝐷2𝐻-0.00260.79541.21870.4026阔叶混交林区域地下𝑀𝑎0.01210.3354类型a叶𝐷2𝐻𝑀=𝑎⋅𝑀𝑎𝑏0.01190.02250.6205茎枝𝐷2𝐻𝑀=𝑎+𝑏⋅𝐷2𝐻0.04870.6156地上𝐷2𝐻𝑀=𝑎+𝑏⋅𝐷2𝐻0.06690.0397地下𝑀𝑎𝑀=𝑎+𝑏⋅𝐷2𝐻0.11690.6181暖温带落叶阔叶整株𝐷2𝐻0.37530.6172林区域地下𝑀𝑎𝑀=𝑎⋅𝑀𝑎𝑏⋅𝜆0.2652类型a叶𝐴𝑐𝑀=𝑎+𝑏⋅𝐷2𝐻0.65610.00280.0778𝑀=𝑎⋅𝑀𝑎𝑏⋅𝜆0.0266𝑀=𝑎⋅𝐴𝑐𝑏⋅𝜆亚热带常绿阔叶0.0294林区域类型a0.80941.43150.0367亚热带常绿阔叶类型b0.79981.1974林区域0.90591.266950植被分区枝干形态器官自变量模型形式模型系数R2abλ茎枝𝐴𝑐𝑀=𝑎+𝑏⋅𝐴𝑐0.02110.42080.5613地上𝐴𝑐𝑀=𝑎+𝑏⋅𝐴𝑐0.02490.50940.6301温带草原区域类型b地下𝑀𝑎𝑀=𝑎⋅𝑀𝑎𝑏⋅𝜆0.72940.96741.3231整株𝐴𝑐𝑀=𝑎+𝑏⋅𝐴𝑐-0.02141.35200.5550地上𝐴𝑐𝑀=𝑎+𝑏⋅𝐴𝑐0.14801.14800.5435温带荒漠区域类型b地下𝐴𝑐𝑀=𝑎+𝑏⋅𝐴𝑐0.08550.40360.5621整株𝐴𝑐𝑀=𝑎+𝑏⋅𝐴𝑐0.23061.56330.6258地上𝐴𝑐𝑀=𝑎+𝑏⋅𝐴𝑐0.13750.20290.3488青藏被高区原域高寒植类型b地下𝑀𝑎𝑀=𝑎+𝑏⋅𝑀𝑎-0.03531.02960.7104𝑀=𝑎+𝑏⋅𝐴𝑐0.22130.46320.4515整株𝐴𝑐注1:根据个体的枝干形态,可将灌木划分为2种类型:类型a为分枝明确、枝干离散可数的灌木;类型b为分枝不明确、枝干成簇成团的灌木。注2:𝑀𝑎为单株地上生物量,单位为千克每株(kg·stem-1);M为对应器官的生物量,单位为千克每株(kg·stem-1);Ac为冠幅投影面积,单位为平方米(m2);D为基径,为地面高度5cm处的树干直径,单位为厘米(cm);H为植株高度,单位为米(m);D2H为基径平方与株高的乘积;a、b为模型系数,λ为幂函数模型中由对数单位(lnM)转换为测量单位(M)时的标准误修正因子;R2表示调整后的决定系数(adj-R2)。51附录B(资料性附录)森林死有机质碳储量变化计算方法B.1森林死有机质碳储量变化的计算本文件采用“储量变化法(StockDifferenceMethod)”计算项目边界内的森林死有机质碳储量在一段时期内的年均变化量:Δ𝐶DOM,𝑡=𝐶DOM,𝑡2−𝐶DOM,𝑡1×44(B.1)𝑡2−𝑡112式中:Δ𝐶DOM,𝑡——项目开始后第t年的森林死有机质碳储量的年变化量,单位为吨二氧化碳当量每年(tCO2e·a-1);𝐶DOM,𝑡1——第t1年时,森林死有机质碳储量,单位为吨碳(tC);𝐶DOM,𝑡2——第t2年时,森林死有机质碳储量,单位为吨碳(tC);44/12——二氧化碳与碳的相对分子质量之比,无量纲;t——项目开始后第t年,无量纲;t1,t2——项目开始后的第t1年和第t2年,单位为年(a),且t1≤t≤t2。B.2森林死有机质碳储量的计算森林死有机质包括枯落物和枯死木。下列方法适用于所有森林类型,包括乔木林、竹林和灌木林:𝐶DOM,𝑡=∑∑𝐴𝑖,𝑗,𝑡×(𝐿𝐼𝑖,𝑗,𝑡×𝐶𝐹LI,𝑖,𝑗+𝐷𝑊𝑖,𝑗,𝑡×𝐶𝐹DW,𝑖,𝑗)(B.2)𝑖𝑗式中:𝐶DOM,𝑡——第t年时,森林死有机质碳储量,单位为吨碳(tC);𝐴𝑖,𝑗,𝑡——第t年时,第i项目碳层树种j的森林面积,单位为公顷(hm2);𝐿𝐼𝑖,𝑗,𝑡——第t年时,第i项目碳层树种j的单位面积枯落物量,单位为吨每公顷(td.m.·hm-2);𝐷𝑊𝑖,𝑗,𝑡——第t年时,第i项目碳层树种j的单位面积枯死木量,单位为吨每公顷(td.m.·hm-2);𝐶𝐹LI,𝑖,𝑗——第i项目碳层树种j的枯落物含碳率,单位为吨碳每吨(tC·(td.m.)-1),缺省值0.37;𝐶𝐹DW,𝑖,𝑗——第i项目碳层树种j的枯死木含碳率,单位为吨碳每吨(tC·(td.m.)-1),缺省值0.37。B.3枯落物量的计算本文件采用缺省值的方法计算森林枯落物量。首先计算森林单位面积地上生物量,再利用枯落物量与地上生物量的比值,计算单位面积枯落物量:𝐿𝐼𝑖,𝑗,𝑡=𝐴𝐺𝐵𝑖,𝑗,𝑡×𝐷𝐹LI,𝑖,𝑗(B.3)式中:52𝐿𝐼𝑖,𝑗,𝑡——第t年时,第i项目碳层树种j的单位面积枯落物量,单位为吨每𝐴𝐺𝐵𝑖,𝑗,𝑡——公顷(td.m.·hm-2);𝐷𝐹LI,𝑖,𝑗——第t年时,第i项目碳层树种j的单位面积地上生物量,单位为吨每公顷(td.m.·hm-2);第i项目碳层树种j的枯落物量与地上生物量的比值,无量纲。B.4枯死木量的计算本文件采用缺省值的方法估计乔木林枯死木量,竹林或灌木林的枯死木量计为0。首先计算乔木林单位面积地上生物量,再利用枯死木量与地上生物量的比值,计算单位面积枯死木量:𝐷𝑊𝑖,𝑗,𝑡=𝐴𝐺𝐵𝑖,𝑗,𝑡×𝐷𝐹DW,𝑖,𝑗(B.4)式中:𝐷𝑊𝑖,𝑗,𝑡——第t年时,第i项目碳层树种j的单位面积枯死木量,单位为吨每公顷(td.m.·hm-2);𝐴𝐺𝐵𝑖,𝑗,𝑡——第t年时,第i项目碳层树种j的单位面积地上生物量,单位为吨每公顷(td.m.·hm-2);𝐷𝐹DW,𝑖,𝑗——第i项目碳层树种j的枯死木量与地上生物量的比值,无量纲。B.5死有机质参数表B.1中国主要乔木林树种的枯落物量与地上生物量的比值(𝐷𝐹LI,%)(基于文献整合分析)森林类型1年-10年11年-20年适用的林龄范围31年-40年≥41年21年-30年针叶林5.275.545.825.42南方地区阔叶林9.676.924.724.35针阔混7.847.586.784.89针叶林6.019.8014.59北方地区阔叶林7.817.697.69针阔混8.98毛竹林6.63其他竹林17.73灌木林16.30注:南方地区包括上海、江苏、浙江、安徽、江西、福建、广东、广西、海南、四川、云南、贵州、西藏、重庆、湖北、湖南等;北方地区包括黑龙江、吉林、辽宁、内蒙古、山西、河北、北京、天津、山东、河南、陕西、青海、甘肃、宁夏、新疆等。表B.2中国主要乔木林的枯死木量与地上生物量的比值(𝐷𝐹DW,%)(基于文献整合分析)森林类型适用的林龄范围1年-10年11年-20年21年-30年≥31年5.121.74针叶林5.305.82阔叶林南方地区针阔混4.603.2853针叶林3.36北方地区阔叶林3.20针阔混3.28注:南方地区包括上海、江苏、浙江、安徽、江西、福建、广东、广西、海南、四川、云南、贵州、西藏、重庆、湖北、湖南等;北方地区包括黑龙江、吉林、辽宁、内蒙古、山西、河北、北京、天津、山东、河南、陕西、青海、甘肃、宁夏、新疆等。54附录C(资料性附录)森林土壤有机碳储量变化计算方法造林活动由于整地扰动土壤,可能会使项目地块的土壤有机碳储量在整地后发生减少;后期随着林木生长、死亡根系和枯落物返还与分解等,土壤有机碳又会逐渐增加,最终趋于稳定。由于土壤有机碳储量及其变化的监测成本较高、监测结果的不确定性大,基于保守性原则和成本有效性原则,项目业主须基于以下假设条件对土壤有机碳储量及其变化量进行计算:a)整地造林之后0-5年,项目地块的土壤有机碳含量逐渐下降,从第6年之后逐渐上升,恢复至整地前的土壤有机碳水平大约需要20年左右;b)整地造林之后第20-40年,项目地块的土壤有机碳含量呈线性增加,且在第40年后土壤有机碳含量达到稳定状态,即不再增长;c)造林碳汇项目基准线情景下土壤有机碳储量变化量计为0。基于上述假设,项目情景下土壤有机碳储量年变化量可采用如下公式计算:Δ𝑆𝑂𝐶𝑡=∑Δ𝑆𝑂𝐶𝑖,𝑡(C.1)(C.2)𝑖44Δ𝑆𝑂𝐶𝑖,𝑡=𝛿𝑆𝑂𝐶×12×𝐴𝑖,𝑡式中:——整地造林后第t年项目边界内土壤有机碳储量的年变化量,单位Δ𝑆𝑂𝐶𝑡为吨二氧化碳当量每年(tCO2e·a-1);——第i项目碳层整地造林后第t年的土壤有机碳储量的年变化量,Δ𝑆𝑂𝐶𝑖,𝑡单位为吨二氧化碳当量每年(tCO2e·a-1);——整地造林后土壤有机碳密度平均年变化率,单位为吨碳每公顷每𝛿𝑆𝑂𝐶年(tC·hm-2·a-1),见表C.1;——二氧化碳与碳的相对分子质量之比,无量纲;44/12——𝐴𝑖,𝑡——第t年时,第i项目碳层的面积,单位为公顷(hm2);i——t项目碳层,i=1,2,3……,无量纲;自项目整地造林后的年数,无量纲。表C.1整地造林后土壤有机碳密度年变化率参考值(tC·hm-2·a-1)(基于文献整合分析)整地造林后的年限常绿阔叶落叶阔叶针叶竹子灌木0-5年-0.40-0.40-0.40-0.40-0.206-20年+0.20+0.15+0.15+0.15+0.1021-40年+0.70+0.40+0.40+0.40+0.10≥41年0000055附录D(资料性附录)森林火烧引起的温室气体排放量计算方法D.1森林火烧引起的温室气体排放量的计算森林火烧包括火灾和人为火烧。森林火灾会引起林木地上生物质和死有机质的燃烧,人为伐除疫病死木并烧除,都可能会造成非CO2温室气体的排放。(D.1)𝐺𝐻𝐺𝑡=𝐺𝐻𝐺FIRE,AGB,𝑡+𝐺𝐻𝐺FIRE,DOM,𝑡+𝐺𝐻𝐺BURN,𝑡式中:𝐺𝐻𝐺𝑡——第t年时,项目边界内由于森林火烧引起的非CO2温室气体排放量,单位为吨二氧化碳当量每年(tCO2e·a-1);𝐺𝐻𝐺FIRE,AGB,𝑡——第t年时,项目边界内由火灾引起的地上生物质燃烧造成的非CO2温室气体排放量,单位为吨二氧化碳当量每年(tCO2e·a-1);𝐺𝐻𝐺FIRE,DOM,𝑡——第t年时,项目边界内由火灾引起的死有机质燃烧造成的非CO2温室气体排放量,单位为吨二氧化碳当量每年(tCO2e·a-1);𝐺𝐻𝐺BURN,𝑡——第t年时,项目边界内由人为火烧造成的非CO2温室气体排放量,单位为吨二氧化碳当量每年(tCO2e·a-1)。D.2火灾引起地上生物质燃烧造成的温室气体排放的计算火灾引起森林地上生物质燃烧造成的非CO2温室气体排放,使用最近一次项目核查时(TV)划分的碳层、各碳层地上生物量数据和燃烧因子进行计算:𝐺𝐻𝐺FIRE,AGB,𝑡=∑[𝐴BURN,𝑖,𝑡×𝐴𝐺𝐵𝑖,𝑇𝑉×𝐶𝑂𝑀𝐹𝑖(D.2)𝑖×(𝐸𝐹CH4×𝐺𝑊𝑃CH4+𝐸𝐹N2O×𝐺𝑊𝑃N2O)]×10−3式中:第t年时,项目边界内由于森林火灾引起地上生物质燃烧造成𝐺𝐻𝐺FIRE,AGB,𝑡——的非CO2温室气体排放量,单位为吨二氧化碳当量每年(tCO2e·a-1);𝐴BURN,𝑖,𝑡——第t年时,第i项目碳层发生燃烧的土地面积,单位为公顷𝐴𝐺𝐵𝑖,𝑇𝑉——(hm2);𝐶𝑂𝑀𝐹𝑖——火灾发生前,项目最近一次核查时第i项目碳层的地上生物𝐸𝐹CH4——量,单位为吨每公顷(td.m.·hm-2);𝐸𝐹N2O——第i项目碳层的燃烧指数,针对植被类型取值,无量纲;𝐺𝑊𝑃CH4——𝐺𝑊𝑃N2O——CH4排放因子,单位为克甲烷每千克(gCH4·(kgd.m.)-1);N2O排放因子,单位为克氧化亚氮每千克(gN2O·(kgd.m.)-1);CH4的全球增温潜势,用于将CH4转换成CO2e,无量纲;N2O的全球增温潜势,用于将N2O转换成CO2e,无量纲;56i——项目碳层,i=1,2,3……;根据第TV年核查时的分层确定,无量纲;t——自项目开始以来的年数,t=1,2,3……,无量纲;TV——自项目开始至项目最近一次核查的时间,无量纲;10-3——将千克转换成吨的常数。第一次核查时,如果核算期内有火灾发生,但不清楚燃烧前的地上生物量,可保守地采用第一次核查时火灾发生所在的同一碳层的平均单位面积地上生物量进行计算。D.3火灾引起死有机质燃烧造成的温室气体排放的计算森林火灾引起死有机质燃烧造成的非CO2温室气体排放,应使用最近一次核查(TV)的碳层死有机质碳储量来计算。第一次核查时,可保守地采用第一次核查时同一碳层的平均单位面积死有机质碳储量进行计算:(D.3)44𝐺𝐻𝐺FIRE,DOM,𝑡=∑[𝐴BURN,𝑖,𝑡×(𝐷𝑊𝑖,𝑇𝑉×𝐶𝐹DW+𝐿𝐼𝑖,𝑇𝑉×𝐶𝐹LI)]×12×0.07𝑖式中:𝐺𝐻𝐺FIRE,DOM,𝑡——第t年时,项目边界内由于森林火灾引起死有机质燃烧造成的非CO2温室气体排放量,单位为吨二氧化碳当量每年(tCO2e·a-1);𝐴BURN,𝑖,𝑡——第t年时,第i项目碳层发生燃烧的土地面积,单位为公顷(hm2);𝐷𝑊𝑖,𝑇𝑉——火灾发生前,项目最近一次核查时第i项目碳层的枯死木单位面积干物质量,单位为吨每公顷(td.m.·hm-2),使用附录B的方法计算;𝐶𝐹DW——枯死木干物质含碳率,单位为吨碳每吨(tC·(td.m.)-1);𝐿𝐼𝑖,𝑇𝑉——火灾发生前,项目最近一次核查时第i项目碳层的枯落物单位面积干物质量,单位为吨每公顷(td.m.·hm-2),使用附录B的方法计算;𝐶𝐹LI——枯落物干物质量含碳率,单位为吨碳每吨(tC·(td.m.)-1);i——项目碳层,i=1,2,3……;根据第Tv年核查时的分层确定,无量纲;t——自项目开始以来的年数,t=1,2,3……,无量纲;Tv——自项目开始至项目最近一次核查的时间,无量纲;——二氧化碳与碳的相对分子质量之比,无量纲;4412——非CO2排放量占碳储量的比例。0.07D.4人为火烧引起的温室气体排放项目期内如果有伐除病原疫木并烧除的活动,则须计算人为火烧活动引起的非CO2温室气体排放。可通过调查人为火烧活动发生的碳层内采伐病原疫木的数量比例(如蓄积量比例或株数比例),使用最近一次项目核查时(TV)划分的相同碳层的平均地上生物量数57据来计算燃烧的地上生物量,结合燃烧因子计算人为火烧活动引起的非CO2温室气体排放:𝐺𝐻𝐺BURN,𝑡=∑[𝐴𝑖,𝑡×𝐴𝐺𝐵𝑖,𝑇𝑉×𝑅BURN,𝑖,𝑡×𝐶𝑂𝑀𝐹𝑖(D.4)𝑖×(𝐸𝐹CH4×𝐺𝑊𝑃CH4+𝐸𝐹N2O×𝐺𝑊𝑃N2O)]×10−3式中:𝐺𝐻𝐺BURN,𝑡——第t年时,项目边界内由于人为火烧引起的非CO2温室气体排放量,单位为吨二氧化碳当量每年(tCO2e·a-1);𝐴𝑖,𝑡——𝐴𝐺𝐵𝑖,𝑇𝑉——第t年时,第i项目碳层的土地面积,单位为公顷(hm2);火灾发生前,项目最近一次核查时第i项目碳层的地上生物量,单位为吨每公顷(td.m.·hm-2),使用附录A的方法计算;𝑅BURN,𝑖,𝑡——第t年时,第i项目碳层内烧除病原疫木的数量占比(如蓄积量比例或株数比例),单位为百分比(%);𝐶𝑂𝑀𝐹𝑖——第i项目碳层的燃烧指数(针对每个植被类型),无量纲;𝐸𝐹CH4——𝐸𝐹N2O——CH4排放因子,单位为克甲烷每千克(gCH4·(kgd.m.)-1);N2O排放因子,单位为克氧化亚氮每千克(gN2O·(kgd.m.)-1);𝐺𝑊𝑃CH4——CH4的全球增温潜势,用于将CH4转换成CO2e,无量纲;𝐺𝑊𝑃N2O——i——N2O的全球增温潜势,用于将N2O转换成CO2e,无量纲;项目碳层,i=1,2,3……;根据第TV年核查时的分层确定,无量纲;t——自项目开始以来的年数,t=1,2,3……,无量纲;TV——自项目开始至项目最近一次核查的时间,无量纲;10-3——将千克转换成吨的常数。第一次核查时,如果核算期内有人为火烧发生,但不清楚燃烧前的地上生物量,可保守地采用第一次核查时火灾发生所在的同一碳层的平均单位面积地上生物量。表D.1不同植被类型燃烧因子参考值森林类型林龄缺省值3年-5年0.46亚热带/热带森林6年-10年0.6711年-17年0.50寒温带森林≥18年0.32温带森林全林龄0.40全林龄0.4558附录E(资料性附录)监测样地数量计算与样地布设方法E.1抽样设计要求项目监测所需的样地数量,可采用如下方法进行计算:第一步:计算样地数量n。如果n≥30,则该步骤所得样地数为n值;如果n<30,则须采用自由度为n-1时的𝑡𝑉𝐴𝐿值,进行第二次迭代计算,得到的n值即为最终的样地数;𝑁×𝑡𝑉2𝐴𝐿×(∑𝑖𝑤𝑖×𝑆𝑖)2(E.1)𝑛=𝑁×𝐸2+𝑡𝑉2𝐴𝐿×∑𝑖(𝑤𝑖×𝑆𝑖2)式中:n——项目边界内计算生物质碳储量所需的监测样地数量,无量纲;N——项目边界内监测样地的抽样总体,N=A/Ap;其中A是项目总面积,Ap是样地面积,无量纲;𝑡𝑉𝐴𝐿——可靠性指标。在一定的可靠性水平下,自由度为无穷(∞)时查t-分布双侧t-分位数表的t值,取值为1.645,无量纲;𝑤𝑖——项目边界内第i项目碳层的面积权重,wi=Ai/A,其中A是项目总面积(hm2),Ai是第i项目碳层的面积(hm2),无量纲;𝑆𝑖——项目边界内第i项目碳层单位面积生物质碳储量估计值的标准差,单位为吨碳每公顷(tC·hm-2);项目设计阶段,采用碳层单位面积生物质碳储量估计值的10%;E——项目单位面积生物质碳储量估计值允许的误差范围,单位为吨碳每公顷(tC·hm-2);项目设计阶段,采用项目单位面积生物质碳储量估计值的10%;i——项目碳层,i=1,2,3……,无量纲。第二步:当抽样面积较大时(抽样面积大于项目面积的5%),按公式(E.1)进行计算获得样地数n之后,按照公式(E.2)对n值进行调整,从而确定最终的样地数(𝑛𝑎):1(E.2)𝑛𝑎=𝑛×1+𝑛⁄𝑁式中:𝑛𝑎——调整后项目边界内计算生物质碳储量所需的监测样地数量,无量纲;n——项目边界内计算生物质碳储量所需的监测样地数量,无量纲;N——项目边界内监测样地的抽样总体,无量纲。第三步:当抽样面积较小时(抽样面积小于项目面积的5%),可采用简化方式计算:(E.3)𝑡𝑉𝐴𝐿22𝑛=(𝐸)×(∑𝑤𝑖×𝑆𝑖)𝑖式中:n——项目边界内计算生物质碳储量所需的监测样地数量,无量纲;𝑡𝑉𝐴𝐿——可靠性指标。在一定的可靠性水平下,自由度为无穷(∞)时查t-分布双侧t-分位数表的t值,取值为1.645,无量纲;59𝑤𝑖——项目边界内第i项目碳层的面积权重,wi=Ai/A,其中A是项目总面积(hm2),Ai是第i项目碳层的面积(hm2),无量纲;𝑆𝑖——项目边界内第i项目碳层单位面积碳储量估计值的标准差,单位为吨碳每公顷(tC·hm-2);项目设计阶段,采用碳层单位面积生物质碳储量估计值的10%;E——项目单位面积生物质碳储量估计值允许的绝对误差限,单位为吨碳每公顷(tC·hm-2);项目设计阶段,采用项目单位面积生物质碳储量估计值的10%;i——项目碳层,i=1,2,3……,无量纲。第四步:分配到各碳层的监测样地数量,采用最优分配法进行计算。分配样地数量不足3个的碳层,最少设置3个样地:𝑛=𝑛×𝑤𝑖×𝑆𝑖𝑖(E.4)∑𝑖(𝑤𝑖×𝑆𝑖)式中:𝑛𝑖——项目边界内第i项目碳层计算生物质碳储量所需的监测样地数量,无量纲;n——项目边界内计算生物质碳储量所需的监测样地数量,无量纲;𝑤𝑖——项目边界内第i项目碳层的面积权重,wi=Ai/A,其中A是项目总面积(hm2),Ai是第i项目碳层的面积(hm2),无量纲;𝑆𝑖——项目边界内第i项目碳层单位面积生物质碳储量标准差的估计值,单位为吨碳每公顷(tC·hm-2);项目设计阶段,采用碳层单位面积生物质碳储量估计值的10%;i——项目碳层,i=1,2,3……,无量纲。E.2样地布设方案样地的空间布设须采用随机起点、系统布点的方法,具体操作流程如下:a)采用GIS等空间工具将每个碳层网格化,每个网格面积大小与监测样地面积大小相同。计算第i个项目碳层的总网格交叉点的数量(𝑁𝑖),将每个网格交叉点按照固定顺序编号,从1、2、3……直到𝑁𝑖;b)在1~𝑁𝑖之间产生一个随机数(如在Excel表格中,使用随机数公式f(x)=ROUND(RAND()(𝑁𝑖),0)产生一个随机数),该随机数代表的网格交叉点编号即为第i项目碳层的第1个监测样地的中心点。第2个样地的中心点等于第1个样地的网格交叉点编号加间隔的网格交叉点数,该间隔数等于第i项目碳层的总网格交叉点数量(𝑁𝑖)除以该碳层样地数量(𝑛𝑖)后取整数;第3个样地中心点的网格交叉点编号等于第2个样地的网格交叉点编号加间隔的网格交叉点数,依此类推。c)将每个确定的网格交叉点作为监测样地的中心点,在GIS等空间分析工具的帮助下,确定每个监测样地的经纬度坐标。坐标以度表示,至少保留6位小数。E.3样地设置根据确定的样地中心点坐标,使用RTK等定位工具找到样地中心点准确地理位置。如果样地边界距离林缘、悬崖等地形地物小于10m,可沿这些地形地物边缘线的垂直方向移动一定距离,最多不超过10m,并记录新的样地中心点坐标。坐标以度表示,至少保留6位小数。60首次监测时,可在样地中心点设置永久性标志,便于后续监测时的位置识别。样地边界除用于测定时识别外,不宜建立永久性标志。样地面积为0.04hm2~0.06hm2,样地形状采用矩形(样地测量闭合差≤0.5%)或圆形。对于在坡地上的样地,须进行坡度校正。在同一个项目中,所有样地的面积和形状应相同,样地内林木和管理方式应与样地外完全一致。对于灌木造林,须在每个样地的四个角和中心点位置,设置5个固定样方,每个固定样方面积不小于4m2。现场记录经纬度坐标(以度表示的坐标至少保留6位小数)、位置(县、乡、村和小地名)、样地名称/编号、样地的形状和面积大小、树种和造林时间等信息。固定样地复位率须达100%。为了便于核查,项目开始前可对样地内的原有林木单独进行标记,在项目计入期内可不进行监测。E.4样地调整项目实施阶段,如果重新调整了碳层划分,或为了满足抽样精度需要额外增加样地,须对碳层内的样地数量和布设进行调整。每个碳层在保留已有样地的基础上,在新碳层内按照上述原则和步骤,补充并布设新的监测样地,以确保抽样精度。61附录F(资料性附录)森林生物质碳储量变化的样地监测方法第一步:样地每木检尺,实测样地内除基准线林木以外所有活立木的胸径(或竹径),根据选择测量树高(或竹高),起测胸径为2cm。对于灌木林可只选择监测样方内灌木株数、单株灌木的基径、分枝数、灌高、灌幅等,或监测样方内的灌木盖度。第二步:采用附录A的方法计算单株林木(或竹子、灌木)的生物量,累加得到样地水平的生物量,再结合生物量含碳率计算样地水平的生物质碳储量,以及各碳层的平均单位面积生物质碳储量。第三步:计算第i项目碳层样本平均数(平均单位面积生物质碳储量的估计值)及其方差:𝑐Biomass,𝑖,𝑡=(∑𝑐Biomass,𝑖,𝑝,𝑡)⁄𝑛𝑖(F.1)𝑝𝑆2=𝑛𝑖×∑𝑝𝑐B2iomass,𝑖,𝑝,𝑡−(∑𝑝𝑐Biomass,𝑖,𝑝,𝑡)2(F.2)𝑐Biomass,𝑖,𝑡𝑛×(𝑛−1)𝑖𝑖式中:𝑐Biomass,𝑖,𝑡——第t年时,第i项目碳层平均单位面积生物质碳储量,单位为吨碳每公顷(tC·hm-2);𝑐Biomass,𝑖,𝑝,𝑡——第t年时,第i项目碳层样地p的单位面积生物质碳储量,单位为吨碳每公顷(tC·hm-2);𝑛𝑖——第i项目碳层的监测样地数,无量纲;𝑆2——第t年时,第i项目碳层平均单位面积生物质碳储量的方差,单𝑐Boimass,𝑖,𝑡位为吨碳每公顷的平方((tC·hm-2)2);第i项目碳层中的样地,p=1,2,3……,𝑛𝑖,无量纲;p——项目碳层,i=1,2,3……,无量纲;i——t——自项目开始以来的年数,t=1,2,3……,无量纲。第四步:计算项目总体平均数估计值(平均单位面积生物质碳储量的估计值)及其方差:𝑐Biomass,𝑡=∑(𝑤𝑖×𝑐Biomass,𝑖,𝑡)(F.3)(F.4)𝑖𝑆2𝑐Biomass,𝑡=∑(𝑤2𝑖×2𝑆)𝑐Biomass,𝑖,𝑡𝑛𝑖𝑖式中:——第t年时,项目边界内的平均单位面积生物质碳储量,单位为吨𝑐Biomass,𝑡——碳每公顷(tC·hm-2);——项目边界内第i项目碳层的面积权重,wi=Ai/A,其中A是项目总𝑤𝑖——面积(hm2),Ai是第i项目碳层的面积(hm2),无量纲;𝑐Biomass,𝑖,𝑡——第t年时,第i项目碳层的平均单位面积生物质碳储量,单位为𝑆2吨碳每公顷(tC·hm-2);第t年时,项目平均单位面积生物质碳储量的方差,单位为吨碳𝑐Biomass,𝑡每公顷的平方((tC·hm-2)2);第t年时,第i项目碳层平均单位面积生物质碳储量的方差,单𝑆262𝑐Biomass,𝑖,𝑡位为吨碳每公顷的平方((tC·hm-2)2);𝑛𝑖——第i项目碳层的监测样地数,无量纲;i——项目碳层,i=1,2,3……,无量纲;t——自项目开始以来的年数,t=1,2,3……,无量纲。第五步:计算项目边界内平均单位面积生物质碳储量的不确定性:𝑡𝑉𝐴𝐿×𝑆𝑐Biomass,𝑡(F.5)𝑢𝑐Biomass,𝑡=𝑐Biomass,𝑡式中:𝑢𝑐Biomass,𝑡——第t年时,项目边界内平均单位面积生物质碳储量的估计值的不确定性,即相对误差限,单位为百分比(%)。要求相对误差不大于10%,即抽样精度不低于90%;𝑡𝑉𝐴𝐿——可靠性指标,自由度等于n-M(其中n是项目边界内样地总数,M是林木生物量计算的碳层数),置信水平为90%,查t-分布双侧分位数表获得,无量纲。如置信水平为90%,自由度为45时,双侧t-分布的t值在Excel电子表中输入“=TINV(0.10,45)”可计算得到t值为𝑆𝑐Biomass,𝑡——1.6794;第t年时,项目边界内平均单位面积生物质碳储量方差的平方根,即𝑐Biomass,𝑡——标准误差,单位为吨碳每公顷(tC·hm-2);第t年时,项目边界内的平均单位面积生物质碳储量,单位为吨碳每t——公顷(tC·hm-2);自项目开始以来的年数,t=1,2,3……,无量纲。第六步:计算第t年项目边界内的生物质总碳储量:𝐶Biomass,𝑡=𝐴×𝑐Biomass,𝑡(F.6)式中:𝐶Biomass,𝑡——第t年时,项目边界内生物质碳储量,单位为吨碳(tC);A——项目边界内各碳层的面积总和,单位为公顷(hm2);𝑐Biomass,𝑡——第t年时,项目边界内的平均单位面积生物质碳储量,单位吨t——碳每公顷(tC·hm-2);自项目开始以来的年数,t=1,2,3……,无量纲。第七步:计算核算期内第t年项目边界内生物质碳储量的年变化量。假设在核算期内,生物质碳储量变化是线性的:Δ𝐶Biomass,𝑡=𝐶Biomass,𝑡2−𝐶Biomass,𝑡1×44(F.7)𝑡2−𝑡112式中:Δ𝐶Biomass,𝑡——核算期内第t年的项目边界内生物质碳储量的年变化量,单位为吨二氧化碳当量每年(tCO2e·a-1);𝐶Biomass,𝑡1——第t1年时,项目边界内生物质碳储量,单位为吨二氧化碳当量每年(tCO2e·a-1);𝐶Biomass,𝑡2——第t2年时,项目边界内生物质碳储量,单位为吨二氧化碳当量每年(tCO2e·a-1);44/12——二氧化碳与碳的相对分子质量之比,无量纲;t1,t2——核算期的间隔时间,单位为年(a);开始于第t1年,结束于第t2年,且t1≤t≤t2;t——自项目开始以来的年数,t=1,2,3……,无量纲。如果乔木造林使用的苗木平均胸径<2cm,可假定造林时的林木生物质碳储量为0。如果乔木造林活动采用了平均胸径≥2cm以上的大苗,须根据造林时各碳层的造林树种、平63均胸径、平均树高、造林密度、造林面积等,采用附录A的方法计算造林时的林木生物质碳储量。64附录G(资料性附录)森林生物质碳储量的机载激光雷达监测方法本附录文件主要适用于乔木(或竹子)造林。在项目监测阶段,项目业主须在减排量核算中明确说明实施过程,并提交实施过程中完整的数据文件,包括但不限于激光雷达点云文件(.las)、点云数据配套的边界矢量数据文件(.shp或.kml)、数据处理代码文件以及相应的模型计算结果文件(.csv或.xlsx)。首次使用机载激光雷达进行森林生物质碳储量监测的操作流程如图G.1所示:图G.1首次使用机载激光雷达监测森林生物质碳储量的操作流程第一步:样本单木选择根据项目树种选择在项目所在地针对每个树种随机选择不少于60株单木作为建模样本。样本单木胸径不小于2cm,按照2cm~3cm的径阶间隔进行划分,每个径阶不少于2株~3株样木,尽可能覆盖所有径阶。所有被选择的样本单木应较容易从激光雷达点云数据中分割。使用RTK等设备对单木坐标进行精确定位(定位误差小于0.3m)并标记编号(k=1,2,653……)。第二步:样本单木生物量实测实测所有样本单木的胸径和(或)树高,并利用生物量方程法或生物量转换因子法计算单木k的全株(或地上)生物量,以此作为单木k的实测全株(或地上)生物量(𝑏𝑘m)。第三步:激光雷达单木特征参数获取参照《CH/T8024-2011机载激光雷达数据获取技术规范》,使用机载激光雷达对所有样木所在范围内的林木进行扫描,获取机载激光雷达数据,进行噪声去除、点云分类(地面点与非地面点)、点云归一化、单木分割、单木特征参数提取等处理。单木特征参数包括但不限于树高、百分位树高、冠幅、激光雷达生物量指数(LiDARBiomassIndex,LBI)等。a)单木位置和树高利用基于栅格的方法或基于点云的方法对机载激光雷达数据进行单木分割,获取每株单木的机载激光点云,并从单木机载激光点云中计算获取单木的位置和树高。单木位置和树高可由最高点的空间坐标来确定:(𝑋Tree,𝑌Tree,𝑍Tree)=(𝑋Highest,𝑌Highest,𝑍Highest)(G.1)式中:(𝑋Tree,𝑌Tree,𝑍Tree)——最终确定的单木三维坐标,单位为米(m);(𝑋Highest,𝑌Highest,𝑍Highest)——单木最高点的三维坐标,单位为米(m)。b)单木冠幅单木点云在X方向和Y方向上垂直投影的平均值,视为冠幅直径,可采用下式计算:𝐶=(𝑋max−𝑋min)+(𝑌max−𝑌min)(G.2)2式中:𝐶——单木的冠幅,单位为米(m);——单木点云在X方向垂直投影的最大值,单位为米(m);𝑋max——单木点云在X方向垂直投影的最小值,单位为米(m);𝑋min——单木点云在Y方向垂直投影的最大值,单位为米(m);𝑌max——单木点云在Y方向垂直投影的最小值,单位为米(m)。𝑌min第四步:单木定位匹配选取统一的参考坐标系(WGS84或CGCS2000)对激光雷达单木定位与每木检尺样本单木定位进行匹配,获取每个样本单木对应的激光雷达特征参数和实测单木全株(或地上)生物量,形成单木特征参数与样地实测全株(或地上)生物量的数据集。在进行该步骤时须同时检验单木分割精度,分割精度应大于等于90%,否则应重新检验该技术在当前场景下的适用性。𝐾l−𝐾m(G.3)𝜔=(1−𝐾m)×100%式中:𝜔——单木分割精度,单位百分比(%);𝐾l——点云单木分割数量,单位为株(stem);66𝐾m——样地实测单木数量,单位为株(stem)。第五步:激光雷达生物量模型构建根据不同林地的实际情况,可选择适当的拟合方法(如最小二乘法)构建激光雷达生物量模型。随机选取从第四步数据集中一定数量的样本作为训练集(如训练集与验证集按3:1的比例),构建激光雷达生物量模型:𝑏𝑘m=𝑓(𝑥𝑘l)𝑚𝑖𝑛𝑖𝑚𝑖𝑧𝑖𝑛𝑔𝑆𝑆𝐸𝑏𝑘l=𝑓(̂𝑥𝑘l)(G.4)→式中:𝑏𝑘m——第k株单木的样地实测全株(或地上)生物量,单位为千克(kg𝑏𝑘l——d.m.);𝑥𝑘l——第k株单木的激光雷达计算的全株(或地上)生物量,单位为千SSE——克(kgd.m.);𝑓(𝑥𝑘l)——第k株单木的激光雷达特征参数,单位为米(m)或无量纲;𝑓(̂𝑥𝑘l)——残差平方和,𝑆𝑆𝐸=∑𝑘(𝑏𝑘m−𝑏𝑘l)2;样地实测全株(或地上)生物量与激光雷达特征参数之间的假设l——函数形式;通过最小化调查样地数据内的SSE,可对𝑓(𝑥𝑘l)求解参数,并将其记为𝑓(̂𝑥𝑘l),单位为千克每株(kgd.m.·stem-1);拟合所得的激光雷达全株(或地上)生物量计算函数;也即对𝑓(𝑥𝑘l)求解参数的结果,单位为千克每株(kgd.m.·stem-1);通过激光雷达直接获取或计算的值,无量纲;m——样地实测直接获取或计算的值,无量纲;k——样本单木,k=1,2,3……,无量纲。此处以基于激光雷达生物量指数和树高的拟合方法为例,建立激光雷达特征参数与林木(或林地)特征参数的关系模型,也可选择采用其他建模方法:a)LBI的计算对任一样本分割单木的激光点云数据从枝下高度(树冠低点)到树冠顶点高度以一定的高度间隔进行分层,并根据每层的点云计算对应层的体积,最终对所有层进行积分,获取单木树冠的激光雷达生物量指数:(G.5)𝐻𝑇𝐿𝐵𝐼=𝑙𝑖𝑚∑{𝑈𝐿(𝐻)×𝜋×[𝑟(𝐻)]2×𝛥𝐻×𝐻}𝛥𝐻→0𝐻=𝐻𝐶式中:𝐿𝐵𝐼——单木的激光雷达生物量指数,无量纲;𝑈𝐿(𝐻)——从解析木的点云数据中提取的单木叶面积的体积密度函数,单位为平方米每立方米(m2·m-3);𝑟(𝐻)——高度为H处的树冠的半径,通过点云单木特征参数提取,单位为米(m);𝛥𝐻——树冠纵向分层的间隔,单位为米(m);𝐻𝐶——枝下高度(树冠低点),单位为米(m);𝐻𝑇——单木树冠顶点高度,单位为米(m)。b)生物量模型的构建结合从机载激光雷达中提取的所有样本单木的激光雷达生物量指数和树高,及野外实测获取的单木全株(或地上)生物量,采用非线性回归方法构建单木生物量模型,公式(G.6)以使用LBI反演地上生物量的方法为例。也可采用其他类似的建模方法,选择不同的激光雷达特征参数和林木(或林地)特征参数,建立其与生物量之间的关系模型。672𝛽(G.6)ln𝐴𝐺𝐵=ln𝑝+𝛽ln𝐻𝑇+𝛼ln𝐿𝐵𝐼式中:——𝐴𝐺𝐵——单木地上生物量,单位为千克每株(kg·stem-1);𝐻𝑇——单木的树冠高度,单位为米(m);𝐿𝐵𝐼——单木的激光雷达生物量指数,无量纲;p——模型参数,无量纲;α——模型参数,无量纲;β模型参数,无量纲。第六步:模型精度检验计算模型决定系数(𝑅2)。对于单一树种,训练集决定系数应不低于0.9,验证集决定系数应不低于0.8;对于多树种混交林,训练集与验证集决定系数均应不低于0.7。若满足精度要求,后续监测可复用首次监测时建立的激光雷达生物量模型进行计算。若不能满足精度要求,则应重新建模,直至达到精度要求。𝑅2=∑𝑘(𝑏𝑘m−𝑏𝑘l)2(G.7)1−∑(𝑏m−𝑏̄m)2𝑘𝑘𝑘式中:𝑅2——模型决定系数,无量纲;𝑏̄𝑘m——样地实测全株(或地上)生物量的平均值,单位为千克(kgd.m.);𝑏𝑘m——第k株单木的样地实测全株(或地上)生物量,单位为千克(kgd.m.);𝑏𝑘l——第k株单木的激光雷达计算的全株(或地上)生物量,单位为千克(kgd.m.);k——样本单木,k=1,2,3……,无量纲;l——通过激光雷达直接获取或计算的值,无量纲;m——样地实测直接获取或计算的值,无量纲。第七步:碳层生物量和碳储量计算若模型满足第六步中的精度要求,则使用第五步得到的激光雷达生物量模型计算第i项目碳层内每株单木的全株(或地上)生物量,形成第i项目碳层内j树种激光雷达单木全株(或地上)生物量数据集𝑏𝑖l,𝑗,1,𝑏𝑖l,𝑗,2⋯⋯,𝑏𝑖l,𝑗,𝑠,再结合各类型树种生物量含碳率(若计算的是地上生物量,则须利用地下生物量与地上生物量的比值转化为全株生物量)将单木生物量转换为单木生物质碳储量,再计算各碳层的生物质碳储量:𝑐𝑖l=∑∑𝑐𝑖l,𝑗,𝑠(G.8)𝑗𝑠式中:第i项目碳层的激光雷达生物质碳储量,单位为吨碳(tC);𝑐𝑖l——𝑐𝑖l,𝑗,𝑠——第i项目碳层第j类树种第s株单木的激光雷达生物质碳储量,单位为吨碳(tC);s——第i项目碳层内第j类树种的单木,s=1,2,3……,单位为株(stem);i——项目碳层,i=1,2,3……,无量纲;j——树种,j=1,2,3……,无量纲;l——通过激光雷达直接获取或计算的值,无量纲。68第八步:项目边界内生物质碳储量计算利用各碳层生物质碳储量,计算项目边界内生物质碳储量:𝑐l=∑𝑐𝑖l(G.9)式中:𝑖𝑐l——项目边界内激光雷达生物质碳储量,单位为吨碳(tC);𝑐𝑖l——第i项目碳层的激光雷达生物质碳储量,单位为吨碳(tC);i——项目碳层,i=1,2,3……,无量纲;l——通过激光雷达直接获取或计算的值,无量纲。第九步:项目边界内生物质碳储量年变化量的计算利用附录A公式(A.1)计算项目边界内森林生物质碳储量的年变化量。69