中国东北地区气象数据空间平稳性检验.pdfVIP

第十二届全国数学地质与地学信息学术研讨会 中国东北地区气象数据的空间平稳性检验· 许家琦1李颜伶2舒红1刘艳3 ‘1．武汉大学测绘遥感信息工程国家重点实验室，湖北省武汉市，430079； 2．青岛勘察涓绘研究院，山东省青岛市，266032；3．中国气象局乌鲁木齐沙漠气象研究所) 摘要：地质统计是以随机场平稳为前提条件，然而在现实中这个前提条件未必成立。因此在统计建模前。 需要检验数据在空间上的平稳性。本文针对中国东北地区的气象数据(日气温，日风速和月降水)，采用 散点图法，趋势线法和邻域圆法探讨了空问数据平稳性检验的问题。 结果日气温数据有南高北低的趋势， 月降水数据有东高西低，南高北低的趋势，日风速数据在空间上平稳。这些结果与实际地理情况捆吻合。 关键词：平稳性：地质统计；散点图法：趋势线法；邻域圆法 O引言 在对空问数据的分析中总是假设空间协方差结构平稳，然而这一假设在实际中常常不成 立。例如降水数据就显示出空间非平稳性‘1·羽，月降水中硫酸浓度也呈现出了非平稳结构【3】． 空间数据的平稳性检验有多种方法。地理加权回归法将空间的概念引入到一般线性回归的回 归参数中，回归系数随着空间位置的变化而变化[41。地理加权回归法检验的是变量之间相关 关系的空间异质性，GWR在已知数据非平稳的情况下分析这种非平稳性，而不是用于数据 检验是否平稳。本文采用的三种方法(三维散点图法，趋势线法和邻域圆法)都是用来检验 属性数据本身的空间异质性。谱方法用局部平稳过程的卷积来表示整个非平稳过程，并用这 个非平稳模型来进行谱密度的参数和非参数估计【5，们。由于每个被视为平稳的局部区域中都 要有足够多的数据来计算空间变异函数然或者空间协方差函数【刀，因此这种方法并不适用于 我国较为稀疏的气象站点数据。三维散点图法拟合的是整个研究区域内属性数据的趋势面； 趋势线法拟合的是投影面上所有属性数据的趋势：虽然邻域圆法也是把空间分割为一些小的 区域来研究，但邻域圆法关注的是局部区域内的二阶矩。均值和方差的计算与经验变异函数 的计算不同，对于参与计算的气象站点没有太多数量上的要求。在非平稳的空间协方差结构 分析中，把不平稳各向异性的地理空间映射到平稳各向同性的频谱空间中进行处理，然后将 处理结果逆变换到地理空间中【8loSG模型计算复杂，物理含义不够直观明确。而本文采用 的三种方法：(1)三维散点图法有着简单直观的优点，通过查看气象数据在三维空间中的数 值大小来初步判断数据是否平稳。(2)在趋势线法中，把气象数据投影到东西和南北两个面 上，在这两个面上分别用多项式法拟合出两条趋势线，从而判断出空间数据在东西方向和南 +通讯作者：舒红 Email：shu_hongelmars．whu．edu．cn 193 第十二届全国数学地质与地学信息学术研讨会 北方向上是否平稳[91。(3)由于随机变量的一阶矩和二阶矩已经足够为大多数实际问题提供 可以接受的近似解决方法【lol，在邻域圆法中通过统计圆内样本数据的均值和方差来检验我 国东北地区气象数据的空间平稳性【¨】。 1研究区域概况 本文采用我国东北地区的气象观测站点数据，包括2006年1月1日的日平均气温数据， 2006年1月1日的日平均风速数据以及2006年5月的月降水数据作为实验数据。气象站点 分布如图1所示。 冬季，受极地大陆气团控制，东北区气候严寒。1月是全年最冷月，全区的月平均气温 均在零下4摄氏度以下。温度梯度之大，为全年之冠。大小兴安岭北段的平均气温在零下 30摄氏度左右，是我国冬季最冷的地方。冬季(以1月为代表)在蒙古冷高压影响下，空 气稳定，风速较小。特别是松花江流域以北地区，因位置最北，常在高压内部控制之下，冬 季为全年风速最小季节。1月份全区平均风速在l-5m／s之间；山区1-3m／s：平原多在3-4m／s。 本区的降水主要集中于夏季。各地降水一般从5月份明显增加，9月份又陡降。降水量自东 南向西北有递减的趋势。降水日多于20天的区域由东北部的山区延伸到了松辽平原的东部， 降水日数的高值区仍是小兴安岭南部和长白山一带112l。 图l中国东北地区气象观测站点分布图 ofChina