VIC分布式水文模型——Frozen Soil Model 冻土模型原理

空空君shy

已于 2023-12-09 21:00:11 修改

阅读量425

点赞数 2

分类专栏： VIC 文章标签：学习方法数学建模

于 2023-12-09 20:49:50 首次发布

本文链接：https://blog.csdn.net/weixin_48366455/article/details/134897528

版权

VIC 专栏收录该内容

5 篇文章 3 订阅

订阅专栏

1 概述

1.1 冻土在水文过程中的作用

冻土会影响水文过程。土壤中的冰会减少渗入地表的降水量和融雪量，足够高的含冰量会使土壤几乎不透水。冻土还能在冬季储存更多的土壤水分。冻土的土壤水分无法排出或蒸发，因此有冻土的地区春季土壤水分含量会明显高于没有冰冻土壤的地区。

1.2 VIC冻土算法概述

VIC的冻土算法分为两种，分别为Liang et al. (1999) 提出的采用土壤热通量方程（soil thermal flux equations）,另一种是Cherkauer and Lettenmaier (1999) 提出的一种更明晰的数值技术，土壤温度在土壤柱中的几个“热节点（thermal nodes）”处求解，节点的数量和位置可由用户定义。目前后者使用更为广泛。

土壤热通量（soil thermal flux）和湿通量（moisture flux）是一个耦合过程，当土壤温度低于0℃时，土壤水分开始冻结。由于土壤内部的压力以及与土壤颗粒的相互作用，即使土壤温度远低于0℃，仍有一部分土壤水分可以保持不冻。冰不会从土层中排出，也会减少渗入土层的新水量。冰还会改变土壤的热特性。与水相比，冰具有更高的导热性和更低的体积比热，且当水结冰时会释放出潜热，使得周围的土壤温度升高。

VIC分两个阶段处理这种耦合关系。

首先，模型求解土壤柱的能量平衡。这涉及到地表能量平衡的求解，地表能量平衡包括来自土壤薄表层的热通量和热储量。对地表温度进行迭代求解，以找到最接近于封闭能量平衡方程的值，每一步还计算出了每个土壤热节点的热通量。

接着，当求解出满足精度的地表温度时，得出的土壤热节点温度将用于确定土壤层的含冰量，进而计算出湿通量。时间步长结束时土层含冰量和含水量将用于估算下一时间步长开始时的土壤热属性。

2 VIC冻土算法

2.1 热通量求解

通过土壤柱的热通量方程如下

$C_s \frac{\partial T}{\partial t}= \frac{\partial}{\partial z}(\kappa \frac{\partial T}{\partial z})+\rho _iL_f (\frac{\partial \theta_i}{\partial_t})$

式中，κ为土壤的导热性（W/m/K），Cs为土壤体积比热（J m-3 K-1)，T为土壤温度（℃），ρi为冰的密度（kg/m3），Lf为潜热（J/kg），θi为土层中的含冰量（m3/m3），t为时间（s），z为深度（m），从地表向下为正方向。式子右边的最后一项表示因冻结而引起的潜热变化，因此仅适用于土壤冻结的情况。

Liang et al. (1999) 的算法为根据地表温度Ts，第一层土壤表面温度T1，土壤柱底层温度Tb来预测土壤柱的热通量，T1和Tb之间的温度分布假设为指数分布，而Ts是作为地表能量平衡的一部分来求解。然而初步测设表明，上述三个温度不足以准确求解土柱中的土壤冻融状态。

Cherkaure and Lettenmaier (1999）的算法则不使用指数温度曲线，而是在用户定义的土壤柱节点（如上图的T1-T6）上明确求解土壤温度，Ts，T1和Tb的定义不变。测试表明，减少节点数量不仅会减少计算时间，还会降低求解精度。VIC目前该算法使用广泛。

对上面方程的求解采用显式有限差分，下边界条件Tb设置为年平均气温，上边界温度或地表温度，作为地表能量平衡计算的一部分进行预测。当数值格式计算新的地表温度时，通过土壤柱计算热通量，以提供最新的地表热通量输入。当数值模式将地表能量平衡误差最小化到允许误差范围内时，数值模式停止计算。

其中地表能量平衡积分基本方程为

$R_n=H+\rho_wL_eE+G+\Delta H_s$