基于GPS和DEMETER卫星数据的地震电离层电子浓度异常变化研究

闫相相

（中国地震局地质研究所，北京100029）

基于GPS和DEMETER卫星数据的地震电离层电子浓度异常变化研究

闫相相

（中国地震局地质研究所，北京100029）

地震电磁现象的观测与研究是认识地球内部物质特性、地震孕育发生过程以及空间信号传播和介质特性变化的重要途径。空间对地观测技术的发展，推动了地震电磁信号在大气层、电离层的传播与耦合机制的研究。从20世纪60年代发现美国阿拉斯加大地震震中区上空出现电离层异常扰动现象以来，有关地震空间电磁、电离层扰动的现象引起了广泛的关注，成为地震电磁学领域一个新的研究热点。研究表明，在地震孕育发生过程中，地下介质由弹性变形进入非弹性形变阶段，原有微裂隙扩展、新生微裂隙发育并定向排列，体积膨胀，水溶液进入和重新分布等物理过程，产生压磁、压电、动电、感应电磁、热磁效应等物理电磁效应的同步和定向激发，从而引起在距离震源一定范围内可观测的明显地震电磁前兆现象。大量的震例和统计分析显示，大地震和火山喷发前几天到几小时内，对应孕震区、火山喷发区上空大气层、电离层存在电磁异常扰动现象。主要包括：①LF／VLF／ELF／ULF信号相位、振幅异常变化；②电离层电子浓度总含量（total electron content，TEC）和电离层F2层临界频率foF2异常；③电离层等离子体参数变化；④高能粒子扰动；⑤热红外辐射异常等。

全球定位系统（Global Positioning System，GPS）作为当前探测电离层的有效手段已经得到广泛应用，通过其获取的电离层TEC则是揭示电离层时空变化规律和特征的重要物理参量。2000年以来，俄罗斯、美国、法国、乌克兰等国先后发射了地震电磁卫星，尤其是法国DEMETER卫星，其科学数据被广泛地应用于地震电离层扰动研究，取得了极大进展。中国的地震电磁监测试验卫星（China Seismo-Electromagnetic Satellite，CSES）计划已经列入“十二五”立项。

本论文依托地震电磁立体观测体系发展思路，配合中国地震电磁卫星计划的推进，开展了基于DEMETER卫星和GPS数据的地震空间电磁电离层扰动现象的研究。本文研究内容主要可分为3大部分：①研究GPS测量电离层TEC的原理，掌握电磁卫星数据处理方法和异常提取技术；开展DEMETER卫星各载荷数据处理及异常提取方法研究。②分别针对汶川地震、日本地震和玉树地震开展了震前电离层扰动的多参数分析，研究和对比了3次强震前电离层异常的特征；分别利用DEMETER卫星数据和GPS数据统计分析了全球7．0级以上强震和中国境内6．0级以上强震期间的电子浓度异常变化规律；针对我国西南地震多发区，开展伴随着多次地震的电离层TEC长时间序列分析，深化对地震电离层异常扰动的研究。③调研和总结目前地震电离层耦合模拟的最新成果，详细阐述孕震区扰动电场（seisomogenic electrostatic field，SEF）的计算过程，分析其对电离层的效应。

本文研究工作获得以下主要认识：

（1）汶川地震的电离层扰动多参数分析结果表明，震前多种参量都出现了一定程度的扰动现象，不同参数在时间和空间上的扰动特征并不一致，但这些参数之间存在相应的联系。

利用GPS TEC、DEMETER卫星和NOAA／AVHRR卫星的观测数据，系统梳理了汶川地震前多种参数（包括电离层TEC、Ne、Ni、Tb、Ti以及高能粒子通量等）的变化特征，探索总结了汶川地震前大气层-电离层响应。结果表明，电离层电子浓度在5月6—10日有连续的负异常现象，5月9日则明显增强；异常的范围主要位于震中南部区域，南北向约1 100～ 1 670 km，东西向约1 600～3 700 km，扰动的峰值距震中约700～900 km。NOAA／AVHRR卫星数据得到汶川地震前5月7日和8日出现了显著的热红外亮温异常，而DEMETER卫星的原位探测则发现5月9日离子温度有明显增高，幅度约12．5%；以上异常区域都位于震中西北，范围约100～300 km。此外，DEMETER卫星IDP探测器高能粒子数据结果显示5月6日夜间100～600 KeV能谱段出现明显增强，增幅达6σ；这与前人得到的异常电场增强时间相一致。研究结果表明，不同参数在时间和空间上的扰动特征并不一致，但这些参数之间存在相应的联系，因此有必要联合地基和天基手段共同观测。

（2）日本地震和玉树地震的电离层电子浓度变化研究结果表明，震前TEC都出现了不同程度的增强现象，而异常同时可能受到太阳或地磁活动的共同作用。

利用IGS（International GNSS Service）提供的TEC数据和日本NICT提供的电离层测高仪数据，分析了2011年3月11日日本MW9．0特大地震（Tohoku-oki earthquake）前上空的电离层变化。研究发现，震前3天即3月8日04～14时（UT）期间，赤道异常区电离层TEC出现了明显的增强现象，异常最大幅度达40TECU，并伴随有南半球磁共轭区增强现象。日本境内的4个电离层测高仪数据分析显示，3月8日4个台站同时出现了foF2增强现象，其中Okinawa和Yamagawa 2个站增幅明显。需要指出，3月6—12日太阳活动水平有所增强，可见太阳和地磁活动对3月8日大幅度的TEC扰动有所贡献。然而，分析表明，仅仅是太阳活动增强无法造成如此强度的电离层扰动，即无论从持续时间、空间分布特征以及异常强度等方面都表明，以上TEC变化跟此次大地震有很大的关系。

同时，综合GPS和DEMETER卫星2种观测手段研究了玉树地震前的电离层等离子体参量变化情况。结果表明，GPS观测的电离层TEC在震前20 h左右有明显的增强现象，TEC增强的区域位于震中南部15°N～30°N范围内，扰动的峰值距震中约1 000 km。TEC增强区域显示为典型的EIA增强现象，与汶川和日本地震前观测到的TEC空间扰动特征较为一致。通过对DEMETER卫星IAP和ISL载荷记录的电离层等离子体参数时序变化和空间纬度变化分析，没有发现明显的跟地震相关的扰动现象。以上也表明了2种观测手段由于测量电子浓度的方式不同，其观测结果存在一定的差别。

（3）DEMETER卫星数据统计全球强震期间电子浓度Ne变化结果表明，在赤道和中低纬地区更容易观测到明显的Ne扰动，震级大小则与扰动幅度基本呈正相关。基于GPS数据全国强震期间TEC变化统计分析结果则表明，TEC异常扰动多集中在震前2～6天，同样以增强为主，且集中在地方时12～20时。

基于法国DEMETER卫星ISL升轨数据，统计分析了2005—2009年全球37个M≥7．0地震前后的电离层电子浓度变化。结果发现，共有19个地震（51%）前观测到了较为明显的电子浓度异常扰动现象，其中大部分表现为异常增强；观测到的电子浓度异常现象一般出现在震前1～5天内，也有部分地震前出现2次及以上扰动现象。分析表明，在赤道和中低纬地区更容易观测到明显的电子浓度扰动，即异常震例基本发生在纬度±40°以内（18次），更是以±20°内居多（13次）。此外，震源深度对电子浓度扰动的影响并不明显，而震级大小则与扰动幅度基本上呈正相关；震前出现的电子浓度异常有时会受到地磁活跃的共同影响，此时的扰动幅度一般较大。

基于GPS数据统计分析了2000—2010年期间中国境内（包括台湾地区）30次MW6．0以上地震的TEC异常扰动情况，结果表明，20次地震（66．7%）前观测到了较为明显的TEC异常扰动，扰动以增强为主；扰动多集中在震前2～6天，而且具有明显的地方时特点，扰动一般集中在地方时12～20时，一般平均持续4～6 h，有的地震可以持续几天。震前电离层TEC异常的空间变化主要表现为对电离层赤道异常双峰（EIA）的影响上，不但会增强或降低峰值位置的电子浓度值，而且会造成驼峰最大值位置向磁赤道方向偏移。

（4）通过对中国西南地震多发区以及“检验区”电离层TEC变化长时间序列分析结果表明，电离层是一个复杂的系统，受到太阳、地磁以及地震活动等因素影响，其扰动具有多源性。

通过中国地壳运动网络提供的GPS观测数据，获取了高精度电离层TEC分布，采用滑动四分位法分析研究了中国西南区域2008年4—10月（太阳和地磁活动平静时段）6次连续MW6．0以上地震期间孕震区电离层TEC长时间变化及其异常分布；并在此基础上利用GIM（Global Ionosphere Maps）数据分析了全球TEC变化特征。鉴于电离层主要受到太阳和地磁等空间天气的影响，我们将TEC变化与太阳EUV（extreme ultraviolet）辐射、行星际磁场（interplanetary magnetic field，IMF）南向分量Bz以及地磁活动指数Dst和Kp进行了比较。研究发现，该时段内电离层TEC异常扰动与太阳和地磁活动有很好的相关性；而除汶川地震外，其他地震前没有发现明显的跟地震相关的TEC异常扰动现象。同时，对比分析了与上述研究区位于同一地磁纬度的“检验区”（30°～50°E，15°～35°N）的GPS TEC随时间变化和异常分布情况，结果显示，TEC异常分布的时空特征与研究区域较为一致。结果表明，电离层是一个复杂的系统，其扰动具有多源性。

（5）通过对地震电离层耦合机理研究初探，认识到电离层TEC在SEF的作用下出现明显变化，在赤道向（longitude equatorward）较为明显。

介绍了电离层电流和电导率张量以及电离层O＋光化学作用的数学物理方程，总结了孕震区扰动电场的数学模型及公式推导。结果表明，夜间扰动电场更容易传播到电离层，这与观测现象相一致。扰动电场向上传播至电离层F层高度，导致电离层电子浓度变化。结果显示，电离层TEC在SEF的作用下出现明显变化，而在赤道向上更为明显，这同样也与观测现象一致。

地震-电离层；全球定位系统（GPS）；DEMETER卫星；电子浓度总含量（TEC）；异常信息；统计分析；孕震区电场

（作者电子信箱，闫相相：ahuyanxiang＠gmail．com）

P315．72；

A；

10．3969／j．issn．0235-4975．2014．07．013