VIC分布式水文模型——参数估计和数据源

目录

1 前言

2 背景参数

2.1 地理数据

 2.2 气象数据

 3 植被参数

4 土壤参数

5 水文参数

6 总结

---

1 前言

如前面的模型原理中提到的，VIC的输入参数中包括背景参数（地理数据和气象强迫数据），植被参数，土壤参数和水文参数。VIC模型建立了模型需要的全球植被和土壤两大参数库，各参数值以网格为单元估计。

参数率定时会对6个水文参数进行率定。下面介绍背景参数，植被参数和土壤参数植被的数据源，还有在参数率定时不同水文参数会对结果带来的作用。

2 背景参数

包括流域的DEM数据和气象数据。

2.1 地理数据

首先我们需要获取覆盖流域的DEM数据，这个数据可以从地理空间数据云获取，链接如下。https://www.gscloud.cn/#page1/2

 进入网页后，选择高级选择，筛选数据集选择DEM数字高程数据，根据需要选择30m或者90m。根据经纬度或者行政区选择流域的位置，根据筛选出来的数据源进行勾选，将所研究流域覆盖即可。

2.2 气象数据

需要流域内及附近每个气象站的坐标以及各自的历时观测数据，需要日观测数据，需包含降雨量，日平均气温或者日最高/最低气温以及风速。有时这些数据记录是不完全的，有的数据项因为各种原因会在某些日子缺数据，或者缺失某个连续时间段的记录，一般这种情况数据中会使用某特定的数字来表示。必要时可以用数据处理工具或者编程处理，如R语言等。

如《中国国家级地面气象站基本气象要素日值数据集》的说明文档中的标注：

台站海拔高度      +100000   当台站海拔高度为估测值时，在估测数据基础上加100000

各要素项              32766   数据缺测或无观测任务

气压日极值          +20000  气压极值取自定时值，在原值上加20000

日最小相对湿度  +300    最小相对湿度取自定时值，在原值上加300

风速         +1000   当风速超过仪器测量上限时，在上限数据基础上加1000

风向         1-17    用数字表示风向方位，17表示静风

                +100    当表示风向为八风向时，在原值上加100

                90X 风向出现X个时，风向数据用个数X表示

                95X 风向至少出现X个时，风向数据用个数X表示

降水量     32700   表示降水"微量"

                32XXX   XXX为纯雾露霜

                31XXX   XXX为雨和雪的总量

                30XXX   XXX为雪量(仅包括雨夹雪，雪暴）

蒸发量     32700   表示蒸发器结冰

                +1000 蒸发器中注入的水全部蒸发，在注入的水量数据基础上加1000

0cm地温 +10000  实际温度（零上）超仪器上限刻度，在上限数据基础上加10000

                -10000  实际温度（零下）超仪器下限刻度，在下限数据基础上减10000

 3 植被参数

在植被参数库中，植被类型按照马里兰大学开发的1km AdVanced Very High Resolution Radiometer (AVHRR)土地覆盖类型分类划分。该分类包含了11种植被参数和2种土地利用参数（裸土和城市建筑）。植被的网格覆盖率由植被在单元网格内的覆盖面积占网格面积的比重大小决定。植被高度、最小气孔阻抗等按植被类型不同估计。

马里兰大学的土地覆盖类型图的下载链接如下🔗：https://daac.ornl.gov/cgi-bin/dsviewer.pl?ds_id=969

按照Calder方法根据植被的平均高度确定粗糙率相关参数。结构阻抗按照Ducoudre等的方法确定。植被叶面指数（LAI）以月尺度变化区分。

4 土壤参数

土壤参数，比如土壤结构和土壤容重密度从5-min Food and Agricultural Organization-United Nations Educational, Scientific, and Cultural Organization  (FAO-UNESCO)的数字土壤地图和World Inventory of Soil Emission Potentials (WISE)土体参数库提取。其余土壤参数，如孔隙率、饱和水力传导度、非饱和水力传导度指数根据Cosby计算方法获得。

可使用的土壤数据的来源很多，如HWSD，1km分辨率的土壤类型数据库等，可从IIASA网站或者国家冰川冻土沙漠科学数据中心下载。链接如下：

http://www.ncdc.ac.cn/portal/metadata/a948627d-4b71-4f68-b1b6-fe02e302af09

 另外，在用冻土模式或能量平衡模式时，需要对土壤气泡压力BUBBLE进行赋值，前面的文章中也写了该值可以用水力传导度expt进行估算，也可以用《Handbook of Hydrology》中的表格进行估计。表格如下：

5 水文参数

即前面文章中提到的需要率定的几个参数，其取值范围也给出了，可参考之前的文章VIC分布式水文模型——模型原理。

• 其中，BI值将入渗水量描述成饱和区域网格面积的函数，其值高说明下渗透水量低，从而增大产流量；
• 当Ds值增大，在低土壤含水量时下层土壤的基流将增大；
• Ws类似于Ds，其值高将提高产生非线性基流反映的土壤含水量，势必延迟洪峰的出现；
• 对于土壤深度D2，D3，其值增大将增加蒸发量，延缓基流的峰值。

对于缺失资料的地区来说，其参数可以进行参数移植，将相似气候和土壤类型的地区的参数移植过去。Tory等对其进行了研究，结果表明，在不同时间尺度上的率定移植对模型模拟结果影响较小，但在空间尺度上的移植差异较大。、

6 总结

上文对VIC模型计算所需准备的参数和需要率定的参数进行了总结。其中的一些描述参考了李致家等人的《现代水文模拟与预报技术》。一些文字表述如有不对之处欢迎指出。