2. 中国地震局地震预测研究所, 北京 100036
2. Institute of Earthquake Science, China Earthquake Administration, Beijing 100036, China
近几年来,探索地震与电离层扰动关系成为研究地震监测热点之一,很多与地震相关的电离层扰动现象的文章陆续发表,得到的主要结论有电离层电子温度、氧离子浓度、等离子体、TEC、F2层临界频率等参量在地震发生前后,地震活动区的上方确实存在异常现象,电离层受影响区域取决于震级,受 影响变化最大的电离层区域与震中在电离层的垂直 投影位置不一致(张红波,2011; Pulinets and Ouzounov, 2011;Namgaladza et al., 2012).但地震电离层扰动的机理解释以及如何区分不同扰动因素对电离层产生的作用,还有待进一步研究.
目前,很多学者的研究工作专注于震例的分析,包括单个地震的震例研究和多个地震的统计分析. 单个震例研究方面,在汶川MS8.0地震(Akhoondzadeh et al., 2010; 余涛等,2009; Liu,2009)、玉树MS7.1地震(熊晶等,2010)、日本MS9.0地震(Akhoondzadeh,2012; 姚宜斌等,2012)等强震前观测到TEC、Ne、NmF2发生异常扰动.多个震例统计分析方面,Liu等(2010)使用15天滑动平均的方法统计分析了2001—2007年发生在台湾地区的150个MS≥5.0地震,结果表明,地震前后几天GPS TEC出现扰动现象.Akhoondzadeh(2010)使用四分位法研究了强震前的TEC变化,结果显示震前TEC的异常幅度依赖于地震震级.Le(2011)使用标准差统计分析了2002—2010年736个MS≥6.0地震,结果表明GPS-TEC震前发生异常扰动,大震和浅源地震的扰动幅度更强.关于汶川地震,目前的研究没有涉及对GPS TEC进行相关分析和功率谱密度分析,对此本文使用相关技术和Pwelch法计算功率谱密度重点分析了汶川地震前后GPS TEC在时域和频域上的变化,同时,使用COSMIC数据分析NmF2在汶川地震前后空间分布变化,以深化认识.
2 数据分析方法
2.1 数据 2008年5月12日14时28分四川省汶川县发生MS8.0特大地震,震中位于103.4°E、31°N,震源深度约14 km.本文TEC数据使用震中周围的DLHA、XNIN、YANC、LHAS、XIAG、LUZH 地基GPS观测台2008年3月1日—5月20日数据,观测站与震中分布如图 1.
COSMIC数据使用汶川地震震前三个月至震后一个月NmF2数据,按kp≤3筛选数据,考察2月12日至6月11日11 ∶ 00—14 ∶ 00(地方时)震中周围区域NmF2空间分布变化.
2.2 分析方法
2.2.1 调和分析方法 调和分析法基于傅里叶级数及傅里叶变换,是将满足一定条件的周期信号s(t)分解为基本波形叠加的方法(韩鹏等,2009;余优等,2009):
使用相关分析方法处理TEC拟合数据.相关分析方法是研究变量之间是否存在某种依存关系,并对具有依存关系的变量探讨其相关程度的一种统计分析方法(Pulients et al., 2004):
使用welch法计算震中周围GPS台站TEC的功率谱密度.该方法可以选择适当的窗函数,并在分段时,允许各段之间有重叠,减小方差,提高数据使用率.窗函数本文选择使用海明窗(Barbe et al., 2010):
关于电离层参量空间分布参考值的计算和地震异常判断的方法有很多种.Liu以COSMIC震前7到22天数据作为参考值,Kakinami(2010)以COSMIC 经验模型值作为参考值研究震前电离层NmF2分布变化.本文使用12阶次球谐函数计算的模型值作为考察地震前后NmF2空间分布变化的 参考值,球谐函数数学公式(章红平,2006)如下:
本文使用调和分析方法拟合数据讨论TEC汶川地震前后变化特征,首先对TEC数据做地方时校正.然后,计算2008年3月1日—5月20日TEC拟合数据,局部显示4月25日—4月30日TEC原始数据和拟合后数据如图 2,可以看出,3次谐波能较好的反映TEC日变形态,且可以剔除掉原始信号中频率相对较高的干扰信号,降低干扰信号的影响.
再使用相关分析方法处理TEC拟合数据讨论TEC在时域上的变化,先计算出每个台站计算时间段的时均值数据表示为 =[(1),(2),(3),(24)],然后计算出每天数据时均值,TECi=[TECi(1),TECi(2),TECi(3),…TECi(24)],i=1,2,…,81,最后求出每个台站TECi 的相关系数.
同时,分析TEC在频域上的变化,使用Welch法计算功率谱密度,窗函数选择海明窗,计算数据采样频率1/900 Hz,考虑TEC具有明显的日变特征,因此,本文重点考察的周期成分是24 h、12 h,窗口 长度选择96×3个点,窗口重叠率50%,数据长度 94×12个点. 3.1.2 COSMIC NmF2数据处理
计算2月12日至6月11日11 ∶ 00—14 ∶ 00(地方时)NmF2数据,30天为间隔,使用基于球谐函数计算的模型值作为参考值,实际观测值作为考察值.首先,按2°×2°划分网格,如果网格内有多个值取均值;然后,计算模型值,再计算δ,δ=(xo-xm)/xm,xm表示模型值,xo表示观测值,估算模型误差的影响,以2008年2月12日—3月12日数据为例,这段时间太阳活动水平较低,计算数据按条件kp≤3进行筛选排除空间活动影响,计算得到平均误差为0.18,最大误差为0.52,最小误差接近于0,误差数据76.7%分布在-0.2~0.2区间内. 3.2 结果分析与讨论
按上述方法处理汶川地震震中周围YANC、XNIN、DLHA、LHAS、LUZH和XIAG台观测的TEC数据,分析TEC相关系数变化,同时给出kp和F107指数如图 3和图 4.
可以看出,各台站的日变与各自的平均TEC值具有很好的一致性,相关系数基本都在0.9以上,而在震前5月9日YANC、XNIN、DLHA、LUZH台开始有显著变化,相关系数降低到0.8以下,地震发生当天及震后5月15日YANC、XNIN、DLHA台相关系数减小到0.6以下,LHAS、XIAG台相关系数也有不同程度的变化,表明震前5月9日开始,各台站平常具有的良好相关性被打破,震前相关系数变化明显,检测到TEC信号在震前存在异常变化,并随地震的发生,震后相关系数异常变化相比震前更加显著,存在明显的同震响应现象.探讨相关系数产生变化的原因,对原始数据进行了分析核实,没有发现明显的干扰,并且数据已被多次应用(丁宗华等,2010;赵莹等,2010),因而数据本身是可靠的.同时,考察了空间活动的影响,地震前后全球磁情指数kp显示并无强烈磁场扰动,太阳活动水平也处于低水平,所以震中周围多个台站平时具有的一致性被打破不太可能是空间活动造成的.各台站TEC在地震前后的异常变化表明在震中周围发生了局部的电磁扰动现象,下面应用Welch法计算功率谱密度讨论TEC在频域上的变化.
时域分析的结果显示,5月9日开始各台站平时具有的良好相关性被打破,所以本文选择5月9日—20日数据计算功率谱密度,同时,使用地磁活动平静时期4月11日—22日数据计算功率谱密度进行对比分析,4月11日—22日这段时间每天ΣKP最大值为23,计算结果如图 5.
可以看出各台站24 h全日波成分功率谱密度幅度值最大,对比平静时期功率谱密度幅度值,YANC台TEC12 h半日波成分在地震发生时间功率谱密度幅度显著增大、24 h全日波成分幅度值降低,LUZH、XIAG和LHAS台TEC24 h全日波、12 h半日波成分在地震发生时间功率谱密度幅度明显降低,DLHA、XNIN台TEC24 h全日波在地震发生时间功率谱密度幅度值降低、12 h半日波幅度值变化不明显.其他周期成分在地震时间和平静时间功率谱密度幅度值变化趋于一致.
各台站TEC值在5月9日—5月20日频域上发生的变化,验证了地震前后TEC确实出现异常扰动.通过对原始数据的核实和排除空间环境影响,我们认为震前TEC在时域和频域上发生的扰动和汶川地震的孕育相关,各台站震前电离层TEC出现的扰动变化存在一定差异,可能受到地震孕育过程电磁信号的传播特性和地震电离层耦合作用的影响.同时,我们进一步讨论NmF2空间分布变化,计算2月12日—6月11日11 ∶ 00—14 ∶ 00(地方时)NmF2观测值相对模型值变化,30天间隔,计算δ,δ主要分布在-0.2~0.2区间,少数分布在-0.5~-0.2和0.2~0.5区间.这段范围可能由模型误差引起,不作为判断电离层扰动的依据.显示经度50—150°E和纬度10—70°N范围内δ分布变化如图 6.
计算结果表明,震前2月12日—5月11日,δ变化没有显著差异,4月12日—5月11日震中东南部区域电离层扰动幅度减小,δ分布在-0.5~-0.3区间,但是可能由模型误差引起,所以不能判定与汶川地震的孕育相关.5月12日—6月11日,在经度60—85°E、纬度26—54°N和经度110—130°E、 纬度范围46—62°N δ变化幅度明显增大,基本分布 在1.2~1.5区间,这段时间太阳活动处于低水平,数据按照地磁指数kp≤3筛选数据,电离层发生的异常扰动不太可能是空间活动造成的,2008年5月12日—6月11日在汶川及邻近区域发生MS≥6.0强余震5次,MS≥5.0余震32次,认为NmF2出现的强扰动现象可能与汶川地区发生的多次地震活动相关,存在电离层同震扰动现象.NmF2的空间分析表明,由于存在模型误差的原因,我们不能判定震前电离层扰动与汶川地震相关,但是,震后发生的电离层强扰动现象可能与地震活动存在一定联系.此外,在处理分析COSMIC星座获得的NmF2数据中还存在一些问题,比如,数据的精度还有待于进一步提高,同时由于数据精度不足,也降低了分析方法的可靠性. 4 讨论和结论
很多学者利用GPS TEC数据和其他电离层参量分析研究了汶川地震前后的电离层扰动现象,主要研究成果表明,汶川地震震前三天即5月9日震中附近的NmF2和TEC等电离层参量发生扰动,异常显著增强,同时其扰动特征在南北磁共轭区同时出现,并具有向磁赤道漂移的趋势(Zhao et al., 2008;Yu et al., 2009;丁宗华等,2010;赵莹等,2010);张学民(2009a,2009b)使用GPS TEC、电离 层测高仪得到的foF2、地电阻率观测和DEMETER 卫星观测的离子密度和温度等数据研究发现,在汶川地震前各种电磁观测技术记录到了明显的异常现象;林剑(2009)基于GPS TEC数据,研究发现在汶川地震前后一个星期孕震区上空电离层连续出现异常扰动;吴云(2011)利用卡尔曼法确认地震前第三天和第二天的显著TEC异常与地震相关;何宇飞(2009)基于DEMETER数据研究发现,震中上空区域不同频率信号对应的信噪比在震前明显降低;朱涛(2011)研究发现,地震发生时段14 kHz以下昼测和夜测电场的平均功率谱密度幅度减小,汶川地 震可能是导致2008年5月9日—5月24日时段10~14 kHz 昼测和10~18 kHz夜测电场平均功率谱密度幅度减小的原因;马新欣(2013)基于COSMIC掩星数据研究发现汶川地震震中邻近区域震前Ne显著变化.目前的研究没有涉及汶川地震前后TEC的相关分析和功率谱密度的分析,本文针对此,研究了GPS TEC的相关系数和功率谱密度变化,并分析了 地震前后NmF2空间分布变化,得出主要研究结论:
(1)从震前5月9日开始各台站相关性有不同程度的变化,YANC、XNIN、DLHA、LUZH台相关系数变化显著,地震发生当天和震后5月15日YANC、XNIN、DLHA台相关系数减小到0.6以下,TEC平常具有良好的相关性被打破,并且,震后相关系数异常变化相比震前更加显著.
(2)YANC、LUZH、XIAG和LHAS台TEC24 h全日波、12 h 半日波成分在地震发生时间功率谱密度幅度值变化明显,DLHA、XNIN台TEC24 h全日波成分在地震发生时间功率谱密度幅度值降低、12 h半日波成分幅度值变化不明显,其他周期成分在地震发生时间和平静时间功率谱密度幅度值变化趋于一致.
(3)2008年5月12日—6月11日在汶川地震邻近区域δ变化幅度明显增大,NmF2发生强扰动,我们认为可能与汶川地区发生的多次地震活动相关,震前电离层的扰动现象,由于存在模型误差的影响,我们不能判定与汶川地震孕育相关.
研究结果显示,TEC和NmF2时空分布变化存在电离层震前扰动和同震响应现象,并且TEC震后的异常变化相比震前更加显著.本文采用的方法已经被应用并发表(Pulinets et al., 2004;朱涛等,2011;马新欣等,2013),因而方法本身是可靠的.与以前的研究结论对比表明,使用的方法不仅可以探测到震前电离层异常信号,还可以提取到同震震后扰动信息,对进行地震短期预报具有参考价值.但是地震-电离层耦合关系是非常复杂的物理和化学过程,对耦合机制目前还没有统一的解释,研究认为,地震产生的电磁效应主要以化学、声学、电磁途径向空中传播,进而引起电离层扰动,主要机制为声重波、空气重力波与垂直电场机制(Freund,2004; Liperovsky et al., 2008; Du et al., 2011; Pulinets et al., 2011).现阶段,主要工作仍集中在震例分析和模型假说研究,鉴于岩石圈-大气层-电离层耦合的复杂性和电离层异常的多样性,若将震前电离层参量应用于实际地震预报中,在耦合机制和有效的提取电离层前兆异常信号方面,还需要做很多工作.
[1] | Akhoondzadeh M, Parrot M, Saradjian M R. 2010. Electron and ion density variations before strong earthquakes (M>6.0) using DEMETER and GPS data. Nat. Hazards Earth Syst. Sci. , 10(1): 7-18. |
[2] | Akhoondzadeh M. 2012. Anomalous TEC variations associated with the powerful Tohoku earthquake of 11 March 2011. Nat. Hazards Earth Syst. Sci. , 12(5): 1453-1462. |
[3] | Barbe K, Pintelon R, Schoukens J. 2010. Welch method revisited: nonparametric power spectrum estimation via circular overlap. Signal Processing, 58(2): 553-565. |
[4] | Du A M, Tsurutani B T, Sun W. 2011. Solar wind energy input during prolonged, intense Northward interplanetary magnetic fields: A new coupling function. Journal of Geophysical Research, 116: A12215, doi: 10. 1029/2011JA016718. |
[5] | Ding Z H, Wu J, Sun S J,et al. 2010. The variation of ionosphere on some days before the Wenchuan Earthquake. Chinese J. Geophys.(in Chinese), 53(1): 30-38. |
[6] | Freund F A, Takeuchi A, Lau B W S, et al. 2004. Stress-induced changes in the electrical conductivity of igneous rocks and the generation of ground currents. Terrestrial Atmospheric and Oceanic Sciences, 15(3): 437-468. |
[7] | Han P, Huang Q H, Xiu J G,et al. 2009. Principal component analysis of geomagnetic diurnal variation associated with earthquakes: case study of the M6.1 Iwate-ken Nairku Hokubu earthquake. Chinese J. Geophys.(in Chinese), 52(6): 1556-1563. |
[8] | Hen Y F, Yang D M,Chen H R,et al. 2009. SNR changes of VLF radio signals detected onboard the DEMETER satellite and their possible relationship to the Wenchuan earthquake. Science in China Series D: Earth Sciences,39(4): 403-412. |
[9] | Kakinami Y S, Liu J Y, Tsai L C, et al. 2010. Ionospheric electron content anomalies detected by a FORMOSAT-3/COSMIC empirical model before and after the Wenchuan Earthquake. International Journal of Remote Sensing, 31(13): 3571-3578. |
[10] | Le H, Liu J Y, Liu L. 2011. A statistical analysis of ionospheric anomalies before 736 M6.0+ earthquakes during 2002-2010. J. Geophys. Res., 116: A02303, doi: 10. 1029/2010JA015781. |
[11] | Lin J, Wu Y, Zhu F Y, et al.. 2009. Wenchuan earthquake ionosphere TEC anomaly detected by GPS. Chinese J. Geophys.(in Chinese), 52(1): 297-300. |
[12] | Liperovsky V A, Pokhotelov O A, Meister C V, et al. 2008. Physical models of coupling in the lithosphere-atmosphere-ionosphere system before earthquakes. Geomagnetism and Aeronomy, 48(6): 795-806. |
[13] | Liu J Y, Chen Y I, Chen C H, et al. 2009. Seismoionospheric GPS total electron content anomalies observed before the 12 May 2008 Mw7.9 Wenchuan earthquake. Journal of Geophysical Research, 114: A04320, doi: 10. 1029/2008JA013698. |
[14] | Liu J Y, Chen C H, Chen Y I, et al. 2010. A statistical study of ionopheric earthquake precursors monitored by using equatorial ionization anomaly of GPS TEC in Taiwan during 2001-2007. Journal of Asian Earth Sciences, 39(1-2): 76-80. |
[15] | Ma X X, Lin Z,Chen H R, et al. 2013. Ionosphere anomaly before the Wenchuan Ms8.0 earthquake detected by COSMIC occultation data. Acta Seismologica Sinica (in Chinese), 35(6): 848-855. |
[16] | Namgaladza A A, Zolotov O V, Karpov M, et al. 2012. Manifestations of the earthquake preparations in the ionosphere total electron content variations. Natural Science, 4(11): 848-855. |
[17] | Pulinets S, Ouzounov D. 2011. Lithosphere-Atmosphere-Ionosphere Coupling (LAIC) model-An unified concept for earthquake precursors validation. Journal of Asian Earth Sciences, 41(4): 371-382. |
[18] | Pulinets S A, Boyarchuk K A, Hegai V V, et al. 2000. Quasielectrostatic model of Atmosphere-Thermosphere-Ionosphere coupling. Adv. Space Res. , 26(8): 1209-1218. |
[19] | Pulinets S A, Gaivoronska T B, Leyva C A, et al. 2004. Correlation analysis technique revealing ionospheric precursors of earthquakes. Nat. Hazards Earth Syst. Sci. , 4(5-6): 697-702. |
[20] | Wu Y, Fu N B, Lin J, et al. 2011. Research on TEC anomalies before Ms8.0 Wenchuan earthquake by using Kalman filtering. Journal of Geodesy and Geodynamics (in Chinese), 31(2): 23-27. |
[21] | Xiong J, Zhou Y,Wu Y, et al. 2010. Ionospheric VTEC anomaly detected by ground-based GNSS before Ms7.1 Yushu earthquake. Journal of Geodesy and Geodynamics (in Chinese), 30(5): 24-27. |
[22] | Yao Y B, Chen P, Wu H, et al. 2012. Analysis of Ionospheric anomalies before the 2011 Ms9.0 Japan earthquake. Chin. Sci. Bull. , 57(5): 355-365. |
[23] | Yu T, Mao T,Wang Y G, et al. 2009. Study of the Ionospheric anomaly before the Wencuan earthquake. Chin. Sci. Bull., 54(4): 493-499. |
[24] | Yu T, Mao T, Wang Y G, et al. 2009. Study of the Ionospheric anomaly before the Wenchuan earthquake. Chinese Science Bulletin, 54(6): 1080-1086. |
[25] | Yu R, Wang W X,Liu L B, et al. 2009. A global Ionospheric TEC perturbation index. Chinese J. Geophys. (in Chinese), 52(9): 2189-2194. |
[26] | Zhang H B, Ou M,Zhen W M, et al. 2011. Ionospheric anomalies analysis of Eastern Japan great earthquake based on comprehensive observed data from the space and the ground. Chinese Journal of Radio Science (in Chinese), 26(3): 486-493. |
[27] | Zhang H P. 2006. Research on monitoring and delay correction of ionosphere using ground-based GPS data in Chinese region[Ph. D. thesis]. Shanghai: Shanghai Astronomical Observatory, Chinese Academy of Sciences. |
[28] | Zhang X M, Li M, Guan H P. 2009a. Anomaly analysis of earth resistivity observations before the Wenchuan earthquake. Earthquake (in Chinese), 29(1): 108-115. |
[29] | Zhang X M, Ding J H, Shen X H. 2009b. Electromagnetic perturbations before Wenchuan M8 earthquake and stereo electromagnetic observation system. Chinese Journal of Radio Science (in Chinese), 24(1): 1-8. |
[30] | Zhao Y, Zhang X H , Liu J N. 2010. Perturbation analysis of the Ionospheric TEC before and after the Wenchuan earthquake. Progress in Geophys. (in Chinese), 25(2): 447-453. |
[31] | Zhu T, Wang L W. 2011. LF electric field anomalies related to Wenchuan earthquake observed by DEMETER satellite. Chinese J. Geophys.(in Chinese), 54(3): 717-727. |
[32] | Zhao B, Wang M, Yu T, et al. 2008. Is an unusual large enhancement of Ionospheric electron density linked with the 2008 great Wenchuan earthquake? Journal of Geophysical Research, 113: A11304, doi: 10. 1029/2008JA013613. |
[33] | 丁宗华, 吴建, 孙树计等. 2010. 汶川大地震前电离层参量的变化特征与分析. 地球物理学报, 53(1): 30-38. |
[34] | 韩鹏, 黄清华, 修济刚. 2009. 地磁日变与地震活动关系的主成分分析——以日本岩手县北部6.1级地震为例. 地球物理学报, 52(6): 1556-1563. |
[35] | 何宇飞, 杨冬梅, 陈化然等. 2009. DEMETER卫星探测到可能与汶川地震有关的地面VLF发射站信号的信噪比变化. 中国科学(D辑), 39(4):403-412. |
[36] | 林剑, 吴云, 祝芙英等. 2009. 基于GPS探测汶川地震电离层TEC的异常. 地球物理学报, 52(1): 297-300. |
[37] | 马新欣, 林湛, 陈化然等. 2013. 基于COSMIC掩星数据探测汶川MS8.0地震震前电离层异常. 地震学报, 35(6): 848-855. |
[38] | 吴云, 付宁波, 林剑等. 2011. 用卡尔曼滤波法分析汶川MS8.0地震TEC异常. 大地测量与地球动力学, 31(2): 23-27. |
[39] | 熊晶, 周义炎, 吴云. 2010. 玉树MS7.1地震前电离层VTEC异常. 大地测量与地球动力学, 30(5): 24-27. |
[40] | 姚宜斌, 陈鹏, 吴寒等. 2012. 2011年3月11日日本地震震前电离层异常变化分析. 科学通报, 57(5): 355-365. |
[41] | 余涛, 毛田, 王云冈等. 2009. 汶川特大地震地震前电离层主要参量变化. 科学通报, 54(4): 493-499. |
[42] | 余优, 万卫星, 刘立波等. 2009. 全球电离层TEC起伏特性分析. 地球物理学报, 52(9): 2189-2194. |
[43] | 张红波, 欧明, 甄卫民等. 2011. 天地基综合观测的东日本大地震电离层异常分析. 电波科学学报, 26(3): 487-493. |
[44] | 章红平. 2006. 基于地基GPS的中国区域电离层监测与延迟改正研究[博士学位论文]. 上海: 中国科学院上海天文台, 33-34. |
[45] | 张学民, 李美, 关华平. 2009a. 汶川8.0级地震前的地电阻率异常分析. 地震, 29(1): 108-115. |
[46] | 张学民, 丁鉴海, 申旭辉等. 2009b. 汶川8级地震前电磁扰动与电磁立体监测体系. 电波科学学报, 24(1): 1-8. |
[47] | 赵莹, 张小红, 刘经南. 2010. 汶川大地震前后电离层电子含量的扰动分析. 地球物理学进展, 25(2): 447-453. |
[48] | 朱涛, 王兰炜. 2011. DEMETER卫星观测到的与汶川地震有关的LF电场异常. 地球物理学报, 54(3): 717-727. |