基于耦合热传导-冰流模型的冰川钻孔温度重建的气候反演算法比较

孙欢, 王宁练, 张泉

Comparison of climate inversion algorithms for glacier borehole temperature reconstruction based on coupled heat conduction-ice flow model

Huan SUN, Ninglian WANG, Quan ZHANG

图4 边界条件分别为 $μ_{1} (t) = 2 s i n ? (2 π t / t_{f})$ 和 $μ_{2} (t) = 2 s i n ? (π t / t_{f})$ 时三种反演算法重建的冰面温度变化［对 $θ_{1} (z)$ 添加 $\pm 0.02$ ℃随机扰动前后的重建结果（a）、（b）；对 $θ_{2} (z)$ 添加 $\pm 0.02$ ℃随机扰动前后的重建结果（c）、（d）］

Fig. 4 Reconstructed changes in glacier surface temperature under the boundary conditions of $μ_{1} (t) = 2 s i n (2 π t / t_{f})$ and $μ_{2} (t) = 2 s i n (π t / t_{f})$ （black line） using the least square method （blue line）， Tikhonov regularization method （red line） and Monte Carlo algorithm （green line）［the input data $θ_{1} (z)$ and $θ_{2} (z)$ without errors，（a） and （c）， the input data $θ_{1} (z)$ and $θ_{2} (z)$ with the random errors $\pm 0.02$ ℃，（b） and （d）］