12号冰川冰内结构特征分析.pdfVIP

第２６卷　第６期 地球科学进展 Ｖｏｌ．２６　Ｎｏ．６ ２０１１年６月 ＡＤＶＡＮＣＥＳＩＮＥＡＲＴＨＳＣＩＥＮＣＥ Ｊｕｎ．，２０１１ 文章编号：１００１８１６６（２０１１）０６０６３１１１ 祁连山老虎沟１２号冰川冰内结构特征分析 １ １ １ １，２ 武　震 ，张世强 ，刘时银 ，杜文涛 （１．中国科学院寒区旱区环境与工程研究所，冰冻圈科学国家重点实验室，甘肃　兰州　７３００００； ２．中国科学院寒区旱区环境与工程研究所冰冻圈科学国家重点实验室，祁连山冰川与 生态环境综合观测研究站，甘肃　兰州　７３００００） 摘　要：利用探地雷达（ＧｒｏｕｎｄＰｅｎｅｔｒａｔｉｏｎＲａｄａｒ，ＧＰＲ）调查冰川冰内结构和冰层厚度是一种监测 冰川变化常用手段。应用麦克斯维（Ｍａｘｗｅｌｌ）方程的二维时域有限差分 （ＦＤＴＤ）方法，通过将模拟 图像与老虎沟１２号冰川的实测图像的对比，分析了雷达剖面中的几种反射特征，如冰裂隙、融洞、 暖冰等。对比结果表明，ＧＰＲ的反射图像会受到诸多因素的影响，除了介电常数，还包括介质的厚 度，形态，粗糙度等因素的影响。探测结果表明，具有极大陆型特征的老虎沟１２号冰川的局部的冰 内温度呈现出较高的态势。同时，ＦＤＴＤ模拟也为我们未来利用ＧＰＲ来监测冰川随气候的变化提 供了很好的分析手段。 关　键　词：探地雷达；老虎沟１２号冰川；时域有限差分；模拟 中图分类号：Ｐ９１４　　　文献标志码：Ａ 译而解译难度取决于地下介质的结构或介质复杂复 １　引　言 杂，以及入射波产生的反射波，绕射波和噪音的叠加 ［１］ 影响。因此，开发辅助方法来提高解译精度就显得 冰川被认为是气候变化敏感的指示器 ，近一 个多世纪以来的升温，已引起包括两极冰盖在内的 十分必要。 ［２］ 由于时域有限差分法（ＦｉｎｉｔｅＤｉｆｆｅｒｅｎｃｅＴｉｍｅ 全球范围陆地冰的强烈退缩 冰川规模变化伴随 着冰内或冰下过程变化，这一变化反过来又会对冰 Ｄｏｍａｉｎ，ＦＤＴＤ）具备简单使用和计算速度快的优 ［１２］ ［１２～１９］ 川波动幅度带来影响，这其中包括消融机制及其与 点 ，其在电磁波模拟中得到广泛应用 ，虽然 ［３］ 有很多人将这种方法应用于实测雷达图像的辅助分 冰内冰下水文特征的关联 ，因而要一定的技术和 方法探测冰内／冰下结构。一些研究发现，伴随冰川 析，但将 ＦＤＴＤ方法应用在冰川研究中的例子还 ［４］ 较少。 退缩，冰内结构和冰川性质正在发生变化 如大陆 性冰川内部含水量呈增加的迹象［５～７］，表明气候变 为了研究典型冰川对气候的响应机理，本文利 暖无论是对冰川规模还是冰体性质都带来了强烈影 用ＦＤＴＤ的２Ｄ模拟对老虎沟１２号冰川的局部冰内 响。ＧＰＲ因对冰有良好的穿透性和易于对冰内和 结构特征进行了模拟验证。通过将模拟图像与实测 冰下介质埋藏深度和类型进行辨别［５，６］，而在冰川 图像的对比，很好地还原了冰内结构特征，为认识冰 研究中已有较长的应用历史［８～１１］。ＧＰＲ应用最关 川内部特征提供了方法基础。在本文中，除