SWAT模型学习——土壤数据制备过程的学习记录

目录

1.源数据下载

 2.数据初步处理

 3.土壤参数整理与计算——土壤属性索引表

3.1 HWSD自带参数对照

3.2 SPAW计算剩余参数

 3.3  土壤侵蚀力因子USLE_K计算

3.4 水文分组HYDGRP计算

---

1.源数据下载

使用HWSD数据集，点击 Download database 和 HWSD Raster，下载得到原始全球的栅格数据及其对应的参数库。

https://www.fao.org/soils-portal/soil-survey/soil-maps-and-databases/harmonized-world-soil-database-v12/en/

 数据下载后即可在ArcGIS中打开查看

 2.数据初步处理

Access打开.mdb数据库，并将HWSD_DATA另存为.xls文件备用

ArcGIS打开.bil的栅格文件

—— 按研究区域进行栅格裁剪

—— 定义投影（Data Management - Projections and Transformations - Define Projections, WGS-1984）

——创建属性表（Data Management - Raster - Raster Properties - Build Raster Attribute Table）

—— 表连接HWSD_DATA数据（直接链接mdb数据库里的，否则数据太大，加载不出来），以Value字段对应MU_GLOBAL字段连接，仅保留匹配值

导出属性表，此时导出的属性表就包含栅格自身的value和count以及连接的属性，将其导出。

完成这一步后，我研究区内的土壤属性值（Value）共108种，种类太多，进一步进行重分类的操作。对土壤属性值进行整理，一共有31种HWSD的亚类，按土壤分组一共有14类，因此按照土壤分组对其进行重分类。每个分组内占比最高的作为最终土壤的类别，采取其属性。进行重分类（Spatial Analyst Tools - Reclass - Reclassify），最终得到14类土壤。

重分类分组参考官方文档

isricu_i9264_001.pdf (fao.org)

 3.土壤参数整理与计算——土壤属性索引表

基于SWAT自带的土壤属性表，查看其属性库SWAT2012里的usersoil属性表，导出，并获得其表头，构建自身的土壤属性表。

3.1 HWSD自带参数对照

各指标之间的对应情况如图所示，其中，SOL_Z1和SOL_Z2分别表示第一层和第二层，对应HWSD中的T_和S_，一一填写参数。

注：HWSD里土壤参数最多到第二层，两层的时候一般第一层土壤表层到底层的深度SOL_Z1为300，第二层SOL_Z2为700

3.2 SPAW计算剩余参数

第一步后，使用SPAW软件计算TEXTURE、SOL_BD(可选)、SOL_AWC、SOL_K。

SPAW设置：

 参数输入并计算：

依次输入四个值Sand(SAND1), Clay(CLAY1), Organic Matter (SOL_CBN1/0.58), Gravel (ROCK1)，第二层的也为对应参数，每个参数输入后按Enter，最后得到结果

需要注意的是：Organic Matter = SOL_CBN1/0.58，有教程提到Salinity输入对结果也有一定的影响，因此输入该参数，Salinity = TECE （HWSD）

 最终得到的结果对应参数分别为：

Texture Class —— TEXTURE （土壤结构），用英文首字母缩写表示，土壤层之间用-连接。

Available Water —— SOL_AWC （土壤有效含水量）（或Field Capacity - Wilting Point）

Sat.Hydraulic Cond —— SOL_K （饱和水力传导系数）

Matric Bulk Density —— SOL_BD（湿容重）（仍采用自带的结果）

依次计算并填入参数库中。

 3.3  土壤侵蚀力因子USLE_K计算

基于公式进行计算：

计算完成，填入参数库。

3.4 水文分组HYDGRP计算

根据土壤中含沙量（SAND，%）分别计算土壤各层的平均颗粒粒径Y，然后计算下渗率X，取其中最小的X作为最小下渗率，考虑两个公式：

Y=SAND/10×0.03+0.002

X=(20*Y)^1.8

当沙粒含量为0时，Y取0.01mm，当沙粒含量为100%时，Y取0.3mm，粘土含量为100%时，Y取0.002mm。

注：这里两个不确定的点：一是Y的计算有的公式是0.002，有的地方公式是0.02，经多方查找，这里使用0.002；二是水文学分组的划分，有提到 “若最小下渗率出现在土层上层深度小于500mm时，则参考正常标准；若最小下渗率出现在土层上层深度500-1000mm时，则将土壤水文单元上调一类，即B调至A；若最小下渗率出现在土层上层深度1000mm之下，则基于1000mm之上的土壤下渗率来划分水文分组。” 这里暂未上调，直接按照最小下渗率结果确定水文学分组。

至此完成所有数据库参数的计算，导入到SWAT2012的数据库中。

4.土壤索引表建立

直接新建.txt文档（或者excel另存为），对应VALUE和NAME，完成土壤数据制备。

 参考：

SWAT模型教程---土地利用、土壤数据、气象数据的处理_hydgrp_DoWayChan的博客-CSDN博客https://blog.csdn.net/weixin_41917608/article/details/114634911世界土壤数据库（HWSD）全球土壤数据下载处理 (qq.com)https://mp.weixin.qq.com/s?__biz=MzU4MTQ5MTM5MA==&mid=2247487185&idx=5&sn=d2f1a66cf02ef1e7c89b3722967600cc&chksm=fd478100ca3008164f8c318dd695bbc1c766c5bad9d44c46ccdd18695042a438aa591a626252（SWAT-3）SWAT土壤数据库的建立 - 知乎 (zhihu.com)https://zhuanlan.zhihu.com/p/442738029HWSD土壤数据库介绍 - 知乎 (zhihu.com)https://zhuanlan.zhihu.com/p/489387024基于HWSD数据集的SWAT模型土壤数据库建立_哔哩哔哩_bilibilihttps://www.bilibili.com/video/BV1nS4y1N7kS/?spm_id_from=333.788.top_right_bar_window_history.content.click&vd_source=df3a04d66b98d6b97562de571e167374SWAT模型土壤数据库建立方法 - 百度文库 (baidu.com)https://wenku.baidu.com/view/b2edeffb270c844769eae009581b6bd97f19bce0.html?_wkts_=1686554709799&bdQuery=%E5%9C%9F%E5%A3%A4%E6%9C%80%E5%B0%8F%E4%B8%8B%E6%B8%97%E7%8E%87%E8%AE%A1%E7%AE%97