陆地生态系统生产总值（GEP）核算技术指南TheTechnicalGuidelineonGrossEcosystemProduct（GEP）（1.0version）生态环境部环境规划院ChineseAcademyofEnvironmentalPlanning中国科学院生态环境研究中心ResearchCenterforEco-EnvironmentalSciences,CAS二〇二〇年九月前言为全面贯彻落实习近平生态文明思想和“绿水青山就是金山银山”理念，以及《中共中央关于全面深化改革若干重大问题的决定》、《关于加快推进生态文明建设的意见》和《生态文明体制改革总体方案》中建立生态效益评估机制、促进人与自然和谐的部署，保障国家和区域生态安全，指导和规范陆地生态系统生产总值核算工作，提高陆地生态系统生产总值实物量与价值量核算的科学性、规范性和可操作性，制定本指南，供地方参考使用。本指南可为生态效益纳入经济社会发展评价体系、完善发展成果考核评价体系提供重要支撑，为建立生态产品实现机制、区域生态补偿、自然资源资产审计、自然资产负债表编制等制度的实施提供科学依据。本指南规定了陆地生态系统生产总值实物量与价值量核算的技术流程、指标体系与核算方法等内容。本指南附录均为资料性附录。本指南由生态环境部环境规划院与中国科学院生态环境研究中心共同编制。本指南主要起草人：欧阳志云、王金南、肖燚、於方、郑华、马国霞，宋昌素、周夏飞、邹梓颖，杨威杉，逯非，周颖，韩宝龙，吴春生，彭菲，宁佳。1目录1总则............................................................11.1适用范围......................................................11.2规范性引用文件................................................11.3术语和定义....................................................12核算程序与内容..................................................42.1核算原则......................................................42.2核算程序......................................................42.3核算内容与指标................................................72.4生态系统服务实物量核算........................................92.5生态系统服务价值量核算.......................................103物质产品实物量和价值量核算.....................................113.1物质产品实物量核算...........................................113.2物质产品价值量核算...........................................134水源涵养实物量和价值量核算......................................144.1水源涵养实物量核算...........................................144.2水源涵养价值量核算...........................................155土壤保持实物量和价值量核算......................................165.1土壤保持实物量核算...........................................165.2土壤保持价值量核算...........................................176防风固沙实物量和价值量核算......................................186.1防风固沙实物量核算...........................................186.2防风固沙价值量核算...........................................197海岸带防护实物量和价值量核算....................................207.1海岸带防护实物量核算.........................................207.2海岸带防护价值量核算.........................................218洪水调蓄实物量和价值量核算......................................218.1洪水调蓄实物量核算...........................................2128.2洪水调蓄价值量核算...........................................259碳固定实物量和价值量核算........................................269.1碳固定实物量核算.............................................269.2碳固定价值量核算.............................................3210氧气提供实物量和价值量核算.....................................3210.1氧气提供实物量核算..........................................3210.2氧气提供价值量核算..........................................3311空气净化实物量和价值量核算.....................................3411.1空气净化实物量核算..........................................3411.2空气净化价值量核算..........................................3512水质净化实物量和价值量核算.....................................3612.1水质净化实物量核算..........................................3612.2水质净化价值量核算..........................................3813气候调节实物量和价值量核算.....................................3913.1气候调节实物量核算..........................................3913.2气候调节价值量核算..........................................4114物种保育实物量和价值量核算.....................................4114.1物种保育实物量核算..........................................4114.2物种保育价值量核算..........................................4215休闲旅游实物量和价值量核算.....................................4315.1休闲旅游实物量核算..........................................4315.2休闲旅游价值量核算..........................................4416景观价值实物量和价值量核算.....................................4516.1景观价值实物量核算..........................................4516.2景观价值价值量核算..........................................4617附则..........................................................4617.1核算数据....................................................4617.2成果表达形式................................................4717.3结果校核....................................................47317.4信息管理....................................................47附录A（资料性附录）生态系统生产总值核算数据来源表................48附录B（资料性附录）陆地生态系统生产总值核算报告编写大纲..........51附录C（资料性附录）生态系统生产总值实物量核算参数参考表..........52附录C.1土壤保持实物量核算参数..................................52附录C.2防风固沙实物量核算参数..................................55附录C.3固碳实物量核算参数......................................61附录C.4空气净化实物量核算参数..................................65附录C.5水质净化实物量核算参数...................................6611总则1.1适用范围本指南适用于省、市、县等行政区域开展生态系统生产总值核算工作。自然地理区域可参考本指南开展生态系统生产总值核算。1.2规范性引用文件本指南内容引用了下列文件中的条款。凡是不注日期的引用文件，其有效版本适用于本指南。《环境经济核算体系——实验生态账户》(2012)GB3095环境空气质量标准GB3838地表水环境质量标准GB/T12343国家基本比例尺地图编绘规范GB/T13923基础地理信息要素分类与代码CH/T9005基础地理信息数据库基本规定GB-T38582森林生态系统服务功能评估规范LY∕T2899湿地生态系统服务评估规范1.3术语和定义以下术语和定义适用于本指南。（1）陆地生态系统（Terrestrialecosystem）陆地生态系统是指地球表面陆地生物及其环境通过能流、物流、信息流形成的功能整体。陆地生态系统包括森林生态系统、草地生态2系统、湿地生态系统、荒漠生态系统、农田生态系统、城市生态系统等类型。（2）湿地生态系统（Wetlandecosystem）湿地生态系统由陆地和水域相互作用而形成的兼顾水域和陆地生态系统特征的自然综合系统，包括陆地所有淡水生态系统、陆地和海洋过渡地带的滨海湿地生态系统和海洋边缘部分咸水、半咸水水域。（3）草地生态系统（Grasslandecosystem）草地生态系统是以饲用植物和食草动物为主体的生物群落与其生存环境共同构成的开放生态系统，包括人工草地生态系统和天然草地生态系统两大类。（4）森林生态系统（Forestecosystem）森林生态系统是以乔木为主体的生物群落与其非生物环境相互作用，并产生能量转换和物质循环的综合系统，包括天然林生态系统和人工林生态系统。（5）荒漠生态系统（Desertecosystem）荒漠生态系统指由超强耐旱生物及其干旱环境所组成的生态系统。（6）农田生态系统（Farmlandecosystem）农田生态系统是指人类在以作物为中心的农田中，利用生物和非生物环境之间以及生物种群之间的相互关系，通过合理的生态结构和高效生态机能，进行能量转化和物质循环，并按人类社会需要进行物质生产的综合体。（7）城市生态系统（Urbanecosystem）城市生态系统是城市居民与其环境相互作用而形成的统一整体，3由自然环境、社会经济和文化科学技术共同组成。（8）生态系统服务（Ecosystemservices）人类从生态系统中得到的惠益，包括物质产品供给、调节服务、文化服务以及支持服务。（9）生态系统生产总值（GrossEcosystemProduct,GEP）生态系统生产总值也可简称为生态产品总值，是指生态系统为人类福祉和经济社会可持续发展提供的各种最终产品与服务价值的总和，主要包括生态系统提供的物质产品、调节服务和文化服务。（10）物质产品（Materialservices）人类从生态系统获取的可在市场交换的各种物质产品，如食物、纤维、木材、药物、装饰材料与其他物质材料。（11）调节服务（Regulatingservices）生态系统提供改善人类生存与生活环境的惠益，如调节气候、涵养水源、保持土壤、调蓄洪水、降解污染物、固定二氧化碳、氧气提供等。（12）文化服务（Culturalservices）人类通过精神感受、知识获取、休闲娱乐和美学体验从生态系统获得的非物质惠益。（13）实物量（Biophysicalvalue）生态系统产品与服务的实物量，如粮食产量、木材生产量、水产品捕捞量、洪水调蓄量、土壤保持量、碳固定量与景点旅游人数等。（14）价值量（MonetaryValue）生态系统产品与服务的货币价值。42核算程序与内容2.1核算原则陆地生态系统生产总值核算应遵循以下原则：科学性原则。根据生态系统与人类福祉的关系构建核算框架与指标体系。以生态系统结构、格局和过程与生态系统服务关系为基础，构建核算方法。实用性原则。根据核算区域生态系统的特点、核算目的和数据可获得性，确定适合本区域的核算内容、指标和方法。系统性原则。生态系统生产总值核算既要考虑为当地人提供的惠益，也需要考虑为其它地区提供的惠益。开放性原则。应根据生态系统核算研究的最新成果，发展和完善核算的指标与方法。2.2核算程序陆地生态系统生产总值核算的主要工作程序包括：根据核算目的，确定生态系统生产总值核算区域范围，明确生态系统类型以及生态系统服务的清单，确定核算模型方法与适用技术参数，开展生态系统生产总值实物量与价值量核算（见图2-1）。5图2-1生态系统生产总值核算工作程序6（1）确定核算的区域范围：根据核算目的，确定生态系统生产总值核算的空间范围。核算区域可以是行政区域，如村、乡、县、市或省，也可以是功能相对完整的生态地理单元，如一片森林、一个湖泊、一片沼泽或不同尺度的流域，以及由不同生态系统类型组合而成的地域单元。（2）明确生态系统类型与分布：调查分析核算区域内的森林、草地、湿地、荒漠、农田、城镇等生态系统类型、面积与分布，绘制生态系统空间分布图。（3）编制生态系统产品与服务清单：根据生态系统类型及核算的用途，如生态补偿、离任审计、生态产品交易，调查核算范围的生态系统服务的种类，编制生态系统产品和服务清单。当核算目标为评估生态保护成效时，可以只核算生态调节服务和生态文化服务价值。（4）收集资料与补充调查：收集开展生态系统生产总值核算所需要的相关文献资料、监测与统计等信息数据以及基础地理与地形图件，开展必要的实地观测调查，进行数据预处理以及参数本地化。（5）开展生态系统产品与服务实物量核算：选择科学合理、符合核算区域特点的实物量核算方法与技术参数，根据确定的核算基准时间，核算各类生态系统产品与服务的实物量。（6）开展生态系统产品与服务价值量核算：根据生态系统产品与服务实物量，运用市场价值法、替代成本法等方法，核算生态系统产品与服务的货币价值；无法获得核算年份价格数据时，利用已有年份数据，按照价格指数进行折算。（7）核算生态系统生产总值：将核算区域范围的生态产品与服务价值加总，得到生态系统生产总值。72.3核算内容与指标2.3.1核算内容生态系统生产总值核算生态系统物质产品价值、调节服务价值和文化服务价值，不包括生态支持服务价值。根据不同的核算目的，核算不同类型的生态系统生产总值：（1）核算生态系统对人类福祉和经济社会发展支撑作用时，核算生态系统的物质产品价值、调节服务价值和文化服务价值之和。GEP=EPV+ERV+ECV（2）核算生态保护成效与生态效益时，核算生态系统的调节服务价值和文化服务价值。GEP=ERV+ECV式中，GEP为生态系统生产总值；EPV为生态系统物质产品价值；ERV为生态系统调节服务价值；ECV为生态系统文化服务价值。2.3.2指标体系生态系统生产总值核算指标体系由物质产品、调节服务和文化服务3大类服务构成，其中：物质产品主要包括农业产品、林业产品、畜牧业产品、渔业产品、生态能源和其他；调节服务主要包括水源涵养、土壤保持、防风固沙、海岸带防护、洪水调蓄、碳固定、氧气提供、空气净化、水质净化、气候调节和物种保育；文化服务主要包括休闲旅游、景观价值，见表2-1、表2-2。8表2-1生态系统生产总值（GEP）核算指标序号一级指标二级指标指标说明1物质产品农业产品从农业生态系统中获得的初级产品，如稻谷、玉米、谷子、豆类、薯类、油料、棉花、麻类、糖类、烟叶、茶叶、药材、蔬菜、水果等。2林业产品林木产品、林产品以及与森林资源相关的初级产品，如木材、竹材、松脂、生漆、油桐籽等。3畜牧业产品利用放牧、圈养或者两者结合的方式，饲养禽畜获得的产品，如牛、羊、猪、家禽、奶类、禽蛋等。4渔业产品利用水域中生物的物质转化功能，通过捕捞、养殖等方式取的水产品，如鱼类、其他水生动物等。5生态能源生态系统中的生物物质及其所含的能量，如沼气、秸秆、薪柴、水能等。6其他用于装饰品的一些产品（例如动物皮毛）和花卉、苗木等。7调节服务水源涵养生态系统通过其结构和过程拦截滞蓄降水，增强土壤下渗，涵养土壤水分和补充地下水、调节河川流量，增加可利用水资源量的功能。8土壤保持生态系统通过其结构与过程保护土壤、降低雨水的侵蚀能力，减少土壤流失的功能。9防风固沙生态系统通过增加土壤抗风能力，降低风力侵蚀和风沙危害的功能。10海岸带防护生态系统减低海浪，避免或减小海堤或海岸侵蚀的功能。11洪水调蓄生态系统通过调节暴雨径流、削减洪峰，减轻洪水危害的功能。12碳固定生态系统吸收二氧化碳合成有机物质，将碳固定在植物和土壤中，降低大气中二氧化碳浓度的功能。13氧气提供生态系统通过光合作用释放出氧气，维持大气氧气浓度稳定的功能。14空气净化生态系统吸收、阻滤大气中的污染物，如SO2、NOx、颗粒物等，降低空气污染浓度，改善空气环境的功能。15水质净化生态系统通过物理和生化过程对水体污染物吸附、降解以及生物吸收等，降低水体污染物浓度、净化水环境的功能。16气候调节生态系统通过植被蒸腾作用和水面蒸发过程吸收能量、降低气温、提高湿度的功能。17物种保育生态系统为珍稀濒危物种提供生存与繁衍场所的作用和价值。18文化服务休闲旅游人类通过精神感受、知识获取、休闲娱乐和美学体验、康养等旅游休闲方式，从生态系统获得的非物质惠益。19景观价值生态系统为人类提供美学体验、精神愉悦，从而提高周边土地、房产价值的功能。9表2-2生态系统生产总值实物量及价值量核算指标体系服务类别核算科目实物量指标价值量指标物质产品农业产品农业产品产量农业产品产值林业产品林业产品产量林业产品产值畜牧业产品畜牧业产品产量畜牧业产品产值渔业产品渔业产品产量渔业产品产值生态能源生态能源总量生态能源产值其他产品装饰观赏资源总量等装饰观赏资源产值等调节服务水源涵养水源涵养量水源涵养价值土壤保持土壤保持量减少泥沙淤积价值减少面源污染价值防风固沙固沙量草地恢复成本海岸带防护海岸带防护面积海岸带防护价值洪水调蓄洪水调蓄量调蓄洪水价值空气净化净化二氧化硫量净化二氧化硫价值净化氮氧化物量净化氮氧化物价值净化颗粒物量净化颗粒物价值水质净化净化COD量净化COD价值净化总氮量净化总氮价值净化总磷量净化总磷价值碳固定固定二氧化碳量碳固定价值氧气提供氧气提供量氧气提供价值气候调节植被蒸腾消耗能量植被蒸腾调节温湿度价值水面蒸发消耗能量水面蒸发调节温湿度价值物种保育珍稀濒危物种数量保护区面积珍稀濒危物种保育价值保护区保育价值文化服务休闲旅游游客总人数游憩康养价值景观价值受益土地与房产面积土地、房产升值2.4生态系统服务实物量核算生态系统生产总值实物量核算包括三大类，即物质产品实物量核算、调节服务实物量核算、文化服务实物量核算。生态系统生产总值实物量核算的核算项目、实物量指标和核算方法见表2-3。10表2-3生态系统生产总值实物量核算方法服务类别核算项目实物量指标核算方法物质产品农业产品农业产品产量统计调查林业产品林业产品产量畜牧业产品畜牧业产品产量渔业产品渔业产品产量生态能源生态能源总量其他装饰观赏资源产量等调节服务水源涵养水源涵养量水量平衡法，水量供给法土壤保持土壤保持量修正通用土壤流失方程（RUSLE）防风固沙固沙量修正风力侵蚀模型（REWQ）海岸带防护海岸带防护面积统计调查洪水调蓄湖泊：可调蓄水量水量储存模型水库：防洪库容沼泽：滞水量空气净化净化二氧化硫量污染物净化模型净化氮氧化物量净化颗粒物量水质净化净化COD量污染物净化模型净化总氮量净化总磷量碳固定固定二氧化碳量固碳机理模型氧气提供氧气提供量释氧机理模型气候调节植被蒸腾消耗能量蒸散模型水面蒸发消耗能量物种保育珍稀濒危物种数量统计调查文化服务休闲旅游游客总人数统计调查景观价值受益土地面积或公众2.5生态系统服务价值量核算在生态系统生产总值实物量核算的基础上，确定各类生态系统服务的价格，核算生态服务价值。具体而言，生态系统生产总值价值量核算中，物质产品价值主要用市场价值法核算，调节服务价值主要用替代成本法进行核算，文化服务价值使用旅行费用法。生态系统生产总值价值量核算的核算项目、价值量指标和核算方法，见表2-4。11表2-4生态系统生产总值价值量核算方法服务类别核算项目价值量指标核算方法物质产品农业产品农业产品产值市场价值法林业产品林业产品产值畜牧业产品畜牧业产品产值渔业产品渔业产品产值生态能源生态能源产值其他装饰观赏资源产值等调节服务水源涵养水源涵养价值替代成本法土壤保持减少泥沙淤积价值替代成本法减少面源污染价值替代成本法防风固沙固沙价值恢复成本法海岸带防护由于防护减少的损失价值替代成本法洪水调蓄调蓄洪水价值影子工程法空气净化净化二氧化硫价值替代成本法净化氮氧化物价值替代成本法净化颗粒物价值替代成本法水质净化净化总氮价值替代成本法净化总磷价值替代成本法净化COD价值替代成本法碳固定固定二氧化碳价值替代成本法氧气提供氧气提供价值替代成本法气候调节植被蒸腾调节温湿度价值替代成本法水面蒸发调节温湿度价值物种保育物种保育价值保育价值法文化服务休闲旅游休闲旅游价值旅行费用法景观价值景观价值享乐价格法3物质产品实物量和价值量核算3.1物质产品实物量核算3.1.1指标内涵物质产品是指人类从生态系统获取的能够在市场交易的产品，满足人类生活、生产与发展的物质需求，包括农业、林业、畜牧业、渔业等产品与生态能源。具体指标见表3-1。12表3-1生态系统的物质产品指标类别项目内容指标农业产品粮食作物谷物水稻、小麦、玉米、谷子、高粱、其他谷物等豆类大豆、绿豆、红小豆薯类马铃薯油料油料花生、油菜籽、向日葵籽、芝麻、胡麻籽棉花棉花棉花麻类麻类黄红麻、亚麻、大麻、苎麻糖类糖类甜菜、甘蔗烟叶烟叶烟叶药材药材药材茶叶茶叶红毛茶、绿毛茶蔬菜蔬菜蔬菜（含菜用瓜）瓜类瓜类西瓜、甜瓜、草莓水果水果香蕉、苹果、梨、葡萄、柑桔、红枣、柿子、菠萝、其他园林水果等林业产品木材木材木材其他林产品橡胶橡胶松脂松脂生漆生漆油桐籽油桐籽油茶籽油茶籽畜牧业产品畜禽畜禽牛、羊、猪、家禽等奶类奶类牛奶禽蛋禽蛋禽蛋动物皮毛羊毛细绵羊毛、半细绵羊毛、羊绒、山羊毛其他畜产品蜂蜜蜂蜜蚕茧柞蚕茧、桑蚕茧渔业产品淡水产品淡水产品鱼类、贝类、虾蟹类、其他海水产品海水产品鱼类、贝类、虾蟹类、藻类、其他生态能源生态能源水能水能发电量薪柴薪柴量秸秆固化产量沼气沼气量其它装饰观赏资源等装饰观赏资源等花卉、苗木等133.1.2核算方法生态系统在一定时间内提供的各类产品的产量可以通过现有的经济核算体系获得，实物产品的产量可以通过统计资料获取。𝐸𝑝𝑟𝑜=∑𝐸𝑖𝑛𝑖=1（1）式中，Epro为产品总产量（t/a）；Ei为第i种产品的产量（t/a）；i为核算区产品种类，i=1，2，3，……，n。3.1.3核算参数及数据来源产量数据来自林业、农业、渔业及统计部门相关资料，以及实地调查。3.2物质产品价值量核算3.2.1定价思路生态系统物质产品价值是指生态系统通过初级生产、次级生产为人类提供农产品、林业产品、畜牧业产品、渔业产品、生态能源等的经济价值。由于生态系统物质产品能够在市场上进行交易，存在相应的市场价格，可以运用市场价值法对生态系统的物质产品服务进行价值核算。3.2.2价值量核算模型𝑉𝑚=∑𝐸𝑖×𝑃𝑖𝑛𝑖=1（2）式中，Vm为生态系统物质产品价值（元/a）；Ei为第i类生态系统产品产量（根据产品的计量单位确定，如kg/a），Pi为第i类生态系统产品的价格（根据产品的计量单位确定，如元/kg）。143.2.3定价参数与数据来源产品产量由实物量核算得到，产品价格从林业、农业、渔业及统计部门获得或根据市场定价获得。4水源涵养实物量和价值量核算4.1水源涵养实物量核算4.1.1指标内涵水源涵养服务是生态系统拦截滞蓄降水，增强土壤下渗、蓄积，涵养土壤水分、调节暴雨径流和补充地下水，增加可利用水资源的功能。水源涵养量大的地区不仅满足核算区内生产生活的水源需求，还持续地向区域外提供水资源。选用水源涵养量，作为生态系统水源涵养实物量的评价指标。4.1.2核算方法水源涵养量核算主要有两种方法，推荐优先采用水量平衡法，技术参数缺失的情况下，可以选择水量供给法。4.1.2.1水量平衡法通过水量平衡方程计算。水量平衡方程是指在一定的时空内，水分在生态系统中保持质量守恒，即生态系统水源涵养量是降水输入与暴雨径流和生态系统自身水分消耗量的差值。𝑄𝑤𝑟=∑𝐴𝑖×(𝑃𝑖−𝑅𝑖−𝐸𝑇𝑖)×10−3𝑛𝑖=1（3）15式中，Qwr为水源涵养量（m3/a）；Pi为产流降雨量（mm/a），Ri为地表径流量（mm/a），ETi为蒸散发量（mm/a）；Ai为i类生态系统的面积（m2）；i为生态系统类型，n为生态系统类型总数。4.1.2.2水量供给法水源涵养量是生态系统为本地区和下游地区提供的水资源总量，包括本地区的用水量和净出境水量。𝑄𝑤𝑟=(𝑈𝑄𝑤−𝑇𝑄𝑤)+(𝐿𝑄𝑤−𝐸𝑄𝑤)×(1−𝛿)（4）式中，Qwr为水源涵养量（m3/a）；UQw为核算区内的用水量（包括工业、生活用水量），TQw为跨流域净调水量，LQw为区域出境水量（m3/a），EQw为区域入境水量（m3/a）；δ为区域产流径流系数。4.1.3核算参数及数据来源核算区域的产流降雨量数据通过气象部门获取，地表径流量、蒸散发量等在技术条件允许情况下进行实测，或从遥感数据及核算区域的相关文献中获取，用水量、区域出入境水量等数据通过水利或统计部门获得。4.2水源涵养价值量核算4.2.1定价思路水源涵养价值主要表现在蓄水保水的经济价值。可运用影子工程法，即模拟建设蓄水量与生态系统水源涵养量相当的水利设施，以建设该水利设施所需要的成本核算水源涵养价值。164.2.2价值量核算模型𝑉𝑤𝑟=𝑄𝑤𝑟×𝐶𝑤𝑒（5）式中，Vwr为水源涵养价值（元/a）；Qwr为核算区内总的水源涵养量（m3/a），Cwe为水资源交易市场价格，当交易市场未建立时，以水库建设的工程及维护成本或水资源影子价格（元/m3）替代，水库单位库容的工程造价及维护成本（元/m3）。4.2.3定价参数与数据来源生态系统水源涵养量由实物量核算得到。水库单位库容的工程造价及维护成本等数据来自发改委、水利等部门发布的工程预算依据，或公开发表的参考文献，并根据价格指数折算得到核算年份的价格。5土壤保持实物量和价值量核算5.1土壤保持实物量核算5.1.1指标内涵土壤保持功能是生态系统（如森林、草地等）通过林冠层、枯落物、根系等各个层次保护土壤、消减降雨侵蚀力，增加土壤抗蚀性，减少土壤流失，保持土壤的功能。选用土壤保持量，即生态系统减少的土壤侵蚀量（用潜在土壤侵蚀量与实际土壤侵蚀量的差值测度）作为生态系统水土保持功能的评价指标。其中，实际土壤侵蚀是指当前地表植被覆盖情形下的土壤侵蚀量，潜在土壤侵蚀是指没有地表植被覆盖情形下可能发生的土壤侵蚀量。175.1.2核算方法土壤保持量核算主要基于修正的通用水土流失方程（RUSLE）计算。𝑄𝑠𝑟=𝑅×𝐾×𝐿×𝑆×(1−𝐶×𝑃)（6）式中，Qsr为土壤保持量（t/a）；R为降雨侵蚀力因子，用多年平均年降雨侵蚀力指数表示；K为土壤可蚀性因子，通常用标准样方上单位降雨侵蚀力所引起的土壤流失量来表示；L为坡长因子（无量纲），S为坡度因子（无量纲），C为植被覆盖和管理因子（无量纲），P为水土保持措施因子（无量纲）。5.1.3核算参数及数据来源降雨侵蚀力因子R、土壤可蚀性因子K、坡长因子L、坡度因子S的算法以及覆盖和管理因子C以及水土保持措施因子P参考附录C.1或查阅相关文献。5.2土壤保持价值量核算5.2.1定价思路生态系统土壤保持价值主要包括减少面源污染和减少泥沙淤积两个方面的价值。生态系统通过保持土壤，减少氮、磷等土壤营养物质进入下游水体（包括河流、湖泊、水库和海湾等），可降低下游水体的面源污染。根据土壤保持量和土壤中氮、磷的含量，运用替代成本法（即污染物处理的成本）核算减少面源污染的价值。18生态系统通过保持土壤，减少水库、河流、湖泊的泥沙淤积，有利于降低干旱、洪涝灾害发生的风险。根据土壤保持量和淤积量，运用替代成本法（即水库清淤工程的费用）核算减少泥沙淤积价值。5.2.2价值量核算模型Vsr=𝑉𝑠𝑑+𝑉𝑑𝑝𝑑（7）Vsd=λ×(𝑄𝑠𝑟/𝜌)×c（8）𝑉𝑑𝑝𝑑=∑Qsr×Ci×𝑃𝑖ni=1（9）式中，Vsr为生态系统土壤保持价值（元/a）；Vsd为减少泥沙淤积价值（元/a）；Vdpd为减少面源污染价值（元/a）；Qsr为土壤保持量（t/a）；c为单位水库清淤工程费用（元/m3）；ρ为土壤容重（t/m3）；λ为泥沙淤积系数；i为土壤中氮、磷等营养物质数量，i=1，2，…，n；Ci为土壤中氮、磷等营养物质的纯含量（%）；𝑃𝑖为处理成本。5.2.3定价参数与数据来源土壤保持量由实物量核算得到。土壤容重、氮、磷、钾含量、单位水库清淤工程费、单位污染物处理成本、肥料价格等数据来源于当地土壤调查、文献、专项调查以及发改委等物价部门。6防风固沙实物量和价值量核算6.1防风固沙实物量核算6.1.1指标内涵防风固沙功能是指生态系统减少因大风导致的土壤流失和风沙危害的功能。在风蚀过程中，植被减少土壤裸露，对土壤形成保护，19减少风蚀输沙量，还可以通过根系固定表层土壤，改良土壤结构，提高土壤抗风蚀的能力，植被还可以通过增加地表粗糙度、阻截等方式降低风速、降低大风风力侵蚀和风沙危害。选用防风固沙量，即通过生态系统减少的风蚀量（潜在风蚀量与实际风蚀量的差值），作为生态系统防风固沙功能的评价指标。6.1.2核算方法防风固沙量核算主要基于修正的风力侵蚀模型（RWEQ）计算。𝑄𝑠𝑓=0.1699×(𝑊𝐹×𝐸𝐹×𝑆𝐶𝐹×𝐾′)1.3711×(1−𝐶1.3711)（10）式中，Qsf为防风固沙量（t/a）；WF为气候侵蚀因子（kg/m）；EF为土壤侵蚀因子；SCF为土壤结皮因子；K'为地表糙度因子；C为植被覆盖因子。6.1.3核算参数及数据来源气候侵蚀因子、地表糙度因子、土壤侵蚀因子、土壤结皮因子、植被覆盖因子参考附录C.2。6.2防风固沙价值量核算6.2.1定价思路根据防风固沙量和土壤沙化盖沙厚度，核算出减少的沙化土地面积；运用恢复成本法，根据单位面积沙化土地治理费用或单位植被恢复成本核算生态系统防风固沙功能的价值。206.2.2价值量核算模型𝑉𝑠𝑓=𝑄𝑠𝑓𝜌∙ℎ×𝑐（11）式中，Vsf为防风固沙价值（元/a）；Qsf为防风固沙量（t/a）；ρ为土壤容重（t/m3）；h为土壤沙化覆沙厚度（m）；c为单位治沙工程的成本（元/m3）或单位植被恢复成本（元/m2）。6.2.3定价参数与数据来源防风固沙量由实物量核算得到，土壤容重来自土壤调查或文献资料，单位治沙工程成本或单位植被恢复成本来自政府相关文件规定。7海岸带防护实物量和价值量核算7.1海岸带防护实物量核算7.1.1指标内涵海岸带防护是指滨海盐沼、红树林、珊瑚礁等生态系统减低海浪，避免或减小海堤或海岸侵蚀的功能。7.1.2核算方法运用滨海盐沼、红树林、珊瑚礁等生态系统防护或替代海堤等防护工程的长度核算。𝐷𝑐𝑙=∑𝐷𝑐𝑙𝑖𝑛1(12)式中，𝐷𝑐𝑙为生态系统防护的海岸带总长度（km/a）；𝐷𝑐𝑙𝑖为i类生态系统防护的海岸带长度（km/a）；i为研究区生态系统类型，i=1，2，…，n，无量纲；n为研究区生态系统类型数量，无量纲。217.1.3核算参数及数据来源海岸带长度从自然资源部门获得，或利用遥感数据分析结合实地调查获得。7.2海岸带防护价值量核算7.2.1定价思路运用替代成本法（即海浪防护工程建设成本）核算滨海盐沼、红树林、珊瑚礁等生态系统防风护堤的价值。7.2.2价值核算模型𝑉𝑐𝑙=∑𝐷𝑐𝑙𝑖𝑛1×𝐶𝑐𝑙𝑖（13）式中，Vcl为海岸带防护价值（元/a）；Dcli为i类生态系统防护的海岸带长度（km）；𝐶𝑐𝑙𝑖为i类生态系统海浪防护工程单位长度建设维护成本（元/km）。7.2.3核算参数及数据来源生态系统防护的海岸带长度由实物量核算得到。防护工程的单位长度建设成本从发改委、海洋等部门获得。8洪水调蓄实物量和价值量核算8.1洪水调蓄实物量核算8.1.1指标内涵洪水调蓄功能是指自然生态系统所特有的生态结构能够吸纳大量的降水和过境水，蓄积洪峰水量，削减并滞后洪峰，以缓解汛期洪22峰造成的威胁和损失的功能。原则上只核算年降雨量大于400mm地区的洪水调蓄价值。选用植被调蓄水量和洪水期滞水量（库塘、湖泊、沼泽）表征生态系统的洪水调蓄能力，即调节洪水的潜在能力。8.1.2核算方法洪水调蓄能力：𝐶𝑓𝑚=𝐶𝑣𝑐+𝐶𝑟𝑐+𝐶𝑙𝑐+𝐶𝑚𝑐（14）式中，Cfm为洪水调蓄量（m3/a）；𝐶𝑣𝑐为植被洪水调蓄量（m3/a）；Crc为库塘洪水调蓄量（m3/a），Clc为湖区洪水调蓄量（m3/a），Cmc为沼泽洪水调蓄量（m3/a）。植被洪水调蓄量：Cvc=∑(Ph-Rfi)×Sivni=1×10-3（15）式中，Cvc为植被洪水调蓄量（m3/a）；Ph为大暴雨产流降雨量（mm）；Rfi为第i种生态系统产生的地表径流量（mm）；Siv为第i种自然植被生态系统的面积（km2）；𝑖为自然植被生态系统类型，i=1,2,…,n；n为自然植被生态系统类型数量，无量纲。库塘洪水调蓄量－方法1：库塘湿地调蓄洪水的能力，通过库塘湿地的实际洪水调蓄库容来计算，库塘湿地的实际洪水调蓄库容按其总库容的35%进行计算。𝐶𝑟𝑐=0.35×𝐶𝑡（16）式中，𝐶𝑟𝑐为水库防洪库容（m3/a）；Ct为水库总库容（m3）。库塘洪水调蓄量－方法2：23在实测数据可以支撑的条件下，库塘湿地调蓄洪水的能力，可通过洪水期水库的进出水总量进行计算。𝐶𝑟𝑐=𝐶𝑖−𝐶𝑜（17）式中，𝐶𝑟𝑐为水库防洪库容（m3/a）；𝐶𝑖为洪水期水库进水总量（m3/a），𝐶𝑜为洪水期水库出水总量（m3/a）。湖泊洪水调蓄量－方法1：湖泊洪水调蓄量的大小与湖泊面积、水深、换水周期，以及流域内的水文过程密切相关。可根据湖泊水文学过程，通过汛期湖泊入湖、出湖流量随时间的变化计算湖泊在某一段时间内洪水调蓄量：𝐶𝑙𝑐=∫(𝑄𝐼−𝑄𝑂)𝑑𝑡(𝑄𝐼>𝑄𝑂)𝑡2𝑡1（18）式中，Clc为湖泊t1-t2时间段内洪水调蓄量（m3）；QI为入湖流量（m3/s），QO为出湖流量（m3/s）。湖泊洪水调蓄量－方法2：考虑区域气候条件的差异，根据《中国湖泊志》，将全国湖泊划分为东部平原、蒙新高原、云贵高原、青藏高原、东北平原与山区5个湖泊区，分区构建湖泊换水次数与补给系数的模型，通过补给系数核算湖泊换水次数；按不同湖泊区，基于湖面面积与湖泊换水次数建立湖泊水量调节量评价模型：东部平原区:𝐶𝑙𝑐=𝑒4.924×𝐴1.128×3.19（19）蒙新高原区:𝐶𝑙𝑐=𝑒5.653×𝐴0.680×0.26（20）24云贵高原区:𝐶𝑙𝑐=𝑒4.904×𝐴0.927×0.36（21）青藏高原区:𝐶𝑙𝑐=𝑒6.636×𝐴0.678×0.14（22）东北平原与山区:𝐶𝑙𝑐=𝑒5.808×𝐴0.866×0.98（23）𝐶𝑙𝑐为湖泊的洪水调蓄量（104m3/a），A为湖泊面积（km2）。沼泽洪水调蓄量：沼泽土壤具有特殊的水文物理性质，草根层和泥炭层孔隙度达72%-93%，饱和持水量达830%-1030%，在汛期可储存大量水分。沼泽湿地像一个巨大的天然蓄水库消纳洪水；同时沼泽湿地植被可拦截径流减缓洪水流速，削减和滞后洪峰。基于沼泽土壤蓄水量和地表滞水量模型核算沼泽湿地洪水调蓄量能力：𝐶𝑚𝑐=𝐶𝑠𝑤𝑠+𝐶𝑠𝑟（24）式中，Cmc为沼泽洪水调蓄量（亿m3/a）；Csws为沼泽土壤蓄水量（亿m3/a），Csr为沼泽地表滞水量（亿m3/a）。沼泽土壤蓄水量的计算公式如下：𝐶𝑠𝑤𝑠=𝑆×ℎ×𝜌×(𝐹−𝐸)×10−2/𝜌𝑤（25）式中，Csws为沼泽土壤蓄水量（亿m3/a）；S为沼泽总面积（km2），h为沼泽湿地土壤蓄水深度（m），ρ为沼泽湿地土壤容重（g/cm3），ρw为水的密度（g/cm3），F为沼泽湿地土壤饱和含水率（无量纲），E为沼泽湿地洪水淹没前的自然含水率（无量纲）。25沼泽地表滞水量的计算公式如下：𝐶𝑠𝑟=𝑆×𝐻×10−2（26）式中，Csr为沼泽地表滞水量（亿m3/a）；S为沼泽湿地总面积（km2），H为沼泽湿地地表滞水高度（m）。8.1.3核算参数及数据来源湖泊进出水量通过水利部门统计资料获得。各类植被面积、湖面面积、沼泽面积来源自然资源、水利或其他相关统计部门。出入湖流量来源于水利、水文监测站点的实际观测数据。暴雨降雨量数据来源于气象部门。暴雨径流量参考相关研究文献或进行实测。8.2洪水调蓄价值量核算8.2.1定价思路运用替代成本法（即水库的建设成本）核算自然生态系统的洪水调蓄价值。8.2.2价值量核算模型𝑉𝑓𝑚=𝐶𝑓𝑚×𝐶𝑤𝑒（27）式中，Vfm为生态系统洪水调蓄价值（元/a）；Cfm为生态系统洪水调蓄量（m3/a）；𝐶𝑤𝑒为水库单位库容的工程造价及维护成本（元/m3）。8.2.3核算参数及数据来源生态系统洪水调蓄量由实物量核算得到。水库单位库容的工程造价及维护成本等数据来自发改委、水利等部门发布的工程预算依据，或公开发表的参考文献。269碳固定实物量和价值量核算9.1碳固定实物量核算9.1.1陆地生态系统固碳9.1.1.1指标内涵生态系统固碳功能是指自然生态系统吸收大气中的二氧化碳（CO2）合成有机质，将碳固定在植物或土壤中的功能。该功能有利于降低大气中二氧化碳浓度，减缓温室效应。生态系统的固碳功能，对降低减排压力具有重要意义。选用固定二氧化碳量作为生态系统固碳功能的评价指标。生态系统固碳量为陆地生态系统固碳和岩溶固碳总量。𝑄𝐶𝑂2=𝑄𝑡𝐶𝑂2+𝑄𝑘𝐶𝑂2（28）9.1.1.2核算方法陆地生态系统固碳功能的计算方法主要有两种，根据数据可得性选择固碳速率法或净生态系统生产力法（NEP法）。方法1：固碳速率法𝑄𝑡𝐶𝑂2=𝑀𝐶𝑂2/𝑀𝐶×(𝐹𝐶𝑆+𝐺𝑆𝐶𝑆+𝑊𝐶𝑆+𝐶𝑆𝐶𝑆)（29）式中，𝑄𝑡𝐶𝑂2为陆地生态系统二氧化碳总固定量（tCO2/a）；FCS为森林（及灌丛）固碳量（tC/a），GSCS为草地固碳量（tC/a），WCS为湿地固碳量（tC/a），CSCS为农田固碳量（tC/a），𝑀𝐶𝑂2𝑀𝐶=44/12，为C转化为CO2的系数。27固碳速率法-森林及灌丛固碳量：𝐹𝐶𝑆=𝐹𝐶𝑆𝑅×𝑆𝐹×（1+𝛽）（30）式中，FCSR为森林及灌丛的固碳速率（tC·ha-1·a-1），SF为森林及灌丛面积（ha），β为森林及灌丛土壤固碳系数。固碳速率法-草地固碳量：由于草地植被每年都会枯落，其固定的碳又返还回大气或者进入土壤中，故不考虑草地植被的固碳量，只考虑草地的土壤固碳量。𝐺𝑆𝐶𝑆=𝐺𝑆𝑅×𝑆𝐺（31）式中，GSR为草地土壤的固碳速率（tC·ha-1·a-1），SG为草地面积（ha）。固碳速率法-湿地固碳量：WCS=∑SCSR𝑖×𝑆𝑊𝑖×10-2ni=1（32）式中，SCSRi为第i类水域湿地的固碳速率（gC·m-2·a-1），𝑆𝑊𝑖为第i类水域湿地的面积（ha），i=1，2，…，n。固碳速率法-农田土壤固碳量：由于农田植被每年都会被收获，其固定的碳又返还回大气或者进入土壤中，故不考虑农田植被的固碳量，只考虑农田的土壤固碳量。𝐶𝑆𝐶𝑆=(𝐵𝑆𝑆+𝑆𝐶𝑆𝑅𝑁+𝑃𝑅×𝑆𝐶𝑆𝑅𝑆)×𝑆𝐶（33）式中，CSCS为农田土壤固碳量（tC/a）；BSS为无固碳措施条件下的农田土壤固碳速率（tC·ha-1·a-1），SCSRN为施用化学氮肥和复合肥的农田土壤固碳速率（tC·ha-1·a-1），SCSRS为秸秆全部还田的农田28土壤固碳速率（tC·ha-1·a-1），PR为农田秸秆还田推广施行率（%），SC为农田面积（ha）。无固碳措施条件下的农田土壤固碳速率：BSS=NSC×BD×H×0.1（34）式中，NSC为无化学肥料和有机肥料施用的情况下，我国农田土壤有机碳的变化，单位g·kg-1·a-1；BD为各省土壤容重，H为土壤厚度。施用化学氮肥、复合肥和秸秆还田的土壤固碳速率：施用化学氮肥的固碳速率：东北农区：SCSRN=1.7385×TNF-104.03（35）华北农区：SCSRN=0.5286×TNF+1.5973（36）西北农区：SCSRN=0.6352×TNF-1.0834（37）南方农区：SCSRN=1.5339×TNF-266.7（38）式中，TNF为单位面积耕地化学氮肥、复合肥总施用量（kgN·ha-1·a-1），计算公式如下：TNF=(NF+CF×0.3)/𝑆𝑃（39）式中，NF和CF为化学氮肥和复合肥施用量（t），𝑆𝑃为耕地面积（ha）。秸秆还田的固碳速率：东北农区SCSRS=40.524×S+340.33（40）华北农区SCSRS=40.607×S+181.9（41）西北农区SCSRS=17.116×S+30.553（42）29南方农区SCSRS=43.548×S+375.1（43）式中，S为单位耕地面积秸秆还田量（t·ha-1·a-1），计算方法如下式所示：S=∑𝐶𝑌j𝑛𝑗=1×𝑆𝐺𝑅𝑗/𝑆𝑃（44）式中，CYj为作物j在当年的产量（t），𝑆𝑃为耕地面积（ha），SGRj为作物j的草谷比。方法2：净生态系统生产力法森林、草地和湿地固碳：净生态系统生产力（NEP）是定量化分析生态系统碳源/汇的重要科学指标，生态系统固碳量可以用NEP衡量。NEP广泛应用于碳循环研究中，NEP可由净初级生产力（NPP）减去异氧呼吸消耗得到，或根据NPP与NEP的相关转换系数换算得到，然后测算出陆地生态系统固定二氧化碳的质量：𝑄𝑡𝐶𝑂2=𝑀𝐶𝑂2/𝑀𝐶×𝑁𝐸𝑃（45）式中，𝑄𝑡𝐶𝑂2为陆地生态系统固碳量（t·CO2/a）；𝑀𝐶𝑂2/𝑀𝐶同上；NEP为净生态系统生产力（tC/a）。其中，净生态系统生产力NEP有两种算法：①由净初级生产力（NPP）减去异氧呼吸消耗得到：𝑁𝐸𝑃=𝑁𝑃𝑃−𝑅𝑆（46）式中，NEP为净生态系统生产力（t·C/a）；NPP为净初级生产力（t·C/a），RS为土壤呼吸消耗碳量（t·C/a）。②按照各省市NEP和NPP的转换系数，根据NPP计算得到NEP：30𝑁𝐸𝑃=𝛼∗𝑁𝑃𝑃∗𝑀𝐶6/𝑀𝐶6𝐻10𝑂5（47）式中，NEP为净生态系统生产力（t·C/a）；α为NEP和NPP的转换系数，NPP为净初级生产力（t·干物质/a），𝑀𝐶6/𝑀𝐶6𝐻10𝑂5=72/162为干物质转化为C的系数。土壤固碳：𝑄𝑆𝐶𝑂2=∑𝐴𝑖×𝑆𝑖𝑛𝑖（48）式中，𝑄𝑆𝐶𝑂2为土壤固碳量（t·CO2/a）；𝐴𝑖为不同生态系统的土壤面积（ha），𝑆𝑖为不同生态系统实测土壤固碳量（t·CO2/ha·a-1）。9.1.1.3核算参数及数据来源净初级生产力NPP、土壤呼吸消耗碳量RS、生物量数据、各类陆地生态系统面积、化学氮肥NF、复合肥施用量CF和作物j在当年的产量CYj等数据来自然资源、林业、农业和统计等部门的遥感数据、统计数据、实地调查或相关文献数据；NPP与NEP转换系数、生物量-碳转换系数、森林及灌丛固碳速率推荐值、森林及灌丛土壤固碳系数、草地土壤固碳速率、湿地固碳速率、各省土壤容重、无化学肥料和有机肥料施用的情况下我国农田土壤有机碳的变化、土壤厚度各省土壤容重以及作物j的草谷比等见附录C.3。9.1.2岩溶固碳9.1.2.1指标内涵岩溶作用，也称喀斯特作用，是指地下水和地表水对可溶性岩石进行以化学作用为主的溶蚀，每形成1molHCO3-需要从大气吸收310.5molCO2，将碳以HCO3-形式存储在水体中。我国广大的岩溶面积以及岩溶形成机制对吸收大气CO2具有一定作用。选用固定CO2量作为岩溶固碳功能的评价指标。岩溶固碳仅在含有岩溶地貌的区域选择核算。9.1.2.2核算方法岩溶地区固碳量可以用以下公式核算：𝑄𝑘𝐶𝑂2=12×𝑆×𝑀×𝐶[HCO3−]×𝑀CO2𝑀HCO3−（49）式中，𝑄𝑘𝐶𝑂2为岩溶作用总固碳量（tCO2/a）；S为岩溶地区面积（km2），M为区域地下水径流模数（107L/(km2·a)），𝐶[HCO3−]为[HCO3-]浓度（g/L）；MCO2为CO2的摩尔质量44，MHCO3-为HCO3-的摩尔质量61。根据相关数据，式（49）可简化为式（50）：𝑄𝑘𝐶𝑂2=0.3607×𝑆×𝑀×𝐶[HCO3−]（50）系数0.3607由CO2摩尔质量（MCO2）44，HCO3-摩尔质量（MHCO3-）61，每形成1摩尔质量HCO3-需要从大气吸收0.5摩尔质量CO2这三个系数确定。9.1.2.3核算参数及数据来源区域岩溶面积数据可从地方自然资源、水文和地质部门获得，区域地下水径流模数以及HCO3-浓度根据地方水文部门数据确定或实测获得，附录C.3供参考。329.2碳固定价值量核算9.2.1定价思路生态系统固碳价值可以采用替代成本法（造林成本法、工业减排成本）与市场价值法（碳交易价格）核算生态系统固碳的价值。9.2.2价值量核算模型VCf=𝑄𝐶𝑂2×CC（51）式中，VCf为生态系统固碳价值（元/a）；𝑄𝐶𝑂2为生态系统固碳总量（tCO2/a）；CC为碳价格（元/t）。9.2.3定价参数与数据来源固碳量由实物量核算得到。单位造林固碳成本、工业碳减排成本、碳交易市场价格参考相关文献，建议采用碳市场交易价格。10氧气提供实物量和价值量核算10.1氧气提供实物量核算10.1.1指标内涵生态系统的释氧功能指植物在光合作用过程中，释放出氧气的功能。这种功能对于维护大气中氧气的稳定，改善人居环境具有重要意义。选用释氧量作为生态系统释氧功能的评价指标。10.1.2核算方法根据光合作用化学方程式可知，植物每生产吸收1molCO2，就会释放1mol氧气。NEP可由净初级生产力（NPP）减去异氧呼吸消耗33得到，或根据NPP与NEP的相关转换系数获得，然后测算出生态系统释放氧气的质量：𝑄𝑜𝑝=𝑀𝑂2/𝑀𝐶𝑂2×𝑄𝐶𝑂2（52）式中，Qop为生态系统释氧量（t氧气/a）；𝑀𝑂2/𝑀𝐶𝑂2=32/44为𝐶𝑂2转化为O2的系数，𝑄𝐶𝑂2为生态系统固碳量（tC/a）。10.1.3核算参数及数据来源生态系统固碳量𝑄𝐶𝑂2来源于9.1。10.2氧气提供价值量核算10.2.1定价思路采用市场价值法（即制氧价格）核算生态系统提供氧气的价值。10.2.2价值量核算模型Vop=𝑄𝑜𝑝×𝐶𝑜（53）式中，Vop为生态系统释氧价值（元/a）；Qop为生态系统氧气释放量（t氧气/a），𝐶𝑜为工业制氧价格（元/t）。10.2.3定价参数与数据来源生态系统释氧量由实物量核算得到，制氧价格根据实际调查获得。3411空气净化实物量和价值量核算11.1空气净化实物量核算11.1.1指标内涵空气净化功能是指生态系统吸收、过滤、阻隔和分解大气污染物（如二氧化硫、氮氧化物、颗粒物等），净化空气污染物，改善大气环境的功能。空气净化功能主要体现在净化污染物和阻滞颗粒物方面。11.1.2核算方法空气净化功能核算依据污染物浓度是否超过环境空气功能区质量标准而选择不同的方法。若污染物浓度未超过环境空气功能区质量标准，则采用方法1核算污染物净化量；若污染物浓度超过环境空气功能区质量标准，则采用方法2核算污染物净化量。环境空气功能区质量标准见附录C.4。空气净化服务－方法1：如果污染物排放量未超过环境空气功能区质量标准，则采用污染物排放量核算实物量。𝑄𝑎𝑝=∑𝑄𝑖𝑛𝑖=1（54）式中，Qap为大气污染物排放总量（kg/a）；Qi为第i类大气污染物排放量（kg/a）；i为污染物类别，i=1，2，…，n，无量纲；n为大气污染物类别的数量，无量纲。空气净化服务－方法2：如果污染物排放量超过环境空气功能区质量标准，则采用生态系统自净能力核算实物量。𝑄𝑎𝑝=∑∑𝑄𝑖𝑗×𝐴𝑖𝑛𝑗=1𝑚i=1（55）35式中，Qap为生态系统空气净化能力（kg/a）；Qij为第i类生态系统第j种大气污染物的单位面积净化量（kg/km2·a），i为生态系统类型，i=1，2，…，m，无量纲。j为大气污染物类别，j=1，2，…，n，无量纲。Ai为第i类生态系统类型面积（km2）；m为生态系统类型的数量，无量纲；n为大气污染物类别的数量，无量纲。11.1.3核算参数及数据来源污染物排放数据从生态环境部门获取；各类生态系统面积来源于自然资源部门；生态系统对污染物的单位面积净化量来源于参考文献或实地监测。11.2空气净化价值量核算11.2.1定价思路生态系统空气净化价值是指生态系统吸收、过滤、阻隔和分解降低大气污染物（如二氧化硫、氮氧化物、颗粒物等），使大气环境得到改善产生的生态效应。采用替代成本法（工业治理大气污染物成本），核算生态系统空气净化价值。11.2.2价值量核算模型二氧化硫、氮氧化物、颗粒物净化价值计算方法：运用二氧化硫、氮氧化物、颗粒物三种污染物空气净化实物量，分别乘以单位二氧化硫、氮氧化物、颗粒物处理的费用，核算空气净化价值。𝑉𝑎=∑𝑄𝑎𝑝𝑖×𝐶𝑖𝑛𝑖=1（56）36式中，Va为生态系统大气环境净化的价值（元/a）；𝑄𝑎𝑝𝑖为第i种大气污染物的净化量（t/a），j为大气污染物类别，i=1，2，…，n，无量纲；𝐶𝑖为第i类大气污染物的治理成本（元/t）。11.2.3定价参数与数据来源污染物净化量由实物量核算得到。单位治理成本采用核算地方印发的排污费征收标准，没有地方标准的，可以参考国家发展和改革委员会发布的《排污费征收标准及计算方法》收费标准或者《中华人民共和国环境保护税法》中的税收标准。12水质净化实物量和价值量核算12.1水质净化实物量核算12.1.1指标内涵水质净化功能是指湖泊、河流、沼泽等水域湿地生态系统吸附、降解、转化水体污染物，净化水环境的功能。水质净化服务价值核算主要是利用监测数据，根据水体生态系统中污染物构成和浓度变化，选取适当的指标对其进行定量化核算。常用指标包括氨氮、COD、总氮、总磷以及部分重金属等。根据我国《地表水环境质量标准》（GB3838-2002）中对水环境质量应控制的项目和限值的规定，选取指标作为生态系统水环境净化功能的评价指标。12.1.2核算方法水质净化功能核算依据污染物浓度是否超过地表水水域环境功能和保护目标而选择不同的方法。若污染物浓度未超过地表水水域环37境功能标准限值，则采用方法1进行核算；若污染物排放浓度超过地表水水域环境功能标准限值，则采用方法2进行核算。地表水水域环境功能标准限值见附录C.5表C.10。水质净化功能－方法1：如果污染物排放量未超过地表水水域环境功能标准限值，则采用污染物排放量核算实物量。𝑄𝑤𝑝=∑𝑃𝑖𝑛𝑖=1（57）式中，Qwp为水体污染物净化量（kg/a）；𝑃𝑖为i类污染物排放量（kg/a），包括总氮、总磷、COD等；i为污染物类别，i=1，2，…，n，无量纲；n为水体污染物类别的数量，无量纲。水质净化功能－方法2：如果污染物排放量超过地表水水域环境功能标准限值，根据生态系统的自净能力核算实物量：𝑄𝑤𝑝=∑∑𝑃𝑖𝑗×𝐴𝑖𝑛𝑗=1𝑚i=1（58）式中，Qwp为污染物净化总量（kg）；𝑃𝑖𝑗为某种生态系统单位面积污染物净化量（kg/km2）；𝐴𝑖为生态系统面积（km2）；i=1，2，…，m，无量纲；m为生态系统类型的数量，无量纲；j为污染物类别，j=1，2，…，n，无量纲；n为水体污染物类别的数量，无量纲。12.1.3核算参数及数据来源污染物排放数据从生态环境部门获取；各类生态系统面积来源于自然资源部门；生态系统对污染物的单位面积净化量来源于参考文献或实地监测。3812.2水质净化价值量核算12.2.1定价思路核算生态系统降解水体污染物、净化水质的价值，运用替代成本法核算生态系统水质净化功能的价值。12.2.2价值量核算模型水质净化价值量计算采用替代成本法，通过工业治理水污染物成本核算生态系统水质净化价值。化学需氧量、氨氮净化价值计算方法：运用化学需氧量、氨氮两种污染物水质净化实物量，分别乘以单位化学需氧量、氨氮处理的费用，核算水体净化价值。𝑉𝑤=∑𝑄𝑤𝑝𝑖𝑛𝑖=1×𝐶𝑖（59）式中，𝑉𝑤为生态系统水质净化的价值（元/a）；𝑄𝑤𝑝𝑖为第i类水污染物的净化量（t/a）；𝐶𝑖为第i类水污染物的单位治理成本（元/t）；i为研究区第i类水体污染物类别，i=1，2，…，n，无量纲。n为研究区水体污染物类别的数量，无量纲。12.2.3定价参数与数据来源污染物净化量由实物量核算得到。COD、氨氮水质等污染物单位治理成本采用核算地方印发的排污费征收标准，没有地方标准的，可以参考国家发展和改革委员会发布的《排污费征收标准及计算方法》收费标准或者《中华人民共和国环境保护税法》中的税收标准。3913气候调节实物量和价值量核算13.1气候调节实物量核算13.1.1指标内涵生态系统气候调节服务是指生态系统通过植被蒸腾作用、水面蒸发过程吸收太阳能，降低气温、增加空气湿度，改善人居环境舒适程度的生态功能。选用生态系统蒸散发过程消耗的能量作为生态系统气候调节服务的评价指标。13.1.2核算方法气候调节服务核算可用实际测量生态系统内外温差、生态系统消耗的太阳能量和生态系统的总蒸散量进行核算，优先选择实际测量方法，其次根据数据可得性选取生态系统消耗的太阳能量方法或生态系统的总蒸散量进行核算。方法一：采用实际测量生态系统内外温差进行实物量转换。𝑄=∑𝛥𝑇𝑖×𝜌𝑐×𝑉𝑛𝑖=1（60）式中，𝑄为吸收的大气热量（J/a）；𝜌𝑐为空气的比热容（J/m3·℃-1）；𝑉为生态系统内空气的体积（m3）；𝛥𝑇𝑖为第𝑖天生态系统内外实测温差（℃）；𝑛为年内日最高温超过26℃的总天数。方法二：采用生态系统消耗的太阳能量作为气候调节的实物量。40𝐶𝑅𝑄=𝐸𝑇𝐸−𝑁𝑅𝐸（61）式中，CRQ为生态系统消耗的太阳能量（J/a）；ETE为森林、草地、灌丛、湿地等生态系统蒸腾作用消耗的太阳能量（J/a）；NRE为森林、草地、湿地等生态系统吸收的太阳净辐射能量（J/a）。方法三：采用生态系统蒸腾蒸发总消耗的能量作为气候调节的实物量。𝐸𝑡𝑡=𝐸𝑝𝑡+𝐸𝑤𝑒（62）𝐸𝑝𝑡=∑𝐸𝑃𝑃𝑖×𝑆𝑖×𝐷3𝑖×106/(3600∗𝑟)（63）𝐸𝑤𝑒=𝐸𝑤×𝑞×103/(3600)+𝐸𝑤×𝑦（64）式中，Ett为生态系统蒸腾蒸发消耗的总能量（kW·h/a）；Ept为生态系统植被蒸腾消耗的能量（kW·h/a）；Ewe为湿地生态系统蒸发消耗的能量（kW·h/a）；EPPi为i类生态系统单位面积蒸腾消耗热量（kJ.m-2d-1）；Si为i类生态系统面积（km2）；D为日最高气温大于26℃天数；r为空调能效比：3.0，无量纲；i为生态系统类型（森林、灌丛、草地）；Ew为蒸发量（m3）；q为挥发潜热，即蒸发1克水所需要的热量（J/g）；𝑦为加湿器将1m3水转化为蒸汽的耗电量（kw·h），仅计算湿度小于45%时的增湿功能。13.1.3核算参数及数据来源水面蒸发量、植被蒸散发量、生态系统面积、单位面积蒸腾耗热量等数据来自气象、自然资源、林业等相关部门和文献资料。4113.2气候调节价值量核算13.2.1定价思路运用替代成本法（即人工调节温度和湿度所需要的耗电量）核算生态系统蒸腾调节温度或湿度价值和水面蒸发调节温度或湿度价值。13.2.2价值量核算模型𝑉𝑡𝑡=𝐸𝑡𝑡×𝑃𝑒（65）式中，Vtt为生态系统气候调节的价值（元/a）；𝐸𝑡𝑡为生态系统调节温度或湿度消耗的总能量（kW·h/a）；𝑃𝑒为当地电价（元/kW·h）。13.2.3定价参数与数据来源生态系统调节温度或湿度所耗能量由实物量核算得到。电价从核算地方发展与改革委员会发布的相关文件或供电部门获取，一般参考工业电价。14物种保育实物量和价值量核算14.1物种保育实物量核算14.1.1指标内涵物种多样性是生物多样性最主要的结构和功能单位，可以为生态系统演替与生物进化提供必须的物种与遗传资源，是人类生存和发展的基础。物种保育服务是指生态系统为珍稀濒危动植物物种提供生存与繁衍的场所，从而对其起到保育作用的作用和价值。4214.1.2核算方法方法1：Shannnon-Weiner指数法把濒危动植物、特有动植物和古树名木的数量纳入到计算中，公式如下：𝐺𝑏𝑖𝑜=𝐴×(1+0.1×∑𝐸𝑚𝑥𝑚=1+0.1×∑𝐵𝑛𝑦𝑛=1+0.1×∑𝑂𝑟𝑧𝑟=1)(66)式中，Gbio物种保育的实物量；Em区域内物种m的濒危分值；Bn区域内物种n的特有值；Or区域内物种r的古树年龄指数；x为计算濒危指数物种数量；y为计算特有种指数物种数量；z为计算古树年龄指数物种数量；A为群落面积（ha）。方法2：保护区保护法采用区域保护区面积进行核算。𝐺𝑏𝑖𝑜𝑝=𝑆(67)式中，Gbiop物种保育的实物量；S为自然保护区面积。14.2物种保育价值量核算14.2.1定价思路运用单位面积保育成本核算物种保育服务的价值量。14.2.2价值量核算模型实物量与物种保育成本相乘得到价值量：方法1：Shannnon-Weiner价值法实物量与物种保育价值相乘得到价值量：𝑉𝑏𝑖𝑜=𝐺𝑏𝑖𝑜×𝑆生（68）43式中，𝑉𝑏𝑖𝑜生物多样性价值（元/a）；𝑆生为单位面积物种保育价值（元/ha·a-1）。方法2：保护区保护价值法𝑉𝑏𝑖𝑜𝑝=𝐺𝑏𝑖𝑜𝑝×𝑆𝑐（69）式中，𝑉𝑏𝑖𝑜𝑝物种保育价值（元/a）；G为物种保育实物量；𝑆𝑐为自然保护区单位面积保育成本（元/ha·a-1）。14.2.3核算参数及数据来源保护区总投入、物种数量参照相关统计、报告和文献资料。不同Shannnon-Weiner指数对应等级划分，以及相应的单位面积物种多样性保护价值参照附录C.6。15休闲旅游实物量和价值量核算15.1休闲旅游实物量核算15.1.1指标内涵人类通过精神感受、知识获取、休闲娱乐和美学体验从生态系统获得的非物质惠益。15.1.2核算方法采用区域内自然景观的年旅游总人次作为文化服务的实物量评价指标。Nt=∑Ntini=1（70）式中，Nt为游客总人数；Nti为第i个旅游区的人数；n为旅游区44个数，i=1，2，…，n。15.1.3核算参数及数据来源自然景观名录、旅游人数与旅客来源通过旅游、园林、统计等部门或问卷调查获取。15.2休闲旅游价值量核算15.2.1定价思路运用旅行费用法核算人们通过休闲旅游活动体验生态系统与自然景观美学价值，并获得知识和精神愉悦的非物质价值。15.2.2价值量核算模型𝑉𝑟=∑𝑁𝑗×𝑇𝐶𝑗𝐽𝑗=1（71）𝑇𝐶𝑗=𝑇𝑗×𝑊𝑗+𝐶𝑗（72）𝐶𝑗=𝐶𝑡𝑐,𝑗+𝐶𝑙𝑓,𝑗+𝐶𝑒𝑓,𝑗（73）式中，𝑉𝑟表示被核算地点的休闲旅游价值（元/a）；𝑁𝑗表示j地到核算地区旅游的总人数（人/a）；𝑗=1,2…，𝐽表示来被核算地点旅游的游客所在区域（区域按距离核算地点的距离划同心圆，如省内、省外等）。𝑇𝐶𝑗表示来自j地的游客的平均旅行成本（元/人）；𝑇𝑗表示来自j地的游客用于旅途和核算旅游地点的平均时间（天/人）；𝑊𝑗表示来自j地的游客的当地平均工资（元/（人·天））；𝐶𝑗表示来自j地的游客花费的平均直接旅行费用（元/人），其中包括游客从j地到核算区域的交通费用𝐶𝑡𝑐,𝑗（元/人）、食宿花费𝐶𝑙𝑓,𝑗（元/人）和门票费用𝐶𝑒𝑓,𝑗（元/人）。4515.2.3定价参数与数据来源自然景观名录、旅游人数通过旅游、园林等部门获取，游客的社会经济特征、旅行费用情况等通过问卷调查获得。16景观价值实物量和价值量核算16.1景观价值实物量核算16.1.1指标内涵生态系统的景观价值是指森林、湖泊、河流、海洋等生态系统可以为其周边的人群提供美学体验、精神愉悦的功能，从而提高周边土地、房产价值。16.1.2核算方法采用能直接从自然生态系统获得景观价值的土地与居住小区房产面积作为景观价值实物量评价指标。Al=∑Alini=1（74）式中，𝐴𝑙为从自然生态系统景观获得升值的土地与居住小区房产总面积（km2/a）；𝐴𝑙𝑖为第i区的房产面积（km2），i=1，2，…，n。16.1.3核算参数及数据来源受益土地与居住区名录及面积通过调查获取。4616.2景观价值价值量核算16.2.1定价思路运用享乐价值法核算生态系统为其周边地区人群提供美学体验、精神愉悦功能的价值。16.2.2价值核算模型𝑉𝑎=𝐴𝑎×𝑃𝑎（75）式中，𝑉𝑎为景观价值（元/a）；𝐴𝑎为受益总面积（km2）；𝑃𝑎为由生态系统带来的单位面积溢价（元/km2·a-1）。16.2.3核算参数及数据来源受益面积由实物量核算得到。生态系统带来的单位面积溢价由实地调研获取。17附则17.1核算数据基础数据是开展生态系统生产总值核算的重要保障，涉及的数据内容包括发改、农业、国土、林业、水务、气象、环保、旅游、电力、统计等部门数据。基础数据获取时，应确保数据的权威性、准确性、时效性和可获得性。根据核算需求，应从数据的完整性、代表性、准确性与可比性等角度进行分析筛选，采用权威部门的遥感、监测、普查、调查、统计等数据，尽可能使用高分辨率空间数据，选取适合核算工作需要的数据资料（参考来源见附录A）。同时，在核算过程中，47应保证所处理数据资料的完整性、合理性和有效性，并对原始数据、拷贝数据及数据分析结果进行校核，数据处理应符合相应规范的要求。为使核算报告准确可信，在核算过程中应对数据的获取与处理进行质量控制，降低不确定性。在核算过程中，要确保使用的模型和收集到的数据能够代表实际情况。17.2成果表达形式核算成果主要包括核算报告、生态系统生产总值核算总表。核算总表用表格形式表达GEP各项指标的实物量与价值量。17.3结果校核核算完成后，应该将核算结果与其他地区已有的核算结果对比，对比不同生态系统服务功能单位面积实物量、价值量的差异，验证结果的合理性。17.4信息管理核算主体应记录并保存下列资料，保存时间不少于5年：（1）核算方法相关信息：选择基于模型的方法时，应保存以下内容：a）获取数据和参数的相关资料；b）不确定性及如何降低不确定性的相关说明。选择基于实地调查的方法时，应保存以下内容：a）有关职能部门出具的测量仪器证明文件；b）连续测量的所有原始数据；c）不确定性及如何降低不确定性的相关说明；d）验证计算，应保留所有基于计算的保存内容。（2）数据质量控制相关记录文件。48附录A（资料性附录）生态系统生产总值核算数据来源表表A.1生态系统生产总值核算数据来源表服务类别评估科目评估指标单位数据来源物质产品农业产品谷物谷物产量、单价、产值-统计部门，农业部门薯类薯类油料油料豆类豆类糖料糖料蔬菜蔬菜水果水果…………（根据评估样点实际产品更新）小计林业产品竹木采伐木材统计部门，林业部门竹材林下产品五倍籽核桃板栗松子……（根据评估样点实际产品更新）小计畜牧业产品畜禽产量牛肉统计部门,畜牧业部门羊肉猪肉禽肉猪出栏数牛出栏数羊出栏数家禽出栏数……（根据评估样点实际产品更新）奶类牛奶……禽蛋禽蛋小计49渔业产品水产品养殖统计部门捕捞……（根据评估样点实际产品更新）小计生态能源生态能源水能统计部门生物质能源（秸秆、薪柴等）小计调节服务水源涵养气候调节土壤保持防风固沙森林面积km²自然资源部门灌丛面积km²草地面积km²湖泊面积km²沼泽面积km²河流面积km²湿地面积km²水田面积km²旱地面积km²园地面积km²多年均降雨量mm气象部门多年均蒸发量mm年降雨量mm年蒸发量mm多年月均温度℃多年月均湿度%年均夏季气温℃土壤容重t/m3自然资源部门水源涵养洪水调蓄土壤保持防风固沙水库单位库容造价元/m3水利部门水务部门农业部门环境部门区域入境水量亿m3区域出境水量亿m3水库清淤工程费用元/m3土壤中氮的纯含量%土壤中磷的纯含量%环境工程降解氮工程费用元/t环境工程降解磷工程费用元/t水库总库容万m3水库防洪库容万m3碳固定氧气提NPP亿t林业部门生物量亿t50供空气净化二氧化硫排放量万t/a生态环境部门发改部门氮氧化物排放量万t/a颗粒物排放量万t/a二氧化硫治理费用元/t氮氧化物治理费用元/t颗粒物治理费用元/t水质净化废水排放总量万t/a工业废水排放量万t/a生活污水排放量万t/a污水处理总量（污水厂）万t/aCOD排放量万t/a污水处理厂处理量t/a总氮排放量万t/a污水处理厂处理量t/a总磷排放量万t/a污水处理厂处理量t/aCOD治理费用元/t氨氮治理费用元/t因生态系统降低水环境治理成本的公司数量个问卷调查统计部门因生态系统降低水环境治理的费用万元气候调节电价元/kW·h电力部门文化服务休闲旅游价值旅游人数万人/年旅游发展改革委员会旅游收入万元/年景观价值受益土地面积万亩自然资源部门、建设部门受益居住区名录51附录B（资料性附录）陆地生态系统生产总值核算报告编写大纲前言介绍核算背景、目的意义、任务来源等1区域概况介绍核算区域地理范围、自然环境状况、经济社会状况、生态环境保护状况等基本情况。2核算目标与原则介绍核算目标、核算原则、核算依据、核算基准年。3核算方法与数据介绍主要核算思路、方法、数据来源与数据处理方法。4生态系统生产总值实物量介绍生态系统生产总值实物量的计算过程与结果，包括物质产品实物量、调节服务实物量和文化服务实物量。5生态系统生产总值价值量介绍生态系统生产总值价值量的计算过程与结果，包括物质产品价值量、调节服务价值量和文化服务价值量。6结论与建议介绍生态系统生产总值核算的结论，提出相关政策建议。7附件包括生态系统生产总值核算过程中相关的技术资料及附表、附图等。52附录C（资料性附录）生态系统生产总值实物量核算参数参考表附录C.1土壤保持实物量核算参数降雨侵蚀力因子R、土壤可蚀性因子K、坡长坡度因子L、S的算法以及植被覆盖和管理因子C以及水土保持措施因子P的计算方法可参照下述方法：（1）降雨侵蚀力因子（R）降雨侵蚀力是土壤侵蚀的驱动因子，与土壤侵蚀强度有直接的关系。降雨侵蚀力R计算可分为EI30经典计算方法和常规气象资料简易算法两类。由于降雨动能E和30min降雨强度I30资料获取难度较大，所以国内外许多学者根据区域性降雨侵蚀特点，建立了基于常规降雨量资料的简易模型。可采用全国日降雨量拟合模型来核算降雨侵蚀力，是基于日降雨量资料的半月降雨侵蚀力模型。其公式如下：𝑀𝑖=α∑𝐷𝑗𝛽𝑘𝑗=1其中，Mi为某半月时段的降雨侵蚀力值(MJ·mm·hm-2·h-1·a-1)；Dj表示半月时段内第j天的侵蚀性日雨量（要求日雨量大于等于12mm，否则以0计算，阈值12mm与中国侵蚀性降雨标准一致；k表示半月时段内的天数，半月时段的划分以每月第15日为界，每月前15天作为一个半月时段，该月剩下部分作为另一个半月时段，将全年依次划分为24个时段。α、β是模型待定参数：𝛽=0.8363+18.144𝑃𝑑12+24.455𝑃𝑦1253𝛼=21.586𝛽−7.1891其中，Pd12表示日雨量12mm以上（包括等于12mm）的日平均雨量；Py12表示日雨量12mm以上(包括12mm)的年平均雨量。（2）土壤可蚀性因子（K）土壤是土壤侵蚀发生的主体，土壤可蚀性是表征土壤性质对侵蚀敏感程度的指标，即在标准单位小区上测得的特定土壤在单位降雨侵蚀力作用下的土壤流失率，推荐采用Nomo图法。根据主要土壤性质，建立了K值与土壤性质之间的诺谟图Nomo模型。其计算公式如下：𝐾=[2.1×10−4(12−𝑂𝑀)𝑀1.14+3.25(𝑆−2)+2.5(𝑃−3)]/100×0.1317其中，K为土壤可蚀性值；OM为土壤有机质含量百分比（%）；M为土壤颗粒级配参数，为美国粒径分级制中（粉粒+极细砂）与（100-粘粒）百分比之积；S为土壤结构系数；P为渗透等级。粒径等级：粘粒为（<0.002mm）；粉粒为（0.002-0.05mm）；极细砂为（0.05-0.1mm）；砂粒为（0.1-2.0mm）。表C.1诺谟图中结构性指数与可渗透性指数的定义结构性指数S含义可渗透性指数P含义1非常坚固（verystructuredorparticulate）1快速（Papid）2很坚固（fairlystructured）2中快速(Moderatetorapid)3较坚固(slightlystructured)3中速(Moderate)4坚固(solid)4中慢速(Moderatetoslow)5慢速(Slow)6极慢(Veryslow)54（3）坡长和坡度因子（L、S）由于坡度和坡长因子相互之间联系较为紧密，因此通常将它们作为一个整体进行考虑。坡长因子是指在其它条件相同的情况下，某一长度的田块破面上的土壤流失量与72.6英尺（标准单位小区的长度）长坡面上的流失量的比值；坡度因子是指在其它条件相同的情况下，某一坡度的田块坡面上的土壤流失量与9%（标准单位小区的坡度）坡度的坡面上流失量的比值。坡度坡长因子的算法可采用建立在以下核心算法为：𝑆={10.8𝑠𝑖𝑛𝜃+0.03𝑡<9%16.8𝑠𝑖𝑛𝜃−0.59%≤𝑡<18%21.91𝑠𝑖𝑛𝜃−0.96𝑡≥18%𝐿=(𝜆22.13)𝑚𝑚=𝛽1+𝛽𝛽=(𝑠𝑖𝑛𝜃/0.089)[3.0×(𝑠𝑖𝑛𝜃)0.8+0.56]其中，L为坡长因子；𝜆为坡长，通过汇流计算得到，m为无量纲常数，S为坡度因子，t为百分比坡度，𝜃为坡度，单位是弧度。（4）植被覆盖和管理因子（C）C因子是指在一定的植被覆盖度和管理措施下，一定面积土地上的土壤流失量与采取连续清耕、休闲处理的相同面积土地上的流失量的比值，为无量纲数，介于0-1之间。要确定C因子的值，需要详细的气候、土地利用、前期作物残留量、土壤湿度等资料，在大尺度研究中，一般难以获取这些资料，且C值的经典算法非常复杂，国内部分学者采用植被覆盖度求解C值，可采用以下C值计算方法：55C={1𝑓=0.10.6508−0.3436𝑙𝑔𝑓0𝑓>78.3%0.1<𝑓≤78.3%上述公式中，植被覆盖度f基于植被指数NDVI数据计算得到，公式如下：𝑓=(𝑁𝐷𝑉𝐼−𝑁𝐷𝑉𝐼𝑠𝑜𝑖𝑙)(𝑁𝐷𝑉𝐼𝑚𝑎𝑥−𝑁𝐷𝑉𝐼𝑠𝑜𝑖𝑙)其中，𝑁𝐷𝑉𝐼𝑠𝑜𝑖𝑙为纯裸土象元的NDVI值；𝑁𝐷𝑉𝐼𝑚𝑎𝑥纯植被象元的NDVI值。附录C.2防风固沙实物量核算参数气候侵蚀因子、地表糙度因子、土壤侵蚀因子、土壤结皮因子、植被覆盖因子的计算方法可参照下述方法：（1）气候侵蚀因子气候侵蚀因子中的风因子和土壤湿度因子可采用从中国气象科学数据共享服务网上下载（http://cdc.cma.gov.cn）的国家台站的日均风速、降水、温度、日照时数，纬度等来计算完成；雪盖因子可采用从中国西部环境与生态科学数据中心(http://westdc.westgis.ac.cn)下载的中国雪深长时间序列数据集来计算雪盖因子。式中，WF为气候因子，kg/m；WS2为2m处风速，m/s；WSt为2m处临界风速（假定为5m/s）；N为风速的观测次数（一般500次）；Nd为试验的天数，d；ρ为空气密度，kg/m3；g为重力加速度，m/s2；SW为土壤湿度因子无量纲；SD为雪覆盖因子。SDSWgNNWS-WSWSWFNid2t12256i)其中空气密度通过下式求得：式中，EL为海拔高度，单位为km，T为绝对温度，单位为开氏度。ii)土壤湿度通过下式求得：Etp可采用SamaniandPessaralkli（1986）的方法，公式如下：式中：SW为土壤湿度因子；ETp为潜在相对蒸发量，mm；R为降雨量，mm；I为灌溉量，mm；Rd为降雨次数和（或）灌溉天数：Nd为天数，d（一般15d）；SR为太阳辐射总量，cal/cm2；DT为平均温度，℃。其中太阳辐射采用方法如下：𝑅𝑠0=(𝑎𝑠+𝑏𝑠𝑛𝑁)𝑅𝑎𝑁=24𝜋∙𝜔𝑠𝑅𝑎=24(60)𝜋∙𝐺𝑠𝑟∙𝑑𝑟∙[𝜔𝑠∙𝑠𝑖𝑛(𝜑)∙𝑠𝑖𝑛(𝛿)+𝑐𝑜𝑠(𝜑)∙𝑐𝑜𝑠(𝛿)∙𝑠𝑖𝑛(𝜔𝑠)]𝑑𝑟=1+0.033∙𝑐𝑜𝑠(2𝜋365∙𝐽)𝛿=0.409∙𝑠𝑖𝑛(2𝜋365∙𝐽−1.39)𝜔𝑠=𝑎𝑟𝑐𝑐𝑜𝑠[−𝑡𝑎𝑛(𝜑)∙𝑡𝑎𝑛(𝛿)]T0.0048EL0.1183EL-1.013348.02pddpETNRIRETSW)8.17(5.580162.0DTSRETp57或𝜔𝑠=𝜋2−𝑎𝑟𝑐𝑡𝑎𝑛[−𝑡𝑎𝑛(𝜑)∙𝑡𝑎𝑛(𝛿)𝑋0.5]𝑋=1−[𝑡𝑎𝑛(𝜑)]2∙[𝑡𝑎𝑛(𝛿)]2𝑋=0.00001𝑖𝑓𝑋≤0其中，Rs：太阳辐射（MJ·m-2·day-1）；Rs0：晴空下的太阳辐射（MJ·m-2·day-1）；as和bs的取值：最好以当地校正结果为准，在无实测校正地区推荐采用as为0.25，bs为0.50的标准值。Ra：地球外辐射（MJ·m-2·day-1）；n：实际日照时数；N：最大可能日照时数；Gsr：太阳常数（0.0820MJ·m-2·min-1）；dr：日地相对距离；ωs：日落时角（rad）；φ：纬度（rad）；δ：太阳倾角（rad）；J：对应于一年中的第几天。iii)雪盖因子通过下式求得：其中，P(snowcover>25.4mm)为计算时段内积雪覆盖深度大于25.4mm的概率。（2）土壤侵蚀因子土壤颗粒分为可蚀性土粒和非可蚀性土粒，粒径大于0.84mm的土粒不易被风蚀，称为非可蚀性颗粒。土壤可蚀性因子则为土壤表层直径小于0.84mm的颗粒的含量。可用以下方程来描述土壤可蚀性因子EF的值：其中，Sa为土壤砂粒含量；Si为土壤粉砂含量；Sa/Cl为土壤砂粒和粘土含量比；OM为有机质含量；CaCO3为碳酸钙含量。25.4mmcoversnowP-1SD10095.059.2/33.017.031.009.293CaCOOMClSaSiSaEF58（3）土壤结皮因子土壤结皮（soilcrust）为土壤颗粒物（特别是黏土、粉砂与有机质颗粒）的胶结作用而在土壤表面生成一层物理、化学和生物性状均较特殊的土壤微层。其中，CL为粘土含量（5.0～39.3），OM为有机质含量（0.18～4.79）。表C.2RWEQ标准数据库中物质含量范围表SaSiClSa/ClOMCaCO3范围（%）5.5～93.60.5～69.55.0-39.3l.2～53.00.18～4.790.0～25.2土壤数据可采用中国西部环境与生态科学数据中心(http://westdc.westgis.ac.cn)提供的1：100万土壤图所附的土壤属性表和空间数据。对于土壤可蚀性和土壤结皮而言，由于我国土壤颗粒分级与美国制不同，RWEQ中的分级使用的是美国制，为此需要先对土壤颗粒含量进行粒径转换，且实测的土壤颗粒含量参数需符合RWEQ标准数据库中的物质含量范围表，当实测值不符合要求时，可以使用RWEQ内嵌的土壤质地资料的输入参数。表C.3RWEQ模型内嵌适用的土壤资料%土壤类型SaSiOMCaCO3砂土9340.31壤质砂土84100.52砂质壤土64260.53砂质粘壤土591313砂质粘土5271322)(021.0)(0066.011OMClSCF59泥土6881.53粉质壤土21671.53壤土41411.53粉质黏壤土105623粉质粘土6472.53粘壤土32342.53黏土202033（4）地表粗糙度因子土壤糙度因子K′取决于自由糙度RR和定向糙度OR，可采用一种滚轴式链条法来测定地表糙度。其基本原理是：两点间直线距离最短，当地表糙度增加时，其地表距离随之增加。于是当一个给定长度为L1的链条放于粗糙的地表时，其水平长度将缩小为L2，L1和L2的差值和地表粗糙程度密切相关，可用下式来计算地表糙度。式中，Crr为任意方向上的地表糙度。这里还涉及到一个垂直于垄向的糙度PR，当用链条法测出RR和PR时，OR就可通过下式来求得。OR=PR-RR土垄糙度因子可通过下式来进行计算，假设土垄呈等腰三角形则：式中：Nr为在长度L2范围内土垄的数量。根据实验结果，当不考虑风向的影响时，土壤糙度因子（K′）与土垄糙度因子（Kr）和随机糙度（RR）的回归方程为：100112LLCrrrrrrCK1K6eK124.04.28.1934.0′60也可采用风向与垄成任意角度时的糙度因子公式：式中，θ为风向与垄平行方向的夹角（角度）。（5）植被覆盖因子该因子用来确定生长植被、枯萎植被（农作物为倒伏残茬）以及农作物的直立残茬对土壤风蚀的影响。综合植被覆盖因子COG为枯萎植被、直立残茬和生长植被覆盖三因子的乘积。i)枯萎植被式中，SLRf为枯萎植被的土壤流失比率；SC为枯萎植被地表覆盖率。ii)直立残茬式中，SLRs为直立残茬土壤流失比率；SA为直立残茬当量面积，是1m2内直立秸秆的个数乘以秸秆直径的平均值（cm）再乘以秸秆高度（cm）。iii)生长植被式中，SLRc为植被覆盖土壤流失比率；cc是土表植被覆盖度。rrrrcCKKReK124.041.286.1934.0′342221074.31033.771.4100.1cR）（sc-0.0438feSLR）（0.6413SA-0.0344seSLR）（0.7366cc-5.614ceSLR61附录C.3固碳实物量核算参数NEP和NPP的转换系数α取值如表C.4所示：表C.4各地理区和各省不同生态系统NEP和NPP的转换系数地理区省份森林灌丛草地湿地华北平均0.3300.2970.1520.197北京0.3610.3460.2810.086天津0.3630.3600.0680.053河北0.3540.3400.2420.103山西0.3760.3350.2180.077内蒙0.3150.1120.1290.208东北平均0.2110.1460.1730.147辽宁0.1720.1690.1810.076吉林0.1390.0570.1250.075黑龙江0.2520.1660.1980.160华东平均0.1490.1160.1300.064上海0.1490.1160.0150.059江苏0.1400.0820.0940.069浙江0.1580.1510.1120.052安徽0.1630.1120.1290.045福建0.1750.1390.1380.080江西0.1120.0850.1180.065山东0.1880.1970.1460.094中南平均0.1200.0990.0980.045河南0.2080.2000.1730.062湖北0.1620.1300.1130.038湖南0.0970.0700.0840.084广东0.1360.0900.1060.031广西0.0850.0750.0630.020海南0.1690.1430.0320.025西南平均0.2090.1420.2260.189重庆0.1760.1830.1700.082四川0.2550.2030.2870.146贵州0.0490.0080.0600.010云南0.2260.1630.1360.080西藏0.2150.0870.2620.30262西北平均0.2040.1040.1660.208陕西0.2220.2390.2330.096甘肃0.2540.1960.1460.150青海0.2390.1910.2780.233宁夏0.2100.1350.1690.056新疆0.0350.1910.0150.148生物量-碳转换系数𝐶𝐶：森林和灌丛的转化系数为0.5，草地的转化系数为0.45。森林及灌丛的固碳速率FCSR由森林清查数据计算获得，如表C.5所示：表C.5各省森林（及灌丛）生态系统固碳速率区域森林(及灌丛)植被固碳速率（tC·ha-1·a-1）森林(及灌丛)土壤固碳速率（tC·ha-1·a-1）2000年2010年2015年2000年2010年2015年全国0.530.750.840.340.480.54北京0.911.000.840.590.650.55天津0.901.342.110.580.861.36河北0.171.120.770.110.720.50山西0.520.851.020.330.550.66内蒙古0.280.461.040.180.300.67辽宁0.240.921.260.150.590.82吉林0.360.260.970.230.160.63黑龙江0.260.750.630.170.480.40上海1.524.161.970.982.691.27江苏0.482.242.600.311.451.68浙江0.360.850.710.230.550.46安徽0.661.111.150.430.710.74福建0.870.801.750.560.521.13江西0.380.920.290.250.590.19山东0.352.831.050.221.830.68河南0.631.941.120.411.250.73湖北0.621.141.260.400.740.81湖南0.541.001.000.350.650.65广东0.290.190.790.190.130.51广西0.560.950.490.360.610.31海南0.260.031.000.170.020.6563重庆1.331.331.240.860.860.80四川0.510.650.630.330.420.41贵州1.161.281.170.750.820.76云南1.020.700.660.660.450.43西藏1.500.460.120.970.300.08陕西0.400.590.840.260.380.54甘肃0.790.890.900.510.570.58青海0.200.750.920.130.490.60宁夏1.041.041.310.670.670.84新疆0.991.001.360.640.640.88森林及灌丛土壤固碳系数β取值0.646。全国草地（除青藏高原外）土壤的固碳速率为0.02tC·ha-1·a-1，青藏高原区域为0.03tC·ha-1·a-1。草地土壤固碳速率GSR见表C.6。表C.6草地土壤固碳速率区域草地土壤固碳速率（tC·ha-1·a-1）区域草地土壤固碳速率（tC·ha-1·a-1）全国0.02河南0.02北京0.02湖北0.02天津0.02湖南0.02河北0.02广东0.02山西0.02广西0.02内蒙古0.06海南0.02辽宁0.02重庆0.02吉林0.02四川0.03黑龙江0.02贵州0.02上海0.02云南0.03江苏0.02西藏0.03浙江0.02陕西0.02安徽0.02甘肃0.02福建0.02青海0.03江西0.02宁夏0.02山东0.02新疆0.03水域湿地的固碳速率SCSRi取值见表C.7：64表C.7湿地固碳速率类型固碳速率（gC·m-2·a-1）湖泊湿地类型东部平原地区湖泊湿地56.67蒙新高原地区湖泊湿地30.26云贵高原地区湖泊湿地20.08青藏高原地区湖泊湿地12.57东北平原与山区湖泊湿地4.49沼泽湿地类型泥炭和苔藓泥炭沼泽24.80腐泥沼泽32.48内陆盐沼67.11沿海滩涂盐沼235.62红树林沼泽444.27无化学肥料和有机肥料施用的情况下，我国农田土壤有机碳的变化NSC取-0.06；各省土壤容重BD取值见表C.8；土壤厚度取20cm。表C.8土壤容重省份容重（g/cm3）省份容重（g/cm3）全国1.250河南1.312北京1.289湖北1.261天津1.309湖南1.241河北1.294广东1.225山西1.261广西1.245内蒙古1.271海南1.244辽宁1.288重庆1.233吉林1.188四川1.240黑龙江1.141贵州1.207上海1.194云南1.230江苏1.253西藏1.225浙江1.177陕西1.274安徽1.241甘肃1.261福建1.297青海1.236江西1.281宁夏1.248山东1.304新疆1.298作物j的草谷比SGRj取值如表C.9所示：65表C.9不同作物的草谷比作物草谷比SGRj作物草谷比SGRj水稻0.623油菜2小麦1.366向日葵2玉米2棉花8.1高粱1甘蔗0.1马铃薯0.5区域地下水径流模数以及HCO3-浓度见表C.10：表C.10区域地下水径流模数以及HCO3-浓度岩溶类型区HCO3-浓度（g/L）M-岩溶水径流模数（107L/(km2·a)）南方岩溶区0.23140.59北方岩溶区0.24520.86青藏高原区0.14519.95埋藏岩溶区0.2250.375附录C.4空气净化实物量核算参数表C.11环境空气污染物浓度限值污染物平均时间年平均浓度限制单位一级二级二氧化硫年平均2060ug/m3二氧化氮年平均4040颗粒物PM10年平均4070颗粒物PM2.5年平均1535注：环境空气功能区分为两类，一类区为自然保护区、风景名胜区和其他需要特殊保护的区域；二类区为居住区、商业交通居民混合区、文化区、工业区和农村地区。一类区适用一类浓度限值，二类区适用二类浓度限值。核算过程中，将核算区域大气污染物监测点位的算术平均值与所在功能区的空气浓度限值进行比较，来确定核算方法。66附录C.5水质净化实物量核算参数表C.12地表水污染物浓度限值污染物I类II类III类IV类V类化学需氧量1515203040氨氮0.150.511.52注：地表水水环境功能分为五类，I类适用于源头水、国家自然保护区；II类适用于集中式生活饮用水地表水源地一级保护区、珍稀水生生物栖息地、鱼虾类产卵场、仔稚幼鱼的索饵场等；III类适用于集中式生活饮用水地表水源地二级保护区、鱼虾类越冬场、洄游通道、水产养殖区等渔业水域及游泳区；IV类适用于一般工业用水区及人体非直接接触的娱乐用水区；V类适用于农业用水区及一般景观要求水域。核算过程中，将核算区域水质监测断面的污染物浓度算术平均值与所在功能区的污染物浓度限值进行比较，来确定核算方法。附录C.6生物多样性保护核算参数表C.13物种濒危指数体系濒危指数濒危等级4极危3濒危2易危1近危表C.14特有种指数体系特有种指数分布范围4仅限于范围不大的山峰或特殊的自然地理环境下分布3仅限于某些较大的自然地理环境下分布的类群，如仅仅分布于较大的海岛、高原、若干个山脉等2仅限于某个大陆分布的分布群1至少在2个大陆都有分布的分布群0世界广布的分类群67表C.15古树年龄指数体系指数等级古树年龄1100-299年2300-499年3≥500年表C.16Shannon-Wiener指数等级划分及其价值量等级Shannon-Wiener指数单价/（元/hm2·a）I指数≥666760II5≤指数≤653410III4≤指数≤540060IV3≤指数≤426700V2≤指数≤313350VI1≤指数≤26680