生态系统服务功能模块之水源涵养计算(一)

嘿嘿灵活的小胖子

已于 2024-06-01 22:21:45 修改

阅读量1.1k

点赞数 23

文章标签： arcgis

于 2024-05-10 19:50:59 首次发布

本文链接：https://blog.csdn.net/m0_53209421/article/details/138678549

版权

前言：（碎碎念：啊啊啊，老早就想通过CSDN写一些东西了，但是本人实在有点懒，坚持不住啊！哭唧唧）最近在研究关于生态系统服务功能方面的内容，借此平台记录一下学习过程。今天学习生态系统服务功能的水源涵养功能，下面让我们开始吧！

水源涵养的定义：生态系统通过其结构和过程拦截滞蓄降水，增强土壤下渗，涵养土壤水分和补充地下水，调节河川流量，增加可利用水资源量的功能。

水源涵养量的计算公式为： 𝑊𝑅=𝑀𝑖𝑛(249/𝑉,1)∗𝑀𝑖𝑛(1,0.9∗𝑇1/3)∗𝑀𝑖𝑛(1,𝐾/300)∗𝑌

其中：WR为水源涵养量（mm）；V为流速系数；K为土壤饱和导水率；Y为产水量；T1为地形指数。

要想计算水源涵养量，需要确定V、K、Y、T1，下面分别对它们是如何确定的进行说明。

①流速系数V：视植被类型而定，通过文献，水田、水体、建设用地均取2012，林地、草地、岩石及裸地分别取200、500、1500。

②土壤饱和导水率K： $K=1.14810^{(-0.6+1.260.01C2-6.40.001*C1)}$

C1、C2分别为土壤黏粒与砂粒含量（%）。对于土壤黏粒与砂粒含量可以通过下载土壤数据获得（可以私信我），后文也会放链接。

③地形因子T1：通过阅读文献发现， 𝑇1=𝑙𝑛(𝑎/𝑡𝑎𝑛𝑏)

a为单位等高线长度的上坡集水区面积，b为局部地段坡度。用arcgis可以计算。

④产水量Y，对于产水量Y的计算通过Invest模型计算，我们需要收集的数据通过下图可以看到（碎碎念：数据太难收集了，孩子差点崩溃，/(ㄒoㄒ)/~~）。我分别对其进行说明。

备注：1、所有空间输入必须具有完全相同的投影坐标系 (线性米单位), 而不是地理坐标系 (以度为单位)。

2、栅格输入可能具有不同的像元大小，并且将对其进行重采样以匹配土地利用/土地覆被栅格的像元大小。因此，所有模型结果的像元大小将与土地利用/土地覆被栅格相同。

3、数据单位是米、厘米之类的都需转化为毫米！！！

①Workspace ：工作空间，建个放结果的文件夹。

②文件后缀：虽然是可选，但是建议写上。假设你要做2020年的某个地区的水源涵养量，就可以在这个地方填2020。

③年降水量：年累积降水量，数据下载链接：中国1km分辨率逐月降水量数据集（1901-2021）- 时空三极环境大数据平台 (tpdc.ac.cn)。

注意！！！这个链接下载的数据单位是0.1mm（需转化为mm）, 数据格式是nc格式，第一次处理nc格式的数据，需要将nc格式转化为tif，两种方法：一种通过arcgis处理，比较繁琐复杂；另外一种通过编程批量处理，这里建议采用python编程批量处理（需要代码的可以私信我哦）。

④蒸散发数据：放个链接：http://data.tpdc.ac.cn/zh-hans/data/8b11da09-1a40-4014-bd3d-2b86e6dccad4/?q= 单位：0.1mm（需转化为mm）

⑤植被根系限制层数据：可用基岩深度代替。

⑥土层内植被可用水量（PAWC） （取值范围0-1）

式中：sand、sily、clay、OM分别表示土壤砂粒、粉粒、黏粒、有机质的含量，单位为%。通过土壤数据获得。基于世界土壤数据库（HWSD）的中国土壤数据集(v1.1) - 国家冰川冻土沙漠科学数据中心 (ncdc.ac.cn)

注意：土壤有机质通过土壤有机碳计算。土壤有机质=土壤有机碳*1.724。

⑦土壤利用数据：我通过GEE下载的。Modies数据集500米分辨率的（我做黄河流域大尺度的）。

⑧植被生理系数表：LULC生物物理参数表。LULC栅格中的所有值都必须在该表中有相应的条目。

1.lucode (integer, required): 你所使用的土地利用数据的每一类土地利用类型的值代码。
2.lulc_veg (integer, required): 表示LULC等级是否为植被的代码，用于AET的目的。除湿地外，所有植被等级均为1，所有其他等级均为0，包括湿地、城市区域、水体等.

3.root_depth (number, units: mm, required): 该LULC的植物最大根系深度。只用于'lulc_veg'值为1的类。
这通常被给出为一种植被类型95%的根系生物量发生的深度。对于不使用一般Budyko曲线的土地用途，不需要根深。在这些情况下，根深度字段将被忽略，并可能被设置为一个值，例如-1，以指示不使用该字段。

4.kc (number, units: unitless, required): LULC的作物系数。用于计算潜在蒸散发，以修正参考蒸散发。

⑨Z参数：通过研究区的水资源公报调整。直接浏览器搜索研究区的政府网站，下载相应的水资源公报。看下面这张图片产水量/降雨量=0.505，0.505就是产水系数，所以计算出来的结果也要保证产水系数接近于0.505，如果invest结果计算出来的产水系数为0.3，那就调高产水量，进而降低Z值（产水量与Z值成反比），再去计算一遍直到跑出来合适的结果。