复杂地形区域平均气温空间插值方法研究.docxVIP

第 ４２卷 第 ６期２０１４年 １２月气象科技ＭＥＴＥＯＲＯＬＯＧＩＣＡＬＳＣＩＥＮＣＥ ＡＮＤ ＴＥＣＨＮＯＬＯＧＹＶｏｌ．４２，Ｎｏ．６ Ｄｅｃ．２０１４复杂地形区域平均气温空间插值方法研究赵平伟１鲁镁２彭贵芬３罗睿１（１ 云南省临沧市气象局 ，临 沧 ６７７０９９；２ 云南省沧源县勐来乡中心校 ，临 沧 ６７７０９９；３ 云南省气象台 ，昆 明 ６５００３４）摘要 采用反距离权重 、克 立格 、径 向基函数法 、三维二次趋势面和地理加权回归５种 方法 ，对 云南１２５个 站１月 、４ 月 、７月 、１０月 和年的 ３０年平均气温进行空间插 值并比较分析发现 ：反 距 离 权 重 、克 立 格 、径 向 基 函 数 法 ３ 种 常 规 插值方法对地形复杂的云南气温空间分布模拟精度不高 ；三维二次趋势面和地理 加权回归模型对云南气温空间分 布模拟较好 ，交叉检验结果显示前者的平均绝对误差（ＭＡＥ）为 ０．４３～１．０２℃ ，均 方根误差（ＲＭＳＥ）为０．６７～１．７７℃ ，其推算的云南气温栅格数据能较好地反映出云南各地气温的分布和差异 ；ＧＷＲ 模型对气温的交叉验证在 ５ 种插值方法中误差最小 ，插 值结果 ＭＡＥ 在 ０．６５ ℃以下 ，ＲＭＳＥ 在 ０．８ ℃以下 ；进 而使用 “地理加权回归模型插值 ＋ 反距离权重残差内插”叠加法对抽取的 １０个检验站平均气温进行插值检验 ，６４％ 的插值结果绝对误差在 ０．５ ℃以 内 ，７４％ 的插值结果相对误差在 ５％ 以内 ，且实测值与插值估算值回归关系决定系数 Ｒ２ 在 ０．９以 上。关键词 气温 空间插值 三维二次趋势面 地理加权回归 复杂地形引言云南地处中国西南边陲，总面积３９．４×１０４ ｋｍ２。 它的地势走向为北高南低 ，西 北最高，东 南最低，由 北向南呈阶 梯 状 下 降，境 内 山 高 谷 深，海 拔 高 差 达 ６６６３．６ ｍ。根据气候带划分，云 南可划分为高原气 候带、温带、北亚热带、中亚热带、南亚热带和北热带 等 ６ 个 气 候 类 型［１］，是 研 究 山 地 气 候 的 理 想 区 域。 目前，气象资源的 空间化已被广泛应用于生态系统 模拟与农业区划中，不少学者致力 于气象要素空间 插值，发现常规的插值模型对地形复杂 、气象站点稀 少地区气象要素空间分布的模拟精度并不高。针对 这一问题，众多学者通过增加相关的辅助地理信息 ， 再对地形复杂的区域进行气象要素空间插值 ，来 提 高复杂地形下的气象要素的插值精度 ：李 正泉［２］等 采用趋势面分析对东北地区 ３０ａ平均降水与 湿度资 料进行空间插值；张 海静［３］等 应用克里格插值法和 建立最低气温与地形的多元线性回归模型 、多 维二 次趋势面模型对青海省无测站地区栅格单元的月最 低气 温 值 进 行 推 算；郭 兆 夏［４］等 通 过 引 入 经 度、纬 度、海拔高度作为 影响因子对陕西及周边地区的年 降水数据分区构建三维二次趋势面模型。王智［５］等 利用统计方法分析新疆多年平均降水量的最佳抽样 间距，建立降 水量与经纬度、高 程、坡 度、ＮＤＶＩ间 的 多元回归方程。此外，吴文玉［６］、仲 嘉亮［７］、李 军［８］、 朱蕾［９］、舒守娟［１０］等人也通 过增加辅助地理信息的 方法对气象要素的空间插值进行深入研究。这些方 法都是尽量将更多影响气温空间变化的因素考虑进 去，如地理经纬度、海拔高度、坡度、坡向、空间变异、 太阳辐射和植被覆盖等。 但是，用于空间插值的气 象要素数据属于横截面数据 ，上述研究多基于最小 二乘的多元线性回归和趋势面方法开展 ，及 直接假 定气象要素横截面单元是同质的 ，而多数情况下，由 于气候资源禀赋、地理位置的不同，同一气象因素在 不同地区可能表现出完全不同的影响关系。采用气 象要素横截面数据与辅助地理信息进行回归研究时 假定截面单元之间的同质 ，使得研究结果由于忽视 了空间效 应 （空 间 相 关性和空间异质性 ）而 可 能 有 偏，进而导致研究 得出的各种结果和推论缺乏应有 的解释力。且已有的文献中研究区像云南这样一省 含多个气候带且受焚风效应 、冷湖效应和屏障效应 等地形影响［１１］明 显 的 复 杂 山 地 气 候 环 境 下 的 气 温ｈｔｔｐ：／／ｗｗｗ．ｑｘｋｊ．ｎｅｔ．ｃｎ气 象科技 作者简介 ：赵 平伟 ，男 ，１９８３年 生 ，学 士 ，工 程师 ，主要从事气象数理统计分析和短期预报工作 ，Ｅｍａｉｌ：７２０３９５６＠ｑｑ．ｃｏｍ收稿日期 ：２０１３年 １２月 ２５日 ；定 稿日期 ：２０１４年 ２月 １２日第 ６期赵平伟等 ：复杂地形区域平均气温空间插值方法研究 １００３空间插值的研究鲜有报道 ；同时，气象要素的空间插 值没有绝对的最