采用黄河源区2012年8月—2013年7月的日尺度地温观测数据,运用含有水相变化的地下水运动的水-热耦合原理,利用SUTRA-ice软件建立了模拟冻土在冻融过程中温度场分布的二维模型,在校正模型参数的基础上,预测了未来情景下黄河源区多年冻土活动层厚度的变化情况。结果表明:黄河源区冻土活动层浅层地温最高为14℃,最低为-9.5℃,分别出现在7月、1月;土壤温度变幅从浅层20 cm到深层150 cm不断减小且存在滞后现象,与20 cm处的冻结日期相比,50、90、120、150 cm处分别滞后9、24、41、51 d;在假定升温率为0.04℃/a的情景下,黄河源区冻土活动层厚度由3.90 m增加到20 a后的5.86 m,增加率为9.8 cm/a。
鉴于黄河源区实测季节冻土最大冻结深度资料极其匮乏,基于1996~2008年黄河源区及其周边气象站点季节冻土数据,分析了季节冻土最大冻结深度与负积温的关系及时空变化规律,建立了最大冻结深度估算公式,对气温空间插值并由最大负积温和高程估算季节冻土最大冻结深度。结果表明,在季节冻土的迅速发展期,冻结下界深度与负积温呈线性相关;最大冻结深度和最大负积温相关性显著;最大负积温可反映季节冻土随时间的变化趋势;黄河源区仅阿尼玛卿山及其周边存在山地多年冻土,季节冻土最大冻结深度由西北向东南方向递减;气温较高年份中东部地区最大冻结深度明显变小,最大冻结深度小的地区对气候变暖更加敏感;1996~2007年间,最大冻结深度随时间推移呈减小趋势。
黄河源区位于青藏高原多年冻土区东北部边缘地带,是季节冻土、岛状多年冻土和在大片连续多年冻土并存地带.多年冻土层在垂向分布上有衔接状和不衔接状两大类.不衔接状又可分为浅埋藏(8m)和双层多年冻土等形式.从20世纪80年代以来,源区气温以0.02℃.a-1增温率持续上升,人类经济活动日益增强,导致冻土呈区域性退化.多年冻土下界普遍升高50~80m,最大季节冻深平均减少了0.12m,浅层地下水温度上升0.5~0.7℃.冻土退化总体趋势是由大片状分布逐渐变为岛状、斑状分布,多年冻土层变薄,冻土面积缩小,融区范围扩大.部分多年冻土岛完全消失变为季节冻土.
近40年来,黄河源区地温长期处于增温状态,多年冻土出现表层融化。冻土退化后,土壤含水率减少,导致寒区径流的变化。所以研究冻土为寒区水资源估算提供了重要的依据。利用土壤水分迁移方程和热传导方程,采用中心差分格式并在一定假设前提下建立起冻土水热耦合迁移数学模型。对黄河源区站点的冻融深度进行模拟,结果表明:计算值和近2 a冻土的上、下限实测值吻合较好。由此表明采用的数学模型和数值方法是合理、可信的,在此基础上讨论了土壤冻融过程的一般规律。
在对不同植被覆盖度(95%,70%~80%,40%~50%和10%)下的土壤水分(θV)和土壤温度(Ts)进行日观测的基础上,研究了冻融过程中植被覆盖变化对土壤水分分布和温度的影响.土壤温度与水分关系的回归分析表明:土壤冻融过程明显受植被覆盖变化的影响,植被覆盖变化还导致土壤水分和温度的耦合变化.使用了一个土壤水分和温度的耦合模型来研究植被覆盖的影响,结果证明了这一方法的有效性.结果表明:土壤水分对土壤温度的变化范围和幅度都有影响,高盖度下的土壤比低盖度土壤持水性强.此外,在冻结过程中,由于水的热容量大于土壤的,高盖度土壤能够抑制土壤温度的降低幅度,高盖度土壤具有较好的绝热功能.对于黄河源区不同植被类型覆盖下季节冻土冻融过程中的土壤温湿空间变化研究有利于为高寒冻土地区冻土和生态环境的保护及合理利用提供科学依据.
黄河源区由于近年来气候变化的影响,打破了高寒植被与冻土环境之间稳定的适应性关系,由此引发了一系列生态环境退化的现象.在黄河源区多年野外工作的基础上,定量分析了冻土与植被之间的关系.研究表明:多年冻土埋深通过影响浅层土壤含水量影响植被生长的,多年冻土的埋深与浅层土壤含水率和植被的覆盖率具有良好的相关性规律.冻土埋深2 m时,冻结层上水水位低、补给量少,冻结层上水水量小,毛细上升高度不能达到植被根系分布的浅层土壤中,植被生长环境干旱化,多数植被生长受限制,这时只有少量根系发达的耐旱植被存活,覆盖率小,一般不超过35%.因此,2 m的多年冻土埋深为“生态冻土埋深”.近20 a来,黄河源区地温长期处于增温状态,多年冻土出现表层融化,形成深埋的或少冰的冻土等现象;部分地带完全融化消失,连续多年冻土变成不连续冻土或岛状冻土.多年冻土退化后,土壤含水量减少,导致植被物种更替、“黑土滩”等退化现象.