变分方法估算西北半干旱区地表热通量适用性研究.pdfVIP

* 变分方法估算西北半干旱区地表热通量的适用性研究 曹帮军，张述文 （兰州大学 大气科学学院/甘肃省干旱气候变化与减灾重点实验室，甘肃 兰州 730000 ） 摘 要 利用 2009 年 6-8 月兰州大学半干旱气候与环境观测站（SACOL ）的地表辐射观测资料，近地层大气温度、 湿度和水平风速的三层梯度观测资料，以及热通量的涡动相关观测资料，比较了变分法估算的地表热通量 与涡动相关观测值间差异，评估了不同相似函数对变分法估算地表热通量的影响。结果表明：变分法估算 的地表热通量与涡动相关观测值在时间变化趋势上符合得很好，但中午前后时段具体数值还存在差异。层 结不稳定时，不同相似函数对变分法估算热通量影响较大，但当不稳定变弱时其影响显著减小；层结稳定 时，不同相似函数对变分法估算热通量影响非常小，与观测误差相比可以忽略。 关键词：变分法；涡动相关法；地表热通量；西北半干旱区；相似函数 1 引言 用变分法估算地表热通量的研究工作已开展不少[1-6]，其基本思想是利用近地层大气温 度、湿度和水平风速的梯度观测资料，基于 Monin-Obukhov 相似性理论[7-8]在适当地表能量 平衡的约束下寻求一类泛函方程的最优解，进而估算出地表热通量。由于早期边界层观测试 验缺乏涡动相关法得到的地表热通量的相关观测资料，已有变分法估算地表热通量的工作只 [9] [10] 是与传统的廓线法 和 Bowen 比能量平衡法 得出的热通量结果进行比较，发现变分法在 计算精度和稳定性方面都优于两类传统方法。目前，尽管涡动相关法观测的湍流通量还不算 完美，但是一般认为该方法得出的脉动观测在经过严格数据处理和质量控制后统计出的地表 热通量的精度是很高的[11-12] 。基于此，本文将它视为地表热通量观测的标准值，用来比较变 分法估算的地表热通量，考察变分法估算值的精度。另一方面，上述不同方法计算地表热通 量的对比研究还没有在西北半干旱区进行；半干旱区介于湿润区和干旱区之间，是一个非常 敏感的地带[13] 。研究该区的地表热通量的变化特征对理解干旱灾害发生的机理和提高半干 旱区的预测水平都具有重要意义。因此，有必要开展变分法在该区估算地表热通量的适用性                                                                基金项目：国家自然科学基金项目资助） 已被《高原气象》杂志接收，近期将发表 已获“兰州大学2011 年创新创业行动计划”二等奖  1    研究，填补以往研究空白，这是本文第一部分研究内容。 此外，考虑到以往变分法估算通量研究中只选取代表某一特定形式的温度、湿度和水平 风速的相似函数，如稳定时只选 Beljaars 和 Holtslag 拟合的关系式[14]，不稳定时只选 Paulson 的公式[15] 。已有研究发现不同相似函数对廓线法计算地表热通量的结果影响很大[16] ，由于 变分法是把廓线法和 Bowen 比能量平衡法在最优估计理论基础上融合在一起，是否也对温 度、湿度和水平风速的相似函数选取敏感？这需要试验研究才能回答。为此，本文在层结不 稳定情况下增加另外四组相似函数，在稳定情形时增加另外三组相似函数形式，定量研究变 分法对不同形式的相似函数敏感程度，完善对变分法估算地表热通量的性能评估，这是本文 第二部分研究内容。 2 资料介绍 本文资料为兰州大学半干旱气候与环境观测站(SACOL)在 2009 年夏季 6、7 和 8 三个月 的近地层大气温度、湿度、水平风速的观测，地表辐射的观测和涡动相关法观测的湍流通量。 该站位于甘肃省榆中县萃英山上(35.57 oN, 104.08 oE) ，站点海拔高度为 1965.8 m，下垫面为 覆盖原生植被的黄土高原地貌。该地区年平均降水量为 381.8 mm，属温带半干旱气候。观 测站设有一个 32 m 高度的边界层观测