1、3附件生态保护红线划定指南环境保护部 国家发展改革委2017 年 5 月4目录1 1 适用范围适用范围.6 62 2 编制依据编制依据.6 63 3 术语和定义术语和定义.8 84 4 划定原则划定原则.9 94.1 科学性原则.94.2 整体性原则.94.3 协调性原则.104.4 动态性原则.105 5 管控要求管控要求.10106 6 划定工作程序划定工作程序.11116.1 制定工作方案和技术方案.116.2 开展划定工作.116.3 各省（区、市）方案上报与审核.116.4 方案批准与发布.127 7 划定技术流程划定技术流程.12127.1 开展科学评估.127.2 校验划定范围.
2、147.3 确定红线边界.157.4 形成划定成果.1657.5 开展勘界定标.168 8 命名与编码命名与编码.18188.1 命名.188.2 编码.189 9 成果要求成果要求.19199.1 文本.199.2 图件.199.3 登记表.209.4 台账数据库.209.5 技术报告.201010 附则附则.2020附录附录 A A 生态系统服务功能重要性评估方法生态系统服务功能重要性评估方法.2121附录附录 B B 生态环境敏感性评估方法生态环境敏感性评估方法.4242附录附录 C C 生态保护红线综合制图生态保护红线综合制图.5151附录附录 D D 生态保护红线汇总表生态保护红线汇
3、总表.5454附录附录 E E 生态保护红线登记表生态保护红线登记表.5555附录附录 F F 生态保护红线划定技术报告编写大纲生态保护红线划定技术报告编写大纲.56566为贯彻中华人民共和国环境保护法 中共中央关于全面深化改革若干重大问题的决定，落实关于划定并严守生态保护红线的若干意见（以下简称若干意见），指导全国生态保护红线划定工作，保障国家生态安全，制定本指南。1 适用范围本指南适用于中华人民共和国陆地国土空间生态保护红线的划定。2 编制依据中华人民共和国环境保护法中华人民共和国国家安全法中华人民共和国水土保持法中华人民共和国土地管理法中华人民共和国水法中华人民共和国草原法中华人民共和国
4、防沙治沙法中华人民共和国森林法中共中央 国务院关于加快推进生态文明建设的意见（中发201512 号）生态文明体制改革总体方案（中发201525 号）关于划定并严守生态保护红线的若干意见（厅字 2017 2 号）国务院关于印发全国主体功能区规划的通知（国发201046 号）国务院关于印发全国国土规划纲要（20162030 年）的通知（国发20173 号）7国务院办公厅关于印发湿地保护修复制度方案的通知（国办201689 号）国务院关于印发“十三五”生态环境保护规划的通知（国发201665 号）国务院关于全国水土保持规划（20152030 年）的批复（国函2015160 号）国务院关于全国重要江河
5、湖泊水功能区划（20112030 年）的批复（国函2011167 号）国务院关于全国林地保护利用规划纲要（20102020 年）的批复（国函201069 号）国务院关于印发全国土地利用总体规划纲要（2006-2020 年）的通知（国发200833 号）关于印发全国土地利用总体规划纲要（2006-2020 年）调整方案的通知（国土资发201667 号）水利部关于印发全国重要饮用水水源地名录（2016 年）的通知（水资源函2016383 号）农业部关于印发全国草原保护建设利用总体规划的通知（农计发200711 号）关于印发全国生态功能区划（修编版）的公告（环境保护部 中国科学院公告 2015 年第
6、 61 号）8GB/T 2260中华人民共和国行政区划代码GB/T 12343国家基本比例尺地图编绘规范GB/T 13923基础地理信息要素分类与代码GB/T21010-2007土地利用现状分类HJ/T 338饮用水水源保护区划分技术规范SL 190土壤侵蚀分类分级标准CH/T 9005基础地理信息数据库基本规定3 术语和定义生态保护红线：指在生态空间范围内具有特殊重要生态功能、必须强制性严格保护的区域，是保障和维护国家生态安全的底线和生命线，通常包括具有重要水源涵养、生物多样性维护、水土保持、防风固沙、海岸生态稳定等功能的生态功能重要区域，以及水土流失、土地沙化、石漠化、盐渍化等生态环境敏感
7、脆弱区域。国土空间：指国家主权与主权权利管辖下的地域空间，是国民生存的场所和环境，包括陆地、陆上水域、内水、领海、领空等。生态空间：指具有自然属性、以提供生态服务或生态产品为主体功能的国土空间，包括森林、草原、湿地、河流、湖泊、滩涂、岸线、海洋、荒地、荒漠、戈壁、冰川、高山冻原、无居民海岛等。重点生态功能区：指生态系统十分重要，关系全国或区域生态安全，需要在国土空间开发中限制进行大规模高强度工业化城镇化开发，以保持并提高生态产品供给能力的区域，主要类型包括水源涵养区、水土保持区、防风固沙区和生物多样性维护区。生态环境敏感脆弱区：指生态系统稳定性差，容易受到外界活动影响而产生生态退化且难以自我修
8、复的区域。9禁止开发区域：指依法设立的各级各类自然文化资源保护区域，以及其他禁止进行工业化城镇化开发、需要特殊保护的重点生态功能区。生态安全：指在国家或区域尺度上，生态系统结构合理、功能完善、格局稳定，并能够为人类生存和经济社会发展持续提供生态服务的状态，是国家安全的重要组成部分。生态安全格局：指由事关国家和区域生态安全的关键性保护地构成的结构完整、功能完备、分布连续的生态空间布局。勘界定标：指对已划定的生态保护红线边界进行实地勘查、测绘，核准拐点坐标，勘定精确界线，设立统一规范的界碑界桩和标识标牌的行为。4 划定原则4.14.1 科学性原则科学性原则以构建国家生态安全格局为目标，采取定量评估
9、与定性判定相结合的方法划定生态保护红线。在资源环境承载能力和国土空间开发适宜性评价的基础上，按生态系统服务功能（以下简称生态功能）重要性、生态环境敏感性识别生态保护红线范围，并落实到国土空间，确保生态保护红线布局合理、落地准确、边界清晰。4.24.2 整体性整体性原则原则统筹考虑自然生态整体性和系统性，结合山脉、河流、地貌单元、植被等自然边界以及生态廊道的连通性，合理划定生态保护红线，应划尽划，避免生境破碎化，加强跨区域间生态保护红线的有序衔接。104.34.3 协调性原则协调性原则建立协调有序的生态保护红线划定工作机制，强化部门联动，上下结合，充分与主体功能区规划、生态功能区划、水功能区划及
10、土地利用现状、城乡发展布局、国家应对气候变化规划等相衔接，与永久基本农田保护红线和城镇开发边界相协调，与经济社会发展需求和当前监管能力相适应，统筹划定生态保护红线。4.44.4 动态性原则动态性原则根据构建国家和区域生态安全格局，提升生态保护能力和生态系统完整性的需要，生态保护红线布局应不断优化和完善，面积只增不减。5 管控要求生态保护红线原则上按禁止开发区域的要求进行管理。严禁不符合主体功能定位的各类开发活动，严禁任意改变用途，确保生态功能不降低、面积不减少、性质不改变。因国家重大基础设施、重大民生保障项目建设等需要调整的，由省级政府组织论证，提出调整方案，经环境保护部、国家发展改革委会同有
11、关部门提出审核意见后，报国务院批准。功能不降低。生态保护红线内的自然生态系统结构保持相对稳定，退化生态系统功能不断改善，质量不断提升。面积不减少。生态保护红线边界保持相对固定，生态保护红线面积只能增加，不能减少。性质不改变。严格实施生态保护红线国土空间用途管制，严禁随意改变用地性质。116 划定工作程序采取自上而下和自下而上相结合的方式划定全国和各省（区、市）生态保护红线。6.16.1 制定工作方案和技术方案制定工作方案和技术方案各省（区、市）依照本指南和生态保护红线划定的总体要求，建立划定工作责任制和协调机制，制定各省（区、市）工作方案和技术方案，明确职责分工，组织专门队伍，有序推进划定工作
12、。6.26.2 开展划定工作开展划定工作环境保护部、国家发展改革委会同有关部门开展国家生态保护红线顶层设计，提出各省(区、市)生态保护红线空间格局和分布建议方案，明确需要保护的湿地、草原、森林等生态系统分布范围，指导各地生态保护红线划定，并做好跨省域生态保护红线的衔接与协调。各省(区、市)依据工作方案和技术方案组织开展划定工作，参照国家生态保护红线空间格局和分布建议方案，结合本地实际情况，形成本行政区生态保护红线划定初步方案（含文本、图件、登记表），征求相关部门和地方政府意见，开展专家论证。经修改完善报省（区、市）人民政府审议同意后，形成生态保护红线划定方案（送审稿）。6.36.3 各省（区、
13、市）方案上报与审核各省（区、市）方案上报与审核各省（区、市）人民政府将生态保护红线划定方案（送审稿）报送环境保护部、国家发展改革委，环境保护部、国家发展改革委会同有关部门组织对各省（区、市）生态保护红线开展技术审核并12提出意见，各省（区、市）人民政府根据意见修改完善后，形成生态保护红线划定方案（报批稿）。6.46.4 方案批准与发布方案批准与发布环境保护部、国家发展改革委会同有关部门将各省（区、市）生态保护红线划定方案报国务院审批后，由各省（区、市）人民政府发布实施。环境保护部、国家发展改革委会同有关部门在各省（区、市）生态保护红线划定方案基础上进行汇总，形成全国生态保护红线划定方案，报国务
14、院同意后，向社会发布。已初步完成或发布实施生态保护红线划定方案的有关省（区、市）按照若干意见和本指南要求，进一步优化调整划定方案，按程序上报审批。7 划定技术流程按照定量与定性相结合的原则，通过科学评估，识别生态保护的重点类型和重要区域，合理划定生态保护红线。7.17.1 开展科学评估开展科学评估在国土空间范围内，按照资源环境承载能力和国土空间开发适宜性评价技术方法，开展生态功能重要性评估和生态环境敏感性评估，确定水源涵养、生物多样性维护、水土保持、防风固沙等生态功能极重要区域及极敏感区域，纳入生态保护红线。科学评估的主要步骤包括：确定基本评估单元、选择评估类型与方法、数据准备、模型运算、评估
15、分级和现场校验。（1 1）确定基本评估单元）确定基本评估单元根据生态评估参数的数据可获取性，统一评估工作精度要求。13原则上评估的基本空间单元应为 250m250m 网格，有条件的地区可进一步提高精度。评估工作运行环境采用地理信息系统软件。（2 2）选择评估类型与方法）选择评估类型与方法根据本地区生态环境特征和主要生态问题，确定生态功能和生态环境敏感性类型，并结合数据条件，选取适宜的评估方法（参见附录 A 和附录 B）。（3 3）数据）数据准备准备根据评估方法，搜集评估所需的各类数据，如基础地理信息数据、土地利用现状及年度调查监测数据、气象观测数据、遥感影像、地表参量、生态系统类型与分布数据等
16、。评估的基础数据类型为栅格数据，非栅格数据应进行预处理，统一转换为便于空间计算的网格化栅格数据。（4 4）模型运算）模型运算根据评估公式，在地理信息系统软件中输入评估所需的各项参数，计算生态系统服务功能重要性和生态环境敏感性指数。（5 5）评估分级）评估分级根据评估结果，将生态功能重要性依次划分为一般重要、重要和极重要 3 个等级，将生态环境敏感性依次划分为一般敏感、敏感和极敏感 3 个等级。（6 6）现场校核）现场校核根据相关规划、区划中重要生态区域空间分布，结合专家知识，综合判断评估结果与实际生态状况的相符性。针对不符合实际情况的评估结果开展现场核查校验与调整，使评估结果趋于合理。147.
17、27.2 校验划定范围校验划定范围根据科学评估结果，将评估得到的生态功能极重要区和生态环境极敏感区进行叠加合并，并与以下保护地进行校验，形成生态保护红线空间叠加图，确保划定范围涵盖国家级和省级禁止开发区域，以及其他有必要严格保护的各类保护地。（1 1）国家级和省级禁止开发区域）国家级和省级禁止开发区域国家公园；自然保护区；森林公园的生态保育区和核心景观区；风景名胜区的核心景区；地质公园的地质遗迹保护区；世界自然遗产的核心区和缓冲区；湿地公园的湿地保育区和恢复重建区；饮用水水源地的一级保护区；水产种质资源保护区的核心区；其他类型禁止开发区的核心保护区域。对于上述禁止开发区域内的不同功能分区，应根
18、据生态评估结果最终确定纳入生态保护红线的具体范围。位于生态空间以外或人文景观类的禁止开发区域，不纳入生态保护红线。（2 2）其他各类保护地）其他各类保护地除上述禁止开发区域以外，各地可结合实际情况，根据生态功能重要性，将有必要实施严格保护的各类保护地纳入生态保护红线范围。主要涵盖：极小种群物种分布的栖息地、国家一级公益林、15重要湿地（含滨海湿地）、国家级水土流失重点预防区、沙化土地封禁保护区、野生植物集中分布地、自然岸线、雪山冰川、高原冻土等重要生态保护地。7.7.3 3 确定红线边界确定红线边界将 7.2 确定的生态保护红线叠加图，通过边界处理、现状与规划衔接、跨区域协调、上下对接等步骤，
19、确定生态保护红线边界。（1 1）边界处理）边界处理采用地理信息系统软件，对叠加图层进行图斑聚合处理，合理扣除独立细小斑块和建设用地、基本农田（综合制图方法参见附录C）。边界调整的底图建议采用第一次全国地理普查数据库或土地利用现状及年度调查监测成果，按照保护需要和开发利用现状，结合以下几类界线勾绘调整生态保护红线边界：自然边界，主要是依据地形地貌或生态系统完整性确定的边界，如林线、雪线、流域分界线，以及生态系统分布界线等；自然保护区、风景名胜区等各类保护地边界；江河、湖库，以及海岸等向陆域（或向海）延伸一定距离的边界；地理国情普查、全国土地调查、森林草原湿地荒漠等自然资源调查等明确的地块边界。（
20、2 2）现状与规划衔接）现状与规划衔接将生态保护红线边界与各类规划、区划空间边界及土地利用现状相衔接，综合分析开发建设与生态保护的关系，结合经济社会发展实际，合理确定开发与保护边界，提高生态保护红线划定16合理性和可行性。（3 3）跨区域协调）跨区域协调根据生态安全格局构建需要，综合考虑区域或流域生态系统完整性，以地形、地貌、植被、河流水系等自然界线为依据，充分与相邻行政区域生态保护红线划定结果进行衔接与协调，开展跨区域技术对接，确保生态保护红线空间连续，实现跨区域生态系统整体保护。（4 4）上下对接）上下对接采取上下结合的方式开展技术对接，广泛征求各市县级政府意见，修改完善后达成一致意见，确
21、定生态保护红线边界。7.47.4 形成划定成果形成划定成果在上述工作基础上，编制生态保护红线划定文本、图件、登记表及技术报告，建立台账数据库，形成生态保护红线划定方案。7.57.5 开展勘界定标开展勘界定标根据划定方案确定的生态保护红线分布图，搜集红线附近原有平面控制点坐标成果、控制点网图，以高清正射影像图、地形图和地籍图等相关资料为辅助，调查生态保护红线各类基础信息，明确红线区块边界走向和实地拐点坐标，详细勘定红线边界。选定界桩位置，完成界桩埋设，测定界桩精确空间坐标，建立界桩数据库，形成生态保护红线勘测定界图。设立统一规范的标识标牌，主要内容包括生态保护红线区块的范围、面积、具体拐点坐标、
22、保护对象、主导生态功能、主要管控措施、责任人、监督管理电话等。生态保护红线划定技术流程参见图 1。17图 1生态保护红线划定技术流程水源涵养生物多样性维护水土保持防风固沙生态功能重要性评估生态保护红线勘界定标生态保护红线勘界定标勘查与测定埋设界桩界碑设立标识标牌形成划定成果形成划定成果文本图件登记表技术报告台账数据库确定红线边界确定红线边界现状与规划衔接跨区域协调上下对接校验划定范围校验划定范围禁止开发区域其他各类保护地生态功能极重要区生态环境极敏感区水土流失盐渍化土地沙化石漠化生态环境敏感性评估科学科学评估评估空间叠加边界处理188 命名与编码8.18.1 命名命名生态保护红线命名采取“自然
23、地理单元+主导生态功能+生态保护红线”的命名方式，如“燕山水源涵养生态保护红线”。落到具体地块，增加生态保护红线所属县级行政区，即“县级行政区+自然地理单元+主导生态功能（或生态环境敏感性）+生态保护红线”的命名方式，如“密云区密云水库水源涵养生态保护红线”，以便以县域为基本单元建立生态保护红线台账系统。8.28.2 编码编码为加强生态保护红线信息化管理，对生态保护红线实行统一编码，采用“行政代码-类型代码-数量代码”的三级编码方式(见表 1)。（1）行政代码以县级行政区为单位，由 6 位阿拉伯数字组成。（2）类型代码由 2 位数字组成，第 1 位表示类型特征，其中，1表示生态功能，2 表示生
24、态环境敏感性。后 1 位表示属性分类，其中，生态功能包括：1-水源涵养，2-生物多样性维护，3-水土保持，4-防风固沙，5-其他生态功能。生态环境敏感性包括：1-水土流失，2-土地沙化，3-石漠化，4-盐渍化，5-其他敏感性。（3）数量代码表示某一类型生态保护红线的地块序号，从 001开始编号。19表表 1 1生态保护红线编码方式生态保护红线编码方式行政代码类型代码数量代码一级编码名称二级编码名称0010021生态功能1水源涵养2生物多样性维护3水土保持4防风固沙5其他生态功能2生态环境敏感性1水土流失2土地沙化3石漠化4盐渍化5其他敏感性9 成果要求生态保护红线划定成果包括文本、图件、登记表
25、、台账数据库、技术报告等。9.19.1 文本文本生态保护红线划定文本内容主要包括：划定生态保护红线的重要性和必要性，指导思想、基本原则和目标，生态保护红线类型与概述，管控措施，生态保护红线汇总表（格式参见附录 D）。9.29.2 图件图件生态保护红线图件数据采用 2000 国家大地坐标系统，高斯-克吕格投影，1985 国家高程基准；省级层面基本比例尺为 1：5 万，县级层面基本比例尺原则上不小于 1：1 万，基础数据不满足要求的可采用 1：5 万。生态保护红线图件应包括但不限于：20（1）生态保护红线分布图（2）生态功能重要性评估系列图（3）生态环境敏感性评估系列图（4）禁止开发区域分布图（5
26、）其他保护地分布图9.39.3 登记表登记表以县级行政区为基本单元，编制生态保护红线登记表。登记表内容主要包括红线区块代码、名称、类型、地理位置、面积、人口数量、生态功能、主要生态环境问题、主要人类活动、管控措施等基本信息（登记表格式参见附录 E）。9.49.4 台账数据库台账数据库以县级行政区为基本单元，构建生态保护红线台账数据库。台账信息主要包括红线区块登记表信息、基础地理信息、气象观测要素、社会经济要素、地面监测要素、遥感影像、地表生态参数、土地权属等。9.59.5 技术报告技术报告以文字报告形式表述生态保护红线划定的主要内容（格式参见附录 F）。10 附则本指南由环境保护部、国家发展改
27、革委负责解释。本指南自发布之日起实施，生态保护红线划定技术指南（环发201556 号）同时废止。21附录 A生态系统服务功能重要性评估方法目前，生态系统服务功能采用的评估方法主要有模型评估法和净初级生产力（NPP）定量指标评估法。其中，模型评估法所需参数较多，对数据需求量较大，准确度较高；定量指标法以 NPP 数据为主，参数较少，操作较为简单，但其适用范围具有地域性。为提高评估结论的准确性以及与实地的相符性，评估方法的参数选取可在评估过程进行适当调整和细化，尽可能采用国内权威的、分辨率更高的基础数据。评估结果还需根据实地观测、调查结果进一步校验。对于全国和各省生态保护红线划定，可使用 NPP
28、定量指标法、模型法及其他常用评估方法。鉴于国家发展改革委在资源环境承载力评估中使用的方法为模型法，为保持评估结果的一致性，建议各地优先使用模型法。A.1A.1 模型评估法模型评估法A.1.1A.1.1 水源涵养功能重要性评估水源涵养功能重要性评估水源涵养是生态系统（如森林、草地等）通过其特有的结构与水相互作用，对降水进行截留、渗透、蓄积，并通过蒸散发实现对水流、水循环的调控，主要表现在缓和地表径流、补充地下水、减缓河流流量的季节波动、滞洪补枯、保证水质等方面。以水源涵养量作为生态系统水源涵养功能的评估指标。A.1.1.1A.1.1.1 评估模型评估模型22采用水量平衡方程来计算水源涵养量，计算
29、公式为：3110jiiiiiTQP RETA式中：TQ为总水源涵养量（m3），Pi为降雨量（mm），Ri为地表径流量（mm），ETi为蒸散发（mm），Ai为i类生态系统面积（km2），i为研究区第i类生态系统类型，j为研究区生态系统类型数。A.1.1.2A.1.1.2 数据准备数据准备（1）数据来源与获取根据上述模型，水源涵养功能重要性评估需收集生态系统类型数据集、气象数据集和蒸散发数据集等，具体信息见表 A1。表表 A1A1水源涵养功能重要性评估数据表水源涵养功能重要性评估数据表名称类型分辨率数据来源生态系统类型数据集矢量-全国生态状况遥感调查与评估成果气象数据集文本-中国气象科学数据共享服
30、务网蒸散发数据集栅格1km国家生态系统观测研究网络科技资源服务系统网站（2）数据预处理降雨量因子：降雨量因子：根据气象数据集处理得到。在 Excel 中计算出区域所有气象站点的多年平均降水量，将这些值根据相同的站点名与ArcGIS 中的站点（点图层）数据相连接（Join）。在 Spatial Analyst工具中选择 Interpolate to Raster 选项，选择相应的插值方法得到降水量因子栅格图。地表径流因子地表径流因子：降雨量乘以地表径流系数获得，计算公式如下：R P 式中：R为地表径流量（mm），P为多年平均降雨量（mm），为平均地表径流系数，如表 A2 所示。23表表 A2A2
31、各类型生态系统地表径流系数均值表各类型生态系统地表径流系数均值表生态系统类型 1生态系统类型 2平均地表径流系数（%）森林常绿阔叶林2.67常绿针叶林3.02针阔混交林2.29落叶阔叶林1.33落叶针叶林0.88稀疏林19.20灌丛常绿阔叶灌丛4.26落叶阔叶灌丛4.17针叶灌丛4.17稀疏灌丛19.20草地草甸8.20草原4.78草丛9.37稀疏草地18.27湿地湿地0.00蒸散发因子：蒸散发因子：根据国家生态系统观测研究网络科技资源服务系统网站提供的产品数据。原始数据空间分辨率为 1km，通过 ArcGIS软件重采样为 250m 空间分辨率，得到蒸散发因子栅格图。生态系统面积因子：生态系统
32、面积因子：根据全国生态状况遥感调查与评估成果中的生态系统类型数据集得到。原始数据为矢量数据，通过 ArcGIS 软件转为 250m 空间分辨率的栅格图。A.1.1.3A.1.1.3 模型运算模型运算将各因子统一成 250m 分辨率的栅格数据，在 ArcGIS 栅格计算器（Spatial AnalystRaster Calculator）中，根据公式计算得到生态系统水源涵养量。24A.1.2A.1.2 水土保持功能重要性评估水土保持功能重要性评估水土保持是生态系统（如森林、草地等）通过其结构与过程减少由于水蚀所导致的土壤侵蚀的作用，是生态系统提供的重要调节服务之一。水土保持功能主要与气候、土壤、
33、地形和植被有关。以水土保持量，即潜在土壤侵蚀量与实际土壤侵蚀量的差值，作为生态系统水土保持功能的评估指标。A.1.2.1A.1.2.1 评估模型评估模型采用修正通用水土流失方程（RUSLE）的水土保持服务模型开展评价，公式如下：1cprAAARKLSC式中，Ac为水土保持量（t/hm2a）；Ap为潜在土壤侵蚀量；Ar为实际土壤侵蚀量；R为降雨侵蚀力因子（MJmm/hm2ha）；K为土壤可蚀性因子（thm2h/hm2MJmm）；L、S为地形因子，L表示坡长因子，S表示坡度因子；C为植被覆盖因子。A.1.2.2A.1.2.2 数据准备数据准备（1）数据来源与获取根据上述模型，水土保持功能重要性评估
34、需收集高程数据集、气象数据集和土壤数据集等数据，具体信息见表 A3。25表表 A3A3水土保持功能重要性评估数据表水土保持功能重要性评估数据表名称类型分辨率数据来源高程数据集栅格30m地理空间数据云网站气象数据集文本中国气象科学数据共享服务网土壤数据集矢量/Excel全国生态环境调查数据库中国 1：100 万土壤数据库（2）数据预处理降雨侵蚀降雨侵蚀力力因子因子R R：是指降雨引发土壤侵蚀的潜在能力，通过多年平均年降雨侵蚀力因子反映，计算公式如下：241kkRR半月式中，R为多年平均年降雨侵蚀力（MJmm/hm2ha）；R半月k为第k个半月的降雨侵蚀力（MJmm/hm2ha）；k为一年的 24
35、 个半月，k1,2,24；i为所用降雨资料的年份，i1,2,n；j为第i年第k个半月侵蚀性降雨日的天数，j1,2,m；Pi,j,k为第i年第k个半月第j个侵蚀性日降雨量（mm），可以根据全国范围内气象站点多年的逐日降雨量资料，通过插值获得；或者直接采用国家气象局的逐日降雨量数据产品。为参数，暖季时0.3937，冷季时0.3101。土壤可蚀性因子土壤可蚀性因子K K：指土壤颗粒被水力分离和搬运的难易程度，主要与土壤质地、有机质含量、土体结构、渗透性等土壤理化性质有关，计算公式如下：_260 013830 515750.1317EPICK-.K0.30 20 3exp 0 0256(1/100)/
36、()1 0.25/exp(3.722.95)1 0.7(1/100)/(1/100)exp 5.512 2.9(1/100)ssiltEPICsiltcsiltsssK.-.m-mmmmorgCorgCorgCmmm式中，KEPIC表示修正前的土壤可蚀性因子，K表示修正后的土壤可蚀性因子，mc、msilt、ms和orgC分别为粘粒（0.002mm）、粉粒（0.002mm0.05 mm）、砂粒（0.05 mm2 mm）和有机碳的百分比含量（%），数据来源于中国 1：100 万土壤数据库。在 Excel 表格中，利用上述公式计算K值，然后以土壤类型图为工作底图，在 ArcGIS 中将K值连接（Jo
37、in）到底图上。利用 Conversion Tools 中矢量转栅格工具，转换成空间分辨率为 250m 的土壤可蚀性因子栅格图。地形因子地形因子L L、S S：L表示坡长因子，S表示坡度因子，是反映地形对土壤侵蚀影响的两个因子。在评估中，可以应用地形起伏度，即地面一定距离范围内最大高差，作为区域土壤侵蚀评估的地形指标。选 择 高 程 数 据 集，在 Spatial Analyst 下 使 用 NeighborhoodStatistics，设置 Statistic Type 为最大值和最小值，即得到高程数据集的最大值和最小值，然后在 Spatial Analyst 下使用栅格计算器 Raster
38、 Calculator，公式为最大值-最小值，获取地形起伏度，即地形因子栅格图。植被覆盖因子植被覆盖因子C C：反映了生态系统对土壤侵蚀的影响，是控制土壤侵蚀的积极因素。水田、湿地、城镇和荒漠参照 N-SPECT 的参数分别赋值为 0、0、0.01 和 0.7，旱地按植被覆盖度换算,计算公式如下：271=0.221-0.595logCc旱式中，C旱为旱地的植被覆盖因子，c1为小数形式的植被覆盖度。其余生态系统类型按不同植被覆盖度进行赋值，如表 A4 所示。表表 A4A4不同生态系统类型植被覆盖因子赋值不同生态系统类型植被覆盖因子赋值生态系统类型植被覆盖度90森林0.10.080.060.020
39、.0040.001灌丛0.40.220.140.0850.040.011草地0.450.240.150.090.0430.011乔木园地0.420.230.140.0890.0420.011灌木园地0.40.220.140.0870.0420.011A.1.2.3A.1.2.3 模型运算模型运算将各因子统一成 250m 分辨率的栅格数据，在 ArcGIS 栅格计算器（Spatial AnalystRaster Calculator）中，根据公式计算得到生态系统水土保持量。A.1.3A.1.3 防风固沙功能重要性评估防风固沙功能重要性评估防风固沙是生态系统（如森林、草地等）通过其结构与过程减少由
40、于风蚀所导致的土壤侵蚀的作用，是生态系统提供的重要调节服务之一。防风固沙功能主要与风速、降雨、温度、土壤、地形和植被等因素密切相关。以防风固沙量（潜在风蚀量与实际风蚀量的差值）作为生态系统防风固沙功能的评估指标。A.1.3.1A.1.3.1 评估模型评估模型采用修正风蚀方程来计算防风固沙量，公式如下：=-LLSR SS潜2820.3711/22=109.8=150.71z sLMAXMAXzSQeSQWFEFSCFKSWFEFSCFK潜潜潜潜潜潜式中，SR为固沙量（t km-2a-1）；SL潜为潜在风力侵蚀量（t km-2a-1）；SL为实际风力侵蚀量（t km-2a-1）；QMAX为最大转移
41、量（kg/m)；Z为最大风蚀出现距离（m)；WF为气候因子（kg/m）；K为地表糙度因子；EF为土壤可蚀因子；SCF为土壤结皮因子；C为植被覆盖因子。A.1.3.2A.1.3.2 数据准备数据准备（1）数据来源与获取根据上述评估模型，防风固沙功能重要性评估需用到遥感数据集、高程数据集、气象数据集和土壤数据集等数据，具体信息见表A5。表表 A A5 5防风固沙功能重要性评估数据表防风固沙功能重要性评估数据表名称类型分辨率数据来源遥感数据集栅格250m美国国家航空航天局（NASA）网站或地理空间数据云网站高程数据集栅格30m地理空间数据云网站气象数据集文本中国气象科学数据共享服务网土壤数据集矢量/
42、Excel全国生态环境调查数据库中国 1：100 万土壤数据库中国地区 Modis 雪盖产品数据集栅格0.05 度寒区旱区科学数据中心2(/)20.3711max2150.71()109.8z sLMAXzSQeSSWFEFSCFKCQWFEFSCFKC29（2）数据预处理气候因子气候因子WFWFSDSWgWfWF式中，WF为气候因子，单位为 kg/m，12 个月WF总和得到多年年均WF；Wf为各月多年平均风力因子，为空气密度，g为重力加速度；在 Excel 中计算出区域所有气象站点的多年平均风力，将这些值根据相同的站点名与 ArcGIS 中的站点（点图层）数据相连接（Join）。在 Spat
43、ial Analyst 工具中选择 Interpolate to Raster选项，选择相应的插值方法得到各月多年平均风力因子栅格图。SW为各月多年平均土壤湿度因子，无量纲；SD为雪盖因子，无量纲。雪盖数据来源于寒区旱区科学数据中心的中国地区 Modis 雪盖产品数据集。土壤可蚀因子土壤可蚀因子EFEF329.090.310.170.33/2.590.95100sasisa clOMCacoEF式中，sa为土壤粗砂含量（0.2 mm2 mm）（%）；si为土壤粉砂含量（%）；cl 为土壤粘粒含量（%）；OM为土壤有机质含量（%）；Caco3为碳酸钙含量（%），可不予考虑。土壤结皮因子土壤结皮因
44、子SCFSCF22)(021.0)(0066.011OMclSCF式中，cl 为土壤粘粒含量（%）；OM为土壤有机质含量（%）。植被覆盖因子植被覆盖因子C C30不同植被类型的防风固沙效果不同，研究将植被分为林地、灌丛、草地、农田、裸地和沙漠六个植被类型，根据不同的系数计算各植被覆盖因子C值：式中，SC为植被覆盖度，计算公式见 A.1.2.2；ai为不同植被类型的系数，分别为：林地 0.1535，草地 0.1151，灌丛 0.0921，裸地0.0768，沙地 0.0658，农田 0.0438。地表糙度因子地表糙度因子KK20.2rHKL式中，Kr为土垄糙度，以 Smith-Carson 方程加
45、以计算，单位 cm；Crr为随机糙度因子，取 0，单位 cm；L为地势起伏参数；H为距离L范围内的海拔高程差，在GIS软件中使用Neighborhood statistics工具计算 DEM 数据相邻单元格地形起伏差值获得。A.1.3.3A.1.3.3 模型运算模型运算将各因子统一成 250m 分辨率的栅格数据，在 ArcGIS 栅格计算器（Spatial AnalystRaster Calculator）中，根据公式计算得到生态系统防风固沙量。A.1.4A.1.4 生物多样性维护功能重要性评估生物多样性维护功能重要性评估生物多样性维护功能是生态系统在维持基因、物种、生态系统多样性发挥的作用，
46、是生态系统提供的最主要功能之一。生物多样性维护功能与珍稀濒危和特有动植物的分布丰富程度密切相关，主a()iSCCe)127.041.286.1(934.0CrrKKrreK31要以国家一、二级保护物种和其他具有重要保护价值的物种（含旗舰物种）作为生物多样性保护功能的评估指标。A.1.4.1A.1.4.1 评估模型评估模型以国家一、二级保护物种和其他具有重要保护价值的物种为保护目标，全面收集区域动植物多样性和环境资源数据，建立物种分布数据库。根据关键物种分布点的环境信息和背景信息，应用物种分布模型（Species Distribution Models，SDMs）量化物种对环境的依赖关系，从而预
47、测任何一点某物种分布的概率，结合关键物种的实际分布范围最终划定确保物种长期存活的保护红线。A.1.4.2A.1.4.2 数据准备数据准备（1）物种分布数据库物种分布数据库是以物种名、经纬度和调查时间为核心信息的数据库。物种名分为中文名和拉丁名两个字段；经纬度字段以度为单位，保留小数点后 5 位数字，并记录数据精度，野外调查中通过GPS 记录的数据精度一般在十几米，根据地名信息（县名、乡镇名、河流、山脉等）匹配的经纬度精度一般在几公里至几十公里；时间字段为年（如 1998，2005 等），记录日期与时间、调查的时间阶段等信息。此外，可增加备注字段（memo），记录任何相关信息（该字段不限长度）。
48、各数据来源（调查人、文献等）应记录在数据库中。数据库软件可采用 MS Access。（2）利用物种分布模型预测物种的分布建立物种分布的数据库后，配合环境变量，即可应用物种分布模型预测物种的分布。32环境变量每个环境变量以 GIS 图层的形式储存和管理，范围覆盖整个中国。对于空间范围超过 1000 万平方公里的分析，精度为 1 平方公里的栅格数据是最常用的。常用的环境变量有：地形地貌变量：海拔、坡度、坡向；地表类型变量：土地覆被类型、植被类型、土壤类型；气候变量：年均温、年均降水量的季节变异、1 月最低温、1 月均温、1 月降水量、7 月最高温、7 月均温、7 月降水量、年均温度变化范围、干燥度
49、、辐射强度；生态指标：植被净初级生产力、NDVI、土层厚度、土壤氮含量、土壤碳含量等;人文指标：GDP、人口密度、道路密度、乡镇密度、河流密度等。物种分布模型常用的物种分布模型主要包括回归模型、分类树和混合大量简单模型的神经网络、随机森林等。其中逻辑斯蒂回归是最为简单、应用最广的模型。机器学习类复杂模型（如随机森林、神经网络、Maxent 等）的预测精度较高，近年来应用较多。逻辑斯蒂回归是广义线性模型（Generalized Linear Models）的一种形式。广义线性模型是一般线性模型的扩展，允许因变量为二项分布、泊松分布等离散型的分布；而一般线性模型要求因变量为连续变量，而且其残差为正
50、态分布。大多数统计软件在运行广义线性模型时，都可以应用逐步回归，依据模型拟合优度的统计量33Akaike Information Criterioin(AIC)或 Bayesian InformationCriteria(BIC)去除冗余变量，减少共线性。随机森林（Random Forest）应用 Breiman 的随机森林算法，通过对大量的分类树的计算来进行分类和回归。随机森林把一组解释变量的值（一个向量）输入森林中的每棵分类树中，每棵树都给出这个向量的分类结果（例如物种存在还是不存在）。随机森林对所有的树的分类结果进行打分，并选择得分最高的分类树。整个算法包括树的生长阶段和投票选择阶段。随
