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

今天介绍用Invest模型计算时产水量时需要的一个栅格数据——PAWC（植物有效含水量）。

植物有效含水量是从一些标准土壤图中得到的分数。它被定义为体积场容量与永久萎蔫点的比值之差。植物有效含水量通常以体积值(mm)表示。要得到分数除以土壤深度。如果PAWC不可用，则需要从重平均土壤质地(%粘土，%沙子，%粉土)和土壤孔隙度的多边形形状文件中获得栅格。我是通过后者方式计算得到的PAWC栅格。计算公式如下：

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

注意：土壤数据中不包含有机质属性字段，土壤有机质需要通过土壤有机碳计算。土壤有机质=土壤有机碳*1.724。

有了土壤数据之后怎么实现PAWC的计算呢？让我们接着往下看。土壤数据带有砂粒、粉粒、黏粒、有机碳等字段，首先需要将它们分别单独提取出来，形成各自的栅格数据，那么怎么提取呢？

1、栅格数据某一字段的提取

Spatial Analyst工具——重分类——查找表

输入栅格输入土壤数据，查找字段选择砂粒、黏粒、粉粒、有机碳字段，就能提取出砂粒、黏粒、粉粒、有机碳 ，使其作为单独栅格土层。

有了有机碳栅格数据后，我们来计算有机质。

2、土壤有机质的计算

 打开栅格计算器：Spatial Analyst工具——地图代数——栅格计算器

土壤有机质=土壤有机碳*1.724。

 3、PAWC的计算

同样地，打开栅格计算器，Spatial Analyst工具——地图代数——栅格计算器

输入以上公式， 54.509-0.132 * "Sand2.tif" - 0.003 * "Sand2.tif" * "Sand2.tif" - 0.055 * "Sily2.tif" - 0.006 * "Sily2.tif" * "Sily2.tif" - 0.738 * "Clay2.tif" + 0.007 * "Clay2.tif" * "Clay2.tif" - 2.688 * "OM2.tif" + 0.501 * "OM2.tif" * "OM2.tif

注意：计算出来的PAWC取值为0-100，单位是%，在invest软件输入的PAWC要求单位为1，取值范围为0-1。所以需要再将上述得到的PAWC通过栅格计算器除以100即可。

OK，今天就分享到这里，我们下期再见！！！