华南电离层ec扰动特征分析.docxVIP

华南电离层ec扰动特征分析 1 基于gps的tec测量 电离层扰动主要是指电离层中离子浓度的变化。由于它直接影响地面和卫星的通信，已成为研究的重点。除了太阳活动、地磁运动和子时线沿线移动的中性风外，电离层的扰动还与中强地震有关。例如，1985年3月27日美国发生的嘉陵地震和1995年1月15日日本发生的千冈地震。2008年5月12日，中国四川四川汶川发生了大地震。根据这些地震案例研究，在中和强震易发阶段，岩石圈和电离层之间存在可能的耦合关系。此外，日本的日常饮食也可能导致电离层扰动。 反映电离层扰动的主要参量有,电离层F2层最大等离子体的频率(f0F2),偶发E层频率(f0Es)以及电离层电子浓度总含量(Total Electron Content,TEC).由于GPS卫星探测电离层技术的应用和发展,基于GPS的TEC测量已成为研究电离层的一个重要方法.特别是在全球低纬赤道区域,白天电离层TEC按纬度分布于地磁赤道,呈现赤道异常(或称Appleton异常),这种赤道异常特征经Zakharenkova等研究,发现在地震前几天会显著变化,其峰值在地震前会向地磁赤道漂移. 华南地区地处低纬赤道区域.位于该区域内的台湾海峡地震最为活跃,由地震引起的电离层扰动也较明显.处于华南地区的广州其地理纬度和经度分别为22.5°—24°N,113°—114.5°E,并建有位于广州五山、番禺、花都、从化、增城、南沙、萝岗的8个地面GPS基准站与中心处理站构成的GPS-TEC监测网,如图1(a)所示.该监测网的空间监测范围可达到10°—35°N,100°—125°E,如图1(b)所示.因此,利用广州GPS-TEC监测网,可以对华南地区的电离层进行研究,探讨电离层扰动的各种可能的机理.需要指出的是,电离层赤道异常的北驼峰(地理纬度25°N左右)位于该监测范围内,该北驼峰对于来自各种扰动源的电场反应强烈.因此,可以通过探测北驼峰附近电离层TEC的变化研究电离层扰动和异常源的关系.表1给出了2009年6月1日至7月30日之间,华南地区及其邻区发生的5个5级以上地震,除此之外,还发生过一次日全食事件.为此,本文利用2009年6至7月广州GPS-TEC监测网提供的TEC数据,初步探讨了华南地区电离层TEC扰动与发生的地震以及日全食之间可能的耦合关系. 2 网格划分和数据统计 广州GPS-TEC监测网提供了2009年5月20日至7月27日(6月30日及7月28至30日之间由于仪器故障等原因数据缺失)的TEC原始数据,其数据采样率为15 s.分析所用数据为VTEC,即垂直TEC,这里定义V为VTEC,T为TEC,α为星下点处卫星的高度角,则有V=Tcosα. 对转换后的VTEC数据继续采样,采样率为30s,并去除卫星倾角小于30度的测量数据,减小由电离层水平梯度引起的测量误差.然后对所采得的VTEC数据进行时序分析和空间异常分布分析.在进行时序分析时,时间分辨率为15 s,即每天选取24×4=96个不同时刻,使用滑动窗口进行数据处理.同理,进行空间异常分布分析时,将图1(b)所示的广州GPS监测网的监测范围按1°×1°的分辨率进行网格划分,并在每个网格点上使用滑动窗口.本文所用的滑动窗口其处理数据来自前15天当中8个站所有监测到的数据.所谓滑动窗口,是动态利用前15天当中8个站所有监测到的数据计算每天96个时刻,以及计算空间分布中各节点上的平均值μ、标准差σ及上、下边界值.上边界值定义为μ+1.34σ,下边界定义为μ-1.34σ.以此为基础去分析每天各时刻或空间分布各节点上的异常值ΔX.如果满足μ-1.34σXμ+1.34σ,则ΔX=0:如果满足Xμ+1.34σ,则ΔX=X-(μ+1.34σ):如果满足Xμ-1,34σ,则ΔX=X-(μ-1.34σ).计算得到的正、负异常用来反映电离层TEC的扰动.这种分析和统计方法参见文献. 由于强太阳活动和强地磁活动能引起比地震更强的TEC扰动,因此,有必要考虑相应时间段内的F10.7和Kp地磁指数分布.图2是2009年6至7月太阳射电流F10.7和地磁Kp指数日分布图(数据来源于NOAA).从图2(a)可以看出,2009年6至7月F10.7分布均处于70 sfu (太阳通量单位),显然这段时期内太阳活动是平静的.图2(b)显示,2009年6月24日、25日和29日三天地磁指数达到4以上,其他时间段地磁Kp指数在1～2之间.而图2(c)指出,2009年7月22日地磁Kp指数达到5以上,其他时间地磁Kp指数在3以下.这说明除2009年7月22日外,6-7月太阳活动和地磁活动是平静的,对电离层TEC扰动影响不大. 3 tec数据分析 3.1 波前tec的扰动 图3给出了华南地区电离层2009年6月15日至7月26日VTEC和ΔV