基于地统计方法新江北部雪深空间插值研究.pdfVIP

第十二届全国数学地质与地学信息学术研讨会 基于地统计方法的新疆北部雪深空间插值研究 胡丹桂u。舒红“。刘艳’ (1．武汉大学测绘遥感信息工程国家重点实验室，湖北武汉，430079；2．武汉职业技术学院计算机技 术与软件工程学院，湖北武汉，43∞74；3．中国气象局乌鲁木齐沙漠气象研究所．新疆乌鲁木齐，踟002) 摘要：基于地统计的插值方法有许多，对于目前众多可以提供气候要素的空间插值方法中，具体探讨了 普通克里金法和协同克里金法。以新疆北部区域42．Oo~49．20N’80．d吒91．daE为试验区，利用北疆地区 1969--2008年中每年的1月、2月、12月气象站的雪深观测数据的40年平均值。采用地质统计学普通克里金和 协同克里金空间内插方法，估计试验区的雪深分布．结果表明，内插结果与本区的实际积雪深度的地理分布 趋势基本一致，雪深最大值主要分布在高山一带，雪深最低值分布在海拔低的盆地一带。同时，针对40年 的平均雪深，引入了高程作为协变量的协同克里金的插值方法比普通克里金的插值方法的均方根误差 RMSE值降低T4．25％，预测值与实测值之问的相关系数提高了6．87％。因此，引入了高程作为协变量的协 同克里金插值方法提高了插值精度． 关键字：地统计学；协同克里金；插值；高程；雪深 l引言 中国西北干早地区重要的水资源之一是新疆地区的积雪，积雪也是河川径流的重要补给 来源，也是农业、畜牧业等经济发展的生命线¨1。新疆地区地形较为复杂，阿尔泰山、天山 地势高耸，两疆盆地地区地势低平，冰川、绿洲、沙漠等地貌交替呈现。本文主要分析天山 以北地区，在天山以北的阿尔泰山地区海拔在3000米以上，温度较低，在迎风坡，由于气流 的爬坡过程中水汽凝结，成为我国积雪分布最多的地区之一，也是积雪异常变化的敏感区12I。 然而由于经济和人力原因，积雪观测站点数量极其有限且分布极为不均匀啪。地统计学起 先在地质、矿业、石油、土壤等方面的应用比较广泛，随着研究的不断深入，它的应用领域 不断拓宽和深入。目前已发展了许多地质统计学方法【4l，如泰森多边形嘲、克里金(Kriging) 内插法【6】、反距离加权平均法ITl、趋势面分析法is]、多项式回归法【91等。利用有限的积雪站 点来估算未知点的积雪。这些方法只是局限于积雪站点的观测值，没有考虑到地理位置、地 形地貌、大气环流等因素对积雪深度及时空分布的影响。然而，天山以北的阿尔泰山区气候 资助项目：地理空间信息工程国家测绘地理信息局重点实验室开放研究基金资助项目。气象观测数据分析 的时空统计方法研究”资助(编号201329)；国家自然科学基金项目。时空交互的统计建模”资助(编号： 作者简介：胡丹桂(1980)，女(汉族)，浙江永康人，讲师，博士研究生，主要研究GIS及时空统计分 析，E-maih danguihu@163．com 通讯作者：舒红(1970．)，男(汉)，湖北武汉人，教授，博士，博导，主要研究时空统计及数据同化， E-maih shu．_hong@、vhu．edu．C11． 215 第十二届全国数学地质与地学信息学术研讨会 垂直梯度变化明显，冬季严寒，山谷中少雪而高山地带大雪。塔城市境内地形较为复杂，北 高南低，北部高山连绵、丘陵起伏，南部沼泽、洼地星罗棋布。因此，在这种地形变化较为 复杂的地区，提高积雪的空间插值精度，高程的作用不容忽视。为此，本文将结合高程，采 用协同克里金方法对积雪的深度进行插值，并对插值结果的精度与普通克里金插值方法得到 的结果进行对比分析与探讨。 2研究区域概况及数据预处理 2．1研究区域 本文研究区域位于新疆北部地区(42．00~49．20N，80．oo~91．00E)，地面观测台站所观测 的资料来自1969,-．2008年新疆地区逐日积雪深度。由于气象台站观测资料不可避免存在缺测， 本文剔除缺测较多的台站。共选用资料较好的29个台站，站点分布稀疏且不均匀。 地理数据选用lkmxlkm的北疆数字高程模型(DEM)数据。如图l所示。 且lon’00F11400’O“￡11700’O“￡9000,0“E 图1试验区内气象站点分布图 2．2数据预处理 本文以塔城地区和阿尔泰地区两个研究区为中心，由于研究区