生态系统服务功能模块之水源涵养关键参数计算(二)

前言：（碎碎念：这两天巨开心，昨天一天计算了三期的产水量，哈哈哈哈，原理搞清楚，数据收集齐就成功了一大半，我们要加油哈）。上一期我们知道了水源涵养是如何进行计算的，这期就其中较为关键的几个参数的计算进行详细说明。

一、地形指数T1

𝑇1=𝑙𝑛(𝑎/𝑡𝑎𝑛𝑏)

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

下面详细介绍操作过程：

1.a与b的计算

①对DEM数据进行填洼

ArcToolbox→Spatial Analyst Tools (Spatial Analyst工具)→Hydrology (水文分析)→Fill (填洼)

 

 ②填洼后的数据进行坡度计算，获得slope

ArcToolbox→Spatial Analyst Tools (Spatial Analyst工具)→Surface (表面分析)→Slope(坡度)

 ③填洼后的数据计算水流方向

ArcToolbox→Spatial Analyst Tools (Spatial Analyst工具)→Hydrology (水文分析)→Flow Direction(流向)

 

④水流方向数据计算汇流累积量

ArcToolbox→Spatial Analyst Tools (Spatial Analyst工具)→Hydrology (水文分析)→Flow Accumulation(流量)

 

⑤计算a

 ArcToolbox→Spatial Analyst Tools (Spatial Analyst工具)→Map Alqebra (地图代数)→Raster Calculator (栅格计算器)

 

公式为：

Con("Flow_accumulation.tif"==0,1,"Flow_accumulation.tif")*250000/Con("Flow_direction.tif"==1,500,Con("Flow_direction.tif"==4,500,Con("Flow_direction.tif"==16,500,Con("Flow_direction.tif"==64,500,Con("Flow_direction.tif"==2,500*SquareRoot(2),Con("Flow_direction.tif"==8,500*SquareRoot(2),Con("Flow_direction.tif"==32,500*SquareRoot(2),Con("Flow_direction.tif"==128,500*SquareRoot(2)))))))))      

其中，250000为栅格大小，即栅格分辨率的平方，需要根据自己的数据像元大小进行修改（本文为500m*500m,因此为250000），且当"Flow Accumulation"为0时没意义，至少应为1。对于流向为1、4、16、64 的栅格，除以栅格尺寸，这里是500；对于流向为2、8、32、128的栅格，除以栅格尺寸的√2倍，这里是√2倍的500。

2、T1的计算

𝑇1=𝑙𝑛(𝑎/𝑡𝑎𝑛𝑏)

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

 ArcToolbox→Spatial Analyst Tools (Spatial Analyst工具)→Map Alqebra (地图代数)→Raster Calculator (栅格计算器)

公式为： Ln(a.tif"/Tan(Con("Fill_slope.tif<=0,0.00001,Con("Fill_slope.tif">0,"Fill_slope.tif"*3.1415926/180))))

同理，求Tan（slope）时需要先将角度化为弧度；且当角度为0时将其自定义为0.00001，否则无意义。

我通过Arcgis运行出来的地形指数T1结果如下：

二、植被生理系数表

根据Invest模型的指导手册、联合国粮农组织（FAO）的《作物蒸散量作物需水量计算指南》以及之前学者的相关研究，考虑干旱区的植物根系长度会普遍大于湿润区的实际情况，设定蒸散系数和植物根系深度参数，属性类别显示为非植被覆盖时，模型将不会使用植物根系深度这一指标，因此将其都设为1，如下表所示：

样表1

 OK，今天的分享就到这里了，有什么问题欢迎大家指出，我们下期再见。