《基于集合卡尔曼滤波和HYDRUS-1D模型的土壤剖面含水量同化试验 EnKF and HYDRUS-1D based data assimilation experiments for improving soil moisture profile prediction》.pdfVIP

水 利 学 报 2012年 11月 SHUILI XUEBAO 第43卷 第 11期 文章编号 ：0559—9350(2012)11-1302—10 基于集合卡尔曼滤波和HYDRUS一1D模型的土壤剖面 含水量同化试验 王 文，刘永伟，寇小华，吕海深 (河海大学 水文水资源与水利工程科学 国家重点实验室 ，江苏 南京 210098) 摘要 ：本文建立 了一个基于集合卡尔曼滤波(EnKF)方法和 HYDRUS一1D模型相结合 的单点土壤水分 同化方案 ， 以期通过 同化表层土壤含水量 的实测值来改进土壤剖面含水量预测精度 。结果表明，该 同化方案能有效改进表层 (4-10cm)及浅层 (50cm)的土壤含水量预测精度 ，对下层 (100～150cm)的土壤含水量预测精度仍有影响，但其影 响有好有坏，而对于深层 (≥200cm)的土壤含水量预测精度几乎没有任何影响；对同化过程中的状态变量及关键 状态参数 的变化分析表明，同化效果的产生是 由于通过数据 同化修正了模型的状态变量及关键参数 ，而随着深度 的增加 ，变量和参数的调整量减小，同化效果也就随之减弱。 关键词 ：集合卡尔曼滤波 ；HYDRUS一1D；数据 同化 ；土壤剖面含水量 中图分类号 ：S152．7 文献标识码 ：A 1 研究背景 土壤含水量是降雨径流过程模拟 中的一个重要状态变量 ，具有高度的时空变异性 。土壤含水量 地面实测只能采集有限地点的数据 ，而且采集地表数据比采集深部数据要更容易。遥感可以实现大 面积观测 ，但仅能提供表层观测数据。因此 ，如何根据地表土壤含水量推求土壤剖面含水量信息 已 成为水文研究 的热点之一。 用表层土壤含水量估计土壤剖面含水量的方法大体有4种…：(1)表层观测值 的统计外推法；(2) 表层观测值与剖面水量平衡模型相结合的方法；(3)辐射传输模型的反演方法；(4)参数化的剖面模 型方法 。其 中，基于观测与模型模拟相结合的方法被认为是最有前途的，而将二者结合最常用的手 段是数据同化技术 (包括卡尔曼滤波等)，很多研究者进行了相关尝试 。例如 ，Walker等 利用卡尔曼 滤波 (KF)方法 ，将表层土壤含水量观测数据 同化到一个简化的土壤水分模型中，进行一维土壤水分 剖面反演 的研究 。Crosson等 也利用KF，将地表通量一土壤水模型及辐射转换模型和微波遥感观测 的与土壤含水量密切相关的亮度温度数据相结合 ，实现三维 的土壤含水量估计。Crow和Wood 应用 集合卡尔曼滤波 (EnKF)将微波遥感观测亮度温度同化到陆面模型中，估计根层土壤水分。Jia等 利 用EnKF，将卫星微波遥感AMSR—E亮度温度数据 同化到土壤水模型中，结果表明不仅可改进表层土 壤含水量估计 ，也可以改进 10～50em深度土壤含水量估计 。Tian等 考虑到模型次网格的非均质性和 土壤水的冻融过程 ，构建了一个双通道路面数据 同化系统(DLDAS)，用于同化AMSR—E卫星亮度温 度 (Th)数据 ，对表层 10em的土壤水分估计有明显改进。 HYDRUS一1D模型 是非饱和流和溶质运移模拟应用最广的模型之一 ，在土壤水分运移模拟中也 有许多研究 。本文建立一个基于EnKF方法和HYDRUS一1D模型相结合 的单点土壤水分 同化方案 ， 收稿 日期 ：2012—01—06 基金项 目：国家 自然科学基金项 目40930635) 作者简介 ：王文(1967一)，男，江苏人 ，教授 ，博士，主要从事环境变化对水文过程的影响及水文遥感研究。 E-mail：W ．W ang@ 126．~om 一 1302 — 万方数据 利用实测多层土壤含水量逐 日时间序列数据进行同化试验 ，分析通过同化表层土壤含水量的实测值 来改进土壤剖面含水量估计的效果 。 2 基于EnKF和HYDRUS一1D的数据同化方案 2．1 HYDRUS一1D模型 HYDRUS一1D模型主要是应用数值解法解决变饱和土壤 中垂 向一维水 、 热及溶质 的运移 问题 。本试验用HYDRUS一1D模