内陆河高寒山区流域水热耦合模型的构建及应用.docxVIP

内陆河高寒山区流域水热耦合模型的构建及应用 中国寒冷面积为417.042。这是中国河流和河流的主要起源和生态环境极其薄弱的地区。在全球变暖的背景下,我国大多数寒区出现了冰川退缩、雪线大幅度上升和冻土退化等现象,从而造成了多数河源地区产流量减少、水土流失严重和植被种类演替等环境问题。冰川水文、积雪水文和冻土水文是寒区水文研究的主要内容。相对而言,冻土水文过程在寒区分布的更为广泛。在寒区流域产汇流过程中,多年冻土活动层和季节性冻土融化层底部是深度不断变化的隔水层,冻结和半冻结土壤和岩层中固态含水量的存在,改变了土壤—岩层的能量收支和平衡,既增加了冻结—融化潜热过程,同时改变了土壤—岩石层的导热系数、热容和总能量,以及感热传导过程。固态含水量的存在,还改变了土壤—岩石层的结构,减小了土壤—岩层的有效孔隙度和实际土壤田间持水量,从而改变了土壤液态水分—土壤水势关系(土壤水分特征曲线),改变了土壤—岩石层的实际水力传导率,最终改变了液态水分的流向、流速、流长和流量。陆地表层是水文过程中水分交换最为频繁、最为复杂的地区,是地—气相互作用的汇集地,是对外界影响最为敏感的地区之一。因此,土壤冻融作用和水热交换,贯穿于整个寒区流域的产流、入渗和蒸散发过程之中,是寒区流域水文循环过程的中心环节。探讨一个寒区流域的水文循环过程机理,最重要的一环就是要进行冻土水文过程观测试验与模拟研究。 我国内陆河流域的水资源主要来自于降水量相对较多的高大山区。在一个内陆河山区流域,冰川和多年积雪是“固体水库”,水源涵养林和高山草甸是“绿色水库”,而冻土则是调节和连接这两个“水库”和河川水系的纽带,同时也通过抑制蒸发和冻融过程中水势变化来涵养水源。内陆河山区海拔高度相对青藏高原腹地要低,多年冻土层厚度相对较浅,季节性冻土比例较大,冻土、冰川和积雪对气温的变化更为敏感,在全球变暖的趋势下,许多内陆河高山地区的多年冻土转变为季节性冻土,甚至在乌鲁木齐1号冰川顶部出现了冰面湖。内陆河高寒山区冻土和积雪的变化将直接引起山区流域生态环境的变化,进而导致流域水文过程的变更,最直接的反映将是流域出山口径流量和径流年内分配的改变。这对水资源极为有限且水量主要来源为山区径流的内陆河流域来说,不论是对流域生产规划,还是对流域生态环境恢复和重建,均会产生较大的影响。 因此,需要在内陆河山区流域开展水文循环过程机理研究,在进行小流域综合寒区水文过程观测试验的同时,还需要进行模型模拟研究,将小流域研究成果推广应用到流域尺度,定量水文循环过程中的各变量,以期使全球变化背景下内陆河山区流域所面临的科学问题明朗化。内陆河山区流域季节冻土比例较大,其综合集成成果基本能够包含我国北方多数河源冻土分布区的水文过程。因此,在内陆河山区流域模型研究方面,应该发展一种能够适用于任何存在季节冻土和多年冻土山区流域的分布式模型。 国外分布式水文模型的研究至今已有近40 a的历史,并得以广泛应用。国内分布式水文模型研究起步较晚,但近年来开展得如火如荼[4,5,6,7,8,9,10,11,12,13]。国内许多学者分别从不同角度描述了当前分布式水文模型的研究进展、展望及存在的问题[14,15,16,17,18,19]。但纵观国内外众多的分布式水文模型,很少有模型涉及土壤冻融作用对流域产汇流过程的影响。鉴于上述分析,本文以黑河干流山区流域为例,构建了一个“内陆河高寒山区流域分布式水热耦合模型(DWHC)”,以期能够服务于内陆河山区流域和其他存在季节冻土和多年冻土的流域水文循环过程的研究。由于模型输入变量较多,在流域气象、土壤和植被参数资料等稀少的情况下,拟主要进行模型原理探讨。 1 寒区水文过程的构建技术方案 内陆河高寒山区流域分布式水热耦合模型(DWHC)主要由气象因子模型、植被截留模型、冰川和积雪融化模型、土壤水热耦合模型、蒸散发模型、产流模型、入渗模型和汇流模型等8个子模型组成(图1)。模型根据大气降水到达地表以后各水文过程发生的先后进行构建,考虑了寒区流域冰川水文、积雪水文和冻土水文过程,其中最关键之处是利用土壤水热耦合模型将流域产流、入渗和蒸散发过程融合成一个整体(图1),弥补了分布式水文模型中缺乏冻土水文过程的问题,从而能够全面定量描述整个寒区水文过程。在大量吸收已有成果的基础上,模型在植被截留、入渗、产流和蒸散发计算方面也有所改进和创新,部分模块还设计了多套可选择方案。其中气象因子模型设计了两套方案,一是以地面观测资料进行空间分布,二是由中尺度大气模式MM5输出结果驱动。有关气象因子模型部分,将分别在“Ⅱ:地面资料驱动结果”和“Ⅲ:MM5嵌套结果”中介绍。下面分别就模型输入变量、模型基本原理和模型参数予以介绍。 2 时空模型分辨率和输入变量 2.1 模型建立区域 模型以数字高程模型(DEM)栅格单元为最