陆地生态系统生产总值（GEP）
核算技术指南
生态环境部环境规划院 中国科学院生态环境研究中心
二〇二〇年九月
目录
前言 ............................................................................................................. 1 1 总则 ......................................................................................................... 2
1.1适用范围 ........................................................................................................................................ 2 1.2规范性引用文件 ............................................................................................................................ 2 1.3术语和定义 .................................................................................................................................... 2
2 核算程序与内容 .................................................................................... 5
2.1 核算原则 ....................................................................................................................................... 5 2.2 核算程序 ....................................................................................................................................... 5 2.3 核算内容与指标 ........................................................................................................................... 8 2.4 生态系统服务实物量核算 ......................................................................................................... 10 2.5 生态系统服务价值量核算 ......................................................................................................... 11
3 物质产品实物量和价值量核算 .......................................................... 13
3.1 物质产品实物量核算 ................................................................................................................. 13 3.2 物质产品价值量核算 ................................................................................................................. 14
4水源涵养实物量和价值量核算 ........................................................... 16
4.1 水源涵养实物量核算 ................................................................................................................. 16 4.2 水源涵养价值量核算 ................................................................................................................. 17
5土壤保持实物量和价值量核算 ........................................................... 19
5.1 土壤保持实物量核算 ................................................................................................................. 19 5.2 土壤保持价值量核算 ................................................................................................................. 20
6防风固沙实物量和价值量核算 ........................................................... 22
6.1 防风固沙实物量核算 ................................................................................................................. 22 6.2 防风固沙价值量核算 ................................................................................................................. 23
7海岸带防护实物量和价值量核算 ....................................................... 24
7.1 海岸带防护实物量核算 ............................................................................................................. 24 7.2 海岸带防护价值量核算 ............................................................................................................. 24
8洪水调蓄实物量和价值量核算 ........................................................... 26
8.1 洪水调蓄实物量核算 ................................................................................................................. 26 8.2 洪水调蓄价值量核算 ................................................................................................................. 29
9碳固定实物量和价值量核算 ............................................................... 31
9.1 碳固定实物量核算 ..................................................................................................................... 31
9.2 碳固定价值量核算 ..................................................................................................................... 37
10氧气提供实物量和价值量核算 ......................................................... 38
10.1氧气提供实物量核算 ................................................................................................................ 38 10.2氧气提供价值量核算 ................................................................................................................ 38
11空气净化实物量和价值量核算 ......................................................... 40
11.1 空气净化实物量核算................................................................................................................ 40 11.2 空气净化价值量核算................................................................................................................ 41
12水质净化实物量和价值量核算 ......................................................... 43
12.1 水质净化实物量核算 ............................................................................................................... 43 12.2 水质净化价值量核算 ............................................................................................................... 44
13气候调节实物量和价值量核算 ......................................................... 46
13.1 气候调节实物量核算 ............................................................................................................... 46 13.2 气候调节价值量核算 ............................................................................................................... 48
14物种保育实物量和价值量核算 ......................................................... 49
14.1物种保育实物量核算 ................................................................................................................ 49 14.2物种保育价值量核算 ................................................................................................................ 50
15休闲旅游实物量和价值量核算 ......................................................... 51
15.1休闲旅游实物量核算 ................................................................................................................ 51 15.2 休闲旅游价值量核算 ............................................................................................................... 51
16景观价值实物量和价值量核算 ......................................................... 53
16.1景观价值实物量核算 ................................................................................................................ 53 16.2景观价值价值量核算 ................................................................................................................ 53
17附则 ...................................................................................................... 55
17.1 核算数据 ................................................................................................................................... 55 17.2 成果表达形式 ........................................................................................................................... 55 17.3 结果校核 ................................................................................................................................... 55 17.4 信息管理 ................................................................................................................................... 56
附录A（资料性附录）生态系统生产总值核算数据来源表 .............. 57 附录B（资料性附录）陆地生态系统生产总值核算报告编写大纲 .. 60 附录C（资料性附录）生态系统生产总值实物量核算参数参考表 .. 61
附录C.1 土壤保持实物量核算参数 ................................................................................................ 61 附录C.2 防风固沙实物量核算参数 ................................................................................................ 附录C.3 固碳实物量核算参数 ........................................................................................................ 69
附录C.4 空气净化实物量核算参数 ................................................................................................ 74 附录C.5水质净化实物量核算参数 ................................................................................................ 74
前言
为全面贯彻落实*生态文明思想和“绿水青山就是金山银山”理念，以及《关于全面深化改革若干重大问题的决定》、《关于加快推进生态文明建设的意见》和《生态文明改革总体方案》中建立生态效益评估机制、促进人与自然和谐的部署，保障国家和区域生态安全，指导和规范陆地生态系统生产总值核算工作，提高陆地生态系统生产总值实物量与价值量核算的科学性、规范性和可操作性，制定本指南，供地方参考使用。
本指南可为生态效益纳入经济社会发展评价体系、完善发展成果考核评价体系提供重要支撑，为建立生态产品实现机制、区域生态补偿、自然资源资产审计、自然资产负债表编制等制度的实施提供科学依据。
本指南规定了陆地生态系统生产总值实物量与价值量核算的技术流程、指标体系与核算方法等内容。
本指南附录均为资料性附录。
本指南由生态环境部环境规划院与中国科学院生态环境研究中心共同编制。
本指南主要起草人：欧阳志云、王金南、肖燚、於方、郑华、马国霞，宋昌素、周夏飞、邹梓颖，杨威杉，逯非，周颖，韩宝龙，吴春生，彭菲，宁佳。
1
1 总则
1.1适用范围
本指南适用于省、市、县等行政区域开展生态系统生产总值核算工作。自然地理区域可参考本指南开展生态系统生产总值核算。 1.2规范性引用文件
本指南内容引用了下列文件中的条款。凡是不注日期的引用文件，其有效版本适用于本指南。
《环境经济核算体系——实验生态账户》(2012) GB 3095 环境空气质量标准 GB 3838 地表水环境质量标准
GB/T 12343 国家基本比例尺地图编绘规范 GB/T 13923 基础地理信息要素分类与代码 CH/T 9005 基础地理信息数据库基本规定 GB-T 38582 森林生态系统服务功能评估规范 LY∕T 29 湿地生态系统服务评估规范 1.3术语和定义
以下术语和定义适用于本指南。
（1）陆地生态系统（Terrestrial ecosystem）
陆地生态系统是指地球表面陆地生物及其环境通过能流、物流、信息流形成的功能整体。陆地生态系统包括森林生态系统、草地生态
2
系统、湿地生态系统、荒漠生态系统、农田生态系统、城市生态系统等类型。
（2）湿地生态系统（Wetland ecosystem）
湿地生态系统由陆地和水域相互作用而形成的兼顾水域和陆地生态系统特征的自然综合系统，包括陆地所有淡水生态系统、陆地和海洋过渡地带的滨海湿地生态系统和海洋边缘部分咸水、半咸水水域。
（3）草地生态系统（Grassland ecosystem）
草地生态系统是以饲用植物和食草动物为主体的生物群落与其生存环境共同构成的开放生态系统，包括人工草地生态系统和天然草地生态系统两大类。
（4）森林生态系统（Forest ecosystem）
森林生态系统是以乔木为主体的生物群落与其非生物环境相互作用，并产生能量转换和物质循环的综合系统，包括天然林生态系统和人工林生态系统。
（5）荒漠生态系统（Desert ecosystem）
荒漠生态系统指由超强耐旱生物及其干旱环境所组成的生态系统。
（6）农田生态系统（Farmland ecosystem）
农田生态系统是指人类在以作物为中心的农田中，利用生物和非生物环境之间以及生物种群之间的相互关系，通过合理的生态结构和高效生态机能，进行能量转化和物质循环，并按人类社会需要进行物质生产的综合体。
（7）城市生态系统（Urban ecosystem）
城市生态系统是城市居民与其环境相互作用而形成的统一整体，
3
由自然环境、社会经济和文化科学技术共同组成。
（8）生态系统服务（Ecosystem services）
人类从生态系统中得到的惠益，包括物质产品供给、调节服务、文化服务以及支持服务。
（9）生态系统生产总值（Gross Ecosystem Product, GEP） 生态系统生产总值也可简称为生态产品总值，是指生态系统为人类福祉和经济社会可持续发展提供的各种最终产品与服务价值的总和，主要包括生态系统提供的物质产品、调节服务和文化服务。
（10）物质产品（Material services）
人类从生态系统获取的可在市场交换的各种物质产品，如食物、纤维、木材、药物、装饰材料与其他物质材料。
（11）调节服务（Regulating services）
生态系统提供改善人类生存与生活环境的惠益，如调节气候、涵养水源、保持土壤、调蓄洪水、降解污染物、固定二氧化碳、氧气提供等。
（12）文化服务（Cultural services）
人类通过精神感受、知识获取、休闲娱乐和美学体验从生态系统获得的非物质惠益。
（13）实物量（Biophysical value）
生态系统产品与服务的实物量，如粮食产量、木材生产量、水产品捕捞量、洪水调蓄量、土壤保持量、碳固定量与景点旅游人数等。
（14）价值量（Monetary Value） 生态系统产品与服务的货币价值。
4
2 核算程序与内容
2.1 核算原则
陆地生态系统生产总值核算应遵循以下原则：
科学性原则。根据生态系统与人类福祉的关系构建核算框架与指标体系。以生态系统结构、格局和过程与生态系统服务关系为基础，构建核算方法。
实用性原则。根据核算区域生态系统的特点、核算目的和数据可获得性，确定适合本区域的核算内容、指标和方法。
系统性原则。生态系统生产总值核算既要考虑为当地人提供的惠益，也需要考虑为其它地区提供的惠益。
开放性原则。应根据生态系统核算研究的最新成果，发展和完善核算的指标与方法。 2.2 核算程序
陆地生态系统生产总值核算的主要工作程序包括：根据核算目的，确定生态系统生产总值核算区域范围，明确生态系统类型以及生态系统服务的清单，确定核算模型方法与适用技术参数，开展生态系统生产总值实物量与价值量核算（见图2-1）。
5
图2-1 生态系统生产总值核算工作程序
6
（1）确定核算的区域范围：根据核算目的，确定生态系统生产总值核算的空间范围。核算区域可以是行政区域，如村、乡、县、市或省，也可以是功能相对完整的生态地理单元，如一片森林、一个湖泊、一片沼泽或不同尺度的流域，以及由不同生态系统类型组合而成的地域单元。
（2）明确生态系统类型与分布：调查分析核算区域内的森林、草地、湿地、荒漠、农田、城镇等生态系统类型、面积与分布，绘制生态系统空间分布图。
（3）编制生态系统产品与服务清单：根据生态系统类型及核算的用途，如生态补偿、离任审计、生态产品交易，调查核算范围的生态系统服务的种类，编制生态系统产品和服务清单。当核算目标为评估生态保护成效时，可以只核算生态调节服务和生态文化服务价值。
（4）收集资料与补充调查：收集开展生态系统生产总值核算所需要的相关文献资料、监测与统计等信息数据以及基础地理与地形图件，开展必要的实地观测调查，进行数据预处理以及参数本地化。
（5）开展生态系统产品与服务实物量核算：选择科学合理、符合核算区域特点的实物量核算方法与技术参数，根据确定的核算基准时间，核算各类生态系统产品与服务的实物量。
（6）开展生态系统产品与服务价值量核算：根据生态系统产品与服务实物量，运用市场价值法、替代成本法等方法，核算生态系统产品与服务的货币价值；无法获得核算年份价格数据时，利用已有年份数据，按照价格指数进行折算。
（7）核算生态系统生产总值：将核算区域范围的生态产品与服务价值加总，得到生态系统生产总值。
7
2.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 调节服务 洪水调蓄 碳固定 生态系统吸收二氧化碳合成有机物质，将碳固定在植物和土壤中，降低大气中二氧化碳浓度的功能。 生态系统通过光合作用释放出氧气，维持大气氧气浓度稳定的功能。 生态系统吸收、阻滤大气中的污染物，如SO2、NOx、颗粒物等，降低空气污染浓度，改善空气环境的功能。 生态系统通过物理和生化过程对水体污染物吸附、降解以及生物吸收等，降低水体污染物浓度、净化水环境的功能。 生态系统通过植被蒸腾作用和水面蒸发过程吸收能量、降低气温、提高湿度的功能。 生态系统为珍稀濒危物种提供生存与繁衍场所的作用和价值。 氧气提供 14 15 16 空气净化 水质净化 气候调节 17 18 文化服务 19 物种保育 休闲旅游 人类通过精神感受、知识获取、休闲娱乐和美学体验、康养等旅游休闲方式，从生态系统获得的非物质惠益。 生态系统为人类提供美学体验、精神愉悦，从而提高周边土地、房产价值的功能。 景观价值
9
表2-2 生态系统生产总值实物量及价值量核算指标体系 服务类别 核算科目 农业产品 林业产品 物质产品 畜牧业产品 渔业产品 生态能源 其他产品 水源涵养 土壤保持 防风固沙 海岸带防护 洪水调蓄 空气净化 调节服务 水质净化 碳固定 氧气提供 气候调节 物种保育 文化服务 休闲旅游 景观价值 实物量指标 农业产品产量 林业产品产量 畜牧业产品产量 渔业产品产量 生态能源总量 装饰观赏资源总量等 水源涵养量 土壤保持量 固沙量 海岸带防护面积 洪水调蓄量 净化二氧化硫量 净化氮氧化物量 净化颗粒物量 净化COD量 净化总氮量 净化总磷量 固定二氧化碳量 氧气提供量 植被蒸腾消耗能量 水面蒸发消耗能量 珍稀濒危物种数量 保护区面积 游客总人数 受益土地与房产面积 价值量指标 农业产品产值 林业产品产值 畜牧业产品产值 渔业产品产值 生态能源产值 装饰观赏资源产值等 水源涵养价值 减少泥沙淤积价值 减少面源污染价值 草地恢复成本 海岸带防护价值 调蓄洪水价值 净化二氧化硫价值 净化氮氧化物价值 净化颗粒物价值 净化COD价值 净化总氮价值 净化总磷价值 碳固定价值 氧气提供价值 植被蒸腾调节温湿度价值 水面蒸发调节温湿度价值 珍稀濒危物种保育价值 保护区保育价值 游憩康养价值 土地、房产升值 2.4 生态系统服务实物量核算
生态系统生产总值实物量核算包括三大类，即物质产品实物量核算、调节服务实物量核算、文化服务实物量核算。生态系统生产总值实物量核算的核算项目、实物量指标和核算方法见表2-3。
10
表2-3 生态系统生产总值实物量核算方法
服务类别 核算项目 农业产品 林业产品 物质产品 畜牧业产品 渔业产品 生态能源 其他 水源涵养 土壤保持 防风固沙 海岸带防护 洪水调蓄 调节服务 空气净化 实物量指标 农业产品产量 林业产品产量 畜牧业产品产量 渔业产品产量 生态能源总量 装饰观赏资源产量等 水源涵养量 土壤保持量 固沙量 海岸带防护面积 湖泊：可调蓄水量 水库：防洪库容 沼泽：滞水量 净化二氧化硫量 净化氮氧化物量 净化颗粒物量 净化COD量 水质净化 碳固定 氧气提供 气候调节 物种保育 文化服务 休闲旅游 景观价值 净化总氮量 净化总磷量 固定二氧化碳量 氧气提供量 植被蒸腾消耗能量 水面蒸发消耗能量 珍稀濒危物种数量 游客总人数 受益土地面积或公众 固碳机理模型 释氧机理模型 蒸散模型 统计调查 统计调查 污染物净化模型 污染物净化模型 水量储存模型 水量平衡法，水量供给法 修正通用土壤流失方程（RUSLE） 修正风力侵蚀模型（REWQ） 统计调查 统计调查 核算方法 2.5 生态系统服务价值量核算
在生态系统生产总值实物量核算的基础上，确定各类生态系统服务的价格，核算生态服务价值。具体而言，生态系统生产总值价值量核算中，物质产品价值主要用市场价值法核算，调节服务价值主要用替代成本法进行核算，文化服务价值使用旅行费用法。生态系统生产
11
总值价值量核算的核算项目、价值量指标和核算方法，见表2-4。
表2-4 生态系统生产总值价值量核算方法 服务类别 核算项目 农业产品 林业产品 物质产品 畜牧业产品 渔业产品 生态能源 其他 水源涵养 土壤保持 防风固沙 海岸带防护 洪水调蓄 空气净化 调节服务 水质净化 碳固定 氧气提供 气候调节 物种保育 文化服务
12
价值量指标 农业产品产值 林业产品产值 畜牧业产品产值 渔业产品产值 生态能源产值 装饰观赏资源产值等 水源涵养价值 减少泥沙淤积价值 减少面源污染价值 固沙价值 由于防护减少的损失价值 调蓄洪水价值 净化二氧化硫价值 净化氮氧化物价值 净化颗粒物价值 净化总氮价值 净化总磷价值 净化COD价值 固定二氧化碳价值 氧气提供价值 植被蒸腾调节温湿度价值 水面蒸发调节温湿度价值 物种保育价值 休闲旅游价值 景观价值 核算方法 市场价值法 替代成本法 替代成本法 替代成本法 恢复成本法 替代成本法 影子工程法 替代成本法 替代成本法 替代成本法 替代成本法 替代成本法 替代成本法 替代成本法 替代成本法 替代成本法 保育价值法 旅行费用法 享乐价格法 休闲旅游 景观价值
3 物质产品实物量和价值量核算
3.1 物质产品实物量核算 3.1.1 指标内涵
物质产品是指人类从生态系统获取的能够在市场交易的产品，满足人类生活、生产与发展的物质需求，包括农业、林业、畜牧业、渔业等产品与生态能源。具体指标见表3-1。
表3-1 生态系统的物质产品指标
类别 项目 内容 谷物 粮食作物 豆类 薯类 油料 棉花 麻类 农业产品 糖类 烟叶 药材 茶叶 蔬菜 瓜类 水果 木材 油料 棉花 麻类 糖类 烟叶 药材 茶叶 蔬菜 瓜类 水果 木材 橡胶 林业产品 松脂 其他林产品 生漆 油桐籽 油茶籽 畜禽 奶类 畜牧业产品 禽蛋 动物皮毛 其他畜产品 畜禽 奶类 禽蛋 羊毛 蜂蜜 蚕茧 指标 水稻、小麦、玉米、谷子、高粱、其他谷物等 大豆、绿豆、红小豆 马铃薯 花生、油菜籽、向日葵籽、芝麻、胡麻籽 棉花 黄红麻、亚麻、大麻、苎麻 甜菜、甘蔗 烟叶 药材 红毛茶、绿毛茶 蔬菜（含菜用瓜） 西瓜、甜瓜、草莓 香蕉、苹果、梨、葡萄、柑桔、红枣、柿子、菠萝、其他园林水果等 木材 橡胶 松脂 生漆 油桐籽 油茶籽 牛、羊、猪、家禽等 牛奶 禽蛋 细绵羊毛、半细绵羊毛、羊绒、山羊毛 蜂蜜 柞蚕茧、桑蚕茧 13
淡水产品 海水产品 淡水产品 海水产品 水能 生态能源 生态能源 薪柴 秸秆 沼气 其它 装饰观赏资源等 装饰观赏资源等 鱼类、贝类、虾蟹类、其他 鱼类、贝类、虾蟹类、藻类、其他 水能发电量 薪柴量 固化产量 沼气量 花卉、苗木等 渔业产品 3.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 定价思路
生态系统物质产品价值是指生态系统通过初级生产、次级生产为人类提供农产品、林业产品、畜牧业产品、渔业产品、生态能源等的经济价值。由于生态系统物质产品能够在市场上进行交易，存在相应的市场价格，可以运用市场价值法对生态系统的物质产品服务进行价值核算。
14
3.2.2 价值量核算模型
𝑛
𝑉𝑚 ∑𝑖<1𝐸𝑖 𝑃𝑖 （2）
式中，Vm为生态系统物质产品价值（元/a）；Ei为第i类生态系统产品产量（根据产品的计量单位确定，如kg/a），Pi为第i类生态系统产品的价格（根据产品的计量单位确定，如元/kg）。 3.2.3 定价参数与数据来源
产品产量由实物量核算得到，产品价格从林业、农业、渔业及统计部门获得或根据市场定价获得。
15
4水源涵养实物量和价值量核算
4.1 水源涵养实物量核算 4.1.1指标内涵
水源涵养服务是生态系统拦截滞蓄降水，增强土壤下渗、蓄积，涵养土壤水分、调节暴雨径流和补充地下水，增加可利用水资源的功能。水源涵养量大的地区不仅满足核算区内生产生活的水源需求，还持续地向区域外提供水资源。选用水源涵养量，作为生态系统水源涵养实物量的评价指标。 4.1.2核算方法
水源涵养量核算主要有两种方法，推荐优先采用水量平衡法，技术参数缺失的情况下，可以选择水量供给法。
4.1.2.1水量平衡法
通过水量平衡方程计算。水量平衡方程是指在一定的时空内，水分在生态系统中保持质量守恒，即生态系统水源涵养量是降水输入与暴雨径流和生态系统自身水分消耗量的差值。
;3
𝑄𝑤𝑟=∑𝑛𝑖<1𝐴𝑖 (𝑃𝑖−𝑅𝑖−𝐸𝑇𝑖) 10 （3）
式中，Qwr为水源涵养量（m3/a）；Pi为产流降雨量（mm/a），Ri
为地表径流量（mm/a），ETi为蒸散发量（mm/a）；Ai为i类生态系统的面积（m2）；i为生态系统类型，n为生态系统类型总数。
16
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 定价思路
水源涵养价值主要表现在蓄水保水的经济价值。可运用影子工程法，即模拟建设蓄水量与生态系统水源涵养量相当的水利设施，以建设该水利设施所需要的成本核算水源涵养价值。 4.2.2 价值量核算模型
𝑉𝑤𝑟 𝑄𝑤𝑟 𝐶𝑤 （5）
式中，Vwr为水源涵养价值（元/a）；Qwr为核算区内总的水源涵养量（m3/a），Cwe为水资源交易市场价格，当交易市场未建立时，以
17
水库建设的工程及维护成本或水资源影子价格（元/m3）替代，水库单位库容的工程造价及维护成本（元/m3）。 4.2.3 定价参数与数据来源
生态系统水源涵养量由实物量核算得到。水库单位库容的工程造价及维护成本等数据来自、水利等部门发布的工程预算依据，或公开发表的参考文献，并根据价格指数折算得到核算年份的价格。
18
5土壤保持实物量和价值量核算
5.1 土壤保持实物量核算 5.1.1 指标内涵
土壤保持功能是生态系统（如森林、草地等）通过林冠层、枯落物、根系等各个层次保护土壤、消减降雨侵蚀力，增加土壤抗蚀性，减少土壤流失，保持土壤的功能。
选用土壤保持量，即生态系统减少的土壤侵蚀量（用潜在土壤侵蚀量与实际土壤侵蚀量的差值测度）作为生态系统水土保持功能的评价指标。其中，实际土壤侵蚀是指当前地表植被覆盖情形下的土壤侵蚀量，潜在土壤侵蚀是指没有地表植被覆盖情形下可能发生的土壤侵蚀量。 5.1.2 核算方法
土壤保持量核算主要基于修正的通用水土流失方程（RUSLE）计算。
𝑄𝑠𝑟=𝑅 𝐾 𝐿 𝑆 (1−𝐶 𝑃) （6）
式中，Qsr为土壤保持量（t/ a）；R为降雨侵蚀力因子，用多年平均年降雨侵蚀力指数表示；K为土壤可蚀性因子，通常用标准样方上单位降雨侵蚀力所引起的土壤流失量来表示；L为坡长因子（无量纲），S为坡度因子（无量纲），C为植被覆盖和管理因子（无量纲），P为水土保持措施因子（无量纲）。
19
5.1.3 核算参数及数据来源
降雨侵蚀力因子R、土壤可蚀性因子K、坡长因子L、坡度因子S的算法以及覆盖和管理因子C以及水土保持措施因子P参考附录C.1或查阅相关文献。 5.2 土壤保持价值量核算 5.2.1 定价思路
生态系统土壤保持价值主要包括减少面源污染和减少泥沙淤积两个方面的价值。
生态系统通过保持土壤，减少氮、磷等土壤营养物质进入下游水体（包括河流、湖泊、水库和海湾等），可降低下游水体的面源污染。根据土壤保持量和土壤中氮、磷的含量，运用替代成本法（即污染物处理的成本）核算减少面源污染的价值。
生态系统通过保持土壤，减少水库、河流、湖泊的泥沙淤积，有利于降低干旱、洪涝灾害发生的风险。根据土壤保持量和淤积量，运用替代成本法（即水库清淤工程的费用）核算减少泥沙淤积价值。 5.2.2 价值量核算模型
𝑉𝑠 +𝑉 𝑝 （7） 𝑄𝑠𝑟 ) （8） 𝑉 𝑝 =∑ 𝑃𝑖 （9）
式中，Vsr为生态系统土壤保持价值（元/a）；Vsd为减少泥沙淤积价值（元/a）；Vdpd为减少面源污染价值（元/a）；Qsr为土壤保持量（t/a）；c为单位水库清淤工程费用（元/m3）；ρ为土壤容重（t/m3）； 为泥沙
20
淤积系数；i为土壤中氮、磷等营养物质数量，i = 1，2，…，n； 为土壤中氮、磷等营养物质的纯含量（%）；𝑃𝑖为处理成本。 5.2.3 定价参数与数据来源
土壤保持量由实物量核算得到。土壤容重、氮、磷、钾含量、单位水库清淤工程费、单位污染物处理成本、肥料价格等数据来源于当地土壤调查、文献、专项调查以及等物价部门。
21
6防风固沙实物量和价值量核算
6.1 防风固沙实物量核算 6.1.1 指标内涵
防风固沙功能是指生态系统减少因大风导致的土壤流失和风沙危害的功能。在风蚀过程中，植被减少土壤裸露，对土壤形成保护，减少风蚀输沙量，还可以通过根系固定表层土壤，改良土壤结构，提高土壤抗风蚀的能力，植被还可以通过增加地表粗糙度、阻截等方式降低风速、降低大风风力侵蚀和风沙危害。
选用防风固沙量，即通过生态系统减少的风蚀量（潜在风蚀量与实际风蚀量的差值），作为生态系统防风固沙功能的评价指标。 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。
22
6.2 防风固沙价值量核算 6.2.1 定价思路
根据防风固沙量和土壤沙化盖沙厚度，核算出减少的沙化土地面积；运用恢复成本法，根据单位面积沙化土地治理费用或单位植被恢复成本核算生态系统防风固沙功能的价值。 6.2.2 价值量核算模型
𝑉𝑠𝑓=
𝑐 （11）
式中，Vsf为防风固沙价值（元/a）；Qsf为防风固沙量（t/a）；ρ为土壤容重（t/m3）；h为土壤沙化覆沙厚度（m）；c为单位治沙工程的成本（元/m3）或单位植被恢复成本（元/m2）。 6.2.3 定价参数与数据来源
防风固沙量由实物量核算得到，土壤容重来自土壤调查或文献资料，单位治沙工程成本或单位植被恢复成本来自相关文件规定。
23
7海岸带防护实物量和价值量核算
7.1 海岸带防护实物量核算 7.1.1指标内涵
海岸带防护是指滨海盐沼、红树林、珊瑚礁等生态系统减低海浪，避免或减小海堤或海岸侵蚀的功能。 7.1.2核算方法
运用滨海盐沼、红树林、珊瑚礁等生态系统防护或替代海堤等防护工程的长度核算。
𝐷𝑐𝑙=∑𝑛1𝐷𝑐𝑙𝑖 (12)
式中，𝐷𝑐𝑙为生态系统防护的海岸带总长度（km/a）；𝐷𝑐𝑙𝑖为i类生态系统防护的海岸带长度（km/a）；i为研究区生态系统类型，i = 1，2，…，n，无量纲；n为研究区生态系统类型数量，无量纲。 7.1.3核算参数及数据来源
海岸带长度从自然资源部门获得，或利用遥感数据分析结合实地调查获得。
7.2 海岸带防护价值量核算 7.2.1定价思路
运用替代成本法（即海浪防护工程建设成本）核算滨海盐沼、红树林、珊瑚礁等生态系统防风护堤的价值。
24
7.2.2价值核算模型
𝑉𝑐𝑙=∑𝑛1𝐷𝑐𝑙𝑖 𝐶𝑐𝑙𝑖 （13）
式中，Vcl为海岸带防护价值（元/a）；Dcli为i类生态系统防护的海岸带长度（km）；𝐶𝑐𝑙𝑖为i类生态系统海浪防护工程单位长度建设维护成本（元/km）。 7.2.3核算参数及数据来源
生态系统防护的海岸带长度由实物量核算得到。防护工程的单位长度建设成本从、海洋等部门获得。
25
8洪水调蓄实物量和价值量核算
8.1 洪水调蓄实物量核算 8.1.1 指标内涵
洪水调蓄功能是指自然生态系统所特有的生态结构能够吸纳大量的降水和过境水，蓄积洪峰水量，削减并滞后洪峰，以缓解汛期洪峰造成的威胁和损失的功能。原则上只核算年降雨量大于400mm地区的洪水调蓄价值。
选用植被调蓄水量和洪水期滞水量（库塘、湖泊、沼泽）表征生态系统的洪水调蓄能力，即调节洪水的潜在能力。 8.1.2 核算方法
洪水调蓄能力：
𝐶𝑓𝑚=𝐶𝑣𝑐+𝐶𝑟𝑐+𝐶𝑙𝑐+𝐶𝑚𝑐 （14）
式中，Cfm为洪水调蓄量（m3/a）；𝐶𝑣𝑐为植被洪水调蓄量（m3/a）；Crc为库塘洪水调蓄量（m3/a），Clc为湖区洪水调蓄量（m3/a），Cmc为沼泽洪水调蓄量（m3/a）。
植被洪水调蓄量：
Cvc=∑(Ph-Rfi)×Siv×10-3 （15）
i=1n
式中， v 为植被洪水调蓄量（m3/a）；Ph为大暴雨产流降雨量（mm）；Rf 为第i种生态系统产生的地表径流量（mm）；S v为第i种自然植被生态系统的面积（km2）；𝑖为自然植被生态系统类型，i=1, 2, …, n；n为自然植被生态系统类型数量，无量纲。
26
库塘洪水调蓄量－方法1：
库塘湿地调蓄洪水的能力，通过库塘湿地的实际洪水调蓄库容来计算，库塘湿地的实际洪水调蓄库容按其总库容的35%进行计算。
𝐶𝑟𝑐=0.35 𝐶𝑡 （16）
式中，𝐶𝑟𝑐为水库防洪库容（m3/a）；Ct为水库总库容（m3）。 库塘洪水调蓄量－方法2：
在实测数据可以支撑的条件下，库塘湿地调蓄洪水的能力，可通过洪水期水库的进出水总量进行计算。
𝐶𝑟𝑐=𝐶𝑖−𝐶𝑜 （17）
式中，𝐶𝑟𝑐为水库防洪库容（m3/a）；𝐶𝑖为洪水期水库进水总量（m3/a），𝐶𝑜为洪水期水库出水总量（m3/a）。
湖泊洪水调蓄量－方法1：
湖泊洪水调蓄量的大小与湖泊面积、水深、换水周期，以及流域内的水文过程密切相关。可根据湖泊水文学过程，通过汛期湖泊入湖、出湖流量随时间的变化计算湖泊在某一段时间内洪水调蓄量：
𝐶𝑙𝑐=
𝑡2
(𝑄𝐼∫𝑡1
−𝑄𝑂)𝑑𝑡(𝑄𝐼>𝑄𝑂) （18）
式中，Clc为湖泊t1-t2时间段内洪水调蓄量（m3）；QI为入湖流量（m3/s），QO为出湖流量（m3/s）。
湖泊洪水调蓄量－方法2：
考虑区域气候条件的差异，根据《中国湖泊志》，将全国湖泊划分为东部平原、蒙新高原、云贵高原、青藏高原、东北平原与山区5个湖泊区，分区构建湖泊换水次数与补给系数的模型，通过补给系数
27
核算湖泊换水次数；按不同湖泊区，基于湖面面积与湖泊换水次数建立湖泊水量调节量评价模型：
东部平原区:
𝐶𝑙𝑐=𝑒4.924 𝐴1.128 3.19 （19）
蒙新高原区:
𝐶𝑙𝑐=𝑒5.653 𝐴0.680 0.26 （20）
云贵高原区:
𝐶𝑙𝑐=𝑒4.904 𝐴0.927 0.36 （21）
青藏高原区:
𝐶𝑙𝑐=𝑒6.636 𝐴0.678 0.14 （22）
东北平原与山区:
𝐶𝑙𝑐=𝑒5.808 𝐴0.866 0.98 （23）
沼泽洪水调蓄量：
沼泽土壤具有特殊的水文物理性质，草根层和泥炭层孔隙度达72%-93%，饱和持水量达830%-1030%，在汛期可储存大量水分。沼泽湿地像一个巨大的天然蓄水库消纳洪水；同时沼泽湿地植被可拦截径流减缓洪水流速，削减和滞后洪峰。基于沼泽土壤蓄水量和地表滞水量模型核算沼泽湿地洪水调蓄量能力：
𝐶𝑚𝑐=𝐶𝑠𝑤𝑠+𝐶𝑠𝑟 （24）
式中，Cmc为沼泽洪水调蓄量（亿m3/a）；Csws为沼泽土壤蓄水量（亿m3/a），Csr为沼泽地表滞水量（亿m3/a）。
沼泽土壤蓄水量的计算公式如下：
28
𝐶𝑠𝑤𝑠=𝑆 (𝐹−𝐸) 10;2 𝑤 （25） 式中，Csws为沼泽土壤蓄水量（亿m3/a）；S为沼泽总面积（km2），h为沼泽湿地土壤蓄水深度（m），ρ为沼泽湿地土壤容重（g/cm3），ρw为水的密度（g/cm3），F为沼泽湿地土壤饱和含水率（无量纲），E为沼泽湿地洪水淹没前的自然含水率（无量纲）。
沼泽地表滞水量的计算公式如下：
𝐶𝑠𝑟=𝑆 𝐻 10;2 （26）
式中，Csr为沼泽地表滞水量（亿m3/a）；S为沼泽湿地总面积（km2），H为沼泽湿地地表滞水高度（m）。 8.1.3 核算参数及数据来源
湖泊进出水量通过水利部门统计资料获得。各类植被面积、湖面面积、沼泽面积来源自然资源、水利或其他相关统计部门。出入湖流量来源于水利、水文监测站点的实际观测数据。暴雨降雨量数据来源于气象部门。暴雨径流量参考相关研究文献或进行实测。 8.2 洪水调蓄价值量核算 8.2.1 定价思路
运用替代成本法（即水库的建设成本）核算自然生态系统的洪水调蓄价值。
8.2.2 价值量核算模型
𝑉𝑓𝑚=𝐶𝑓𝑚 𝐶𝑤 （27）
29
式中，Vfm为生态系统洪水调蓄价值（元/a）；Cfm为生态系统洪水调蓄量（m3/a）；𝐶𝑤 为水库单位库容的工程造价及维护成本（元/ m3）。 8.2.3 核算参数及数据来源
生态系统洪水调蓄量由实物量核算得到。水库单位库容的工程造价及维护成本等数据来自、水利等部门发布的工程预算依据，或公开发表的参考文献。
30
9碳固定实物量和价值量核算
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），C转化为CO2的系数。
31
𝑀𝐶𝑂2𝑀𝐶
=44/12，为
固碳速率法-森林及灌丛固碳量：
𝐹𝐶𝑆=𝐹𝐶𝑆𝑅 𝑆𝐹 （1+𝛽） （30）
式中，FCSR为森林及灌丛的固碳速率（tC·ha-1·a-1），SF为森林及灌丛面积（ha），β为森林及灌丛土壤固碳系数。
固碳速率法-草地固碳量：
由于草地植被每年都会枯落，其固定的碳又返还回大气或者进入土壤中，故不考虑草地植被的固碳量，只考虑草地的土壤固碳量。
𝐺𝑆𝐶𝑆=𝐺𝑆𝑅 𝑆𝐺 （31）
式中，GSR为草地土壤的固碳速率（tC·ha-1·a-1），SG为草地面积（ha）。
固碳速率法-湿地固碳量：
-2
S=∑ S SR×𝑆𝑊×10 （32） 𝑖𝑖 =1
式中，SCSRi为第i类水域湿地的固碳速率（g C·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为秸秆全部还田的农田
32
土壤固碳速率（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）
33
南方农区 SCSRS = 43.548×S + 375.1 （43）
式中，S为单位耕地面积秸秆还田量（t·ha-1·a-1），计算方法如下式所示：
S=∑𝑛j 𝑆𝐺𝑅𝑗 𝑆𝑃 （44） 𝑗<1𝐶𝑌
式中，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：
34
𝑁𝐸𝑃=𝛼∗𝑁𝑃𝑃∗𝑀 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-需要从大气吸收
35
0.5mol CO2，将碳以HCO3-形式存储在水体中。我国广大的岩溶面积以及岩溶形成机制对吸收大气CO2具有一定作用。选用固定CO2量作为岩溶固碳功能的评价指标。
岩溶固碳仅在含有岩溶地貌的区域选择核算。 9.1.2.2核算方法
岩溶地区固碳量可以用以下公式核算：
𝑄𝑘 𝑂2= 𝑆 𝑀 𝐶
21
−][HCO3
𝑀CO2𝑀HCO−
3
（49）
式中，𝑄𝑘 𝑂2为岩溶作用总固碳量（tCO2/a）；S为岩溶地区面积
-（km2），M为区域地下水径流模数（107L/(km2·a)），𝐶[HCO−为[HCO]3]3
浓度（g/L）；MCO2为CO2的摩尔质量44，MHCO3-为HCO3-的摩尔质量61。
根据相关数据，式（49）可简化为式（50）：
𝑄𝑘 𝑂2=0.3607 𝑆 𝑀 𝐶[HCO− （50） 3]
系数0.3607由CO2摩尔质量（MCO2）44，HCO3-摩尔质量（MHCO3-）61，每形成1摩尔质量HCO3-需要从大气吸收0.5摩尔质量CO2这三个系数确定。
9.1.2.3核算参数及数据来源
区域岩溶面积数据可从地方自然资源、水文和地质部门获得，区域地下水径流模数以及HCO3-浓度根据地方水文部门数据确定或实测获得，附录C.3供参考。
36
9.2 碳固定价值量核算 9.2.1定价思路
生态系统固碳价值可以采用替代成本法（造林成本法、工业减排成本）与市场价值法（碳交易价格）核算生态系统固碳的价值。 9.2.2价值量核算模型
f 𝑄 𝑂2 （51）
式中，VCf为生态系统固碳价值（元/a）；𝑄 𝑂2为生态系统固碳总量（tCO2/a）；CC为碳价格（元/t）。 9.2.3定价参数与数据来源
固碳量由实物量核算得到。单位造林固碳成本、工业碳减排成本、碳交易市场价格参考相关文献，建议采用碳市场交易价格。
37
10氧气提供实物量和价值量核算
10.1氧气提供实物量核算 10.1.1指标内涵
生态系统的释氧功能指植物在光合作用过程中，释放出氧气的功能。这种功能对于维护大气中氧气的稳定，改善人居环境具有重要意义。选用释氧量作为生态系统释氧功能的评价指标。 10.1.2 核算方法
根据光合作用化学方程式可知，植物每生产吸收1molCO2，就会释放1mol氧气。NEP可由净初级生产力（NPP）减去异氧呼吸消耗得到，或根据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定价思路
采用市场价值法（即制氧价格）核算生态系统提供氧气的价值。
38
10.2.2价值量核算模型
𝑄𝑜𝑝 𝐶𝑜 （53）
式中，Vop为生态系统释氧价值（元/a）；Qop为生态系统氧气释放量（t氧气/a），𝐶𝑜为工业制氧价格（元/t）。 10.2.3定价参数与数据来源
生态系统释氧量由实物量核算得到，制氧价格根据实际调查获得。
39
11空气净化实物量和价值量核算
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∑𝑗<1𝑄𝑖𝑗 𝐴𝑖 （55）
40
式中，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）
41
式中，Va为生态系统大气环境净化的价值（元/a）；𝑄𝑎𝑝𝑖为第i种大气污染物的净化量（t/a），j 为大气污染物类别，i= 1，2，…，n，无量纲； 𝐶𝑖为第i类大气污染物的治理成本（元/t）。 11.2.3 定价参数与数据来源
污染物净化量由实物量核算得到。单位治理成本采用核算地方印发的排污费征收标准，没有地方标准的，可以参考国家发展和改革委员会发布的《排污费征收标准及计算方法》收费标准或者《中华人民共和国环境保护税法》中的税收标准。
42
12水质净化实物量和价值量核算
12.1 水质净化实物量核算 12.1.1 指标内涵
水质净化功能是指湖泊、河流、沼泽等水域湿地生态系统吸附、降解、转化水体污染物，净化水环境的功能。
水质净化服务价值核算主要是利用监测数据，根据水体生态系统中污染物构成和浓度变化，选取适当的指标对其进行定量化核算。常用指标包括氨氮、COD、总氮、总磷以及部分重金属等。根据我国《地表水环境质量标准》（GB3838-2002）中对水环境质量应控制的项目和限值的规定，选取指标作为生态系统水环境净化功能的评价指标。 12.1.2 核算方法
水质净化功能核算依据污染物浓度是否超过地表水水域环境功能和保护目标而选择不同的方法。若污染物浓度未超过地表水水域环境功能标准限值，则采用方法1进行核算；若污染物排放浓度超过地表水水域环境功能标准限值，则采用方法2进行核算。地表水水域环境功能标准限值见附录C.5表C.10。
水质净化功能－方法1：如果污染物排放量未超过地表水水域环境功能标准限值，则采用污染物排放量核算实物量。
𝑄𝑤𝑝=∑𝑛𝑖<1𝑃𝑖 （57）
式中，Qwp为水体污染物净化量（kg/a）；𝑃𝑖为i类污染物排放量（kg/a），包括总氮、总磷、COD等；i 为污染物类别，i = 1，2，…，n，无量纲；n为水体污染物类别的数量，无量纲。
43
水质净化功能－方法2：如果污染物排放量超过地表水水域环境功能标准限值，根据生态系统的自净能力核算实物量：
𝑛
𝑄𝑤𝑝=∑𝑚 <1∑𝑗<1𝑃𝑖𝑗 𝐴𝑖 （58）
式中，Qwp为污染物净化总量（kg）；𝑃𝑖𝑗为某种生态系统单位面积污染物净化量（kg/km2）；𝐴𝑖为生态系统面积（km2）；i = 1，2，…，m，无量纲；m为生态系统类型的数量，无量纲；j为污染物类别，j = 1，2，…，n，无量纲；n为水体污染物类别的数量，无量纲。 12.1.3 核算参数及数据来源
污染物排放数据从生态环境部门获取；各类生态系统面积来源于自然资源部门；生态系统对污染物的单位面积净化量来源于参考文献或实地监测。
12.2 水质净化价值量核算 12.2.1 定价思路
核算生态系统降解水体污染物、净化水质的价值，运用替代成本法核算生态系统水质净化功能的价值。 12.2.2 价值量核算模型
水质净化价值量计算采用替代成本法，通过工业治理水污染物成本核算生态系统水质净化价值。化学需氧量、氨氮净化价值计算方法：运用化学需氧量、氨氮两种污染物水质净化实物量，分别乘以单位化学需氧量、氨氮处理的费用，核算水体净化价值。
𝑛
𝑉𝑤=∑𝑖<1𝑄𝑤𝑝𝑖 𝐶𝑖 （59）
44
式中，𝑉；𝑄𝑤𝑝𝑖为第i类水污𝑤为生态系统水质净化的价值（元/a）染物的净化量（t/a）；𝐶𝑖为第i类水污染物的单位治理成本（元/t）；i为研究区第i类水体污染物类别，i = 1，2，…，n，无量纲。n为研究区水体污染物类别的数量，无量纲。 12.2.3 定价参数与数据来源
污染物净化量由实物量核算得到。COD、氨氮水质等污染物单位治理成本采用核算地方印发的排污费征收标准，没有地方标准的，可以参考国家发展和改革委员会发布的《排污费征收标准及计算方法》收费标准或者《中华人民共和国环境保护税法》中的税收标准。
45
13气候调节实物量和价值量核算
13.1 气候调节实物量核算 13.1.1 指标内涵
生态系统气候调节服务是指生态系统通过植被蒸腾作用、水面蒸发过程吸收太阳能，降低气温、增加空气湿度，改善人居环境舒适程度的生态功能。
选用生态系统蒸散发过程消耗的能量作为生态系统气候调节服务的评价指标。 13.1.2 核算方法
气候调节服务核算可用实际测量生态系统内外温差、生态系统消耗的太阳能量和生态系统的总蒸散量进行核算，优先选择实际测量方法，其次根据数据可得性选取生态系统消耗的太阳能量方法或生态系统的总蒸散量进行核算。
方法一：采用实际测量生态系统内外温差进行实物量转换。
𝑄=∑𝑛𝑖<1𝛥𝑇𝑖 𝑐 𝑉 （60）
式中，𝑄为吸收的大气热量（J/a）； 𝑐为空气的比热容（J/m3·℃-1）；𝑉为生态系统内空气的体积（m3）；𝛥𝑇𝑖为第𝑖天生态系统内外实测温差（℃）； 𝑛为年内日最高温超过26℃的总天数。
方法二：采用生态系统消耗的太阳能量作为气候调节的实物量。
46
𝐶𝑅𝑄=𝐸𝑇𝐸−𝑁𝑅𝐸 （61）
式中，CRQ为生态系统消耗的太阳能量（J/a）；ETE为森林、草地、灌丛、湿地等生态系统蒸腾作用消耗的太阳能量（J/a）；NRE为森林、草地、湿地等生态系统吸收的太阳净辐射能量（J/a）。
方法三：采用生态系统蒸腾蒸发总消耗的能量作为气候调节的实物量。
𝐸𝑡𝑡=𝐸𝑝𝑡+𝐸𝑤 （62）
6
𝐸𝑝𝑡=∑3𝑖𝐸𝑃𝑃𝑖 𝑆𝑖 𝐷 10 (3600∗𝑟) （63）
𝐸𝑤 =𝐸𝑤 𝑞 103 (3600)+𝐸𝑤 𝑦 （）
式中，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 核算参数及数据来源
水面蒸发量、植被蒸散发量、生态系统面积、单位面积蒸腾耗热量等数据来自气象、自然资源、林业等相关部门和文献资料。
47
13.2 气候调节价值量核算 13.2.1 定价思路
运用替代成本法（即人工调节温度和湿度所需要的耗电量）核算生态系统蒸腾调节温度或湿度价值和水面蒸发调节温度或湿度价值。 13.2.2 价值量核算模型
𝑉𝑡𝑡=𝐸𝑡𝑡 𝑃 （65）
式中，Vtt为生态系统气候调节的价值（元/a）；𝐸𝑡𝑡为生态系统调节温度或湿度消耗的总能量（kW·h/a）；𝑃 为当地电价（元/ kW·h）。 13.2.3 定价参数与数据来源
生态系统调节温度或湿度所耗能量由实物量核算得到。电价从核算地方发展与改革委员会发布的相关文件或供电部门获取，一般参考工业电价。
48
14物种保育实物量和价值量核算
14.1物种保育实物量核算 14.1.1指标内涵
物种多样性是生物多样性最主要的结构和功能单位，可以为生态系统演替与生物进化提供必须的物种与遗传资源，是人类生存和发展的基础。物种保育服务是指生态系统为珍稀濒危动植物物种提供生存与繁衍的场所，从而对其起到保育作用的作用和价值。 14.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为自然保护区面积。
49
14.2物种保育价值量核算 14.2.1定价思路
运用单位面积保育成本核算物种保育服务的价值量。 14.2.2价值量核算模型
实物量与物种保育成本相乘得到价值量： 方法1：Shannnon-Weiner价值法 实物量与物种保育价值相乘得到价值量：
𝑉 𝑖𝑜=𝐺 𝑖𝑜 𝑆生 （68）
式中，𝑉 𝑖𝑜生物多样性价值（元/a）；𝑆生为单位面积物种保育价值（元/ha·a-1）。
方法2：保护区保护价值法
𝑉 𝑖𝑜𝑝=𝐺 𝑖𝑜𝑝 𝑆𝑐 （69）
式中，𝑉 𝑖𝑜𝑝 物种保育价值（元/a）；G为物种保育实物量； 𝑆𝑐 为自然保护区单位面积保育成本（元/ha·a-1）。 14.2.3核算参数及数据来源
保护区总投入、物种数量参照相关统计、报告和文献资料。不同Shannnon-Weiner指数对应等级划分，以及相应的单位面积物种多样性保护价值参照附录C.6。
50
15休闲旅游实物量和价值量核算
15.1休闲旅游实物量核算 15.1.1 指标内涵
人类通过精神感受、知识获取、休闲娱乐和美学体验从生态系统获得的非物质惠益。 15.1.2 核算方法
采用区域内自然景观的年旅游总人次作为文化服务的实物量评价指标。
Nt ∑ Nt （70）
式中，Nt为游客总人数；Nt 为第i个旅游区的人数；n为旅游区个数，i=1，2，…，n。 15.1.3 核算参数及数据来源
自然景观名录、旅游人数与旅客来源通过旅游、园林、统计等部门或问卷调查获取。 15.2 休闲旅游价值量核算 15.2.1 定价思路
运用旅行费用法核算人们通过休闲旅游活动体验生态系统与自然景观美学价值，并获得知识和精神愉悦的非物质价值。
51
15.2.2 价值量核算模型
𝑉𝑟=∑𝑗<1𝑁𝑗 𝑇𝐶𝑗 （71） 𝑇𝐶𝑗=𝑇𝑗 𝑊𝑗+𝐶𝑗 （72） 𝐶𝑗=𝐶𝑡𝑐 𝑗:𝐶𝑙𝑓 𝑗:𝐶 𝑓 𝑗 （73）
式中，𝑉𝑟表示被核算地点的休闲旅游价值（元/a）；𝑁𝑗表示j地到核算地区旅游的总人数（人/a）；𝑗=1 2 ，𝐽表示来被核算地点旅游的游客所在区域（区域按距离核算地点的距离划同心圆，如省内、省外等）。𝑇𝐶𝑗表示来自j地的游客的平均旅行成本（元/人）；𝑇𝑗表示来自j地的游客用于旅途和核算旅游地点的平均时间（天/人）；𝑊𝑗表示来自j地的游客的当地平均工资（元/（人·天））；𝐶𝑗表示来自j地的游客花费的平均直接旅行费用（元/人），其中包括游客从j地到核算区域的交通费用𝐶𝑡𝑐 𝑗（元/人）、食宿花费𝐶𝑙𝑓 𝑗（元/人）和门票费用𝐶 𝑓 𝑗（元/人）。 15.2.3 定价参数与数据来源
自然景观名录、旅游人数通过旅游、园林等部门获取，游客的社会经济特征、旅行费用情况等通过问卷调查获得。
𝐽
52
16景观价值实物量和价值量核算
16.1景观价值实物量核算 16.1.1指标内涵
生态系统的景观价值是指森林、湖泊、河流、海洋等生态系统可以为其周边的人群提供美学体验、精神愉悦的功能，从而提高周边土地、房产价值。 16.1.2核算方法
采用能直接从自然生态系统获得景观价值的土地与居住小区房产面积作为景观价值实物量评价指标。
Al=∑ =1Al （74）
式中，𝐴𝑙为从自然生态系统景观获得升值的土地与居住小区房产总面积（km2/a）；𝐴𝑙𝑖为第i区的房产面积（km2），i = 1，2，…，n。 16.1.3核算参数及数据来源
受益土地与居住区名录及面积通过调查获取。 16.2景观价值价值量核算 16.2.1定价思路
运用享乐价值法核算生态系统为其周边地区人群提供美学体验、精神愉悦功能的价值。
53
16.2.2价值核算模型
𝑉𝑎=𝐴𝑎 𝑃𝑎 （75）
2
式中，𝑉为景观价值（元/a）；𝐴为受益总面积（km）；𝑃𝑎为由𝑎𝑎
生态系统带来的单位面积溢价（元/km2·a-1）。 16.2.3核算参数及数据来源
受益面积由实物量核算得到。生态系统带来的单位面积溢价由实地调研获取。
54
17附则
17.1 核算数据
基础数据是开展生态系统生产总值核算的重要保障，涉及的数据内容包括发改、农业、国土、林业、水务、气象、环保、旅游、电力、统计等部门数据。基础数据获取时，应确保数据的权威性、准确性、时效性和可获得性。根据核算需求，应从数据的完整性、代表性、准确性与可比性等角度进行分析筛选，采用权威部门的遥感、监测、普查、调查、统计等数据，尽可能使用高分辨率空间数据，选取适合核算工作需要的数据资料（参考来源见附录A）。同时，在核算过程中，应保证所处理数据资料的完整性、合理性和有效性，并对原始数据、拷贝数据及数据分析结果进行校核，数据处理应符合相应规范的要求。
为使核算报告准确可信，在核算过程中应对数据的获取与处理进行质量控制，降低不确定性。在核算过程中，要确保使用的模型和收集到的数据能够代表实际情况。 17.2 成果表达形式
核算成果主要包括核算报告、生态系统生产总值核算总表。 核算总表用表格形式表达GEP各项指标的实物量、单价与价值量。
17.3 结果校核
核算完成后，应该将核算结果与其他地区已有的核算结果对比，对比不同生态系统服务功能单位面积实物量、价值量的差异，验证结果的合理性。
55
17.4 信息管理
核算主体应记录并保存下列资料，保存时间不少于5年： （1）核算方法相关信息：
选择基于模型的方法时，应保存以下内容：a）获取数据和参数的相关资料；b）不确定性及如何降低不确定性的相关说明。
选择基于实地调查的方法时，应保存以下内容：a）有关职能部门出具的测量仪器证明文件；b）连续测量的所有原始数据；c）不确定性及如何降低不确定性的相关说明；d）验证计算，应保留所有基于计算的保存内容。
（2）数据质量控制相关记录文件。
56
附录A（资料性附录）生态系统生产总值核算数据来源表
表A.1 生态系统生产总值核算数据来源表
服务类别 评估科目 谷物 薯类 油料 豆类 农业产品 糖料 蔬菜 水果 …… 小计 竹木采伐 木材 竹材 五倍籽 核桃 林业产品 林下产品 板栗 松子 物质产品 …… （根据评估样点实际产品更新） 评估指标 谷物 薯类 油料 豆类 糖料 蔬菜 水果 …… 单位 数据来源 统计部门，农业部门 （根据评估样点实际产品更新） 统计部门，林业部门 产量、单价、产值 - 小计 牛肉 羊肉 猪肉 禽肉 畜禽产量 畜牧业 产品 猪出栏数 牛出栏数 羊出栏数 家禽出栏数 …… （根据评估样点实际产品更新） 统计部门,畜牧业部门 奶类 禽蛋 小计 渔业产品 水产品 牛奶 …… 禽蛋 养殖 捕捞 统计部门 57
…… （根据评估样点实际产品更新） 小计 水能 生态能源 生态能源 生物质能源（秸秆、薪柴等） 小计 森林面积 灌丛面积 草地面积 湖泊面积 沼泽面积 河流面积 湿地面积 水源涵养 气候调节 土壤保持 防风固沙 水田面积 旱地面积 园地面积 多年均降雨量 多年均蒸发量 年降雨量 年蒸发量 多年月均温度 多年月均湿度 调节服务 年均夏季气温 土壤容重 水库单位库容造价 区域入境水量 区域出境水量 水源涵养 洪水调蓄 土壤保持 防风固沙 水库清淤工程费用 土壤中氮的纯含量 土壤中磷的纯含量 环境工程降解氮工程费用 环境工程降解磷工程费用 水库总库容 水库防洪库容 碳固定 氧气提供 空气净化 NPP 生物量 二氧化硫排放量 氮氧化物排放量 颗粒物排放量
58
统计部门 km² km² km² km² km² km² km² km² km² km² mm mm mm mm ℃ % ℃ t/m3 元/m3 亿m3 亿m3 元/m3 % % 元/t 元/t 万m3 万m3 亿t 亿t 万t/a 万t/a 万t/a 林业部门 生态环境部门 发改部门 水利部门 水务部门 农业部门 环境部门 自然资源部门 气象部门 自然资源部门
二氧化硫治理费用 氮氧化物治理费用 颗粒物治理费用 废水排放总量 工业废水排放量 生活污水排放量 污水处理总量（污水厂） COD 水质净化 总氮 总磷 排放量 污水处理厂处理量 排放量 污水处理厂处理量 排放量 污水处理厂处理量 COD治理费用 氨氮治理费用 因生态系统降低水环境治理成本的公司数量 因生态系统降低水环境治理的费用 气候调节 休闲旅游价值 景观价值 电价 旅游人数 旅游收入 受益土地面积 受益居住区名录 元/t 元/t 元/t 万t/a 万t/a 万t/a 万t/a 万t/a t/a 万t/a t/a 万t/a t/a 元/t 元/t 个 万元 元/kW·h 问卷调查 统计部门 电力部门 文化服务 万人/年 旅游发展改万元/年 革委员会 万亩 自然资源部门、建设部门
59
附录B（资料性附录）陆地生态系统生产总值核算报告编写大纲 前言
介绍核算背景、目的意义、任务来源等 1 区域概况
介绍核算区域地理范围、自然环境状况、经济社会状况、生态环境保护状况等基本情况。 2 核算目标与原则
介绍核算目标、核算原则、核算依据、核算基准年。 3 核算方法与数据
介绍主要核算思路、方法、数据来源与数据处理方法。 4 生态系统生产总值实物量
介绍生态系统生产总值实物量的计算过程与结果，包括物质产品实物量、调节服务实物量和文化服务实物量。 5 生态系统生产总值价值量
介绍生态系统生产总值价值量的计算过程与结果，包括物质产品价值量、调节服务价值量和文化服务价值量。 6 结论与建议
介绍生态系统生产总值核算的结论，提出相关建议。 7 附件
包括生态系统生产总值核算过程中相关的技术资料及附表、附图等。
60
附录C（资料性附录）生态系统生产总值实物量核算参数参考表 附录C.1 土壤保持实物量核算参数
降雨侵蚀力因子R、土壤可蚀性因子K、坡长坡度因子L、S的算法以及植被覆盖和管理因子C以及水土保持措施因子P的计算方法可参照下述方法：
（1）降雨侵蚀力因子（R）
降雨侵蚀力是土壤侵蚀的驱动因子，与土壤侵蚀强度有直接的关系。降雨侵蚀力R计算可分为EI30经典计算方法和常规气象资料简易算法两类。由于降雨动能E和30 min降雨强度I30资料获取难度较大，所以国内外许多学者根据区域性降雨侵蚀特点，建立了基于常规降雨量资料的简易模型。可采用全国日降雨量拟合模型来核算降雨侵蚀力，是基于日降雨量资料的半月降雨侵蚀力模型。其公式如下：
𝑘
𝑀𝑖=α∑𝐷𝑗𝛽
𝑗<1
其中，Mi 为某半月时段的降雨侵蚀力值 (MJ·mm·hm-2·h-1·a-1) ；Dj表示半月时段内第j天的侵蚀性日雨量（要求日雨量大于等于12 mm，否则以0计算，阈值12 mm与中国侵蚀性降雨标准一致；k表示半月时段内的天数，半月时段的划分以每月第15日为界，每月前15天作为一个半月时段，该月剩下部分作为另一个半月时段，将全年依次划分为24个时段。
α、β是模型待定参数：
𝛽=0.8363+
18.144𝑃 1224.455
+ 𝑃𝑦12
61
𝛼=21.586 𝛽;7.11
其中，Pd12表示日雨量12 mm以上（包括等于12 mm）的日平均雨量；Py12表示日雨量12 mm以上(包括12 mm)的年平均雨量。
（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.002 mm）；粉粒为（0.002-0.05 mm）；极细砂为（0.05-0.1 mm）；砂粒为（0.1-2.0 mm）。
表C.1诺谟图中结构性指数与可渗透性指数的定义
结构性指数S 1 2 3 4 含义 非常坚固（very structured or particulate） 很坚固（fairly structured） 较坚固(slightly structured) 坚固(solid) 可渗透性指数P 1 2 3 4 5 6 含义 快速（Papid） 中快速(Moderate to rapid) 中速(Moderate) 中慢速(Moderate to slow) 慢速(Slow) 极慢(Very slow)
62
（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.0)
𝛽= [3.0 (𝑠𝑖𝑛𝜃)0.8+0.56]其中，L为坡长因子； 𝜆 为坡长，通过汇流计算得到，m为无量纲常数，S为坡度因子，t为百分比坡度，𝜃为坡度，单位是弧度。
（4）植被覆盖和管理因子（C）
C因子是指在一定的植被覆盖度和管理措施下，一定面积土地上的土壤流失量与采取连续清耕、休闲处理的相同面积土地上的流失量的比值，为无量纲数，介于0-1之间。要确定C因子的值，需要详细的气候、土地利用、前期作物残留量、土壤湿度等资料，在大尺度研究中，一般难以获取这些资料，且C值的经典算法非常复杂，国内部分学者采用植被覆盖度求解C值，可采用以下C值计算方法：
63
1 𝑓=0.1
C={0.6508−0.3436 𝑙𝑔𝑓 0.1 <𝑓≤78.3%
0 𝑓>78.3%
上述公式中，植被覆盖度 f 基于植被指数NDVI数据计算得到，公式如下：
(𝑁𝐷𝑉𝐼−𝑁𝐷𝑉𝐼𝑠𝑜𝑖𝑙)
𝑓=
(𝑁𝐷𝑉𝐼𝑚𝑎𝑥−𝑁𝐷𝑉𝐼𝑠𝑜𝑖𝑙)其中，𝑁𝐷𝑉𝐼𝑠𝑜𝑖𝑙为纯裸土象元的NDVI值；𝑁𝐷𝑉𝐼𝑚𝑎𝑥纯植被象元的NDVI值。
附录C.2 防风固沙实物量核算参数
气候侵蚀因子、地表糙度因子、土壤侵蚀因子、土壤结皮因子、植被覆盖因子的计算方法可参照下述方法：
（1）气候侵蚀因子
气候侵蚀因子中的风因子和土壤湿度因子可采用从中国气象科学数据共享服务网上下载（http://cdc.cma.gov.cn）的国家台站的日均风速、降水、温度、日照时数，纬度等来计算完成；雪盖因子可采用从中国西部环境与生态科学数据中心(http://westdc.westgis.ac.cn)下载的中国雪深长时间序列数据集来计算雪盖因子。
WFWSWS2i1N2-WStNd2NgSWSD
式中，WF为气候因子，kg/m；WS2为2 m处风速，m/s；WSt为2 m处临界风速（假定为5 m/s）；N为风速的观测次数（一般500次）；Nd为试验的天数，d；ρ为空气密度，kg/m3；g为重力加速度，m/s2；SW为土壤湿度因子无量纲；SD为雪覆盖因子。
i) 其中空气密度通过下式求得：
1.013-0.1183EL0.0048EL2348.0 T式中，EL为海拔高度，单位为km，T为绝对温度，单位为开氏度。
ii) 土壤湿度通过下式求得：
ETpRISWETpRdNd
Etp可采用Samani and Pessaralkli（1986）的方法，公式如下：
ETp0.0162SR(DT17.8)58.5
式中：SW为土壤湿度因子；ETp为潜在相对蒸发量，mm；R为降雨量，mm；I为灌溉量，mm；Rd为降雨次数和（或）灌溉天数：Nd为天数，d（一般15d）；SR为太阳辐射总量，cal/cm2；DT为平均温度，℃。
其中太阳辐射采用方法如下：
𝑅𝑠0
𝑛
=(𝑎𝑠+𝑏𝑠)𝑅𝑎
𝑁24𝑁= 𝜔𝑠
𝜋24(60)𝑅𝑎= 𝐺𝑠𝑟 𝑑𝑟 [𝜔𝑠 𝑠𝑖𝑛(𝜑) 𝑠𝑖𝑛(𝛿)+𝑐𝑜𝑠(𝜑) 𝑐𝑜𝑠(𝛿) 𝑠𝑖𝑛(𝜔𝑠)]
𝜋2𝜋
𝑑𝑟=1+0.033 𝑐𝑜𝑠( 𝐽)
3652𝜋
𝛿=0.409 𝑠𝑖𝑛( 𝐽−1.39)
365𝜔𝑠=𝑎𝑟𝑐𝑐𝑜𝑠[−𝑡𝑎𝑛(𝜑) 𝑡𝑎𝑛(𝛿)]
65
或 𝜔𝑠=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.0820 MJ·m-2·min-1）；dr：日地相对距离；ωs：日落时角（rad）；φ：纬度（rad）；δ：太阳倾角（rad）；J：对应于一年中的第几天。
iii) 雪盖因子通过下式求得：
SD1-Psnow cover25.4mm
其中，P(snow cover>25.4 mm)为计算时段内积雪覆盖深度大于25.4mm的概率。
（2）土壤侵蚀因子
土壤颗粒分为可蚀性土粒和非可蚀性土粒，粒径大于0.84mm的土粒不易被风蚀，称为非可蚀性颗粒。土壤可蚀性因子则为土壤表层直径小于0.84 mm的颗粒的含量。
可用以下方程来描述土壤可蚀性因子EF的值：
EF29.090.31Sa0.17Si0.33Sa/Cl2.59OM0.95CaCO3
100其中，Sa为土壤砂粒含量；Si为土壤粉砂含量；Sa/Cl为土壤砂粒和粘土含量比；OM为有机质含量；CaCO3为碳酸钙含量。
66
（3）土壤结皮因子
土壤结皮（soil crust）为土壤颗粒物（特别是黏土、粉砂与有机质颗粒）的胶结作用而在土壤表面生成一层物理、化学和生物性状均较特殊的土壤微层。
SCF110.0066(Cl)20.021(OM)2
其中，CL为粘土含量（5.0～39.3），OM为有机质含量（0.18～4.79）。
表C.2 RWEQ标准数据库中物质含量范围表
Sa 范围（%） 5.5～93.6 Si 0.5～69.5 Cl 5.0-39.3 Sa/Cl l.2～53.0 OM 0.18～4.79 CaCO3 0.0～25.2 土壤数据可采用中国西部环境与生态科学数据中心(http://westdc.westgis.ac.cn)提供的1：100万土壤图所附的土壤属性表和空间数据。对于土壤可蚀性和土壤结皮而言，由于我国土壤颗粒分级与美国制不同，RWEQ中的分级使用的是美国制，为此需要先对土壤颗粒含量进行粒径转换，且实测的土壤颗粒含量参数需符合RWEQ标准数据库中的物质含量范围表，当实测值不符合要求时，可以使用RWEQ内嵌的土壤质地资料的输入参数。
表C.3 RWEQ模型内嵌适用的土壤资料 % 土壤类型 砂土 壤质砂土 砂质壤土 砂质粘壤土 砂质粘土 泥土 粉质壤土
Sa 93 84 59 52 6 21 67
Si 4 10 26 13 7 88 67 OM 0.3 0.5 0.5 1 1 1.5 1.5 CaCO3 1 2 3 3 3 3 3
壤土 粉质黏壤土 粉质粘土 粘壤土 黏土 41 10 6 32 20 41 56 47 34 20 1.5 2 2.5 2.5 3 3 3 3 3 3 （4）地表粗糙度因子
土壤糙度因子K′取决于自由糙度RR和定向糙度OR，可采用一种滚轴式链条法来测定地表糙度。其基本原理是：两点间直线距离最短，当地表糙度增加时，其地表距离随之增加。于是当一个给定长度为L1的链条放于粗糙的地表时，其水平长度将缩小为L2，L1和L2的差值和地表粗糙程度密切相关，可用下式来计算地表糙度。
L2Crr1L100
1式中，Crr为任意方向上的地表糙度。这里还涉及到一个垂直于垄向的糙度PR，当用链条法测出RR和PR时，OR就可通过下式来求得。
OR=PR-RR
土垄糙度因子可通过下式来进行计算，假设土垄呈等腰三角形则：
式中：Nr为在长度L2范围内土垄的数量。根据实验结果，当不考虑风向的影响时，土壤糙度因子（K′）与土垄糙度因子（Kr）和随机糙度（RR）的回归方程为：
K′e1.86K2.41Krr0.9340.124Crr
也可采用风向与垄成任意角度时的糙度因子公式：
68
K′eR1.86K2.41Kcrr0.9340.124C
rrRc1.01024.717.3310223.741043
式中，θ为风向与垄平行方向的夹角（角度）。 （5）植被覆盖因子
该因子用来确定生长植被、枯萎植被（农作物为倒伏残茬）以及农作物的直立残茬对土壤风蚀的影响。综合植被覆盖因子COG为枯萎植被、直立残茬和生长植被覆盖三因子的乘积。
i)枯萎植被
（sc） SLRfe-0.0438式中，SLRf为枯萎植被的土壤流失比率；SC为枯萎植被地表覆盖率。
ii)直立残茬
（SASLRse-0.03440.13）
式中，SLRs为直立残茬土壤流失比率；SA为直立残茬当量面积，是1m2内直立秸秆的个数乘以秸秆直径的平均值（cm）再乘以秸秆高度（cm）。
iii)生长植被
（ccSLRce-5.6140.7366）
式中，SLRc为植被覆盖土壤流失比率；cc 是土表植被覆盖度。 附录C.3 固碳实物量核算参数
NEP和NPP的转换系数α取值如表C.4所示：
表C.4 各地理区和各省不同生态系统NEP和NPP的转换系数
69
地理区 省份 平均 北京 天津 河北 山西 内蒙 平均 辽宁 吉林 黑龙江 平均 上海 江苏 浙江 安徽 福建 江西 山东 平均 河南 湖北 森林 0.330 0.361 0.363 0.354 0.376 0.315 0.211 0.172 0.139 0.252 0.149 0.149 0.140 0.158 0.163 0.175 0.112 0.188 0.120 0.208 0.162 0.097 0.136 0.085 0.169 0.209 0.176 0.255 0.049 0.226 0.215 0.204 0.222 0.254 0.239 0.210 0.035 灌丛 0.297 0.346 0.360 0.340 0.335 0.112 0.146 0.169 0.057 0.166 0.116 0.116 0.082 0.151 0.112 0.139 0.085 0.197 0.099 0.200 0.130 0.070 0.090 0.075 0.143 0.142 0.183 0.203 0.008 0.163 0.087 0.104 0.239 0.196 0.191 0.135 0.191 草地 0.152 0.281 0.068 0.242 0.218 0.129 0.173 0.181 0.125 0.198 0.130 0.015 0.094 0.112 0.129 0.138 0.118 0.146 0.098 0.173 0.113 0.084 0.106 0.063 0.032 0.226 0.170 0.287 0.060 0.136 0.262 0.166 0.233 0.146 0.278 0.169 0.015 湿地 0.197 0.086 0.053 0.103 0.077 0.208 0.147 0.076 0.075 0.160 0.0 0.059 0.069 0.052 0.045 0.080 0.065 0.094 0.045 0.062 0.038 0.084 0.031 0.020 0.025 0.1 0.082 0.146 0.010 0.080 0.302 0.208 0.096 0.150 0.233 0.056 0.148 华北 东北 华东 中南 湖南 广东 广西 海南 平均 重庆 四川 贵州 云南 平均 陕西 甘肃 青海 宁夏 西南 西北 生物量-碳转换系数𝐶 ：森林和灌丛的转化系数为0.5，草地的转化系数为0.45。
70
森林及灌丛的固碳速率FCSR由森林清查数据计算获得，如表C.5所示：
表C.5 各省森林（及灌丛）生态系统固碳速率
区域 全国 北京 天津 河北 山西 内蒙古 辽宁 吉林 黑龙江 上海 江苏 浙江 安徽 福建 江西 山东 河南 湖北 湖南 广东 广西 海南 重庆 四川 贵州 云南 陕西 甘肃 青海 宁夏 森林(及灌丛)植被固碳速率（tC·ha-1·a-1） 2000年 0.53 0.91 0.90 0.17 0.52 0.28 0.24 0.36 0.26 1.52 0.48 0.36 0.66 0.87 0.38 0.35 0.63 0.62 0.54 0.29 0.56 0.26 1.33 0.51 1.16 1.02 1.50 0.40 0.79 0.20 1.04 0.99 2010年 0.75 1.00 1.34 1.12 0.85 0.46 0.92 0.26 0.75 4.16 2.24 0.85 1.11 0.80 0.92 2.83 1.94 1.14 1.00 0.19 0.95 0.03 1.33 0.65 1.28 0.70 0.46 0.59 0. 0.75 1.04 1.00 2015年 0.84 0.84 2.11 0.77 1.02 1.04 1.26 0.97 0.63 1.97 2.60 0.71 1.15 1.75 0.29 1.05 1.12 1.26 1.00 0.79 0.49 1.00 1.24 0.63 1.17 0.66 0.12 0.84 0.90 0.92 1.31 1.36 森林(及灌丛)土壤固碳速率 （tC·ha-1·a-1） 2000年 0.34 0.59 0.58 0.11 0.33 0.18 0.15 0.23 0.17 0.98 0.31 0.23 0.43 0.56 0.25 0.22 0.41 0.40 0.35 0.19 0.36 0.17 0.86 0.33 0.75 0.66 0.97 0.26 0.51 0.13 0.67 0. 2010年 0.48 0.65 0.86 0.72 0.55 0.30 0.59 0.16 0.48 2.69 1.45 0.55 0.71 0.52 0.59 1.83 1.25 0.74 0.65 0.13 0.61 0.02 0.86 0.42 0.82 0.45 0.30 0.38 0.57 0.49 0.67 0. 2015年 0.54 0.55 1.36 0.50 0.66 0.67 0.82 0.63 0.40 1.27 1.68 0.46 0.74 1.13 0.19 0.68 0.73 0.81 0.65 0.51 0.31 0.65 0.80 0.41 0.76 0.43 0.08 0.54 0.58 0.60 0.84 0.88 森林及灌丛土壤固碳系数β取值0.6。
71
全国草地（除青藏高原外）土壤的固碳速率为0.02 tC·ha-1·a-1，青藏高原区域为0.03 tC·ha-1·a-1。草地土壤固碳速率GSR见表C.6。
表C.6 草地土壤固碳速率
区域 全国 北京 天津 河北 山西 内蒙古 辽宁 吉林 黑龙江 上海 江苏 浙江 安徽 福建 江西 山东 草地土壤固碳速率 （tC·ha-1·a-1） 0.02 0.02 0.02 0.02 0.02 0.06 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 区域 河南 湖北 湖南 广东 广西 海南 重庆 四川 贵州 云南 陕西 甘肃 青海 宁夏 草地土壤固碳速率 （tC·ha-1·a-1） 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.03 0.02 0.03 0.03 0.02 0.02 0.03 0.02 0.03 水域湿地的固碳速率SCSRi取值见表C.7：
表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 无化学肥料和有机肥料施用的情况下，我国农田土壤有机碳的变
72
化NSC取-0.06；各省土壤容重BD取值见表C.8；土壤厚度取20cm。
表C.8 土壤容重
省份 全国 北京 天津 河北 山西 内蒙古 辽宁 吉林 黑龙江 上海 江苏 浙江 安徽 福建 江西 山东 容重（g/cm3） 1.250 1.2 1.309 1.294 1.261 1.271 1.288 1.188 1.141 1.194 1.253 1.177 1.241 1.297 1.281 1.304 省份 河南 湖北 湖南 广东 广西 海南 重庆 四川 贵州 云南 陕西 甘肃 青海 宁夏 容重（g/cm3） 1.312 1.261 1.241 1.225 1.245 1.244 1.233 1.240 1.207 1.230 1.225 1.274 1.261 1.236 1.248 1.298 作物j的草谷比SGRj取值如表C.9所示：
表C.9 不同作物的草谷比
作物 水稻 小麦 玉米 高粱 马铃薯 草谷比SGRj 0.623 1.366 2 1 0.5 作物 油菜 向日葵 棉花 甘蔗 草谷比SGRj 2 2 8.1 0.1 区域地下水径流模数以及HCO3-浓度见表C.10：
表C.10 区域地下水径流模数以及HCO3- 浓度
岩溶类型区 南方岩溶区 北方岩溶区 青藏高原区 埋藏岩溶区
HCO3-浓度（g/L） 0.231 0.245 0.145 0.225 M-岩溶水径流模数（107L/(km2·a)） 40.59 20.86 19.95 0.375 73
附录C.4 空气净化实物量核算参数
表C.11 环境空气污染物浓度限值
污染物 二氧化硫 二氧化氮 颗粒物PM10 颗粒物PM2.5 平均时间 年平均 年平均 年平均 年平均 年平均浓度 一级 20 40 40 15 二级 60 40 70 35 ug/m3 单位 注：环境空气功能区分为两类，一类区为自然保护区、风景名胜区和其他需要特殊保护的区域；二类区为居住区、商业交通居民混合区、文化区、工业区和农村地区。一类区适用一类浓度限值，二类区适用二类浓度限值。核算过程中，将核算区域大气污染物监测点位的算术平均值与所在功能区的空气浓度限值进行比较，来确定核算方法。
附录C.5水质净化实物量核算参数
表C.12 地表水污染物浓度限值
污染物 化学需氧量 氨氮 I类 15 0.15 II类 15 0.5 III类 20 1 IV类 30 1.5 V类 40 2 注：地表水水环境功能分为五类，I类适用于源头水、国家自然保护区；II类适用于集中式生活饮用水地表水源地一级保护区、珍稀水生生物栖息地、鱼虾类产卵场、仔稚幼鱼的索饵场等；III类适用于集中式生活饮用水地表水源地二级保护区、鱼虾类越冬场、洄游通道、水产养殖区等渔业水域及游泳区；IV类适用于一般工业用水区及人体非直接接触的娱乐用水区；V类适用于农业用水区及一般景观要求水域。核算过程中，将核算区域水质监测断面的污染物浓度算术平均值与所在功能区的污染物浓度限值进行比较，来确定核算方法。
附录C.6生物多样性保护核算参数
表C.13 物种濒危指数体系
濒危指数 4 3 2 濒危等级 极危 濒危 易危 74
1 近危 表C.14 特有种指数体系
特有种指数 4 3 2 1 0 分布范围 仅限于范围不大的山峰或特殊的自然地理环境下分布 仅限于某些较大的自然地理环境下分布的类群，如仅仅分布于较大的海岛、高原、若干个山脉等 仅限于某个分布的分布群 至少在2个都有分布的分布群 世界广布的分类群 表C.15 古树年龄指数体系
指数等级 1 2 3 古树年龄 100-299年 300-499年 ≥500年 表C.16 Shannon-Wiener 指数等级划分及其价值量
等级 I II III IV V VI Shannon-Wiener 指数 指数≥6 5≤指数≤6 4≤指数≤5 3≤指数≤4 2≤指数≤3 1≤指数≤2 单价/（元/hm2·a） 66760 53410 40060 26700 13350 6680
75
因篇幅问题不能全部显示，请点此查看更多更全内容
Copyright © 2019- huatuo9.cn 版权所有 赣ICP备2023008801号-1
违法及侵权请联系：TEL:199 18 7713 E-MAIL:2724546146@qq.com
本站由北京市万商天勤律师事务所王兴未律师提供法律服务