VIC模型简介

声明：内容版权归官网，参考https://vic.readthedocs.io/en/master/，本人搬运工一枚，仅限本人学习笔记

---

一 前言

VIC-5模型与VIC-4模型相比，在版本的基础架构方面有很多的改进。VIC模型可同时对水循环过程中的能量平衡和水量平衡进行模拟，弥补了传统水文模型对能量过程描述的不足。

二 概述

1.主要特征

VIC 模型是大尺度、半分布式的水文模型。它与其它的陆地表面模型有一些共同之处，同时也与全球环流模型相耦合。
①陆地表面被模拟成平坦均一的网格单元
.次网格单元的异质性（例如高程、土地覆盖）是通过统计分布来处理。
②输入是时间分辨率优于日均的气象时序驱动数据（例如降水、大气温度、风速、辐射等）
③陆地表面的地-气通量、水和能量平衡的模拟步长是日尺度或亚日尺度
④进入网格单元的水只能通过大气
.网格单元之间的non-channel流动是忽略的
–在网格单元中到达河川水网的地表径流和地下径流部分假设远大于网格单元之间交换部分
.一旦水进入河川水网中，假期它只会留在水网中不会回流到土壤中

这最后一点对VIC模型的应用会产生几种影响：
①每个网格单元都是独立模拟，网格单元之间是没有交互
②径流的汇流是用单独的模型进行模拟

2.土地覆盖和土壤

Fig1.VIC土地覆盖瓦片和土壤柱，与主要的水分和能量通量

2.1 土地覆盖

①可以将网格单元划分成任意数量的“瓦片”，每个瓦片对应于特定部分的覆盖单元（例如常青松叶林、草地等）（Fig1）
②地理位置或是土地覆盖类型的配置没有考虑；VIC中将多个同样覆盖类型成1个瓦片
③每个网格单元中的1个土地覆盖瓦片可以包含1个湖或湿地流域
④每个网格单元的通量和储量平均值是通过以每个瓦片的通量和储量以面积加权进行统计，写入输出文件
⑤对于给定的瓦片，是使用jarvis方式植被气孔响应来计算蒸腾
⑥相对于地表，植被冠层的能量平衡是单独考虑
⑦土壤蒸发和不同植物间风速和辐射衰减的占比，植被区和非植被区分数通过植被面积比率fv控制。

Fig2
⑧支持在驱动文件中可选输入逐日尺度的LAI,albedo和植被区比率，来替代特定的植被库或植被参数文件的月尺度气候参数
⑨计算潜在蒸散发对潜在蒸腾和潜在土壤蒸发使用面积加权计算。对于当前植被潜在蒸腾，使用建筑阻抗和LAI来计算缺少土壤湿度、水汽压差，温度或日晒限制情况下的冠层阻抗（冠层截流？）

2.2 土壤

①任意数量的土壤层，但通常为3层（Fig1）
②由VIC参数来控制下渗进最顶层
③在根植区的土壤层可因蒸散发耗散湿度
④土壤最上层的蒸发耗散湿度通过ARNO构建
⑤最上层到最下层的水流由重力驱动
⑥ARNO基流构建底层排水
⑦每一个土地覆盖瓦片有它独自的土壤温度廓线
⑧冻土模拟
⑨在土地覆盖瓦片内模拟土壤温度空间异质性
⑩模拟特定多年冻土过程，例如超额地面冰雪融化
⑪用土壤湿度和土壤质地的函数来计算潜水面深度
⑫用Farouki来解释有机质土壤的热特性

3.雪模型

VIC将雪分为如下几个类型：地面积雪、植被冠层上的雪、湖冰上的雪。（Fig3）VIC的雪模型有如下特征：
①地面积雪分为2层；最上面部分的积雪层来单独处理雪面的能量平衡
②考虑局部雪覆盖
③考虑风吹雪的升华

Fig3 VIC雪模型

4.气象

4.1气象输入数据

VIC需要如下气象驱动变量：
①降水
②空气温度
③风速
④长波辐射
⑤短波辐射
⑥大气压强
⑦水汽压强
这些驱动量必须以模型运行所要求的的时间步长来提供（例如 SNOW_STEPS_PER_DAY or MODEL_STEPS_PER_DAY）。通常，VIC用户使用MECLIM包来生成和分解驱动变量。https://vic.readthedocs.io/en/master/Documentation/References/#other-historical-references中了解更多信息。

4.2非气象输入数据

VIC可以读取逐日时序土地覆盖信息例如albedo，LAI，和植被冠层覆盖比例作为驱动变量。https://vic.readthedocs.io/en/master/Documentation/References/#other-historical-references

4.3高程波段

VIC可以考虑降水的空间异质性，由于暴雨前锋线/本地对流或是地形异质性。这里我们根据高程数据来考虑地形因素影响（Fig4）。这主要是对山区雪产生更准确的估计。这个功能主要是通过全局参数文件中的 SNOW_BAND 选项来控制。主要特征：
①可以对网格单元进行细分为任意数量的高程波段，来解释网格单元地形差异。
②针对每个高程波段，气象驱动参数以网格单元平均高程失效
③地理位置或是高程波段的配置没有考虑；VIC将相同高程区域划为1个波段块
④以面积为权重，将每个高程区的流量和储量进行平均写进输出文件
⑤特定高程波段值可单独写进输出文件

Fig4 VIC雪（高程）波段

5.冻土构成

5.1土壤热方程

VIC可以使用近似的土壤温度廓线或者有限差方程来考虑土壤冰部分内容。https://vic.readthedocs.io/en/master/Documentation/References/#other-historical-references
①Liang et al（1999）。设置全局参数文件中 QUICK_FLUX = TRUE; 默认值 FULL_ENERGY = TRUE 和 FROZEN_SOIL = FALSE
②Cherkauer et al（1999）。全局参数文件设置 QUICK_FLUX = FALSE;默认值是FROZEN_SOIL = TRUE
.默认情况下，有限差分公式使用隐格式
.默认情况下，有限差分公式的节点以指数分布

Fig5 VIC冻土构成

5.2多年冻土

这里在 QUICK_FLUX = FALSE 和 FROZEN_SOIL = TRUE情况下使用，例如构成形式https://vic.readthedocs.io/en/master/Documentation/References/#other-historical-references
.全局参数文件选项： IMPLICIT:使用隐格式来解决土壤热廓线，这是默认框架。
.全局参数文件选项：EXP_TRANS:使用指数阶段间隔（近土壤表层密集阶段间隔，深度间隔稀疏节点）（Fig6）

Fig6 VIC多年冻土增强功能

5.3空间霜冻的温度异质性

①细节详见https://vic.readthedocs.io/en/master/Documentation/References/#other-historical-references
②土壤温度在均值附近线性分布（Fig7）
③当平均温度低于冰点，允许土壤中一些水分运动

Fig7 VIC空间霜冻构成

6.动态湖和湿地模型

湖/湿地模型https://vic.readthedocs.io/en/master/Documentation/References/#other-historical-references来处理网格单元内的表面积水。每个网格单元允许湖/湿地系统包含一个土地覆盖瓦片。这里lake指任一表面蓄水，包括永久湖和季节性植被陆地。湖面积可随时间变化，作为蓄水量和地形的函数。在这上下文中，湿地是指土地覆盖瓦片曝露的部分。

6.1湖模型

①多层湖模型（Fig8），基于https://vic.readthedocs.io/en/master/Documentation/References/#other-references，https://vic.readthedocs.io/en/master/Documentation/References/#other-references，https://vic.readthedocs.io/en/master/Documentation/References/#other-references
②能量平衡模型
③混合，辐射衰减，可变冰覆盖
④动态湖面积是蓄容量的函数
.这个关系需要在输入参数中明确
.对于网格单元中只包含一个湖，蓄容量-面积的关系==湖深（+附近地形）
.对于网格单元中包含多个湖情况，模型认为单一的复合湖；这一实例中蓄容量-面积关系！=湖深（这种情况下的算法还在开发中）
.湖体可以直接与河道网络相连接。湖水可以接收同一网格单元中周边高地的径流，和上游网格单元的河流。应用示例https://vic.readthedocs.io/en/master/Documentation/References/#selected-application-references
.湖体排泄直接进入河道网络。湖体流出包括：
–河道径流（模拟成流过宽顶堰）地下径流（如果湿地干燥，将会流入湿地）
.为了运行湖模型，用户需要在全局参数文件中设置LAKES到合适的湖参数文件。网格单元中不包含湖，可以再湖参数文件中指示。
.为了关闭湖模型，用户需要设置LAKES为FALSE，或者删除全局参数文件中任一涉及LAKES的地方。

Fig8 VIC湖模型构成

6.2湿地模型

①在网格单元中，湖和湿地存在于单独“瓦片”
②动态湿地面积=tile 面积 - 动态湖面积
③允许湿地季节性洪水作为湖的增长和衰减
④湿地含水量/能量流的计算与高地瓦片的计算相似
⑤湿地土壤比高地的土壤更加湿润，由于频繁的洪水和湖体水入侵

Fig9 VIC动态湖/湿地交互

7.碳循环处理

①VIC模型可以选择性的模拟光合作用、自养呼吸和异养呼吸，https://vic.readthedocs.io/en/master/Documentation/References/#other-historical-references
②植被物候不是动态的。生物量没有被模拟。当碳循环关闭时，LAI和其它植被特征以相同的方式处理
③有三个土壤碳储池：litter（存留时间2.86年），中间层（存留时间33.3年），slow（滞留时间1000年）
④上一年度总的净初级生产量（NPP；光合作用-自养呼吸）在当前年度以稳定的通量（1/365*（上年度总NPP））加到litter pool
⑤三层土壤碳池的呼吸量与（池中碳总量）e（-residence time）（Lloyd-Taylor 温度相关）*（土壤湿度函数）
⑥各碳池进入大气的CO2是连续量。对于liiter池，剩余的呼吸碳进入到intermediate和slow。中间层进入slow pool。对于slow pool所有的呼吸碳全部进入大气。

8.汇流

汇流模拟使用单独的模型，独立于地面模拟。通常使用Lohmann等提出的汇流模型。原始的汇流模型已不再支持。RVIC版本模型可以获得：http://rvic.readthedocs.org/