Electromagnetic phenomena before 2008 Wenchuan MS8.0 and 2010 Yushu MS7.1 earthquakes
-
摘要: 应用全局小波能谱法对2008年四川汶川MS8.0地震、 2010年青海玉树MS7.1地震前地电、 地磁场变化情况进行了分析研究. 结果表明: ① 在汶川地震震中周围的台站观测到震前地电、 地磁场全局小波能谱值增大的现象; ② 青藏高原东北缘的代乾等3个地电、 地磁台站在玉树地震前均发生显著的小波能谱增大的现象, 震后恢复; ③ 对于上述两次大震周边地区的不同台站的不同测道, 在震前多次出现谱值时间上同步增大或减少的现象; ④ 距玉树地震震中距离相同的山丹和古丰地电场台站, 北南、 北西测向长极距在地震前后能谱值的变化几乎一致. 基于上述分析, 初步研究认为震源孕育激发的电磁辐射是造成震前电磁异常现象的主要原因.Abstract: This paper studied the changes of the geoelectric and geomagnetic fields before the 2008 Wenchuan MS8.0 earthquake in Sichuan Province and the 2010 Yushu MS7.1 earthquake in Qinghai Province using the wavelet power spectrum analysis. The results showed that: ① The increased wavelet power value of the geoelectric and geomagnetic fields were observed by the geoelectric stations and the geoelectromagnetic stations near the epicenter before the Wenchuan earthquake. ② The wavelet power value of the geoelectric and geomagentic fields significantly increased in the TIA station and so on in the northeastern margin of Qinghai-Tibet Plateau before the Yushu earthquake, which recovered after the earthquake. ③ The phenomena that the wavelet power value synchronously increased/decreased were recorded more than once at the different channels of different stations near the above two earthquakes. ④ The changes of the wavelet power value were almost same observed by the NS and NW channels of Sandan station and Gufeng station, which are at the same distance from the epicenter of Yushu earthquake. Based on the above analysis, we drew the conclusion preliminarily that the geoelectromagnetic anomalies were mainly caused by the electromagnetic radiation stimulated by the focal development.
-
引言
410 km和660 km地幔间断面(以下简称410和660)是地幔转换带的上下界面,分别区隔了上、下地幔,其存在形态和性质研究对于理解地球内部物质运移具有重要意义.一般认为,410是橄榄石到瓦兹利石(α→β)的相变面,660则是林伍德石到钙钛矿和镁方铁矿(γ→pv+mw)的相变面(如,Ringwood,1991; Li et al,2000).410和660的相变深度对温度变化敏感(Collier et al,2001),例如,在冷俯冲板块区域,相变深度分别变浅和加深,地幔转换带的厚度随之增大(Bina,Helffrich,1994; Niu et al,2005);在洋中脊或热点区域,相变深度分别加深和变浅,地幔转换带的厚度随之减小(Deuss,2007).此外,地幔矿物中铝、铁成分以及水的介入会改变间断面附近的速度结构和密度分布(Weidner,Wang,1998; Deon et al,2011).
作为西北太平洋地区典型的洋陆俯冲区,日本海俯冲区是研究受俯冲影响的地幔物质组成及速度结构的天然实验室.地震层析成像结果(如,Huang,Zhao,2006; Fukao,Obayashi,2013)显示西太平洋板块在日本海沟—日本海西缘呈西向俯冲的特征,而在朝鲜半岛、中国东北等地区近水平停滞于地幔转换带底部而形成滞留板块(Huang,Zhao,2006; Fukao et al,2009).冷俯冲板块的存在影响了日本海俯冲区周边区域410和660的形态特征.接收函数研究结果显示:在日本列岛地区410出现窄的抬升(10—30 km),而660出现宽的下沉(20—50 km)(如,Tonegawa et al,2005; Niu et al,2005); 在中国东部和朝鲜半岛地区410和660的变化幅度稍小(Gao et al,2010);在中国东北地区660下沉幅度达35 km(Li,Yuan,2003),局部存在可能与石榴子石相变相关的多重间断面(Ai et al,2003).
三重震相方法是探测地球内部速度结构的有力工具(Tajima,Grand,1995).基于密集台网的三重震相波形拟合能有效地揭示间断面附近的精细结构(如,Wang et al,2006; Chu et al,2012).利用P波和S波三重震相对东北亚地区地幔转换带速度结构的研究,多数采用深源地震波形资料揭示该区域660的多样性,即660出现了不同程度的下沉(15—70 km)(如,Tajima,Grand,1995; Wang et al,2006; Zhang et al,2012)和梯度间断面(30—70 km)(如,Wang,Niu,2010; 叶玲玲,李娟,2012)等.因此,利用中源地震的三重震相资料可同时约束410和660的形态特征,更好地揭示地幔转换带的速度结构.
合适震级的地震可以提供相对简单的震源时间函数,密集且相对均匀分布的地震台能够提供高质量的地震波形资料.本文基于中国数字地震台网(China Digital Seismograph Network,简写为CDSN)记录的日本北海道地区2011年10月21日发生的一次中源地震(Mb6.0)的宽频带波形资料(图 1),利用P波和S波三重震相研究日本海俯冲区地幔转换带的速度结构,并以此探讨太平洋俯冲板块对区域结构的影响以及与板块俯冲相关的地球动力学过程.
图 1 本文所用地震位置和台站分布图震源球标出了2011年10月21日中源地震的震中位置,黑色和红色三角形分别表示P波和S波三重震相所用台站,黑色和红色圆点分别为P波和S波射线拐点的地表投影;Ⅰ和Ⅱ为以方位角15°为间隔所划分的两个子区域;灰色线为板块边界(引自DeMets et al,1994),黑色虚线为和达-贝尼奥夫带等深线(引自Gudmundsson,Sambridge,1998)Figure 1. Location of the earthquake and stations used in the study The beach ball represents the epicentral location of the intermediate-depth earthquake on 21 October 2011. Black and red triangles denote the stations for the P- and S-wave triplications,respectively. Black and red dots are the surface projections of P and S-wave turning points,respectively. I and II label two sub-regions divided by the azimuthal interval of 15°. Gray lines indicate the plate boundaries(DeMets et al,1994). Contours shown in black dashed lines represent the Wadati-Benioff Zone(Gudmundsson,Sambridge,1998)1. 数据和方法
1.1 三重震相方法
在地震波传播过程中,地震射线在地球内部速度间断面附近会出现三重震相现象.本文基于IASP91模型(Kennett,Engdahl,1991),利用反射率法(Wang,1999)计算了中源地震(震源深度为188.0 km)的射线路径和理论地震图,如图 2所示.可以看到,地震波在向下传播过程中会遇到410和660两个速度间断面,因此会出现双三重震相现象,相应的折合理论走时曲线呈“双回折”形态.以410附近的三重震相为例,AB震相的斜率对410之上的速度梯度变化敏感,CM则对410之下的速度梯度变化敏感;在AB与CM震相交点O所对应的震中距Δ=15.8°处,台站所接收的410上方和410下方的回折波到时相同(图 2b),交点O对应的震中距对410深度的变化较为敏感.三重震相的射线路径差别主要集中在间断面附近,其相对到时和波形振幅变化能有效地约束间断面附近的速度结构(叶玲玲,李娟,2012; 李国辉等,2014).三重震相通过地震射线拐点的采样可对地球深部速度结构进行约束,采样区域和拐点深度取决于震源深度和震中距的变化(Wang,Niu,2010; 眭怡,周元泽,2015).由于实际资料中三重震相各分支的到时拾取较为困难,因此对观测资料进行波形拟合可有效地揭示地幔转换带的精细结构(如,Wang,Chen,2009; Li et al,2013).
图 2 中源地震P波三重震相原理图(a)三重震相射线路径,包含了410上方的回折波(AB)、410之上的反射波(BC)、410下方的回折波(CD)、660之上的反射波(DE)和660下方的回折波(EF);(b)三重震相的折合走时及波形图,图中字母A-F, M,O和O′是对走时曲线的标记,Δ表示震中距Figure 2. Schematic illustration of the P-wavetriplication for an intermediate-depth earthquake (a)The ray paths for the triplication,including the turning wave above the 410 km discontinuity(hereafter as 410) (AB),the reflection on the 410(BC),the turning wave below the 410(CD),the reflection on the 660 km discontinuity (hereafter as 660)(DE)and the turning wave below the 660(EF);(b)The reduced travel times and synthetic waveforms for the triplication. The traveltime curves are marked with the symbols A-F,M,O and O′,and Δ represents epicentral distance本文使用延时法(Buland,Chapman,1983)计算理论走时曲线,并利用反射率法(Wang,1999)合成理论地震图;基于多次正演试错测试,得到最适P波和S波速度模型,并通过纵横波波速比推测影响地幔物质组成和造成速度异常的因素.
1.2 数据收集及处理
本文选取了中国数字地震台网记录的发生于日本北海道地区一次Mb6.0中源地震的宽频带波形资料,地震数据引自国家测震台网数据备份中心(郑秀芬等,2009).该地震的震源参数引自国际地震中心目录(International Seismological Centre,2012),具体列于表 1;合成波形所采用的震源机制解参数引自全球质心矩张量解(global centroid-moment-tensor,简写为GCMT)(Dziewonski et al,1981; Ekström et al,2012).
表 1 地震事件震源参数表Table 1. Source parameters of the earthquake event年-月-日 时:分:秒 震源深度*/km 北纬/° 东经/° Mb 2011-10-21 08:02:37.62 188.0 43.8729 142.5315 6.0 *经本文重定位的震源深度 震源深度的不确定性会影响三重震相的形态,因此需要对震源深度进行重新定位.本文从美国地震学研究联合会(Incorporated Research Institutions for Seismology,简写为IRIS)数据管理中心下载了震中距处于30°—95°内的14个台网/台阵的宽频带垂向分量波形数据,如全球地震台网(Global Seismographic Network)、美国国家地震台网(United States National Seismic Network)等,并从中挑选出69道信噪比较高的波形,如图 3a所示.考虑到震源机制解(图 3b)会影响不同方位区域的P波初动方向,在手动读取震相P和pP初至时将观测波形(图 3c)与合成理论波形进行了对比分析.基于IASP91模型(Kennett,Engdahl,1991)和CRUST2.0模型(Bassin et al,2000)计算震相pP和P的理论到时. 当扰动震源深度使pP-P的理论走时差与观测走时差的均方根达到最小值(图 3d)时,经重定位后的震源深度为188.0 km.
图 3 地震事件震源深度重定位(a)重定位所用台站分布;(b)震源球图解(Dziewonski et al,1981; Ekström et al,2012);(c)图(a)中台站所记录到的观测波形;(d)pP-P相对时差的残差均方根随震源扰动深度的变化Figure 3. The focal depth relocation of the earthquake event(a)Distribution of the stations used in the relocation;(b)The focal mechanism(Dziewonski et al,1981;Ekström et al,2012);(c)Observed waveforms of the stations in Fig.(a);(d)Root mean square of time residuals of pP-P relative intervals versus focal depth perturbation本文首先对中国数字地震台网记录的原始地震资料进行去均值、去线性趋势以及去仪器响应处理,针对P波和S波分别采用0.05—1.0 Hz和0.02—0.5 Hz带通滤波处理,最后将速度记录通过矩形积分转化为位移记录.考虑到区域速度结构可能存在横向变化,本文以15°方位角间隔将台站划分为I和II两个子区域(图 1),地震射线拐点的投影主要集中在日本海俯冲区且呈现与俯冲区近垂直分布的特征,射线拐点的深度范围约为208—791 km,可较好地约束地幔转换带的速度结构.
2. 结果
本文所用地震与410和660相关的双三重震相出现在10°—30°震中距范围内,I和II子区域的P波、S波波形拟合如图 4和图 5所示.从I和II子区域观测P波资料(图 4a和4c)可以看出:410附近的三重震相呈“窄AOC区”和“窄BOM区”的特征,即OA与OC震相、OB与OM震相之间的相对时差明显变小;AB震相的斜率在震中距处于15°—19.6°范围内时变大,且其终止震中距(尖端B,19.6°)大于IASP91模型的理论预测值(18.8°),这表明410之上存在一高速层;OC与OA震相相对时差明显变小,OC震相较预测震相提前,且两震相交点O对应的震中距(15.3°)小于预测震中距(15.8°),这表明410出现了一定程度的上升.由图 4a和4c还可看到:660附近的三重震相则呈现“宽MO′E区”和“窄DO′F区”的特征,即O′M与O′E震相之间的相对时差变大,O′D与O′F震相之间的相对时差变小;MD震相的斜率在震中距处于23.0°—28.4°范围内时变大,且MD震相的终止震中距(尖端D,28.4°)大于理论预测值(26.3°),这表明在660之上存在一高速层;O′M与O′E震相之间的相对时差变大,O′E震相的起始震中距(尖端E,18.4°)大于理论预测值(16.2°),且两震相交点O′对应的震中距(23.4°)大于理论预测值(21.9°),这表明660出现了一定程度的下沉.
图 4 Ⅰ和Ⅱ子区域P波波形拟合图(a)和(b)分别为Ⅰ子区域观测和合成P波波形图,图(c)和(d)分别为Ⅱ子区域观测和合成P波波形图,图(e)为P波速度模型IASP91(黑线)和JSP(红线).图(b)和(d)中合成波形为基于JSP模型合成,图(a)-(d)中黑色和红色折合走时曲线分别由IASP91和JSP模型计算所得Figure 4. The P waveform fitting for sub-regions Ⅰand ⅡFigs.(a)and(b)show observed and synthetic P waveforms for the sub-region Ⅰ,respectively,and Figs.(c)and(d)show those for the sub-region Ⅱ,respectively. Fig.(e)shows the P-wave velocity models IASP91(black line)and JSP(red line). Synthetic waveforms in(b)and(d)are calculated with the JSP model,and in(a)-(d),black and red reduced traveltime curves are calculated with the IASP91 and JSP models,respectively图 5 Ⅰ和Ⅱ子区域S波波形拟合图(a)和(b)分别为Ⅰ子区域观测和合成S波波形图,图(c)和(d)分别为Ⅱ子区域观测和合成S波波形图, 图(e)为S波速度模型IASP91(黑色线条)and JSS(红色线条). 图(b)和(d)中合成波形是由JSS模型合成,图(a)-(d)中黑色和红色折合走时曲线分别是由IASP91和JSS模型计算所得Figure 5. The S waveform fitting for sub-regions Ⅰand ⅡFigs.(a)and(b)show observed and synthetic S waveforms for the sub-region Ⅰ,respectively,and Figs.(c)and(d)show those for the sub-region Ⅱ,respectively. Fig.(e)shows the S-wave velocity models IASP91(black line)and JSS(red line). Synthetic waveforms in(b)and(d)are calculated with the JSS model,and in(a)-(d),black and red reduced traveltime curves are calculated with the IASP91 and JSS model,respectively通过波形拟合(图 4b和4d)获得的区域P波最适速度模型JSP(Japan Sea P-wave model)(图 4e)显示:410之上存在厚度为85 km的高速层,高速异常最大值为1.1%;410出现了10 km的小幅抬升,其附近的速度跳变为3.3%,小于IASP91模型的3.7%;660之上存在厚度为125 km的高速层,高速异常最大值为1.5%;660出现了25 km的下沉,其下方存在厚度为65 km的低速异常,速度减低约1.0%;660附近的速度跳变为3.7%,小于IASP91模型的5.8%.
需要指出的是,三重震相的一些变化特征可能与多个因素相关,因此,界定各因素的作用程度时需考虑与之相关的三重震相其它变化特征和多次的正演测试分析.例如:410上方存在的高速层和410的抬升均会影响OA与OC震相的相对时差(图 4a);相较该变化特征,AB震相的斜率及OB与OM震相的相对时差对410上方存在高速层这一因素更为敏感,AB与CM震相交点O所对应的震中距对410抬升这一因素也更为敏感.
通过拟合观测S波三重震相(图 5a-d)得到了区域S波最适速度模型JSS(Japan Sea S-wave model)(图 5e),与JSP模型相比,二者的一致性整体上较强.JSS模型显示:410之上存在厚度为185 km的高速层,速度升高1.4%;410出现了10 km的小幅抬升,其附近的速度跳变量为3.1%,小于IASP91模型的4.1%;660之上存在厚度为125 km的高速层,高速异常最大值为1.0%;660出现了25 km的下沉,其下方存在厚度为65 km的低速异常,速度减低约1.0%;660附近的速度跳变为5.4%,小于IASP91模型的6.25%.
通过区域最适P波和S波速度模型可得到纵横波波速比vP/vS,其结果显示:该值在210—400 km深度范围内偏低,约为1.827,小于IASP91模型对应深度范围内的平均值1.844,平均降低约0.92%;而在560—685 km深度范围内数值偏高,约为1.831,大于IASP91模型对应深度范围内的平均值1.823,平均升高约0.44%.
3. 讨论与结论
3.1 410附近的速度结构
本文探测到日本海俯冲区410上方P波和S波均存在高速层,将410之上的P波和S波地震射线反射点分别投影至全球层析成像P波速度模型GAP_P4(Obayashi et al,2013)和S波速度模型S40RTS(Ritsema et al,2011)不同深度的水平图像,并对其分布进行对比分析.图 6a和6b分别给出了410之上的P波和S波射线的反射点分布,可以看到,P波和S波射线的反射点均位于地震波高速异常区,应为西太平洋俯冲板块.Wang(2014)关于vP和vP/vS联合反演结果显示,在日本东北地区,太平洋俯冲板块为一高vP、高vS和低vP/vS的异常区;Xia等(2008)的独立走时层析成像结果也显示出太平洋俯冲板块低泊松比的特征.本文探测到在210—400 km深度范围内低vP/vS(约1.827)的特征,体现了海洋板块低泊松比的特征,与层析成像反演结果相一致.
图 6 410和660之上的P波和S波反射点分布图图(a)和(c)中黑色圆点分别为410和660之上的P波反射点,背景图像分别为GAP_P4模型(Obayashi et al,2013)在400 km和685 km深度的水平图像;图(b)和(d)中红色圆点分别为410和660之上的S波反射点,背景图像分别为S40RTS模型(Ritsema et al,2011)在400 km和685 km深度的水平图像Figure 6. The location of P- and S-wave reflected points on the 410 and 660In Figs.(a)and(c),black dots are the P-wave reflected points on the 410 and 660,respectively; background images are the horizontal images from the GAP_P4 model(Obayashi et al,2013)at the depth of 400 km and 685 km,respectively. In Figs.(b)and(d),red dots are the S-wave reflected points on the 410 and 660,respectively; background images are the horizontal images from the S40RTS model(Ritsema et al,2011)at the depth of 400 km and 685 km,respectively410附近的克拉伯龙(Clapeyron)斜率为2.7—2.9 MPa/K(Bina,Helffrich,1994),10 km的抬升则大概对应于120—129 K的理论低温异常.Kawakatsu和Yoshioka(2011)由数值模拟实验得到的日本海俯冲区温度剖面显示,在400 km深度处俯冲板块内部最冷处约为600 K,板块周围的地幔温度约为1200 K,两者温差约为600 K.本文所换算的地幔低温异常小于数值模拟的实验结果,偏低的低温异常值反映了冷俯冲板块所导致的平均地幔低温效应.
部分研究显示在西太平洋部分地区410上方存在低速层或低速异常体,例如黄海和中国东部地区(Revenaugh,Sipkin,1994)、扬子克拉通(李国辉等,2014)、东海地区(眭怡,周元泽,2015)等,这可能与水及其它挥发分物质的介入引起橄榄石出现部分熔融有关;南千岛地区(Wang,Chen,2009)和日本本州地区(Obayashi et al,2006)的结果则显示可能与俯冲板块下方的下地幔上升热源有关.本文中410之上的地震射线拐点集中于太平洋俯冲板块上(图 5a和5b)且410出现了10 km的抬升,这表明近冷俯冲板块区域的地幔低温效应较为明显,而远离俯冲板块的区域由于受挥发分物质介入或者热物质上涌的影响,410之上可能会出现低速层或低速异常体.
3.2 660附近的速度结构
本文研究结果显示:在560—685 km深度范围内存在高速层,P波和S波高速异常的最大值分别为1.5%和1.0%,660附近的地震射线反射点集中在日本海西缘、朝鲜半岛北部以及中国东北等地区(图 6c和6d);太平洋俯冲板块在此拐点集中区域由俯冲下降的形态过渡至近水平停滞于地幔转换带底部(Zhao et al,2007; Fukao et al,2009).本文探测到660之上的高速层厚度为125 km,明显大于俯冲大洋板块的厚度(76—77 km)(周春银等,2010),这可能与俯冲或滞留板块在地幔转换带内发生堆积和内部变形的过程有关(Fukao et al,2009; Tang et al,2014).
660附近的克拉伯龙斜率为-2.8 MPa/K(Hirose,2002),25 km下沉所对应的理论低温异常约为313 K,大于410附近的低温异常(120—129 K),这可能是由于滞留板块在地幔转换带内与周围地幔的接触面积较大所致.
本文所探测到的纵横波波速比vP/vS在560—685 km深度范围内显示为高值,约1.831.由于S波波速较P波波速对水的介入更为敏感,水的存在可降低S波波速且使vP/vS值增大(Smyth et al,2004; Wang,Chen,2009).Jacobsen等(2004)的高温高压实验结果显示,含水1.0wt%的Fo11林伍德石的vP/vS值为1.740±0.040,而不含水的Fo10林伍德石的vP/vS值为1.706±0.002.假定vP/vS值与林伍德石中的含水量之间呈线性关系,则可推算在560—685 km深度附近地幔转换带约含水0.24wt%.
在660—750 km深度附近存在地震波低速异常,P波和S波波速均降低约1.0%.Zhao等(2007)的层析成像结果显示,在太平洋板块下方存在低速异常,并认为这是由于受太平洋板块俯冲以及部分板块碎片崩塌进入下地幔的影响,在板块下方出现下地幔热物质上涌所致.利用700 km深度附近的地幔热膨胀系数lnvP/T=-0.28×10-4/K(Cammarano et al,2003),低速异常区附近的高温异常约为355 K.此外,俯冲板块可在下地幔顶部通过超水相B分解而释放一定量的水(Ohtani,2005),俯冲或滞留板块下方所释放的水可压入下地幔顶部(Fukao et al,2009).因此,本文推测660下方的低速异常可能与下地幔顶部的热物质上涌有关,俯冲板块深脱水所产生的少量水可能对此低速异常也有所贡献.关于下地幔顶部低速异常的来源,尚需进一步研究和探讨.
中国地震局地球物理研究所国家数字测震台网数据备份中心(doi: 10.7914/SN/CB)和美国地震学研究联合会数据管理中心为本研究提供了地震波形数据,审稿人提出了宝贵的修改意见和建议,作者在此一并表示感谢.
-
图 3 2008年汶川MS8.0地震前后SHN台站地磁场、 地电场全局小波能谱图
(a) 地磁场X分量; (b) 地磁场Y分量; (c) 地电场E1E2测道; (d) 地电场N1S1测道
Figure 3. Wavelet power spectrum of geomagnetic field and geoelectric field at the station SHN before and after the Wenchuan MS8.0 earthquake
(a) X component of geomagnetic field; (b) Y component of geomagnetic field; (c) Measuring channel E1E2 of geoelectric field; (d) Measuring channel N1S1 of electric field
图 5 2008年汶川MS8.0地震前后山丹台地电场全局小波能谱图
(a) NS测道(长极距); (b) EW测道(长极距); (c) NW测道(长极距); (d) NS测道(短极距); (e) EW测道(短极距); (f) NW测道(短极距)
Figure 5. Wavelet power spectrum of geoelectric field at the station Shandan before and after Wenchuan MS8.0 earthquake
(a) NS measuring channel (long distance); (b) EW measuring channel (long distance); (c) NW measuring channel (long distance); (d) NS measuring channel (short distance);(e) EW measuring channel (short distance); (f) NW measuring channel (short distance)
图 6 2008年汶川MS8.0地震前后武都台地电场全局小波能谱图
(a) NS测道(长极距); (b) EW测道(长极距); (c) NW测道(长极距); (d) NS测道(短极距); (e) EW测道(短极距); (f) NW测道(短极距)
Figure 6. Wavelet power spectrum of geoelectric field at the station Wudu before and after Wenchuan MS8.0 earthquake
(a) NS measuring channel (long distance); (b) EW measuring channel (long distance);(c) NW measuring channel (long distance); (d) NS measuring channel (short distance); (e) EW measuring channel (short distance); (f) NW measuring channel (short distance)
图 7 2010年玉树MS7.1地震前后TIA台地电场全局小波能谱图
(a) N1S1测道; (b) N1E1测道; (c) S2S3测道; (d) S2E2测道; (e) N1S3测道; (f) WE2测道
Figure 7. Wavelet power spectrum of geoelectric field at the station TIA before and after Yushu MS7.1 earthquake
(a) N1S1 measuring channel; (b) N1E1 measuring channel; (c) S2S3 measuring channel; (d) S2E2 measuring channel; (e) N1S3 measuring channel; (f) WE2 measuring channel
图 8 玉树MS7.1地震前后山丹台地电场全局小波能谱图
(a) NS测道(长极距);(b) EW测道(长极距); (c) NW测道(长极距); (d) NS测道(短极距); (e) EW测道(短极距); (f) NW测道(短极距)
Figure 8. Wavelet power spectrum of geoelectric field at the station Shandan before and after Yushu MS7.1 earthquake
(a) NS measuring channel (long distance); (b) EW measuring channel (long distance); (c) NW measuring channel (long distance); (d) NS measuring channel (short distance); (e) EW measuring channel (short distance); (f) NW measuring channel (short distance)
图 9 玉树MS7.1地震前后古丰台地电场全局小波能谱图
(a) NS测道(长极距); (b) EW测道(长极距); (c) NW测道(长极距); (d) NS测道(短极距); (e) EW测道(短极距); (f) NW测道(短极距)
Figure 9. Wavelet power spectrum of geoelectric field at the station Gufeng before and after Yushu earthquake
(a) NS measuring channel (long distance); (b) EW measuring channel (long distance); (c) NW measuring channel (long distance); (d) NS measuring channel (short distance); (e) EW measuring channel (short distance); (f) NW measuring channel (short distance)
图 10 2010年4月7—15日山丹台(a)、 古丰台(b)地电场变化曲线和山丹台(c)、 古丰台(d)小波变换滤波后曲线
Figure 10. The curves of geoelectric field changes for the stations Shandan (a) and Gufeng (b) as well as the filtered curves of geoelectric field for the stations Shandan (c) and Gufeng (d) using the wavelet transform from April 7 to 15,2010
表 1 2008年3—6月四川汶川周围发生的MS≥4.0地震
Table 1 MS≥4.0 earthquakes around Wenchuan from March to June of 2008
-
安张辉, 杜学彬, 范莹莹, 刘君, 谭大成, 陈军营, 解滔. 2011. 汶川MS8.0级大震前天基与陆基电场资料联合应用研究[J]. 地球物理学报, 54 (11): 2876-2884. An Z H, Du X B, Fan Y Y, Liu J, Tan D C, Chen J Y, Xie T. 2011. A study of the electric field before the Wenchuan MS8.0 earthquake of 2008 using both space-based and ground-based observational data[J]. Chinese Journal of Geophysics, 54 (11): 2876-2884 (in Chinese).
杜学彬, 刘耀炜, 倪明康. 1997. 强震地下水动态中短临异常空间特征研究[J]. 地震学报, 19 (4): 426-433. Du X B, Liu Y W, Ni M K. 1997. On the spatial characteristic of the short-term and imminent anomalies underground water behaviors before strong earthquake[J]. Acta Seismologica Sinica, 19 (4): 426-433 (in Chinese).
杜学彬. 2010. 在地震预报中的两类视电阻率变化[J]. 中国科学: D辑, 40 (10): 1321-1330. Du X B. 2010. Two types of changes in apparent resistivity in earthquake prediction[J]. Science in China: Serises D, 40 (10): 1321-1330 (in Chinese).
范莹莹, 杜学彬, Jacques Zlotnicki, 谭大成, 刘君, 安张辉, 陈军营, 郑国磊, 解滔. 2010. 汶川MS8.0大震前的电磁现象[J]. 地球物理学报, 53 (12): 2887-2898. Fan Y Y, Du X B, Zlotnicki J, Tan D C, Liu J, An Z H, Chen J Y, Zheng G L, Xie T. 2010. The electromagnetic phenomena before the MS8.0 Wenchuan earthquake[J]. Chinese Journal of Geophysics, 53 (12): 2887-2898 (in Chinese).
高曙德, 汤吉, 杜学彬, 刘小凤, 苏永刚, 陈彦平, 狄国荣, 梅东林, 詹艳, 王立凤. 2010. 汶川8.0地震前后电磁场的变化特征[J]. 地球物理学报, 53 (3): 512-525. Gao S D, Tang J, Du X B, Liu X F, Su Y G, Chen Y P, Di G R, Mei D L, Zhan Y, Wang L F. 2010. The change characteristics of electromagnetic field before to after Wenchuan MS8.0 earthquake[J]. Chinese Journal of Geophysics, 53 (3): 512-525 (in Chinese).
郭自强, 周大庄, 施行觉, 马福胜, 席道瑛, 程纯杰, 周志文. 1988. 岩石破裂中的电子发射[J]. 地球物理学报, 31 (5): 566-571 Guo Z Q, Zhou D Z, Shi X J, Ma F S, Xi D Y, Cheng C J, Zhou Z W. 1988. Electron emission during rock fracture[J]. Chinese Journal of Geophysics, 31 (5): 566-571 (in Chinese).
李宁, 杜学彬, 谭大诚. 2007. 松山观测台地震短临电磁现象[J]. 地震, 27 (增刊): 103-111. Li N, Du X B, Tan D C. 2007. Imminent electro-magnetic phenomenon related to earthquakes recorded at the Songshan station[J]. Earthquake, 27 (Suppl): 103-111 (in Chinese).
李志海, 马宏生, 曲延军. 2009. 2008年3月21日新疆于田7.3级地震发震构造与震前地震活动特征研究[J]. 中国地震, 25 (2): 199-205. Li Z H, Ma H S, Qu Y J. 2009. Study on seismogenic structure and seismic activity characteristics before the Yutian MS7.3 earthquake on March 21, 2008, Xinjiang[J]. Earthquake Research in China, 25 (2): 199-205 (in Chinese).
刘君, 杜学彬, Jacques Zlotnicki, 范莹莹, 安张辉, 解滔, 郑国磊, 谭大成, 陈军营. 2011. 几次大震前的地面和空间电磁场变化[J]. 地球物理学报, 54 (11): 2885-2897. Liu J, Du X B, Zlotnicki J, Fan Y Y, An Z H, Xie T, Zheng G L, Tan D C, Chen J Y. 2011. Changes of the electric and magnetic fields on the ground and in the ionosphere before and after several great earthquakes[J]. Chinese Journal of Geophysics, 54 (11): 2885-2897 (in Chinese).
卢华夏, 贾东, 王良书, 王胜利. 2008. 关于汶川地震发震机制[J]. 高校地质学报, 14 (2): 133-138. Lu H X, Jia D, Wang L S, Wang S L. 2008. On the triggering mechanics of Wenchuan earthquake[J]. Geological Journal of China Universities, 14 (2): 133-138 (in Chinese).
吕品姬, 赵斌, 陈志遥, 张燕, 李正媛. 2011. 小波分解-STFT方法在地形变观测数据中的应用[J]. 大地测量与地球动力学, 31 (5): 136-140. Lü P J, Zhao B, Chen Z Y, Zhang Y, Li Z Y. 2011. Application of wavelet-decomposition and STFT method in continuous deformation observation analysis[J]. Journal of Geodesy and Geodynamics, 31 (5): 136-140 (in Chinese).
马瑾, 马胜利, 刘力强. 1995. 地震前异常的阶段性及其空间分布特征[J]. 地震地质, 17 (4): 363-371. Ma J, Ma S L, Liu L Q. 1995. The stages of anomalies before an earthquake and the characteristics of their spatial distribution[J]. Seismology and Geology, 17 (4): 363-371 (in Chinese).
马宗晋. 1980. 华北地壳的多(应力集中)点场与地震[J]. 地震地质, 2 (1): 39-47. Ma Z J. 1980. Multipoints (concentrated) stress field and earthquake in North China crust[J]. Seismology and Geology, 2 (1): 39-47 (in Chinese).
钱复业, 赵璧如, 钱卫, 赵健, 何世根, 张洪魁, 李世愚, 李绍坤, 严谷良, 汪成民, 孙振凯, 张东宁, 卢军, 张平, 杨国军, 孙加林, 郭纯生, 唐宇雄, 徐建明, 夏坤涛, 鞠航, 殷邦红, 黎明, 杨东生, 起卫罗, 和泰名, 关华平, 赵玉林. 2009. 汶川8.0级地震HRT波地震短临波动前兆及HRT波地震短临预测方法[J]. 中国科学: D辑, 39 (1): 11-23. Qian F Y, Zhao B R, Qian W, Zhao J, He S G, Zhang H K, Li S Y, Li S K, Yan G L, Wang C M, Sun Z K, Zhang D N, Lu J, Zhang P, Yang G J, Sun J L, Guo C S, Tang Y X, Xu J M, Xia K T, Ju H, Yin B H, Li M, Yang D S, Qi W L, He T M, Guan H P, Zhao Y L. 2009. Impending HRT wave precursors to the Wenchuan MS8.0 earthquake and methods of earthquake impending prediction by using HRT wave[J]. Science in China : SerisesD, 39 (1): 11-23 (in Chinese).
钱书清, 任克新, 吕智. 1996. 伴随岩石破裂的VLF、 MF、 HF和VHF电磁辐射特性的研究[J]. 地震学报, 18 (3): 346-351. Qian S Q, Ren K X, Lü Z. 1996. Experimental study on the characteristic of VLF, MF, HF and VHF seismo-electromagnetic radiation with rock burst[J]. Acta Seismologica Sinica, 18 (3): 346-351 (in Chinese).
钱书清, 吕智, 任克新. 1998. 地震电磁辐射前兆不同步现象物理机制的实验研究[J]. 地震学报, 20 (5): 535-540. Qian S Q, Lü Z, Ren K X. 1998. Experimental study on the mechanism of non-synchronism of seismo-electromagnetic radiation precursors[J]. Acta Seismologica Sinica, 20 (5): 535-540 (in Chinese).
汤吉, 詹艳, 王立凤, 董泽义, 赵国泽, 徐建朗. 2010. 汶川地震强余震的电磁同震效应[J]. 地球物理学报, 53 (3): 526-534. Tang J, Zhan Y, Wang L F, Dong Z Y, Zhao G Z, Xu J L. 2010. Electromagnetic coseismic effect associated with aftershock of Wenchuan MS8.0 earthquake[J]. Chinese Journal of Geophysics, 53 (3): 526-534 (in Chinese).
解滔, 杜学彬, 刘君, 范莹莹, 安张辉, 陈军营, 谭大成. 2013. 汶川MS8.0、 海地MS7.0地震电磁信号小波能谱分析[J]. 地震学报, 35 (1): 61-71. Xie T, Du X B, Liu J, Fan Y Y, An Z H, Chen J Y, Tan D C. 2013. Wavelet power spectrum analysis of the electromagnetic signals of Wenchuan MS8.0 and Haiti MS7.0 earthquakes[J]. Acta Seismologica Sinica, 35 (1): 61-71 (in Chinese).
曾中超, 张蓓, 方广有, 王东峰, 阴和俊. 2009. 利用DEMETER卫星数据分析汶川地震前的电离层异常[J]. 地球物理学报, 52 (1): 11-19. Zeng Z C, Zhang B, Fang G Y, Wang D F, Yin H J. 2009. The analysis of ionospheric variations before Wenchuan earthquake with DEMETER data[J]. Chinese Journal of Geophysics, 52 (1): 11-19 (in Chinese).
张元生, 郭晓, 钟美娇, 沈文荣, 李稳, 何斌. 2010. 汶川地震卫星热红外亮温变化[J]. 科学通报, 55 (10): 904-910. Zhang Y S, Guo X, Zhong M J, Shen W R, Li W, He B. 2010. The change of IR-brightness temperature about Wenchuan earthquake[J]. Chinese Science Bulletin, 55 (10): 904-910 (in Chinese).
赵和云, 阮爱国, 扬荣, 梁子斌, 韩德胜. 2002. 天祝大地电场异常与1995年永登5.8、 1996年天祝5.4级地震[J]. 西北地震学报, 24 (1): 56-64. Zhao H Y, Ruan A G, Yang R, Liang Z B, Han D S. 2002. The geoelectric field anormaly of Tianzhu and relationship with Yongdeng earthquake (MS5.8) in 1995 and Tianzhu earthquake (MS5.4) in 1996[J]. Northwestern Seismological Journal, 24 (1): 56-64 (in Chinese).
郑兆苾, 张军, 何康. 2007. 地电场数据的功率谱和小波分析[J]. 地震, 4 (1): 160-169. Zheng Z B, Zhang J, He K. 2007. Power spectrum and wavelet analysis of geoelectric data[J]. Earthquake, 4 (1): 160- 169 (in Chinese).
Краев А П(著). 1951. 张可迁, 陈培光, 张志诚, 顾燕庭, 汪盛辉(译). 1954. 地电原理[M]. 北京: 地质出版社: 100- 236. Краев А П. 1951. Geoelectrics Principle[M].Beijing: Geological Publishing House: 100-236 (in Chinese).
Adler P, Mouël L, Zlotnicki J. 1999. Electrokinetic and magnetic fields generated by flow through a fractured zone: A sensitivity study for La Fournaise volcano[J]. Geophys Res Lett, 26 (6): 795-798.
Du A M, Huang Q H, Yang S F. 2002. Epicenter location by abnormal ULF electromagnetic emissions[J]. Geophys Res Lett, 29 (10): 1455. doi:10.1029/2001GL013616.
Enomoto Y, Hashimoto H. 1992. Transient electrical activity accompanying rock under indentation loading[J]. Tectonophysics, 211 (1/4): 337-344.
Gonzalez W D, Joselyn J A, Kamide Y. 1994. What is a geomagnetic storm?[J]. J Geophys Res, 99 (A4): 5771-5792.
Hattori K, Akinaga Y, Hayakawa M. 2002. ULF magnetic anomaly preceding the 1997 Kagoshima earthquakes[C]//Hayakawa M, Molchanov O A eds. Seismo-Electromagnetics: Lithosphere-Atmosphere-Ionosphere Coupling. Tokyo: Terr Sc Pub: 19-28.
Huang Q H, Ikeya M. 1998. Seismic electromagnetic signals (SEMS) explained by a simulation experiment using electromagnetic waves[J]. Phys Earth Planet Inter, 109 (3/4): 107-114.
Huang Q H. 2005. Controlled analogue experiments on propagation of seismic electromagnetic signals[J]. Chinese Science Bulletin, 50 (17): 1956-1961.
Ishido T, Pritchett J W. 1999. Numerical simulation of electrokinetic potentials associated with subsurface fluid flow[J]. J Geophys Res, 104 (B7): 15247-15259.
Michel S, Zlotnicki J. 1998. Self-potential and magnetic surveying of La Fournaise volcano (Réunion Island): Correlations with faulting, fluid circulation, and eruption[J]. J Geophys Res, 103 (B8): 17845-17857.
Varotsos P, Lazaridou M. 1991. Latest aspects of earthquake prediction in Greece based on seismic electric signals[J]. Tectonophysics, 188 (3/4): 321-347.
Yamada I, Masuda K, Mizutani H. 1989. Electromagnetic and acoustic emission associated with rock fracture[J]. Phys Earth Planet Inter, 57 (1/2): 157-168.
-
期刊类型引用(3)
1. 苏慧,崔清辉,李国辉,尹迪,周元泽. 基于波形拟合的日本南部及琉球岛弧下方滞留板片的形态结构研究. 地震学报. 2025(01): 37-53 . 本站查看
2. 崔冉,周元泽. 地震波三重震相在地幔转换区结构研究中的应用. 地震地磁观测与研究. 2024(05): 11-19 . 百度学术
3. 苏慧,魏荣强,周元泽,崔清辉,李国辉. 东北亚边缘地区地幔过渡带内滞留太平洋板片上界面的三重震相研究. 地球物理学报. 2023(06): 2431-2444 . 百度学术
其他类型引用(4)