摘要
冰川内部温度与冰面温度变化过程密切相关,因此可由实测冰温-深度剖面重建冰面温度变化.耦合热传导-冰流物理模型及其相关反演算法是基于冰温-深度剖面重建冰面温度变化的理论基础和关键.通过构造理想条件下的数值实验,分析比较了最小二乘法、Tikhonov正则化方法和蒙特卡罗算法重建百年尺度冰面温度变化的实验结果.同时,以我国藏北高原腹地的马兰冰川钻孔资料为例,结合气象资料构造真实钻孔数据模拟实验.分别给出了三种反演算法重建10 a和40 a尺度的冰面温度变化结果,并讨论了不同算法的优劣性和适用性.两组模拟实验结果均表明:Tikhonov正则化方法是目前求解该反问题的最优方法.与其他算法相比,Tikhonov正则化方法在各个时间尺度的重建结果与真实冰面温度变化吻合更好.另外,分别对重建百年和40 a尺度的冰温-深度剖面添加±0.02℃和±0.01℃的随机误差,研究发现Tikhonov正则化方法能够有效降低噪声干扰,相较其他算法得到的解更加稳定.这在一定程度上解决了该反问题中解的稳定性问题.
基金项目
国家自然科学基金(42130516)
中国科学院战略性先导科技专项(XDA20060201)
青藏高原综合科学考察研究项目(第二次)(2019QZKK020102)
NETL
NSTL
维普
万方数据