一种基于土壤温湿资料计算地表土壤热通量的温度预报校正法.pdfVIP

中国科学 D 辑 :地球科学 2008 年 第 38 卷 第 2 期 : 243 ~ 250 SCIENCE IN CHINA PRESS 一种基于土壤温湿资料计算地表土壤热通量的 温度预报校正法 ①* ② 阳坤 , 王介民 ① 中国科学院青藏高原研究所, 北京 100085; ② 中国科学院寒区旱区环境与工程研究所, 兰州 730000 * E-mail: yangk@ 收稿日期: 2007-10-30; 接受日期: 2007-12-20 中国科学院引进国外杰出人才项目资助 摘要 分析大气边界层观测站地表能量平衡需要估计其地表土壤热通量. 发展了一种由多 关键词 层土壤温度和湿度观测资料估算土壤热通量的新方法. 该方法首先求解一维热扩散方程得到 土壤热通量 土壤温度的基本廓线, 然后校正所求温度廓线与观测值的偏差, 最后积分温度廓线得到土壤 热传导系数 各层的热通量. 与众多的方法不同, 该方法不需要事先给定不易准确测量和推求的热传导(或 温度修正 热流板 热扩散)系数值. 通过与实测资料对比、模型合成数据试验、以及敏感性分析等, 表明该方法的 计算结果稳定可靠, 对土壤表层数厘米深度内有无观测资料也不敏感. 此外, 指出热流板可准 确测量热通量的方向和相位, 但所测通量值的误差常较大. 地表土壤热通量是地表能量平衡的重要分量; 国等[6]先后比较了这类计算方法, 认为数值离散法和 估计该通量的大小是几乎所有的大气边界层观测站 谐波法的精度较高. 范新岗等 [7]和Gao等 [8,9] 同时考 能量平衡分析都要涉及的问题. 由于土壤温度在表 虑热扩散和热对流机制也提出了新的计算方法. 这 层变化剧烈, 热流变化主要发生在表层, 土壤热通量 些方法或者假设土壤的热传导系数或扩散系数不随 一般随土壤深度近似呈指数衰减 [1]. 土壤热通量可由 土壤深度变化, 因而与实际情况不符; 或者确认土壤 热流板测量, 但除极少数情况外 [2], 热流板通常有一 的垂向非均匀性, 但又在假设各层温度日变化为正 定埋深, 其测量值不能代表地表的土壤热通量. 目前 弦曲线下由观测资料的峰值及其时相确定土壤热扩 的边界层观测站、自动气象站、以及生态网络站等, 散系数等, 计算过程误差较大. 有的经验方法是通过 大多有数层的土壤湿度和温度观测资料. 基于热扩 建立与净辐射关系来近似求地表热通量 [10], 但有关 散方程通过积分这些温湿资料来计算地表热通量仍 参数因地区而异, 且所得土壤热通量在相位上有一 是估计热通量的基本方法; 如Tanaka等 [3]应用此法计 定误差. 关于土壤热通量计算方面的简单回顾可参 算了青藏高原安多站的地表热通量, 李超等 [4]导出了 考张立杰等 [11]. 利用整层的地温信息来计算任一层的热通量的积分 本文旨在基于热传导方程发展一种估计地表热 算法. 还有一类计算方法是基于土壤热通量和温度 通量的新方法, 不需要事先给定土壤的热传导系数 的相位随土壤深度的加深而延迟或者其振幅随土壤 (或热扩散系数), 但能较准确地估算热通量. 需要指 深度的加深而指数衰减的原理等. Horton等 [5]和莫兴 出的是, 土壤深层热通量与地表热通量相差几个量