【双评价学习笔记】生态—水土保持功能重要性
没想到吧，我又更新了
难道是日更的节奏
今天的内容是生态保护重要性评价中的生态系统服务功能重要性评价中的水土保持功能重要性（请一口气读完）。鉴于之前的文章有质疑我的方法的，我这边申明一下，整个双评价系列笔记，都是在7月版本的双评价指南的基础上的糊弄，因此生态部分也还是沿用7月版的内容。其中，指南中说“关于生态系统服务功能重要性评价，评价方法参见《生态保护红线划定指南》（试行）（环办生态〔2017〕48 号）附录中的模型评估法”。在《生态保护红线划定指南》中，写了两种评价方法，分别是模型法和NPP法，由于模型法需要的数据获取过于困难，我这边就都是用NPP法了（包括生物多样性）。想要了解模型法的小伙伴，可以直接网上去找《生态保护红线划定指南》看看。
第三篇 生态保护重要性评价
1
生态系统服务功能重要性评价
1.1 水源涵养功能重要性评价
1.2 水土保持功能重要性评价
1.2.1 评估模型
以生态系统水土保持服务能力指数作为评估指标，计算公式为：
式中：Spro为水土保持服务能力指数，NPPmean为多年植被净初级生产力平均值， Fslo为坡度因子，K为土壤可蚀性因子。
1.2.2 单因子计算
Fslo（坡度因子）：坡度可以直接使用dem生成，之前的农业适宜性和城镇适宜性的文章里面都有提到，因此这边直接使用之前的数据结果。
NPPmean（多年植被净初级生产力平均值）：NPP可以从网络获取（上一篇文章分享了），在获得NPP数据后可以轻易的计算出平均值。下图为京津冀地区的历年平均NPP。
K（土壤可蚀性因子）：指土壤颗粒被水力分离和搬运的难易程度，主要与土壤质地、有机质含量、土体结构、渗透性等土壤理化性质有关，计算公式如下：
式中，KEPIC表示修正前的土壤可蚀性因子，K 表示修正后的土壤可蚀性因子，mc、msilt、ms和orgC 分别为粘粒（＜0.002mm）、粉粒（0.002 mm～0.05 mm）、砂粒（0.05 mm～2 mm）和有机碳的百分比含量（%）。
别看这两个公式好像很复杂，但其实带入mc、msilt、ms和orgC四个数据字段就可以计算得到最后的结果了。这四个字段可以在上一篇文章中分享的“全球土壤数据集”中找到，字段名分别为：T_CLAY: Real (顶层粘土含量) (单位：% wt.)、T_SILT: Real (表层粉沙粒含量)(单位：% wt.)、T_SAND: Real (顶层沙含量) (单位：% wt.)和T_OC: Real (顶层有机碳含量) (单位：% weight)。
在带入四个数据后，可以轻松得到（读书时参考答案里面的“可以轻易得到”的语句终于可以用了）K值，研究区域的K值计算结果如下图所示：
K值（土壤可蚀性因子）计算结果
1.2.3 因子归一化
将各因子归一化处理，使用最大值最小值法将数据归一化到0-1之间。归一化处理可以使用【模糊隶属度】工具。（上一篇文章中提到了）
1.2.4 模型计算
将归一化的各因子，按照前文提到的计算公式使用【栅格计算器】进行计算，就可以得到生态系统水土保持服务能力指数。
今天没有用到新的数据，需要本篇文章里提到的数据的小伙伴，可以看看前一篇文章。
往期回顾：
土地资源评价（分享土壤质地数据、地形高程数据和研究边界数据）戳下面：
水资源评价、气候评价（分享长时间年值、近几年日值气候数据）戳下面：
环境评价、灾害评价中的干旱部分（分享1950年以来的日值降水、平均气温、最高气温、最低气温数据）戳下面：
灾害评价中的雨涝、高温、大风灾害部分（2010年-2017年的全国风速风向数据）戳下面：
生态评价戳下面：
农业集成评价戳下面：
T_T对上面文章的推倒重来戳下面（分享30m地形高程数据）：
城镇建设适宜性评价从这一篇开始，首先是城镇土地资源评价：
城镇适宜性评价中的水资源评价戳下面：
城镇适宜性评价中的气候评价戳下面(分享湿度和温度数据）：
城镇适宜性评价中的区位优势度评价戳下面（分享路网、铁路、高速出入口、火车站、机场等数据）：
城镇适宜性评价中的环境评价戳下面（其实只有大气环境容量==）：
城镇建设集成评价戳下面：
生态重要性评价中的水源涵养功能重要性评价戳下面：
<END>
整理数据不易
熬夜写文秃头
喜欢这篇文章的话
记得帮我们分享点赞
特别特别喜欢的话
也可以赞赏啦
如需转载文章、投稿或者合作
请直接在公众号后台联系
使用关键词搜索历史文章请点击【阅读原文】