竹子造林碳汇项目方法学（版本号V01）2013年10月编制说明竹林作为一种特殊的植被类型，是我国重要的森林类型之一。我国竹子资源十分丰富，是世界上竹类分布最广、资源最丰富的国家，在竹子栽培、利用等方面具有悠久的历史，被誉为“竹子王国”。竹子造林是我国重要的造林类型之一，而现有的CDM造林再造林项目方法学不适于竹子造林。本方法学在传统CDM造林再造林方法学的基础上，增加了竹产品碳库；提供了可供项目参与方选择的新的基线情景识别和额外性论证程序；提供竹子造林碳计量方法。目录第I部分.来源、定义和适用条件11.来源12.规范性引用文件13.定义14.适用条件3第II部分.基线和碳计量方法51.项目边界52.土地合格性53.碳库和温室气体排放源选择64.计入期选择75.基线情景识别和额外性论证76.碳层划分87.基线碳汇量88.项目碳汇量99.泄漏1910.项目减排量19第III部分.监测程序201.项目实施监测202.抽样设计和碳层划分213.精度控制和校正254.不需监测的数据和参数（采用的缺省值或一次性测定值）255.监测的数据和参数356.参考文献39附件:竹子生物量方程411第I部分.来源、定义和适用条件1.来源本方法学参考了下述CDM执行理事会批准的程序、方法学工具和指南：•CDM造林再造林项目活动基线情景识别和额外性论证的组合工具(V01,EB35)；•CDM造林再造林项目活动监测样地数量的计算工具(V02.1.0,EB58)；•CDM造林再造林项目活动导致的生物质燃烧引起的非CO2温室气体排放的估算工具(V04.0.0,EB65)；•CDM造林再造林项目活动导致的土壤有机碳储量变化的估算工具(V01.1.0,EB60)；•CDM造林再造林项目活动林木和灌木碳储量及其变化的估算工具(V03.0.0,EB70)；•CDM造林再造林项目活动枯死木和枯落物碳储量及其变化的估算工具(V02.0.0,EB67)；2.规范性引用文件除参考上述CDM执行理事会批准的最新版本的程序、方法学工具和指南外，下列文件及其更新版本对于本方法学的应用是必不可少的：•中华人民共和国《温室气体自愿减排交易管理暂行办法》（发改气候[2012]1668号）•国家林业局.碳汇造林技术规定（试行）（办造字[2010]84号）•中华人民共和国国家标准：造林技术规程（GB/T15776-2006）•中华人民共和国林业行业标准：毛竹林丰产技术（LY1059-92）•国家森林资源连续清查主要技术规定（林资发[2004]25号）3.定义本方法学及其应用采用下述定义：竹林：是指连续面积不小于1亩、郁闭度不低于20%、成竹竹秆高度不低于2米、竹秆胸径（或眉径）不小于2厘米的以竹类为主的植物群落。竹林是中国森林的一种类型。2小竹丛：是指成竹竹秆高度低于2米或竹秆胸径（或眉径）小于2厘米的任何竹类植物群落。小竹丛不属于森林范畴。大径散生竹林：指成竹竹秆高度大于6米、竹秆胸径（或眉径）大于5厘米的单轴散生型竹林。大径丛生竹林：指成竹竹秆高度大于6米、竹秆胸径（或眉径）大于5厘米的合轴丛生型竹林。小径散生竹林：指成竹竹秆高度大于6米、竹秆胸径（或眉径）2~5厘米的单轴散生竹林。小径丛生竹林：指成竹竹秆高度大于6米、竹秆胸径（或眉径）2~5厘米的合轴丛生竹林。复轴混生型竹林：指成竹竹秆高度大于6米、竹秆胸径（或眉径）大于5厘米的单轴和合轴混生的竹林。立竹度：指单位面积内正常生长的竹子（病死竹、倒伏竹除外）的数量。土壤扰动：是指导致土壤有机碳降低的活动，如整地、松土、翻耕、挖树桩（根）或竹篼等。基线情景：指在没有拟议的竹子造林项目活动时，最能合理地代表项目边界内土地利用和管理未来的可能情景。项目情景：指拟议的竹子造林项目活动下的土地利用和管理情景。碳库：包括地上生物量、地下生物量、枯落物、枯死木和土壤有机质。地上生物量：土壤层以上以干重表示的本本植被（包括竹类）活体的生物量，包括干、桩、枝、皮、种子、花、果和叶等。地下生物量：所有本本植被（包括竹类）活根的生物量。由于细根（直径≤1~2mm）通常很难从土壤有机成分或枯落物中区分出来，因此通常不包括该部分。枯落物：土壤层以上、直径小于5厘米、处于不同分解状态的所有死生物量，包括凋落物、腐殖质，以及不能从经验上从地下生物量中区分出来的活细根（直径≤1~2mm）。枯死木：枯落物以外的所有死生物量，包括枯立木、枯倒木以及直径大于或等于5厘米的枯枝、死根和树桩。土壤有机质：一定深度内（通常为100cm）矿质土和有机土（包括泥炭土）中的有机质，包括不能从经验上从地下生物量中区分出来的活细根。泄漏：指由拟议的竹子造林项目活动引起的、发生在项目边界之外的、可测量的温室气体源排放的增加。计入期：指项目情景相对于基线情景产生额外的温室气体减排量的时间区间。基线碳汇量：指在基线情景下，项目边界内碳库中碳储量变化之和。项目碳汇量：指在项目情景下，项目边界内所选碳库中碳储量变化量，减去由拟议的竹子造林项目活动引起的项目边界内温室气体排放的增加量。3项目减排量：指竹子造林项目活动引起净温室气体减排量，其大小等于项目碳汇量，减去基线碳汇量，再减去泄漏量。额外性：指拟议的竹子造林项目活动产生的项目碳汇量高于基线情景下的基线碳汇量的情形。这种额外的碳汇量在没有拟议的竹子造林项目活动时是不会产生的。4.适用条件本方法学适用于采用竹子进行造林的项目活动，其具体适用条件包括：(a)项目地不属于湿地。(b)如果项目地属下列情况之一，竹子造林或营林过程中对土壤的扰动不超过地表面积的10%：1)土壤为有机土；2)符合下列条件的草地：管理方式有机碳输入改良草地----中度放牧下的可持续利用，至少存在一种改良措施（施肥、草种改良、灌溉）未退化草地----非退化或可持续管理的草地，未实施改良措施中度退化草地----过牧或中度退化，相对于未退化草地，生产力较低，未实施改良措施高输入：实施了除改良草地的措施外的其他改良措施3)符合下列条件的农地：土地利用耕作方式有机碳输入短期作为农地、休（弃）耕地（休耕、弃耕期短于20年，或其他已生长多年生草本植物的闲置农地）。全耕----充分翻耕或频繁（年内）耕作导致强烈土壤扰动。在种植期地表残体盖度低于30%。高输入，且施用粪肥：在中等碳输入的农作系统中定期施用动物粪肥，使碳输入显著增加。减耕----对土壤的扰动较低（通常耕作深度浅，不充分翻耕）。在种植期地表残体盖度通常大于30%。高输入，且施用粪肥：在中等碳输入的农作系统中定期施用动物粪肥，使碳输入显著增加。高输入，不施用粪肥：在中等碳输入的农作系统中，由于种植产生大量作物残体的作物、使用绿肥、种植覆盖作物、植被休耕、灌溉、一年生作物轮作中频繁使用多年生草本等措施，使作物残体输入量明显增大。免耕----播种前不经初耕，仅中等输入：一年生谷类作物残体全部4在播种带上有最低限度的土壤扰动。一般使用杀虫剂控制杂草。返还农地。如果残体被收获，则补施有机肥（如粪肥），也包括施用矿质肥料或轮作固氮作物。高输入，且施用粪肥（同上）高输入，不施用粪肥（同上）长期农耕地（连续耕作20年以上，以一年生作物为主）免耕----播种前不经初耕，仅在播种带上有最低限度的土壤扰动。一般使用杀虫剂控制杂草。高输入，且施用粪肥（同上）(c)项目地适宜竹子生长，种植的竹子最低高度能达到2米，且竹秆胸径（或眉径）至少可达到2厘米，地块连续面积不小于1亩，郁闭度不小于0.20。(d)项目活动不采取烧除的林地清理方式（炼山），对土壤的扰动符合水土保持要求，如沿等高线进行整地，不采用全垦的整地方式。(e)项目活动不清除原有的散生林木。5第II部分.基线和碳计量方法1.项目边界竹子造林项目活动的“项目边界”是指，由拥有土地所有权或使用权的项目参与方实施的竹子造林项目活动的地理范围，也包括以竹子造林活动的产品为原材料生产的竹产品的使用地点。一个竹子造林项目活动可在若干个不同的地块上进行，但每个地块应有特定的地理边界，该边界不包括位于两个或多个地块之间的土地。项目边界包括事前项目边界和事后项目边界。事前项目边界是在项目设计和开发阶段确定的项目边界，是拟实施竹子造林项目活动的地理边界。事前项目边界可采用下述方法之一确定：(a)采用全球定位系统（GPS）、北斗卫星导航系统（Compass）或其他卫星导航系统，进行单点定位或差分技术直接测定项目地块边界的拐点坐标，单点定位误差不大于5米。(b)利用高分辨率的地理空间数据（如卫星影像、航片）、森林分布图、林相图等，在地理信息系统（GIS）辅助下直接读取项目地块的边界坐标。(c)使用大比例尺地形图（比例尺不小于1:10000）进行现场勾绘，结合GPS、Compass等定位系统进行精度控制。事后项目边界是在项目监测时确定的、项目核查时核实的、实际实施的项目活动的边界。事后项目边界可采用上述方法(a)或(b)进行，面积测定允许误差小于5%。在项目审定和核查时，项目参与方须提交地理信息系统（GIS）产出的项目边界的矢量图形文件（.shp文件）。在项目审定时，项目参与方须提供占项目活动总面积三分之二或以上的项目参与方的土地所有权或使用权的证据。在首次核查时，项目参与方须提供所有项目地块的土地所有权或使用权的证据，如县（含县）级以上人民政府核发的土地权属证书或其他有效的证明材料。2.土地合格性项目参与方须采用下述程序证明项目边界内的土地合格性1：(a)提供透明的信息证明，在项目开始时，项目边界内的土地符合下列所有条件:(i)植被状况不符合我国政府定义森林的阈值标准，即植被状况不同时满足下列所有条件：(1)郁闭度≥0.20，(2)树高≥2米，(3)面积≥1亩，(4)如果为竹类，竹秆胸径（或眉径）≥2厘米；(ii)如果地块上有天然或人工幼树，其继续生长不会达到我国政府定义森林的阈值标准；(iii)项目地块不属于因采伐或自然干扰而产生的临时的无林地（迹地）。1基于“证明CDM造林再造林项目活动土地合格性的程序(V01.0,EB35)”修改而来。6(b)提供透明的信息证明，自2005年2月16日起，项目活动所涉每个地块上的植被状况符合上述(a)(i)的条件。(c)为证明上述(a)和(b)，项目参与方须提供下列证据之一，以根据我国政府确定的森林定义标准，区分有林地和无林地，以及可能的土地利用方式的变化：(i)经过地面验证的高分辨率的地理空间数据（如卫星影像、航片）；(ii)森林分布图、林相图或其他林业调查规划空间数据；(iii)土地权属证或其他可用于证明的书面文件。如果没有上述(i)~(iii)的资料，项目参与方须呈交通过参与式乡村评估（PRA）方法获得的书面证据。3.碳库和温室气体排放源选择在项目边界内考虑的碳库如表1。本方法学对项目边界内的温室气体源排放源的选择如表2。表1竹子造林项目活动的碳库选择碳库考虑或不考虑理由或解释地上生物量考虑项目活动的主要碳库地下生物量考虑项目活动的主要碳库枯死木不考虑与基线情景相比该碳库不会降低，因此可保守地忽略不计。枯落物考虑或不考虑与基线情景相比该碳库会增加，但也可保守地选择不考虑该碳库。如果选择该碳库，则项目活动不允许移除地表枯落物。土壤有机碳考虑或不考虑与基线情景相比该碳库会增加，但也可保守地选择不考虑该碳库。如果选择该碳库，则项目活动不允许移除地表枯落物。竹产品考虑或不考虑与基线情景相比该碳库会增加，但也可保守地选择不考虑该碳库。表2项目边界内的温室气体排放源的选择排放源气体考虑或不考虑理由或解释木本植物（包括竹类）生物质燃烧CO2不考虑该CO2排放已在碳储量变化中考虑CH4考虑林地清理、整地或竹林经营过程中由于木本植被（包括竹子）生物质燃烧可引起显著的CH4排放N2O考虑林地清理、整地或竹林经营过程中由于木本植被（包括竹子）生物质燃烧可引起显著的N2O排放化石燃料燃烧CO2CH4N2O不考虑潜在排放量很小，可忽略不计7施肥N2O不考虑潜在排放量很小，可忽略不计4.计入期选择项目参与方须清晰地说明项目的开始日期、计入期和项目期，并解释选择该日期的理由。项目开始日期是指为种植竹子而开始的林地清理和整地的日期。项目开始日期不应早于2005年2月16日。如果项目开始日期早于向国家气候变化主管部门提交备案的日期，项目参与方须提供透明和可核实的证据，证明温室气体减排是项目活动最初的主要目的。这些证据须是发生在项目开始日之前的、官方的或有法律效力的文件。项目期是指实施项目活动的时间区间。计入期是指项目活动相对于基线情景产生额外温室气体减排量的时间区间，计入期的起止日期应与项目期相同。计入期按国家气候变化主管部门规定的方式确定，在颁布相关规定以前，计入期最短为20年，最长不超过30年。5.基线情景识别和额外性论证项目参与方须使用最新版“CDM造林再造林项目活动基线情景识别和额外性论证的组合工具”，来识别竹子造林项目活动的基线情景和论证项目活动的额外性。在使用该工具时，步骤0（STEP0）中的1999年12月31日替换为2005年2月16日。项目参与方也可选用下述“三重测试”程序来识别竹子造林项目活动的的基线情景并论证其额外性2：2.1.符合法律法规的要求项目参与方须证明发生在项目边界内的所有项目活动不会违反任何现有的法律、法规、规章以及其他强制性规定和技术标准。既包括国家级的法律法规和规章以及技术标准，也包括适用的省级和地方的规章以及技术标准。尚未通过的法律或规章则无须考虑。2.2.普遍性做法项目参与方须证明拟议的项目活动不是普遍性做法。如果没有与拟议的项目活动相类似的造林项目活动，该拟议的项目活动就被认为不是普遍性做法，其基线情景则为历史的或现有的土地利用情景。类似的造林项目活动指在项目所在区域、类似的社会经济和生态条件下、普遍实施的与拟议的项目活动相类似的造林活动，包括那些由具有可比性的实体或机构（如大公司、小公司、国家政府项目、地方政府项目等）实施的造林项目活动和那些具有可比性的地理范围、地理位置、环境条件、社会经济条件、制度框架以及投资环境下的造林项目活动，也包括2005年2月16日以前制定的土地利用规划方案。如果项目参与方无法证明拟议的项目活动不是普遍性做法，或者存在与拟议的项目活动相类似的造林项目活动（即拟议的项目活动属于普遍做法），项目参与方须通过下文2.3节的障碍分析，来确定拟议项目的基线情景并证明拟议的项目的额外性。2基于“熊猫标准农林业及其他土地利用行业细则”中的“三重测试”程序修改而来。8项目活动一旦被认定不是普遍性做法，即被认定为在其计入期内具有额外性，并可略去进行下文2.3节的障碍分析。2.3.实施障碍如果拟议的项目活动属于普遍性做法，项目参与方仍可通过实施障碍分析来确定项目活动的基线情景并证明项目活动的额外性，例如由于项目参与方面临相关的障碍，阻碍其在项目区实施通常做法或原有的土地利用规划方案，使得基线情景为维持原有的土地利用方式。实施障碍是指任何可能阻止项目活动开展的因素。项目参与方至少需要对下列三种障碍之一进行评估：财务障碍、技术障碍或机构障碍。项目参与方可以证明存在多种障碍，但只要证明一种障碍存在即可。•财务障碍可以包括高成本、有限的资金，或者在没有项目活动温室气体减排量收益时，内部收益率低于项目参与方预期能接受的最低收益率。如果采用财务障碍测试，项目参与方须提供可靠的定量分析的证据，如净现金流和内部收益率测算，以及相关批准文件等书面材料。•技术障碍包括缺少必需的材料（如种植材料），缺乏有技能的和接受过良好培训的劳动力，缺少法律、传统、市场条件和实践措施等相关知识，缺少实践经验等。•机构障碍包括对技术实施的制度性排斥，技术实施能力不足，管理层缺乏共识等。6.碳层划分如果项目区自然和社会经济条件以及项目活动差异较大，须对项目区进行分层，以提高在一定可靠性下的监测和估计的精度，并降低监测成本。碳层划分包括基线碳层划分和项目碳层划分。基线碳层划分的目的是为了分别基线碳层确定基线情景和估计基线碳汇量。项目参与方可根据项目边界内地块上的主要植被状况（如散生木（竹）盖度和年龄、灌木植被（包括小竹丛）的种类和盖度）和土地利用类型（农地、宜林荒山等）来划分基线碳层。项目碳层划分包括事前项目碳层划分和事后项目碳层划分。事前项目碳层用于项目碳汇量的事前估计，主要根据竹子造林和竹林经营管理计划来划分。事后项目碳层用于项目碳汇量的事后估计，主要根据竹子造林和竹林经营管理实际发生的情况来划分。但是，无论是事前分层还是事后分层，多个竹种可合并为一个碳层，不同时间营造的竹林也可合并为一层，关键是看其是否具有近似的碳储量、相同的计量参数（如生物量生长速率、地下生物量与地上生物量之比、含碳率等）和生物量异速生长方程等，其目的是降低层内变异性，增加层间变异性，从而降低在一定精度要求下所需监测的样地数量。如果发生自然或人为干扰（如火灾、毁林）导致项目的异质性增加，在每次监测和核查时的事后分层调整时均须考虑这些因素的影响。项目参与方可使用项目开始时和发生干扰时的卫星影像来进行碳层划分。7.基线碳汇量根据本方法学的适用条件，基线碳汇量可假定为零，即0,=ΔtBSLC：9t1,2,3,…t竹子造林项目活动开始后的年数（年）8.项目碳汇量项目碳汇量是指在拟议的竹子造林项目活动的情景下，项目边界内所选碳库中碳储量变化量，减去由竹子造林项目活动引起的温室气体排放的增加量，采用下式计算：tEtPtACTUALGHGCC,,,−Δ=Δ(1)式中:tACTUALC,Δ第t年项目碳汇量（tCO2-e.a-1）tPC,Δ第t年项目边界内所选碳库中碳储量年变化量（tCO2-e.a-1）tEGHG,第t年项目活动引起的温室气体排放的年增加量（tCO2-e.a-1）t1,2,3,…t竹子造林项目活动开始后的年数（年）采用下述公式计算项目边界内所选碳库中碳储量的年变化量tPROJHWPtALSOCtPROJLItPROJSHRUBtPROJBAMBOOtPCCCCCC,_,_,_,_,_,Δ+Δ+Δ+Δ+Δ=Δ(2)式中:tPC,Δ第t年项目边界内所选碳库中碳储量的年变化量（tCO2-e.a-1）tPROJBAMBOOC,_Δ项目情景下，第t年项目边界内营造的竹林生物质碳储量的年变化量（tCO2-e.a-1）tPROJSHRUBC,_Δ项目情景下，第t年项目边界内灌木生物质碳储量的年变化量（tCO2-e.a-1）。针对本方法学，灌木包括小竹丛。tPROJLIC,_Δ项目情景下，第t年项目边界内枯落物碳储量的年变化量（tCO2-e.a-1）。对于集约经营的竹林，枯落物碳储量的变化为零。tALSOCC,_Δ项目情景下，第t年项目边界内土壤有机碳储量的年变化（tCO2-e.a-1）。对于集约经营的竹林，土壤有机碳储量的变化为零。tPROJHWPC,_Δ项目情景下，第t年收获的竹材生产的竹产品中碳储量的变化量（tCO2-e.a-1）t1,2,3,…t竹子造林项目活动开始后的年数（年）8.1竹林生物质碳储量变化量（tPROJBAMBOOC,_Δ）的事前估算10tBBPROJBAMBOOtABPROJBAMBOOtPROJBAMBOOCCC,,_,,_,_Δ+Δ=Δ（3）式中：tABPROJBAMBOOC,,_Δ第t年项目边界内营造的竹林地上生物质碳储量的年变化量（tCO2-e.a-1）tBBPROJBAMBOOC,,_Δ第t年项目边界内营造的竹林地下生物质碳储量的年变化量（tCO2-e.a-1）8.1.1地上生物质碳储量变化竹林生长发育分为竹林发育成林阶段（大径散生竹林1-9年，小径散生竹林1-5年，丛生竹1-5年，混生竹1-6年）和竹林成林稳定阶段（大径散生竹林从第10年开始，小径散生竹林从第6年开始，丛生竹第6年开始；混生竹从第7年开始）。达到竹林成林稳定阶段后，由于择伐或自然枯损以及新竹的生长，竹林地上生物量达到动态平衡状态。对于事前估计，根据可获得的数据情况，可从下列方法中选择其中一种方法估算发育成林阶段的地上生物质碳储量的变化。方法I：如果有拟营造的竹林单位面积生物量随竹林年龄变化的相关方程，则可直接用该方程估算造林后各年度的生物质碳储量和碳储量变化量，直到进入竹林成林稳定阶段为止。此后，假定竹林地上生物质碳储量变化量为零。方法II：根据拟营造的竹林的平均胸径、平均高度与竹林年龄的相关方程，再结合单株生物量方程计算平均单株地上生物量，即：)(,ajDBHjtfDBH=(4)),(,DBHtfHajHj=(5)),(,,jjABjABBAMBOOHDBHfB=(6)式中：jDBH发育到at时，竹林平均胸径（厘米)jH发育到at时，竹林平均高（米)),(jjABHDBHf竹种（组）j的地上生物量方程（生物量与直径（胸径、眉径、地径）和竹高的相关方程）jABBAMBOOB,,竹种（组）j的平均单株地上生物量（kg.单株-1）at竹林年龄（年）j竹种或竹种组11然后结合立竹度与竹林年龄的相关方程，计算单位面积竹林地上生物质碳储量：3,,,,,101244,,−⋅⋅⋅⋅⋅=BjjtjtjABBAMBOOBAMBOOCFMNBCtjAB(7))(,,ajNtjtfN=(8)式中：tjABBAMBOOC,,单位面积竹林地上生物质碳储量（tCO2-e.hm-2）tjABBAMBOOB,,,,平均单株地上生物量（Kg.株-1）tjN,对散生或混生竹种，为每公顷立竹度（株.hm-2)；对丛生竹种，为平均每丛的株数（株.丛-1）BiCF,竹子含碳率（tC(td.m.)-1)jM对散生竹取值为1；对丛生竹，为每公顷丛数（丛.hm-2)t项目开始后的年数（年）j竹种或竹种组at竹林年龄（年）；atta−=，其中a为造林发生的年份则营造的竹林地上生物质碳储量的变化为：⎪⎩⎪⎨⎧=Δ≤−⋅>−∑∑jmequilibriuatjiABBAMBOOtjiABBAMBOOtjiBamboojmequilibriuaTtCCATtijtABPROJBAMBOOC,1,,,,,,,,,,,,)(0,,_当当(9)式中:tABPROJBAMBOOC,,_Δ第t年项目边界内营造的竹林地上生物质碳储量的年变化量（tCO2-e.a-1）tjiBambooA,,,第t年第i碳层j竹种（组）的面积（hm2）tjiABBAMBOOC,,,,第t年第i碳层j竹种（组）单位面积竹林地上生物质碳储量（tCO2-e.hm-2）1,,,,−tjiABBAMBOOC第（t-1）年时，第i碳层j竹种（组）单位面积竹林地上生物质碳储量（tCO2-e.hm-2）jmequilibriuT,第j竹种（组）竹林到达成林稳定阶段所需的时间（年）at林龄（年）；atta−=，其中a为造林发生的年份t项目开始后的年数（年）方法III：如果没有上述方法I和方法II所需数据，可采用在达到成林稳定前，按平均生长速率计算，即：12⎪⎩⎪⎨⎧=Δ≤•⋅>∑∑jmequilibriuajmequilibriujimequilibriuABBAMBOOtjiBamboojmequilibriuaTtTCATtijtABPROJBAMBOOC,,,,,,,,,0,,_当当(10)式中:tABPROJBAMBOOC,,_Δ第t年项目情景下竹类地上生物质碳储量的年变化量（tCO2-e.a-1）tjiBambooA,,,第t年第i碳层j竹种（组）的面积（hm2）jmequilibriuABBAMBOOC,,j竹种到达成林稳定阶段时的单位面积地上生物质碳储量（tCO2.hm-2）at林龄（年）；atta−=，其中a为造林发生的年份jmequilibriuT,j竹种到达成林稳定阶段所需的时间（年）t1,2,3,…t竹子造林项目活动开始后的年数（年）8.1.2地下生物质碳储量变化竹林地下生物量通常随着竹林年龄的增加而增加。竹子择伐经营时，通常只移除地上部分（竹秆、竹枝、竹叶），而地下部分（竹蔸、竹根和竹鞭）仍留存于林地中，因此即使竹林到达成林稳定年限后，其地下生物量碳储量通常还会继续增加。竹林地下生物质碳储量的变化可通过竹林地下生物量与地上生物量之比和地上生物质碳储量变化计算，即：tjiBAMBOOtjtjiABPROJBAMBOOtjijtjiABPROJBAMBOOtBBPROJBAMBOOARCRCCaa,,,1,1,,,,_,,,,,_,,_)(⋅⋅−⋅=Δ−−∑∑(11)式中:tBBPROJBAMBOOC,,_Δ第t年项目边界内营造的竹林地下生物质碳储量的年变化量（tCO2-e.a-1）tjiBambooA,,,第t年第i碳层j竹种（组）的面积（hm2）tjiABBAMBOOC,,,,第t年第i碳层j竹种（组）单位面积竹林地上生物质碳储量（tCO2-e.hm-2）1,,,,−tjiABBAMBOOC第（t-1）年时，第i碳层j竹种（组）单位面积竹林地上生物质碳储量（tCO2-e.hm-2）atjR,j竹种（组）在竹林年龄为at时的地下生物量与地上生物量之比1,−atjRj竹种（组）在竹林年龄为)1(−at时的地下生物量与地上生物量之比at林龄（年）；atta−=，其中a为造林发生的年份t项目开始后的年数（年）13如果项目参与方没有竹子地下生物量与地上生物量之比随竹林年龄变化的相关关系，则可假定地下生物量与地上生物量之比为常数。在这种情况下，当竹林到达成林稳定阶段后，地上和地下生物质碳储量的变化均为零。8.2灌木生物质碳储量的变化量（tPROJSHRUBC,_Δ）3对于事前计量，可假定灌木生物质碳储量变化为零。对事后监测和计量。根据灌木盖度对项目边界内的灌木生物量进行分层，并估算每层灌木生物量的碳储量。假定一段时间内（第t1至t2年）灌木生物量的变化是线性的，基线灌木生物质碳储量的年变化量（∆CSHRUB_PROJ,,t）计算如下：∑⎟⎟⎠⎞⎜⎜⎝⎛−−=ΔitiPROJSHRUBtiPROJSHRUBtPROJSHRUBttCCC12,,_,,_,_12(12)式中：tPROJSHRUBC,_Δ第t年项目情景下灌木生物质碳储量的年变化量（tCO2-e·a-1）tiPROJSHRUBC,,_第t年i灌木碳层灌木生物质碳储量（tCO2-e）i1,2,3,……灌木碳层t1,2,3,……自项目开始以来的年数t1,t2项目开始以后的第t1年和第t2年，且t1≤t≤t2采用下式计算第t年i灌木碳层内灌木生物质碳储量：1244)1(,,_,,_,,_⋅⋅⋅+⋅=stiPROJSHRUBstiPROJSHRUBtiPROJSHRUBCFARBC（13）式中：tiPROJSHRUBB,,_第t年i灌木碳层灌木的平均每公顷地上生物量（td.m.hm-2）sR灌木的地下生物量与地上生物量之比（无量纲）tiPROJSHRUBA,,_第t年i灌木碳层的面积(hm-2)sCF灌木生物量中的含碳率（tC(t.d.m.)-1），缺省值为0.47i1,2,3,……灌木碳层t1,2,3,……自项目开始以来的年数44/12将C转换为CO2的分子量比值3参考“CDM造林再造林项目活动林木和灌木碳储量及其变化的估算工具”14灌木平均每公顷生物量的估算方法如下：z灌木盖度<5%时，灌木平均每公顷生物量视为0；z灌木盖度≥5%时，按下列方式进行估算：tiPROJSHRUBFORESTSFtiPROJSHRUBCCBBDRB,,_,,_⋅⋅=（14）式中：SFBDR灌木盖度为1.0时的每公顷灌木生物量与拟议项目所在地区完全郁闭森林每公顷地上生物量之比（无量纲）FORESTB拟议项目所在地区完全郁闭人工林平均每公顷地上生物量（td.m.hm-2）tiPROJSHRUBCC,,_第t年i灌木碳层的灌木盖度，以小数表示（如盖度为10%，则tiPROJSHRUBCC,,_=0.10）（无量纲）i1,2,3,……灌木碳层t1,2,3,……自项目开始以来的年数8.3枯落物碳储量的变化量（tPROJLIC,_Δ）4假定一段时间内枯落物碳储量的年变化量为线性，则一段时间内枯落物碳储量的平均年变化量采用下式计算：∑⎟⎟⎠⎞⎜⎜⎝⎛−−=ΔitiPROJLItiPROJLItPROJLIttCCC12,,_,,_,_12（15）LItiPROJBAMBOOtiPROJLIDFCC⋅=,,_,,,（16）式中：tPROJLIC,_Δ第t年项目情景下枯落物碳储量的年变化量（tCO2-e·a-1）tiPROJLIC,,_第t年i项目碳层的枯落物碳储量（tCO2-e）tiPROJBAMBOOC,,_第t年第i项目碳层的竹林生物质碳储量（tCO2-e）LIDF项目所在地区竹林枯落物碳储量与其活生物质碳储量的比值（无量纲）t1,t2项目开始以后的第t1年和第t2年，且t1≤t≤t2i1,2,3,…，项目的项目碳层4参考“CDM造林再造林项目活动枯死木和枯落物碳储量及其变化的估算工具(V02.0.0,EB67)”15对于集约经营的竹林，枯落物碳储量的变化为零，即0,_=ΔtPROJLIC。8.4土壤有机碳储量的变化量（tALSOCC,_Δ）在估算土壤有机碳储量变化时，本方法学基于以下假设：z项目整地和造林活动在同一年进行；z项目活动的实施将使项目地块的土壤有机碳含量从项目开始前的初始水平提高到相当于天然森林植被下土壤有机碳含量的稳态水平，大约需要20年时间；z从造林活动开始后的20年间，项目情景下土壤有机碳储量的增加线性的；首先确定项目开始前各项目碳层土壤有机碳含量初始值（SOCINITIAL,i）。项目业主或其他项目参与方可以通过国家规定的标准操作程序直接测定项目开始前各碳层的SOCINITIAL,i；也可以采用下列方法估算项目开始前各碳层的SOCINITIAL,i：,,,,,=∗∗∗INTITIALiREFiLUiMGiINiSOCSOCfff（17）式中：SOCINITIAL,i项目开始时，第i项目碳层的土壤有机碳储量（tC.hm-2）SOCREF,i与第i项目碳层具有相似气候、土壤条件的当地自然植被（如当地未退化的、未利用土地上的自然植被）下土壤有机碳储量的参考值（tC.hm-2）fLU,i第i项目碳层与基线土地利用方式相关的碳储量变化因子（无量纲）fMG,i第i项目碳层与基线管理模式相关的碳储量变化因子（无量纲）fIN,i第i项目碳层与基线有机碳输入类型（如农作物秸杆还田、施用肥料）相关的碳储量变化因子（无量纲）i1,2,3,…，项目的林木分层SOCREF,i、fLU,i、fMG,i和fIN,i的取值，可参考本方法学中的参数表。如果选取其它不同的数值，须提供透明和可核实的信息来证明。确定第i项目碳层的造林时间（即由于整地发生土壤扰动的时间，tPREP,i）。对于项目开始以后的第t年，如果：zt≤tPREP,i，则第t年时第i项目碳层的土壤有机碳储量的年变化（dSOCt,i）为0；ztPREP,i<t≤tPREP,I+20，则：16,,,20−=REFiINITIALitiSOCSOCdSOC（18）式中：dSOCt,i第t年i项目碳层的土壤有机碳储量的年变化率（tC.hm-2.a-1）SOCREF,i与项目第i项目碳层具有相似气候、土壤条件的当地自然植被（如当地未退化的、未利用土地上的自然植被）下土壤有机碳储量的参考值（tC.hm-2）SOCINITIAL,i项目开始时，第i层的土壤有机碳库碳储量（tC.hm-2）i1,2,3,…，项目碳层20假定项目地块的土壤有机碳含量从初始水平提高到相当于当地自然植被下土壤有机碳含量的稳态水平需要20年时间由于本方法学采用了基于碳储量变化因子的估算方法。考虑到其精度的不确定性和内在局限性，实际计算过程中土壤有机碳储量的年变化率不超过0.8tC.hm-2.a-1，即：如果dSOCt,i>0.8tC.hm-2.a-1，则dSOCt,i=0.8tC.hm-2.a-1（19）第t年时，项目所有碳层的土壤有机碳储量变化采用下式计算：(),,,144112=Δ=∗∗∗∑ALttitiiSOCAdSOCa（20）式中：∆SOCAL,t第t年时项目情景下土壤有机碳储量的年变化量（tCO2-e·a-1）dSOCt,i第t年i项目碳层的土壤有机碳储量年变化率（tC.hm-2.a-1）At,i第t年i项目碳层的土地面积(hm2)i1,2,3,…，项目碳层t1,2,3,…，项目开始以后的时间1a1年对于集约经营的竹林，土壤有机碳储量的变化为零，即∆SOCAL,t=0。8.5收获竹产品的碳储量变化（tPROJHWPC,_Δ）如果竹子造林项目活动有择伐情况发生，择伐的部分竹材中的碳将以竹产品的形式储存一定时间，而不是立即排放到大气中。对于散生竹类人工林，择伐通常从造林后第8~9年开始，对于丛生竹类人工林，择伐通常从造林后第4~5年开始。我国竹材除传统上用于农业生产和生活工具外，目前主要用于生产竹材人造板，包括竹编胶合板、竹材胶合板、竹材层压板、竹席竹帘胶合板、竹材纤维板和竹材刨花板等，产品17广泛应用于我国的汽车、火车、建筑、集装箱等工业部门。竹木复合人造板和造纸也是当前利用的一种趋势。竹林到达成林稳定阶段后，收获竹材生产的竹产品（HWP）中的碳将是主要的碳汇来源。本方法学假定HWP碳储量的长期变化，等于在产品生产后30年仍在使用和进入垃圾填埋的HWP中的碳量，而其他部分则假定在生产竹产品时立即排放，采用下述公式计算5：tytytytStemBAMBOOtPROJHWPOFBUCC⋅⋅=Δ∑,,,_(21)))2ln(30(tyLTtyeOF⋅−=(22)式中：tPROJHWPC,_Δ第t年项目产生的竹产品碳储量的年变化量（tCO2-e.a-1）tStemBAMBOOC,,第t年项目采伐的竹秆生物质碳储量（tCO2-e）。如果采伐的竹子是以竹秆鲜重计，则应将鲜重通过含水率换算成干重，然后转化为CO2的量。如果采伐利用整株竹子（包括枝和叶），则为地上生物量中的碳储量。tyBU竹子采伐用于ty类竹产品的利用率（%），即竹产品生物量占采伐收获量的百分比。tyOF根据IPCC一阶指数衰减函数确定的、ty类竹产品在生产后30年仍在使用或进入垃圾填埋的比例（无量纲）。ty竹产品种类tyLTty类竹产品的使用寿命(年)t1,2,3,…t竹子造林项目活动开始后的年数（年）8.4项目边界内温室气体排放的估计对于项目事前估计，由于无法预测项目边界内的火灾发生情况，因此可以不考虑森林火灾造成的项目边界内温室气体排放，即GHGE,t=0。对于项目事后估计，由于竹子造林项目活动引起的项目边界内的温室气体排放的增加为6：:tLIEtBAMBOOEtEGHGGHGGHG,,,,,+=(23)式中:EGHG第t年由于竹子造林项目活动的实施引起的项目边界内温室气体排放的增加量(tCO2-e.a-1)5根据2006IPCC国家温室气体清单指南中的一阶衰减函数修改而来。6参考“CDM造林再造林项目活动导致的生物质燃烧引起的非CO2温室气体排放的估算工具”18tBAMBOOEGHG,,第t年项目边界内火灾导致的竹子地上生物质燃烧引起的非CO2温室气体排放的增加量(tCO2-e.a-1)。tLIEGHG,,第t年项目边界内火灾导致的竹林枯落物燃烧引起的非CO2温室气体排放的增加量(tCO2-e.a-1)。t1,2,3,…t竹子造林项目活动开始后的年数(年)火灾引起竹林地上生物质燃烧造成的非CO2温室气体排放，使用最近一次项目核查时各碳层竹林地上生物量数据和燃烧因子进行计算。第一次核查时，无论自然或人为原因引起竹林火灾，其非CO2温室气体排放量都假定为0。()29825001.024,,,,,_⋅+⋅⋅⋅⋅⋅=∑ONCHitiBAMBOOtiburntBAMBOOEEFEFCOMFbAGHGlast(24)式中：tBAMBOOEGHG,,第t年项目边界内火灾导致的竹子地上生物质燃烧引起的非CO2温室气体排放的增加量(tCO2-e.a-1)。tiburnA,,第t年i项目碳层发生火烧的面积（hm2）lasttiBAMBOOb,,火灾发生前，项目最近一次核查时第i项目碳层的竹子地上生物量（td.m·hm-2），详见8.1.1节。如果只是发生地表火，即竹子地上生物量未被燃烧，则lasttiBAMBOOb,,设定为0COMF竹林燃烧系数（无量纲）4CHEFCH4排放因子（gCH4·(kg燃烧的干物质)-1）ONEF2N2O排放因子（gN2O·(kg燃烧的干物质)-1）25CH4的全球增温潜势，用于将CH4转换成CO2当量298N2O的全球增温潜势，用于将N2O转换成CO2当量i1,2,3……第i项目碳层t1,2,3……项目开始以后的年数；年（a）0.001将kg转换成t的常数森林火灾引起枯落物燃烧造成的非CO2温室气体排放，应使用最近一次核查的枯落碳储量来计算。第一次核查时由于火灾导致枯落物燃烧引起的非CO2温室气体排放量设定为0，之后核查时的非CO2温室气体排放量采用下式计算：∑⋅⋅⋅=itiLItiburntLIElastCAGHG,,,,,_124407.0(25)式中：19tLIEGHG,,第t年项目边界内火灾导致的竹林枯落物燃烧引起的非CO2温室气体排放的增加量(tCO2-e.a-1)。tiburnA,,第t年i项目碳层发生火烧的面积（hm2）lasttiLIC,,火灾发生前，项目最近一次核查时第i项目碳层的枯落物单位面积碳储量，使用第8.3节的方法计算（tCO2-e.hm-2）i1,2,3……第i项目碳层t1,2,3……项目开始以后的年数；年（a）0.07常数，非CO2排放量占CO2排放量的比例9.泄漏根据本方法学的适用条件，项目实施可能引起的项目前农业活动的转移，以及项目活动中使用运输工具和燃油机械造成的排放，均可忽略不计。因此0=tLK10.项目减排量竹子造林项目活动引起的项目减排量等于项目碳汇量，减去基线碳汇量，再减去泄漏量，即：ttBSLtACTUALtCHINAARLKCCC−Δ−Δ=Δ,,,_（26）式中:tCHINAARC,−Δ第t年项目减排量（tCO2-e.a-1）tACTUALC,Δ第t年项目碳汇量（tCO2-e.a-1）ΔtBSLC,第t年基线碳汇量（tCO2-e.a-1）tLK第t竹子造林项目活动引起的泄漏量(tCO2-e.a-1)t1,2,3,…t竹子造林项目活动开始后的年数(年)20第III部分.监测程序除非下面的监测变量表中另有要求，所有数据，包括本方法学所用工具中要求的监测项，均须按相关标准进行全面监测和测定。监测过程中收集的所有数据都须以电子版和纸质方式存档，直到计入期结束后至少两年。1.项目实施监测1.1基线碳汇量的监测基线碳汇量在事前确定，计入期内不再对其进行监测。1.2项目边界的监测z采用全球定位系统（GPS）、北斗卫星导航系统（Compass）或其他卫星导航系统，进行单点定位或差分技术直接测定项目地块边界的拐点坐标。也可利用高分辨率的地理空间数据（如卫星影像、航片），在地理信息系统（GIS）辅助下直接读取项目地块的边界坐标。在监测报告中说明使用的坐标系，使用仪器设备的精度；z检查实际边界坐标是否与竹子造林项目设计文件中描述的边界一致；z如果实际边界位于项目设计文件描述的边界之外，则位于项目设计文件确定的边界外的部分将不计入项目边界中；z将测定的拐点坐标或项目边界输入地理信息系统，计算项目地块及各碳层的面积；z在计入期内须对项目边界进行定期监测，如果项目边界发生任何变化，例如发生毁林，应测定毁林的地理坐标和面积，并在下次核查中予以说明。毁林部分地块将移出项目边界之外，在以后不再进行监测。同样，如果某些地块由于某种原因造林失败，并代之以其他土地利用方式，这些地块也可移出项目边界外，且不再进行监测和核查。但是已移出项目边界的地块，在以后不能再纳入项目边界内。而且，如果移出项目边界的地块以前进行过核查，其前期经核查的碳储量应保持不变，纳入碳储量变化的计算中。1.3竹林营造林活动的监测z林地清理和整地的监测：时间、地点（边界）、面积、清理和整地的方式和规格；z造林和幼林管护活动监测：造林和管护的方式、时间、地点、面积、竹种等；z竹林经营管理监测：择伐、松土、除草、施肥等活动的时间和地点；z确保竹子造林和营林各项活动符合本方法学的适用条件。项目参与方须在项目文件中描述，项目的营造林活动及其监测，符合中国竹子营造林相关的技术要求和森林资源清查的技术规范。项目参与方在其监测活动中须制定21标准操作程序（SOP）及质量保证和质量控制程序（QA/QC），包括野外数据的采集、数据记录、管理和存档。最好是采用国家森林资源清查或IPCC指南中的标准操作程序。2.抽样设计和碳层划分2.1碳层更新由于下述原因，每次监测时须对事前或上一次监测划分的碳层进行更新：•实际的竹子造林活动（如造林时间和竹种配置）可能与项目设计发生偏离；•计入期内可能发生无法预计的干扰（如林火），从而增加碳层内的变异性；•竹林经营管理活动（如择伐、施肥、翻耕）活动影响了项目碳层内的均一性；•发生土地利用变化（项目地转化为其他土地利用方式）；•过去的监测发现层内碳储量和碳储量变化的变异性：可将变异性太大的碳层细分为两个或多个碳层，或者将碳储量和碳储量变化及其变异性相近的两个或多个碳层合并为一个碳层；•某些事前或前一次监测划分的碳层可能不复存在。2.2抽样设计竹林生物质碳储量必须基于监测样地的测定。枯落物、土壤有机碳不需进行大量野外测定，可直接采用相关工具进行计算。项目参与方须基于固定样地的连续测定方法，采用碳储量变化法，测定和估计竹林生物质碳库中碳储量的变化。项目监测所需的样地数量，可以采用如下方法进行计算：(1)根据公式（27）计算。如果得到n≥30，则最终的样地数即为n值；如果n<30，则需要采用自由度为n-1时的t值，运用公式（27）进行第二次迭代计算，得到的n值即为最终的样地数；22222VALiiiVALiiiNtwsnNEtws∗∗∗=∗+∗∗⎛⎞⎜⎟⎝⎠∑∑（27）式中：n项目边界内估算生物质碳储量所需的监测样地数量；无量纲N项目边界内监测样地的抽样总体，N=A/Ap，其中A是项目总面积（hm2），Ap是样地面积（一般为0.0667hm2）（无量纲）tVAL可靠性指标。在一定的可靠性水平下，自由度为无穷（∞）时查t分布双侧t分位数表的t值（无量纲）wi项目边界内第i项目碳层的面积权重，wi=Ai/A，其中A是项目总面积（hm2），Ai是第i碳层的面积（hm2）（无量纲）si项目边界内第i项目碳层生物质碳储量估计值的标准差（tC.hm-2）22E项目生物质碳储量估计值允许的误差范围（即置信区间的一半），在每一碳层内用si表示(tC.hm-2)i1,2,3……项目碳层i(2)当抽样面积较大时（抽样面积大于项目面积的5％），按公式(26)进行计算获得样地数n之后，按公式(28)对n值进行调整，从而确定最终的样地数（na）：11/annnN=∗+（28）式中：na调整后项目边界内估算生物质碳储量所需的监测样地数量（无量纲）n项目边界内估算生物质碳储量所需的监测样地数量（无量纲）N项目边界内监测样地的抽样总体（无量纲）(3)当抽样面积较小时（抽样面积小于项目面积的5％），可以采用简化公式（29）计算：22VALiiitnwsE⎛⎞⎛⎞=∗∗⎜⎟⎜⎟⎝⎠⎝⎠∑（29）式中：n=项目边界内估算生物质碳储量所需的监测样地数量；无量纲tVAL=可靠性指标。在一定的可靠性水平下，自由度为无穷（∞）时查t分布双侧t分位数表的t值；无量纲wi=项目边界内第i碳层的面积权重；无量纲si=项目边界内第i碳层生物质碳储量估计值的标准差；tC·hm-2E=项目生物质碳储量估计值允许的误差范围（即置信区间的一半），在每一碳层内用si表示；tC·hm-2i=1,2,3……项目碳层i(4)分配到各层的监测样地数量，采用最优分配法按公式（30）进行计算：iiiiiiwsnnws∗=∗∗∑（30）式中：ni项目边界内第i碳层估算生物质碳储量所需的监测样地数量（无量纲）n项目边界内估算生物质碳储量所需的监测样地数量（无量纲）wi项目边界内第i碳层的面积权重（无量纲）si项目边界内第i碳层生物质碳储量估计值的标准差（tC.hm-2）i1,2,3……项目碳层i23在各项目碳层内，样地的空间分配采用随机起点、系统布点的布设方案。项目参与方须确定首次监测和核查的时间以及间隔期。监测和核查的间隔期为3~10年。2.3竹林碳储量变化测定竹林生物质碳储量的测定和计算步骤如下7：第一步：测定样地内所有竹秆的胸径（DBH）和高度（H）。该胸径也可以用眉径或基径替代，取决于采用的生物量方程中变量的含义。项目边界内原有的散生木不包括在测定的范围内。第二步：利用单株地上生物量方程（),(jjABHDBHf）计算每株竹子的地上生物量。然后可通过地下生物量与地上生物量之比计算全株竹子的生物量（地上生物量乘以)1(,atjR+）。如果没有地上生物量方程，可用全株生物量方程直接计算全株竹子的生物量。将单株竹子生物量累积到样地的生物量和碳储量。所选用的生物量方程须根据“证明估算CDM造林再造林项目活动林木地上生物量的异速生长方程的适用性的程序”证明其适用性。第三步：计算碳层单位面积平均碳储量pinptipBAMBOOtiBAMBOOAnCCi⋅=∑=1,,,,,(31))1(21,,,12,,,2,−×⎟⎟⎠⎞⎜⎜⎝⎛−×=∑∑==iinptipBAMBOOnptipBAMBOOitinnCCnsii(32)式中:tiBAMBOOC,,第t年时，i项目碳层单位面积竹林生物质碳储量（tCO2.hm-2）tipBAMBOOC,,,第t年时，i项目碳层p样地单位面积竹林生物质碳储量（tCO2.hm-2）ini项目碳层的样地数量2,tis第t年时，i项目碳层单位面积竹林生物质碳储量的方差(tCO2.hm-2)2pA样地面积（hm2）第四步：计算项目边界内单位面积竹林生物质碳储量及其方差:∑=×=MitiBAMBOOitBAMBOOCwC1,,,(33)7参考“CDM造林再造林项目活动林木和灌木碳储量及其变化的估算工具”24itiMiiCnswstBAMBOO2,122,×=∑=(34)式中:tBAMBOOC,第t年时，项目边界内单位面积竹林生物质碳储量（tCO2.hm-2）iw碳层i在项目总面积中的面积权重（无量纲）tiBAMBOOC,,第t年时，i项目碳层单位面积竹林生物质碳储量（tCO2.hm-2）2,tBAMBOOCs第t年时，项目单位面积竹林生物质碳储量的方差(tCO2.hm-2)22,tis第t年时，i项目碳层单位面积竹林生物质碳储量的方差(tCO2.hm-2)2ini项目碳层的样地数量M项目碳层数量第五步：计算项目边界内竹林生物质碳储量:tBAMBOOtPROJBAMBOOCAC,,,×=(35)式中:tPROJBAMBOOC,,第t年时，项目边界内竹林生物质碳储量（tCO2）A项目总面积（hm2）tBAMBOOC,第t年时，项目边界内单位面积竹林生物质碳储量（tCO2.hm-2）第六步：计算项目单位面积竹林生物质碳储量的不确定性：tBAMBOOCVALtBAMBOOCstUNCtBAMBOO,,,×=(36)式中:tBAMBOOUNC,以抽样调查的相对误差限（%）表示的项目单位面积竹林生物质碳储量的不确定性（%）VALt可靠性指标：通过危险率（1-置信度）和自由度（n-M）查t分布的双侧分位数表，其中n为项目样地总数，M为项目碳层数量。例如:置信度90%，自由度为45时的可靠性指标可在excel中用“=TINV(0.10,45)”8计算得到1.6794。tBAMBOOCs,项目单位面积竹林生物质碳储量的方差的平方根，即平均值的标准误（tCO2.hm-2）。8在EXCEL2010中采用了T.INV()，而不是TINV()。253.精度控制和校正本方法学仅要求对营造的竹子生物量的监测精度进行控制，要求达到90%可靠性水平下90%的精度。如果不确定性%10,>tBAMBOOUNC，项目参与方可通过增加样地数量，从而使测定结果达到精度要求，也可以选择下述打折的方法。())1(1221,,,,,,,DRCCCtPROJBAMBOOtPROJBAMBOOttPROJBAMBOO−×−=Δ(37)式中:21,,,ttPROJBAMBOOCΔ时间区间t1–t2内竹林生物质碳储量的变化量(tCO2-e)1,,tPROJBAMBOOC时间为t1时竹林生物质碳储量(tCO2-e)2,,tPROJBAMBOOC时间为t2时竹林生物质碳储量(tCO2-e)DR根据项目的不确定性确定的调减因子（%）调减因子表不确定性BAMBOOUNC,（%）DR（%）9()012,,,,>−tPROJBAMBOOtPROJBAMBOOCC()012,,,,<−tPROJBAMBOOtPROJBAMBOOCC小于或等于10%0%0%大于10%小于20%6%-6%大于20%小于30%11%-11%大于或等于30%增加监测样地数量4.不需监测的数据和参数（采用的缺省值或一次性测定值）数据/参数)(,ajDBHtf单位:cm应用的公式编号:（4）描述:在竹林发育成林阶段，竹种（组）j的平均胸径与竹林年龄的相关方程9根据AR-ACM0003调整。26数据源:数据源优先选择次序为：(a)现有的、当地的或相似生态条件下的基于竹种或竹种组的数据；(b)从类似竹种组中选择；(c)根据缺省方程计算：大径散生竹（毛竹）：320056.00452.0572.02000.5aaatttDBH⋅−⋅+⋅+=大径丛生竹（麻竹）：atDBH⋅+=1039603.1960772.1测定步骤（如果有）不适用说明:由于有的文献是用胸径，有的用眉径，也有少量用地径，因此，这里的胸径可以用眉径或地径代替，取决于方程所用的变量的实际含义。数据/参数),(,DBHtfajH单位:m应用的公式编号:（5）描述:在竹林发育成林阶段，竹种（组）j的平均高度与竹林年龄和（或）平均胸径的相关方程数据源:数据源优先选择次序为：(a)现有的、当地的或相似生态条件下的基于竹种或竹种组的数据；(b)从类似竹种组中选择；(c)根据缺省方程计算：大径散生竹（毛竹）：20465.06426.15702.0DBHDBHH⋅−⋅+=大径丛生竹（麻竹）：atH12233144.006452891.01⋅+=测定步骤（如果有）不适用说明:由于有的文献是用胸径，有的用眉径，也有少量用地径，因此，这里的胸径可以用眉径或地径代替，取决于方程所用的变量的实际含义。数据/参数),(jjABHDBHf单位:Kgd.m.株–1应用的公式编号:(6)描述:竹种（组）j的地上生物量方程，即单株地上生物量与胸径和竹高的相关方程27数据源:数据源优先选择次序为：(a)现有的、当地的或相似生态条件下的基于竹种或竹种组的数据；(b)从附件中选择适合的竹类生物量方程。测定步骤（如果有）不适用说明:在选择竹子生物量方程时，须充分考虑竹子类型（丛生、散生、混生、大径竹和小径竹）。这里的胸径可以用眉径或地径代替，取决于方程所用的变量的实际含义。数据/参数BjCF,单位:tC(td.m.)-1应用的公式编号:（7）描述:竹子含碳率数据源:数据源优先选择次序为：(a)当地分别竹种或竹种组的数据(b)省级分别竹种或竹种组的数据(如省级温室气体清单)(c)国家级分别竹种或竹种组的数据(如国家温室气体清单)(d)缺省值0.50测定步骤（如果有）不适用说明:数据/参数)(,ajNtf单位:株.hm-2或株.丛-1应用的公式编号:（8）描述:竹种（组）j的立竹度与竹林年龄的相关方程数据源:数据源优先选择次序为：(a)现有的、当地的或相似生态条件下的基于竹种或竹种组的数据；(b)从类似竹种组中选择；(c)毛竹缺省数据：林龄1234567立竹度(株.hm-2)3756148901,4552,1752,3352,550测定步骤（如果有）不适用说明:28数据/参数atjR,单位:无量纲应用的公式编号:（11）描述:j竹种（组）在竹林年龄为at时的地下生物量与地上生物量之比）数据源:数据源优先选择次序为：(a)现有的、当地的或相似生态条件下的基于竹种或竹种组的数据；(b)国家级基于竹种或竹种组和竹林年龄的数据(如森林资源清查或国家温室气体清单编制中的数据)；(c)如果没有上述相关的数据源可用，则可假定地下生物量与地上生物量之比为常数（不随竹林年龄at而发生变化）并可从下表中选择缺省值：竹子类型代表竹种平均值样本数标准差大径散生竹刚竹属（毛竹）0.605500.071刚竹属(毛环竹)0.688160.023大径丛生竹箣竹属（绿竹）1.127270.112其他大径丛生竹1.376210.330小径散生竹刚竹属（雷竹）0.200750.034小径丛生竹所有小径丛生竹0.632140.153复轴混生型竹寒竹属（方竹）0.389400.051寒竹属（筇竹）0.8161650.097其他混生竹0.928140.162数据源：来自竹林生物量文献的数据库。测定步骤（如果有）不适用说明:推荐优先使用基于竹林年龄的地下生物量与地上生物量之比，直到竹林地下生物量达到稳定为止。数据/参数SCF单位:tC(td.m.)-1应用的公式编号:（13）描述:灌木含碳率29数据源:数据源优先选择次序为：(a)当地分别灌木种（组）的数据；(b)省级分别灌木种（组）的数据(如省级温室气体清单)；(c)国家级分别灌木种（组）的数据(如国家级温室气体清单)；(d)可采用IPCC缺省值：0.47tC(td.m.)-1测定步骤（如果有）不适用说明:数据/参数SR单位:无量纲应用的公式编号:（13）描述:灌木和杂竹丛的地下生物量与地上生物量之比数据源:数据源优先选择次序为：(a)现有的、当地的或相似生态条件下的基于树种或树种组的数据；(b)省级基于灌木种的数据(如森林资源清查或国家温室气体清单编制中的数据)；(c)国家级基于灌木种的数据(如森林资源清查或国家温室气体清单编制中的数据)。如果没有上述数据源的数据，可采用IPCC缺省值0.4。但对于杂竹丛，如果没有数据，可采用上述竹林的缺省值。测定步骤（如果有）不适用说明:数据/参数SFBDR单位:无量纲应用的公式编号:（14）描述:灌木盖度为1.0时的每公顷灌木生物量与拟议项目所在地区完全郁闭森林每公顷地上生物量之比数据源:数据源优先选择次序为：(a)现有的、当地的或相似生态条件下的数据(b)国家级的数据(如森林资源清查或国家温室气体清单编制中的数据)(c)使用缺省值0.10测定步骤（如果有）不适用说明:数据/参数FORESTB单位:td.m.hm-2应用的公式编号:（14）描述:拟议项目所在地区完全郁闭人工林平均每公顷地上生物量30数据源:数据源优先选择次序为：(a)现有的、当地的或相似生态条件下的数据(b)国家级数据(如森林资源清查或国家温室气体清单编制中的数据)测定步骤（如果有）不适用说明:数据/参数DFLI单位:%应用的公式编号:（16）描述:枯落物碳储量与竹林生物质碳储量之比数据源:数据源优先选择次序为：(a)现有的、当地的或相似生态条件下的基于竹种或竹种组的数据；(b)从类似竹种组中选择；(c)从下表中选择缺省值：森林类型平均值样本量标准差散生竹5.28%130.932%丛生竹6.25%110.840%数据来源：基于发表的生物量和枯落物文献建立的数据库。测定步骤（如果有）不适用说明:数据/参数SOCREF数据单位tC.hm-2应用的公式编号(17)、（18）描述与项目第i项目碳层具有相似气候、土壤条件的当地自然植被（如当地未退化的、未利用土地上的自然植被）下土壤有机碳储量的参考值数据源数据源优先顺序：(a)公开出版的与项目区条件相似的数据；(b)相关国家资源调查数据(如土壤普查、森林资源清查或温室气体国家清单)；(c)从下表中选择缺省值：矿质土壤的土壤有机碳库碳储量缺省值(tCha-1，深度0-30cm)气候地区高活性黏土(a)低活性黏土(b)沙质土(c)火山灰土(d)热带，干燥38353150热带，湿润65473970热带，湿润446066130热带，山地88633480(a)高活性黏质土壤（HAC）为轻度到中度风化土壤，由2：1硅酸盐黏土矿物组31成（见世界土壤资源参比基础（WRB）的分类，包括薄层土,变性土，栗钙土，黑钙土，黑土，淋溶土，高活性强酸土，漂白红砂土，碱土，钙积土，石膏土，暗色土，雏形土，粗骨土；在美国农业部（USDA）的分类包括软土、变性土、高基淋溶土、旱成土、始成土）；(b)低活性黏质土壤为高度风化的土壤，由黏土矿质与非结晶铁铝氧化物按1：1的比例组成（见世界土壤资源参比基础（WRB）的分类，包括低活性强酸土、低活性淋溶土、黏绨土、铁铝土；美国农业部的分类包括老成土、氧化土和酸性淋溶土）(c)包括标准结构中砂土比例>70%且粘土比例<8%的所有土壤（与土壤分类无关，在世界土壤资源参比基础（WRB）中包括砂性土，在美国农业部（USDA）的分类中包括砂新成土）；(d)来自火山灰及其同分异构矿质的土壤（在世界土壤资源参比基础（WRB）中为火山灰土，在美国农业部（USDA）的分类为火山灰土）。测定步骤不适用说明数据/参数fLU,i数据单位无量纲应用的公式编号（17）描述第i项目碳层与基线土地利用方式相关的碳储量变化因子数据源数据源优先顺序：(a)公开出版的与项目区条件相似的数据；(b)国家有关资源调查数据(如土壤普查、森林资源清查或温室气体国家清单)；(c)从下表中选择：不同农田土地利用相关的碳储量变化因子（20年内的净效应）土地利用温度条件湿度条件因子值说明长期耕种热带干燥0.58连续耕种20年以上的农田湿润0.48热带山地不适用0.64短期耕种（﹤20年）或闲置（﹤5年）热带干燥0.93连续耕种时间不足20年的农田、和（或）在最近20年的任意时间段内闲置时间少于5年的农田湿润0.82热带山地不适用0.88草地管理碳储量相对变化因子（20年内的净效应）土地利用气候区因子值说明全部全部1.00所有永久草地的土地利用因子值为1测定步骤不适用说明32数据/参数fMG,i数据单位无量纲应用的公式编号（17）描述第i项目碳层与基线管理模式相关的碳储量变化因子数据源数据源优先顺序：（a）现有的、公开发表的、当地的或相似生态条件下的数据；（b）从下表中选择：不同农田管理措施相关的碳储量变化因子（20年内的净效应）土地管理温度条件湿度条件因子值说明及标准全耕全部干燥和湿润/湿地1.00充分翻耕或在一年内频繁耕作导致强烈土壤扰动。在种植期地表残体盖度低于30%少耕热带干燥1.09对土壤的扰动较低（通常耕作深度浅，不充分翻耕）。在种植期地表残体盖度通常大于30%湿润/湿地1.15热带山地1.09草地管理碳储量相对变化因子（20年内的净效应）土地管理气候区因子值说明中等退化草地热带0.97过牧或中度退化，相对于未退化草地，生产力较低，且未实施改良措施热带山地0.96严重退化全部0.70测定步骤不适用说明数据/参数iINf,数据单位无量纲应用的公式编号（17）描述第i项目碳层与基线有机碳输入类型（如农作物秸杆还田、施用肥料）相关的碳储量变化因子数据源数据源优先顺序：（a）现有的、公开发表的、当地的或相似生态条件下的数据；（b）省级或国家水平的适用于项目实施区的数据；（c）从下表中选择缺值：不同农田输入相关的碳储量变化因子（20年内的净效应）33输入水平温度条件湿度条件因子值说明及标准低热带干燥0.95收集去除或燃烧地表残体(如秸秆焚烧）；或频繁裸地休耕；或农作物残体较少（如蔬菜，烟草，棉花等）；或不施矿物肥料、不种植固氮作物等湿润/湿地0.92热带山地0.94中全部干燥和湿润/湿润1.00所有作物残留都返回到田地里。若残留物被移除则添加有机质(如粪肥)。另外，施加矿质肥料或轮作固氮作物。高，不施肥热带干燥1.04通过额外措施（如种植残体较多的农作物、施用绿肥、种植覆盖作物、休耕、灌溉、一年生作物轮作中频繁种植多年生草本植物，但不施有机肥），使作物残体的碳输入量显著增加。湿润1.11热带山地1.08草地管理碳储量相对变化因子（20年内的净效应）输入水平气候区因子值说明中等退化草地热带0.97过牧或中度退化，相对于未退化草地，生产力较低，且未实施改良措施热带山地0.96严重退化全部0.70测定步骤不适用说明数据/参数tyBU单位:%应用的公式编号:（21）描述:竹子采伐用于ty类竹产品利用率数据源:数据源优先选择次序为：(a)当地基于竹产品种类和竹种的数据；(b)国家级基于竹产品种类和竹种的数据；(c)使用如下缺省值：•竹材层压板：50%•竹胶合板：40%•竹地板：20%•纸和纸产品：90%•其他（如生产和生活工具）：50%数据来源：王小青等（2002）34测定步骤（如果有）不适用说明:数据/参数tyLT单位:年应用的公式编号:(22)描述:ty类竹产品的使用寿命数据源:数据源优先选择次序为：(a)公开出版的适于当地条件和产品类型的文献数据；(b)国家级基于木（竹）产品的数据;；(c)从下表选择缺省数据：木（竹）产品类型tyLT建筑50家具30矿柱15车船12包装用材8纸和纸板3锯材30人造板20薪材1数据源:参考下列文献确定：a)IPCCLULUCF优良做法指南；b)COP17关于《京都议定书》第二承诺期LULUCF的决议；c)白彦锋(2010)。测定步骤（如果有）N/A说明:数据/参数COMF数据单位无量纲应用的公式编号（24）描述竹林燃烧系数数据源数据来源的选择应遵循如下顺序：（a）项目实施区当地的调查数据；（b）相邻地区相似条件下的调查数据；（c）国家水平的适用于项目实施区的数据；（d）使用缺省值：0.67测定步骤说明35数据/参数EFCH4数据单位gCH4·(kg燃烧的干物质)-1应用的公式编号（24）描述CH4排放因子数据源数据来源的选择应遵循如下顺序：（a）项目实施区当地的调查数据；（b）相邻地区相似条件下的调查数据；（c）省级或国家水平的适用于项目实施区的数据；（d）采用缺省值：6.8测定步骤说明数据/参数EFN2O数据单位gN2O.(kg燃烧的干物质)-1应用的公式编号(24)描述N2O排放因子数据源应对数据来源进行选择，具体选择顺序如下：（a）项目实施区当地的调查数据；（b）相邻地区相似条件下的调查数据；（c）省级或国家水平的适用于项目实施区的数据；（d）采用缺省值：0.26。测定步骤说明5.监测的数据和参数数据/参数tjiBambooA,,,单位:hm2应用的公式编号:（9）、（10）、(11)、（20）、描述:第t年第i碳层j竹种（组）的面积数据源:野外测定测定步骤采用国家森林资源调查、规划设计调查或作业设计调查使用的标准操作程序(SOP)，如果没有，可采用公开出版的相关技术手册或IPCCGPGLULUCF2003中描述的SOP。其边界数据最好易于输入GIS。监测频率:首次核查后每3-10年一次QA/QC程序:采用国家森林资源调查、规划设计调查或作业设计调查使用的质量保证和质量控制(QA/QC)程序，如果没有，可采用公开出版的相关技术手册或IPCCGPGLULUCF2003中描述的QA/QC程序。说明36数据/参数pA单位:hm2应用的公式编号:（27）、（28）、（31）描述:样地的面积数据源:野外测定测定步骤采用国家森林资源清查使用的标准操作程序(SOP)，如果没有，可采用公开出版的相关技术手册或IPCCGPGLULUCF2003中描述的SOP。监测频率:首次核查开始每3-10年一次QA/QC程序:采用国家森林资源清查使用的质量保证和质量控制(QA/QC)程序，如果没有，可采用公开出版的相关技术手册或IPCCGPGLULUCF2003中描述的QA/QC程序。说明:样地位置应用GPS或Compass记录且在图上标出。数据/参数DBH单位:cm应用的公式编号:描述:通常为胸径，也可以是眉径或地径数据源:野外样地测定。测定步骤采用国家森林资源清查使用的标准操作程序(SOP)，如果没有，可采用公开出版的相关技术手册或IPCCGPGLULUCF2003中描述的SOP。监测频率:首次核查开始每3-10年一次QA/QC程序:采用国家森林资源清查使用的质量保证和质量控制(QA/QC)程序，如果没有，可采用公开出版的相关技术手册或IPCCGPGLULUCF2003中描述的QA/QC程序。说明:数据/参数H单位:m应用的公式编号:描述:竹秆高度数据源:野外样地测定。测定步骤采用国家森林资源清查使用的标准操作程序(SOP)，如果没有，可采用公开出版的相关技术手册或IPCCGPGLULUCF2003中描述的SOP。监测频率:首次核查开始每3-10年一次QA/QC程序:采用国家森林资源清查使用的质量保证和质量控制(QA/QC)程序，如果没有，可采用公开出版的相关技术手册或IPCCGPGLULUCF2003中描述的QA/QC程序。说明:高度可以是全竹秆高，也可以是其他高度，取决于方程中使用的高度含义37数据/参数tiPROJSHRUBA,,_单位:hm2应用的公式编号:(13)描述:t年时i碳层灌木（包括杂竹丛）的面积数据源:野外测定测定步骤采用国家森林资源清查使用的标准操作程序(SOP)，如果没有，可采用公开出版的相关技术手册或IPCCGPGLULUCF2003中描述的SOP。监测频率:首次核查开始每3-10年一次QA/QC程序:采用国家森林资源清查使用的质量保证和质量控制(QA/QC)程序，如果没有，可采用公开出版的相关技术手册或IPCCGPGLULUCF2003中描述的QA/QC程序。说明:数据/参数tiPROJSHRUBCC,,_单位:无量纲应用的公式编号:(14)描述:t年时，第i碳层的灌木（包括杂竹丛）盖度数据源:野外测定测定步骤（如果有）考虑到灌木生物量相对于林木生物量较小，可采用简化的方法测定灌木盖度，如样线方法、目测方法频率首次核查开始每3-10年一次QA/QC程序采用国家森林资源清查使用的标准操作程序(SOP)，如果没有，可采用公开出版的相关技术手册或IPCCGPGLULUCF2003中描述的SOP。说明当基线情景下存在周期性采伐和炼山时，采用0.5做为灌木的平均盖度。如果选用不同的数值，须提供透明和可核实的信息来证明数据/参数tStemBAMBOOC,,单位:tCO2-e应用的公式编号:（21）描述:t年时，项目采伐的竹秆生物质碳储量数据源:调查测定测定步骤（如果有）如果采伐的竹子是以竹秆鲜重计，则应将鲜重通过含水率换算成干重，然后转化为CO2的量。如果采伐利用整株竹子（包括枝和叶），则为地上生物质碳储量。频率每次采伐时QA/QC程序采用国家森林资源调查规划设计使用的标准操作程序(SOP)说明38数据/参数ty单位:无量纲应用的公式编号:（21）、（22）描述:从竹子造林项目采伐的竹子利用产品类型数据源:调查测定测定步骤（如果有）•对于社区采伐，采用PRA的方法调查其采伐的竹子的用途、销售去向，调查样本不少于所涉社区户数的10%。同时跟踪调查所销售竹子的用途和产品种类及其比例•对于企业为主的采伐，记录销售去向，并跟踪调查所销售竹子的用途和产品种类及其比例频率社区采伐每年一次；企业采伐每采伐一次监测一次QA/QC程序说明数据/参数ABURN,i,t数据单位ha应用的公式编号(24)、（25）描述第t年第i层发生火灾的面积；数据源野外测量或者遥感监测测定步骤对发生火灾的区域边界进行定位，可采用地面GPS定位或是通过遥感数据反演监测频率每次森林火灾发生时均须测量QA/AC采用国家森林资源调查使用的质量保证和质量控制(QA/QC)程序说明396.参考文献[1]何亚平，费世民，蒋俊明，陈秀明，余英，唐森强，朱维双.2007.长宁毛竹和苦竹有机碳空间分布格局.四川林业科技,28(5):10-14[2]陈辉，洪伟，兰斌，郑郁善，何东进.1998.闽北毛竹生物量与生产力的研究.林业科学,34(Sp.1):60-64[3]巫启新.贵州毛竹林类型与林分结构的研究.1994.竹子研究汇刊,1983.2(1):112-124[4]聂道平.1994.毛竹林结构的动态特性.林业科学,30(3):201-208.[5]郑郁善，洪伟.1998.毛竹经营学.厦门:厦门大学出版社[6]周国模.2006.毛竹林生态系统中碳储量、固定及其分配与分布的研究.浙江大学博士学位论文.[7]陈双林,吴柏林,吴明,等.2004.退化低丘红壤新造毛竹林地上部分生物量的研究.江西农业大学学报,26(4):527-531.[8]黎曦,鲍雪芳,王福升.2007.赣南毛竹生物量研究.安徽林业科技.[9]徐道旺，陈少红，杨金满.2004.毛环竹笋用林生物量结构调查分析.福建林业科技,31(1):67-70.[10]郑郁善，陈希英，方承，梁鸿集.1997.台湾桂笔生物产量模型研究.福建林学院学报，1（1）：52-55[11]郑郁善,梁鸿燊.1998.台湾桂竹各器官生物量模型研究.竹子研究汇刊,17(1):37-41.[12]梁鸿燊,陈学魁.1998.麻竹单株生物量模型研究.福建林学院学报,18(3):260-262.[13]郑郁善，梁鸿集，游兴早.1997.绿竹生物量模型研究.竹子研究汇刊，16（4）：43-46[14]苏文会,顾小平,关凤英,等.2006.大木竹种群生物量结构及其回归模型.南京林业大学学报（自然科学版）,30(5):51-54.[15]付建生,董文渊,韩梅,等.2007.撑绿竹不同径阶的生物量结构分析.林业科技开发,21(5):47-49.[16]徐小军,周国模,杜华强,等.2011.基于LandsatTM数据估算雷竹林地上生物量.林业科学,47(9):1-6.[17]黄宗安,郑明生,张居文,等.2000.石竹各器官生物量回归模型研究.福建林业科技,27(3):35-37[18]潘红丽,李迈和,田雨,等.2010.卧龙自然保护区油竹子形态学特征及地上部生物量对海拔梯度的响应.四川林业科技,31(3):30-36[19]秦自生,马焱,马恒银.1990.拐棍竹生物生产量的预测模型.四川师范学院学报,11(2):98-102[20]郭子武,李迎春,杨清平,等.2009.花吊丝竹立竹构件与生物量关系的研究.热带亚热带植物学报,17(6):543-548.[21]郑郁善，陈明阳，林金国，赵荣军.1999.肿节少穗竹各器官生物量模型研究.福建林学院学院,18(2):159-162[22]魏泽长，武大宇，王希荣，张德惠，刘年贵.1986.水竹人工林生物量结构的研究.植物生态学与地植物学学报,10(3):190-198.[23]杨前宇,谢锦忠,张玮,等.2011.椽竹各器官生物量模型.浙江农林大学学报,28(3):519-526.[24]郑容妹,郑郁善,闽锋,等.2003.苦竹生物量模型的研究.福建林学院学报,23(1):61-64.[25]林新春,方伟,俞建新,等.2004.苦竹各器官生物量模型.浙江林学院学报,21(2):168-171.[26]郑金双,曹永慧,肖书平,等.2001.茶秆竹生物量模型研究.竹子研究汇刊,4020(4).[27]刘庆，钟章成.1996.斑苦竹无性系种群生物量结构与动态研究.竹类研究,(1):51-56[28]顾大形,陈双林,郭子武,等.2011.四季竹立竹地上现存生物量分配及其与构件因子关系,林业科学研究,24(4):495-499.41附件:竹子生物量方程竹林类型竹种方程（Kg干重/株）建模地点文献大径散生竹刚竹属(毛竹)9.7697+DBH2.3434一DBH0.35132=总W四川长宁[1]1.6579DBH3864951.0=地上W2.3131-0.5805HDBH1341642.0=总W闽北[2]22.34376759DBH10872076.0=干W70.85133807DBH79406626.0=枝叶W12.02978185DBH28040806.0=总W黔北[3]933.03054.02.0811340.1DBH0925.0DBHNW+=地上江西大岗山[4]2.3049DBH1574.0=地上W[5]浙江[6]W地上=-11.497+3.0465DBH+0.1117DBH2江西、浙南[7]W地上=0.04504749281DBH2.2890229H0.28643528赣南[8]刚竹属(毛环竹)4396.20.6278DBH014467.0HW=地上2214.20.59736DBH22128.0HW=总[9]刚竹属(台湾桂竹)8990.11639.0DBHW⋅=地上9756.11718.0DBHW⋅=总福建东部[10]W地上=0.00152DBH2.4094H-0.3028W总=0.000721DBH2.8382H-0.3078福建东部[11]大径丛生竹牡竹属(麻竹)W地上=0.540093DBH1.9305W地上=0.172139DBH1.5684H0.3916福建、海南[12]绿竹属(绿竹)244206.2197169.0DBHW⋅=地上488312.0687679.1194103.0HDBHW⋅=地上福建[13]簕竹属(大木竹)W地上=0.4524DBH2.0347W总=0.5122DBH2.0391浙南[14]撑绿竹W地上=3.11219+0.03232DBH2HHDBHW200033789.055698.3+=总云南水富[15]小径散生竹刚竹属(雷竹)W地上=0.1939DBH1.5654浙江西北部[16]刚竹属(石竹)W总=0.0302DBH2.4123H0.6262单福建尤溪[17]小径丛生箭竹属(油竹子W地上=0.020DBH2H+0.029.12四川[18]772.3)028.0148.0(747.787555.52.771++=TTDW地上42箭竹属(拐棍竹)W地上=0.0719183DBH0.822738H四川[19]牡竹属(花吊丝竹)1龄竹：W地上=-5.45421+1.46011DBH+0.29207H2龄竹：W地上=-3.34805+1.94950DBH+0.13412H3龄竹：W地上=-2.95277+1.84698DBH4龄竹：W地上=-1.45958+1.15918DBH福建华安[20]少穗竹属(肿节少穗竹)1.7668DBH1888.0=秆W1.2135DBH0633.0=枝W1.1853DBH0722.0=叶W1.3836DBH3626.0=总W福建[21]刚竹属(水竹)1.5373DBH6439.0=地上W1.4117DBH3841.0=总W安徽舒城[22]簕竹属(椽竹)W地上=-7445.916+3925.48DBH+45.439DBH2-96.666DBH3W总=-7360.122+3933.155DBH+41.158DBH2-93.171DBH3福建建瓯[23]复轴混生竹大明竹属(苦竹)0.0914+DBH0.0027+DBH0.26682=总W四川长宁[1]W地上=0.1173DBH0.8254H1.1605W总=0.127DBH0.8621H1.1756福建尤溪[24]W地上=432.4468-479.3075DBH+422.8285DBH2W总=0.1910380DBH1.1986H0.2962浙江杭州[25]矢竹属(茶秆竹)W总=1.0493DBH1.4861福建明溪[26]大明竹属(斑苦竹))(DBH0054.02180.02HW+=地上)(DBH0091.08378.02HW+=总重庆缙云山[27]少穗竹属(四季竹)1龄竹：W地上=-0.073+0.51DBH-0.30DBH2+0.112DBH32龄竹：W地上=-0.515+2.37DBH-2.4DBH2+0.9108DBH3浙江西北[28]注：W总：全株总生物量；W地上：地上生物量；DBH:竹子胸径或眉径（cm）；H：竹高(m)；T：竹龄(a)；N：立竹株数(株)