基于landsa影像的中国西部干旱区冰川变化及其对气候变化的响应.docxVIP

基于landsa影像的中国西部干旱区冰川变化及其对气候变化的响应 1 博格达峰地区冰川变化对气候的响应 作为自然资源的一部分，缓冲区是中国西部干旱地区农业和畜牧业生存和发展的重要水源，在区域生态环境和经济发展中发挥着重要作用。在全球变暖的大背景下, 特别是21世纪以来, 中国内陆大部分地区冰川呈加速退缩趋势, 昆仑山 (上官东辉等, 2004) 、北天山 (何毅, 杨太保, 田洪阵等, 2013) 、阿尔泰山 (白金中等, 2012) 、阿拉套山 (何毅, 杨太保, 杜鹃等, 2013) 、祁连山 (Wang et al, 2011) 冰川等均处于加速退缩状态。气温和降水的变化会影响冰川的进退, 因此, 深入全面了解冰川区气温和降水的时空变化情况, 已成为探究冰川变化对气候响应机制及其预测冰川变化的首要条件。 博格达峰是东部天山第一高峰, 两侧盆地冬季受强大的内蒙古冷高压控制, 寒冷而干燥。峰顶冰川积雪终年不化, 温度常在冰点以下, 有少量的水也都冻结。夏季高空盛行西风环流, 是典型的中纬西风带大陆型气候;气温年较差和日较差大, 降水量少而变率大, 且集中分布在夏季 (伍光和, 上田丰等, 1983) , 占全年降水的77%左右。博格达山脉是最早引起中外学者注意并进行考察的冰川区之一 (伍光和, 张顺英等, 1983) , 其冰川融水是白羊河流域、吐鲁番—哈密盆地和柴窝堡盆地可靠的补给水源。该区的现代冰川主要位于南坡的古班博格达河、黑沟和阿克苏河、北坡的三工河、四工河、甘河子沟及白杨河等河流源头。近10年来, 冰川的快速萎缩已经影响到吐鲁番坎儿井的水量 (李开明等, 2011) , 因此深入了解近年博格达峰地区气候变化时空特征及冰川变化的响应特征具有重要的理论和实践意义。基于此, 本文选用天山东部博格达峰地区周边共3个站点1960-2013年气温和降水数据, 从趋势、突变、周期等视角分析近半个世纪以来博格达峰地区气候变化的基本特征, 并结合1972、1990、1999和2013年4期冰川数据, 探讨该区冰川进退对气候变化的响应特征, 以及天山东部复杂的气候变化与冰川波动的变化特征, 为深入探究该区大陆型冰川变化对气候变化的响应机理奠定基础。研究区示意图如图1。 2 数据和方法 2.1 数据来源 2.1.1 中国地面气候资料 气象数据来源于中国气象科学数据共享服务网 (/) 的中国地面气候资料数据集。为对研究区作100 km缓冲区分析, 选取了缓冲区内具有代表性的达坂城、奇台和乌鲁木齐3个气象站点1960-2013年的年均气温、年均降水数据, 具体见表1。 2.1.2 遥感数据 2.2 气候变化分析 2.2.1 显著性及突变情况 采用Mann-Kendall (M-K) 法研究气温和降水量变化趋势的显著性及突变情况 (易湘生等, 2011) , 其公式如下: 式中:Zc为秩相关系数;τ为检验统计变量;p为数据序列所有对偶观测值 (Rj, Ri, ji) 中RjRi出现的次数;n为序列长度;σr2为方差。 2.2.2 小波函数的获取 小波分析是一种时频多分辨率分析方法, 是时间和频率的区域变换, 能有效的从信号中提取信息, 这里简单采用Morlet能谱小波函数, 具体方法如下: 小波变换系数为: 式中:ξ (t′, a) 是小波系数;f (t) 是时间序列;Ψ*是Ψ的共轭函数 (魏凤英, 2007) ;k为系数;Φ为相位;a为小波变换的尺度因子 (伸缩因子) 。 2.3 变化检测方法 2.3.1 冰的人工提取 利用遥感影像提取冰川边界有监督分类、非监督分类、雪盖指数和比值阈值等多种方法 (王高峰等, 2010) 。Raup等 (2007) 研究表明, 人工解译是目前提取冰川信息最为精确的方法。刘宇硕等 (2013) 进行了Landsat影像在冰川面积信息提取的精度评价, 也得出目视解译最为精确的结论。由于本文研究区冰川分布较少, 故采用了人工目视解译的方法提取冰川边界, 结合冰川编目数据和Google Earth三维图像, 目的是为了冰川边界提取判读更准确。对提取的冰川边界转换成Albers等积投影, 最后利用Arc GIS软件统计面积及其变化。运用更高分辨率影像人工提取冰川面积验证Landsat影像的精度是较常用的方法 (Paul, 2000;Bolch et al, 2012) 。Paul (2000) 对Spot影像人工提取的冰川面积和TM影像得到的面积进行对比分析发现, 误差小于1%。本文中运用此种方法, 下载2013年8月Google Earth影像 (分辨率为9.5 m) 进行验证, 发现误差为±2%。 2.3.2 川面积、比值 冰川面积变化速率 (percentages of area changes, PAC)