第二次青藏高原综合科学考察研究
肖敏轩, 赵林, 张钰鑫, 高骏强, 邹德富, 胡国杰, 赵建婷, 李智斌, 赵拥华, 肖瑶, 刘广岳, 杜二计, 刘淼
活动层内部的冻融锋面是冻融过程中冻结土层与融化土层的分界面,其上下土层的水热参数有着显著差异。在陆面过程模式中准确描述冻融锋面的移动过程将有助于提高其对多年冻土水热过程的模拟能力。本研究首先将Noah-MP陆面过程模式的模拟深度扩展到20 m,并将原模式的4层土层增加到19层土层,同时引入前人的有机质方案和植被根系方案,然后在此基础上,通过耦合Stefan方法以加强模式对冻融锋面的模拟能力,进而探究耦合Stefan方法的Noah-MP模式对西大滩多年冻土站点水热过程的模拟效果。研究中设置了不耦合Stefan方法的CTL控制试验和耦合Stefan方法的STE对照试验来分别模拟西大滩多年冻土站点2012年0~20 m的土壤温度与土壤液态含水量,模拟结果用站点0~3.2 m内10个深度的日均土壤温度、土壤液态水含量监测数据以及3 m、6 m和10 m的年均地温监测数据来做验证。研究结果表明,由土壤温度模拟值插值得到的冻融锋面(0 ℃等温线)有明显阶梯状特征,最大冻融深度与实测相比偏大。耦合Stefan方法增强了Noah-MP模式模拟冻融锋面的能力,使得模式能够基于Stefan方法较好地模拟出冻融锋面的变化趋势和最大深度。同时,改进后的模式整体改善了对土壤温度的模拟效果,使得0~3.2 m各土层土壤温度的平均RMSE降至0.89 ℃,减小44%;平均MBE降至-0.13 ℃,减小86%;模拟的3~20 m年均地温与实测数据更为接近。改进后的模式对土壤液态含水量的模拟水平也有一定改善,其中,模拟的0~3.2 m各土层的土壤液态含水量的平均RMSE降至0.06 m3·m-3,减小33%;平均MBE降至-0.01 m3·m-3,减小67%。模式还较好描述了活动层20 cm、40 cm、80 cm、120 cm土壤的融化时间。可以看出,在陆面过程模式中耦合能够较好模拟冻融锋面移动过程的Stefan方法可较大程度提高模式的模拟能力,是陆面过程模式改进的有效途径之一,本研究的结果可为多年冻土区陆面过程模式改进提供参考。