参数准备-土壤参数

Hanjiang_soil.txt：

0.1°栅格重采样及其地理信息：

VIC要求的坐标都是经纬度坐标，因此Arcgis处理数据应统一转换为地理坐标而无需投影坐标。

1 ：包括研究区的SRTM90m影响重采样为0.1°；
2 ： 获取更大范围研究区矢量范围，Cartography Tools -> Masking Tools -> Feature Outline Masks

3 ：用该shp掩膜提取大于实际研究区的rasters；从而保证几乎完全覆盖研究区；
4 ：裁剪后raster to points；（由于通过DEM获取，所以属性表里存有海拔信息和经纬度信息。对于空的信息表，可以：open the attribute table of points.shp，add filed recoding ’long‘ and ‘lat’（经纬度），右键字段选择计算几何，分别得到每个位点的经纬度信息）

5. .dbf文件转为txt或xlxs
6. Matlab

土壤属性（HWSD数据集+SPAW小软件获取）

必须填充输入文件的所有列，但只有在激活热量平衡或冻土模型时才需要定义某些参数。只有当FULL_ENERGY为真时，红色的变量才需要有有意义的值，否则它们将被忽略(但是占位符值-9999必须仍然存在于该字段中，除非文档明确指出该字段是可选的)。可以通过 ## 来添加标题行。

Bulk_density

第一种方式：
HWSD土壤数据处理流程：

1）打开栅格影像（直接打开，无法查看属性表）——右键——属性——拉伸由百分比截断改成无——唯一值——回到栅格影像——右键——打开属性表即可看到value、count属性
2）打开属性数据——HWSD.dbf——HWSD-DATA.dbf——arcgis打开后——toolbox——转换工具——表转excel——转出的excel记为hhhh——arcgis打开hhhh
3）栅格影像——右键——连接——选择栅格影像的value——连接hhhh的mu_global——检验连接可能会waring保留字（没关系）——连接——打开栅格影像属性表，可看到hhhh中的所有字段
4）影像展示——栅格影像右键——属性——符号系统——唯一值——value设置为sym90或者pH等属性字段——符号系统后字段——主显示字段——选择sym90或者pH等属性字段

注：HWSD数据分为栅格和属性数据两部分，其中栅格数据中显示的value值为属性数据中mu_global属性字段。

Topsoil（0-30cm)
Subsoil（30-100cm）

土壤容重顶薄层和第二层使用Topsoil数据，第三层使用Subsoil数据。

通过HWSD判断水体：519/577/689/690/730/731/763/769/770/807

第二种方式：基于南水北调调查项目自己插值：（此处选择IDW）

Topsoil(调查点0-20cm左右）
Subsoil（调查点20- ）

导入xy坐标点————选择IDW方式——Field处选择bulk_density数据——通过调整各项参数使RSM最小——出图后先双击在extent出选择图幅范围与DJK_kuoda相同——然后右键data-转栅格（分辨率0.1）：环境里第三个选择图幅范围与DJK_kuoda一致——掩膜提取——栅格转点——点数据输出

expt

exp指数参数，来源于Campbell方程的指数，=3+2/λ（或者3+2b），值不应小于3.

Wcr_FRACT

该参数作为土壤在临界点的土壤湿度（最大土壤湿度的比例，最大土壤湿度=孔隙度 * 土层厚度），它是由于土壤蒸发，饱和导水率开始下降时的土壤湿度。数值上根据不同的土壤质地，表示为田间持水量的70%。（此处的田间持水量为-33KPa状态下的含水量） ——官网定义

Wpwp_FRACT

该参数作为土壤在萎焉点的土壤湿度（最大土壤湿度的比例，最大土壤湿度=孔隙度 * 土层厚度），萎焉点表示为1500KPa下的土壤含水量，为每一土壤质地的近似值。——官网定义

Ksat

每一土层的饱和水力导度。
（研究区发现：亚表层该值整体大于表层，原因？）

上述变量从SPAW——swc_excel公式中获得：

init_moistInitial

init_moistInitial moisture content of each layer can be set at any reasonable value. One approach is to use fractional soil moisture content (expressed as a fraction of the maximum soil moisture; max. soil moisture = porosity * layer depth) at the critical point, Wcr, which can be computed for each layer as a depth in millimeters by multiplying Wcr by the thickness of the layer in meters, and then multiplying by 1000.

avg_T

土壤热阻尼深度(dp该深度下温度保持恒定）下的温度，通常假定为年均温，该温度作为所有计算土柱热通量的底界。
（气象数据求每一个栅格年均温获得）