歐陽(yáng)新艷 申旭輝

1） 中國(guó)北京100036中國(guó)地震局地震預(yù)測(cè)研究所 2） 中國(guó)北京100871北京大學(xué)地球與空間科學(xué)學(xué)院3） 中國(guó)北京100085中國(guó)地震局地殼應(yīng)力研究所

?

DEMETER衛(wèi)星觀測(cè)的超低頻電場(chǎng)擾動(dòng)預(yù)處理方法及震例應(yīng)用分析

1） 中國(guó)北京100036中國(guó)地震局地震預(yù)測(cè)研究所 2） 中國(guó)北京100871北京大學(xué)地球與空間科學(xué)學(xué)院3） 中國(guó)北京100085中國(guó)地震局地殼應(yīng)力研究所

基于法國(guó)DEMETER衛(wèi)星觀測(cè)的超低頻（ULF）電場(chǎng)波形， 本文提出了更簡(jiǎn)單易實(shí)施的ULF電場(chǎng)擾動(dòng)預(yù)處理方法． 利用該方法， 對(duì)文獻(xiàn)中已報(bào)道的2010年2月27日智利MS8.8地震前的ULF電場(chǎng)擾動(dòng)進(jìn)行了更全面的分析． 結(jié)果表明： 已報(bào)道的20°S—40°S地磁緯度帶內(nèi)的ULF電場(chǎng)擾動(dòng)并非僅出現(xiàn)在最靠近震中的半軌道， 而是出現(xiàn)在較大的經(jīng)度范圍內(nèi)． 大范圍ULF電場(chǎng)擾動(dòng)的區(qū)域正好位于南大西洋異常區(qū)， 很可能受到該區(qū)域異常磁場(chǎng)的影響． 南大西洋異常中心區(qū)域的ULF電場(chǎng)擾動(dòng)表現(xiàn)出共軛特征， 即在南北地磁緯度20°—40°范圍內(nèi)ULF電場(chǎng)均出現(xiàn)約5 mV/m的擾動(dòng)， 北緯地區(qū)的ULF電場(chǎng)擾動(dòng)幅度相比南緯地區(qū)更小． 而在南大西洋異常中心區(qū)域外的ULF電場(chǎng)擾動(dòng)則僅在20°S—40°S地磁緯度帶內(nèi)出現(xiàn)， 在其共軛區(qū)并未觀測(cè)到ULF電場(chǎng)擾動(dòng)． 這可能是由于南大西洋異常中心區(qū)域外的ULF電場(chǎng)擾動(dòng)幅度比中心區(qū)域更小所致， 在電離層中由于碰撞效應(yīng)更易被衰減， 因此ULF電場(chǎng)擾動(dòng)無法傳播至其共軛區(qū)． 本文的方法和結(jié)果可為后續(xù)更全面地分析衛(wèi)星在電離層高度記錄的ULF電場(chǎng)擾動(dòng)提供參考．

DEMETER衛(wèi)星 超低頻（ULF）電場(chǎng)擾動(dòng) 南大西洋異常區(qū)

引言

早在20世紀(jì)80年代， 研究人員便從常規(guī)的地面觀測(cè)中獲得了與地震相關(guān)的電磁輻射現(xiàn)象（Gokhbergetal， 1982）． 由于地震產(chǎn)生的電磁擾動(dòng)具有短臨特征， 使得這一現(xiàn)象受到極大的關(guān)注（Parrot， 1995）． 其中， 超低頻（ultra low frequency, 簡(jiǎn)寫為ULF）（0.01—10 Hz）電磁擾動(dòng)被認(rèn)為是最具前景的地震前兆之一（Hayakawaetal， 2007）． 地面ULF磁場(chǎng)觀測(cè)已記錄到多次地震前的ULF擾動(dòng)（Fraser-Smithetal， 1990; Molchanovetal， 1992; Hayakawaetal， 1996） ．

近些年來， 隨著空間技術(shù)的發(fā)展， 研究人員利用衛(wèi)星資料監(jiān)測(cè)到了越來越多與地震相關(guān)的電磁擾動(dòng)現(xiàn)象（Chmyrevetal， 1989; Molchanovetal， 1993; Parrot， 1994; Němecetal， 2008）． 尤其是DEMETER衛(wèi)星發(fā)射以來， 利用該衛(wèi)星觀測(cè)資料開展地震電磁、 電離層擾動(dòng)的研究報(bào)道很多（歐陽(yáng)新艷等， 2008， 2011； 劉靜等， 2011； 顏蕊等， 2013）． 大多數(shù)研究集中討論了DEMETER衛(wèi)星探測(cè)的極低頻/甚低頻（ELF/VLF）頻段的擾動(dòng)（Bhattacharyaetal， 2007; 何宇飛等， 2009; 曾中超等， 2009; 張學(xué)民等， 2009a， b; Akhoondzadehetal， 2010; Onishietal， 2011; Zhangetal， 2011， 2012a; 朱濤， 王蘭煒， 2011; 澤仁志瑪?shù)龋?2012）． 在利用DEMETER衛(wèi)星觀測(cè)的ULF電場(chǎng)數(shù)據(jù)開展該頻段的擾動(dòng)研究方面， 近年來的研究報(bào)道逐漸增多． 2010年1月12日M7.0海地地震前一個(gè)月， Athanasiou等（2011）發(fā)現(xiàn)DEMETER衛(wèi)星探測(cè)的ULF頻段Ez分量的能量增加了約360%． 汶川地震前， 利用DEMETER衛(wèi)星觀測(cè)的ULF電場(chǎng)數(shù)據(jù)， Zhang等（2012b）發(fā)現(xiàn)ULF電場(chǎng)擾動(dòng)幅度約為3—5 mV/m， 頻率在0.5 Hz以下． 同樣在分析汶川地震前較長(zhǎng)時(shí)段的ULF電場(chǎng)數(shù)據(jù)時(shí)， Walker等（2013）發(fā)現(xiàn)DEMETER衛(wèi)星通過震中上空附近區(qū)域時(shí)記錄到的ULF波動(dòng)事件增加. 頻譜分析表明ULF波動(dòng)為地磁脈動(dòng)事件， 可能是由于震中附近產(chǎn)生的重力波影響所致． Zhang等（2014）選擇印尼和智利兩個(gè)地震區(qū)域中27個(gè)地震事件， 利用ULF電場(chǎng)觀測(cè)進(jìn)行統(tǒng)計(jì)研究， 發(fā)現(xiàn)ULF電場(chǎng)的擾動(dòng)基本小于10 mV/m， 擾動(dòng)位于震中上空．

綜上， 關(guān)于ULF擾動(dòng)的研究均針對(duì)單個(gè)地震事件或小區(qū)域內(nèi)地震事件的研究， 對(duì)于DEMETER衛(wèi)星觀測(cè)的ULF電場(chǎng)擾動(dòng)的全球分布及其特征缺乏全面的認(rèn)識(shí)． 本文旨在提出DEMETER衛(wèi)星觀測(cè)的ULF電場(chǎng)擾動(dòng)預(yù)處理方法， 并以2010年2月27日智利MS8.8地震為例基于該方法分析全球ULF電場(chǎng)擾動(dòng)特征， 以便為后續(xù)空間ULF擾動(dòng)研究提供參考．

1 DEMETER衛(wèi)星及ULF電場(chǎng)觀測(cè)數(shù)據(jù)

法國(guó)DEMETER衛(wèi)星的主要科學(xué)目標(biāo)是研究由地震電磁效應(yīng)引起的電離層擾動(dòng)以及由人類活動(dòng)（如地面電力系統(tǒng)的諧頻輻射， 地面甚低頻發(fā)射機(jī)發(fā)射的VLF波動(dòng)等）引起的電離層擾動(dòng)（Cussacetal， 2006; Parrotetal， 2006）． DEMETER衛(wèi)星于2004年6月發(fā)射進(jìn)入太陽(yáng)同步軌道， 軌道高度為710 km； 2005年12月軌道高度下降為660 km． DEMETER衛(wèi)星的探測(cè)數(shù)據(jù)以半軌道文件組織， 上行的半軌越過赤道的時(shí)間對(duì)應(yīng)于夜間22:30 LT， 下行的半軌對(duì)應(yīng)于白天10:30 LT．

為了研究與地震事件相關(guān)的電離層中的電磁擾動(dòng)， DEMETER衛(wèi)星上裝載了電場(chǎng)探測(cè)儀． 電場(chǎng)探測(cè)儀由4個(gè)球形的傳感器組成， 安裝于4個(gè)伸桿的末端， 衛(wèi)星上電場(chǎng)探測(cè)數(shù)據(jù)的信號(hào)處理由BANT模塊完成（Berthelieretal， 2006）． 電場(chǎng)探測(cè)儀提供4個(gè)頻段的觀測(cè)數(shù)據(jù)， 分別為DC/ULF（0—15 Hz）， ELF（15 Hz—1 kHz）， VLF（15 Hz—17.4 kHz）和高頻段（10 kHz—3.175 MHz）． 根據(jù)衛(wèi)星的巡查（survey）和詳查（burst）模式， 地面接收到不同頻段的電場(chǎng)數(shù)據(jù)為波形或功率譜數(shù)據(jù)． DC/ULF頻段可提供巡查和詳查模式下4通道電勢(shì)的波形數(shù)據(jù)， 通過兩兩傳感器的電壓相減得到電勢(shì)差進(jìn)而可得到3個(gè)方向的電場(chǎng)數(shù)據(jù)． DC/ULF頻段采樣率為39.0625 Hz， 以16位進(jìn)行數(shù)字化， 因而電勢(shì)測(cè)量的分辨率約為0.3 mV， 對(duì)應(yīng)于電場(chǎng)的分辨率約為40 μV/m （Berthelieretal， 2006）．

2 ULF電場(chǎng)波形擾動(dòng)預(yù)處理方法

由于電場(chǎng)探測(cè)儀球形傳感器的表面功函數(shù)不一致而導(dǎo)致不同傳感器的表面電勢(shì)存在差異， 這是引起DC/ULF電場(chǎng)觀測(cè)誤差的主要原因． 對(duì)于自旋穩(wěn)定衛(wèi)星， 兩個(gè)傳感器之間的接觸電勢(shì)差異表現(xiàn)為準(zhǔn)周期信號(hào)， 較容易消除． 而DEMETER衛(wèi)星是三軸穩(wěn)定衛(wèi)星， 很難修正兩個(gè)傳感器之間的接觸電勢(shì)差異． 由此引起大尺度準(zhǔn)DC電場(chǎng)觀測(cè)的絕對(duì)誤差范圍從幾毫伏每米至十毫伏每米． 不過， 由于各個(gè)傳感器的表面電勢(shì)在時(shí)間上是相對(duì)穩(wěn)定的， 小尺度準(zhǔn)DC電場(chǎng)作為疊加在大尺度電場(chǎng)上的擾動(dòng)可以準(zhǔn)確測(cè)量出來（Berthelieretal， 2006）． 緩慢變化的電場(chǎng)信號(hào)主要來源于V×B感應(yīng)電場(chǎng)， 這是由于衛(wèi)星橫跨地球磁力線運(yùn)動(dòng)產(chǎn)生的．

本文提出了更為簡(jiǎn)單實(shí)用的去除DC/ULF頻段趨勢(shì)波形的方法． 由于在空間物理和地震電磁擾動(dòng)研究中， 更關(guān)心的是電場(chǎng)擾動(dòng)而并非空間中電場(chǎng)的絕對(duì)數(shù)值， 因此我們直接對(duì)原始波形進(jìn)行滑動(dòng)平均得到沿此半軌道的趨勢(shì)波形， 利用原始波形減去趨勢(shì)波形， 即得到去除了V×B大尺度感應(yīng)電場(chǎng)而留下小尺度擾動(dòng)電場(chǎng)的波形． 滑動(dòng)平均方法可表示為

（1）

圖1 DEMETER衛(wèi)星DC/ULF頻段Ex， Ey和Ez分量波形示例（2010年2月17日301241半軌）

Δei=ei-Ei,

（2）

式中，ei為原始波形序列，Ei為通過滑動(dòng)平均后得到的趨勢(shì)波形， Δei為擾動(dòng)波形序列. 對(duì)于兩端的取值， 如滑動(dòng)平均窗口取為5， 則E1=e1,E2=（e1+e2+e3）/3；E3=（e1+e2+e3+e4+e5）/5；E4=（e2+e3+e4+e5+e6）/5等． 滑動(dòng)平均窗口（即2n+1）的選取將直接影響數(shù)據(jù)的平滑效果， 若2n+1取值較大， 則平滑作用較大， 有利于抑制高頻變化； 若2n+1取值較小， 則不利于抑制高頻變化從而較難得到大趨勢(shì)性的波形． 實(shí)際中可根據(jù)數(shù)據(jù)的變化情況合理地選擇滑動(dòng)平均窗口的大?。?本文中我們選取的滑動(dòng)平均窗口為2049點(diǎn)， 該采樣數(shù)據(jù)長(zhǎng)度持續(xù)時(shí)間約為52 s， 即滑動(dòng)平均后得到的趨勢(shì)波形反映的是大于52 s的周期變化． 該趨勢(shì)波形相對(duì)于V×B電場(chǎng)仍然反映了非常短周期的變化， 因而更有利于突出電場(chǎng)的小擾動(dòng)變化． 圖1給出了2010年2月17日一條半軌道的示例. 可以看到: 圖1a中的Ex和Ez分量在南北磁緯20°—40°之間均出現(xiàn)了擾動(dòng)信號(hào)． 由于沿半軌道電場(chǎng)變化范圍很大， 兩個(gè)分量的擾動(dòng)相比整個(gè)半軌道電場(chǎng)變化非常微弱， 擾動(dòng)形態(tài)也不清晰； 圖1b中各分量的波形趨勢(shì)與圖1a中各分量的變化一致． 趨勢(shì)波形是否反映原始波形的細(xì)節(jié)取決于滑動(dòng)平均采用的波形數(shù)據(jù)長(zhǎng)度． 圖1c為由圖1a與圖1b對(duì)應(yīng)分量相減后得到的擾動(dòng)波形， 可以看到：Ex，Ey和Ez分量同步出現(xiàn)了電場(chǎng)擾動(dòng)， 其范圍和形態(tài)可以清晰地從圖中識(shí)別， 在磁緯10°S—40°S之間出現(xiàn)了振幅約5 mV/m的電場(chǎng)擾動(dòng)； 在其共軛區(qū)， 約20°N—40°N區(qū)域出現(xiàn)了振幅略小于5 mV/m的擾動(dòng)． 由此示例結(jié)果可知， 本文采用的ULF電場(chǎng)波形擾動(dòng)預(yù)處理方法簡(jiǎn)單易實(shí)施， 并且不會(huì)引入新的誤差． 處理后的擾動(dòng)波形更有助于開展后續(xù)ULF電場(chǎng)擾動(dòng)分析．

3 ULF電場(chǎng)擾動(dòng)分析

2010年2月27日6時(shí)34分（UT）在靠近智利海岸發(fā)生了MS8.8地震， 震源位置為（35.8°S， 72.7°W）， 震源深度為33 km （中國(guó)地震臺(tái)網(wǎng)中心， 2010）． 多篇文獻(xiàn)報(bào)道了該地震前DEMETER衛(wèi)星觀測(cè)到了ULF電場(chǎng)擾動(dòng) （P?aetal， 2011; Zhangetal， 2014）．

圖2 2010年2月17日（智利地震前10天）全球的上行半軌道半軌道號(hào)標(biāo)示于各半軌的右側(cè)， 星形為智利地震的震中位置

圖3 2010年2月17日（智利地震前10天）夜間301271半軌DC/ULF頻段Ex， Ey和Ez分量的原始波形（a）和去除趨勢(shì)后的電場(chǎng)擾動(dòng)波形（b）

圖4 2010年2月26日（智利地震前1天）全球的上行半軌道半軌道號(hào)標(biāo)示于各半軌的右側(cè)， 星形為智利地震的震中位置

針對(duì)這次智利地震， Zhang等（2014）同樣報(bào)道了該地震前一天（2010年2月26日）上行半軌302561在20°S—40°S觀測(cè)到了ULF電場(chǎng)擾動(dòng)． 基于未去除V×B感應(yīng)電場(chǎng)的原始ULF電場(chǎng)波形， Zhang等（2014）得到Ex和Ez分量最大的相對(duì)變化幅度達(dá)到15 mV/m． 同樣， 為了檢驗(yàn)這一半軌出現(xiàn)的ULF電場(chǎng)擾動(dòng)是否與智利地震相關(guān)， 本文同樣掃描了2010年2月26日全球的上行半軌． 圖4給出了這一天所有半軌道的分布． 可以看到， 302561半軌是距震中最近的軌道． 查看全部半軌道的擾動(dòng)波形后， 仍然發(fā)現(xiàn)并非只有302561半軌在20°S—40°S緯度帶內(nèi)出現(xiàn)ULF電場(chǎng)擾動(dòng)， 302571和302581半軌在相同緯度帶亦觀測(cè)到了ULF電場(chǎng)擾動(dòng)．

目前文獻(xiàn)中只報(bào)道了智利地震震前10天和震前1天出現(xiàn)的ULF電場(chǎng)擾動(dòng)． 通過上述分析發(fā)現(xiàn)， ULF異常擾動(dòng)并非只出現(xiàn)在最靠近震中的半軌道． 為了更充分地了解ULF電場(chǎng)擾動(dòng)分布的區(qū)域和特征， 本文進(jìn)一步分析了智利地震前2天至震前9天（即2010年2月18日至2月25日）全球上行半軌道觀測(cè)的ULF電場(chǎng)擾動(dòng)情況． 通過全面分析文獻(xiàn)中未報(bào)道的ULF異常擾動(dòng)的8天中的ULF電場(chǎng)擾動(dòng)波形， 發(fā)現(xiàn)這8天在南大西洋異常區(qū)同樣觀測(cè)到了大范圍的ULF電場(chǎng)擾動(dòng)． 該期間Kp指數(shù)均小于3+， 表明地磁活動(dòng)平靜（World Data Center for Geomagnetism， Kyoto， 2010）． 在這8天中觀測(cè)到的ULF電場(chǎng)擾動(dòng)的分布特征為： 擾動(dòng)出現(xiàn)在20°S—40°S及其共軛緯度帶， 擾動(dòng)的半軌道集中在210°E—330°E范圍內(nèi)．

4 討論與結(jié)論

結(jié)合文獻(xiàn)中報(bào)道的2010年智利MS8.8地震前ULF電場(chǎng)擾動(dòng)的情況， 本文更全面地檢驗(yàn)了觀測(cè)到ULF電場(chǎng)擾動(dòng)的區(qū)域， 發(fā)現(xiàn)并非只在最靠近震中的半軌道出現(xiàn)ULF電場(chǎng)擾動(dòng)現(xiàn)象． 20°S—40°S緯度帶內(nèi)出現(xiàn)ULF電場(chǎng)擾動(dòng)的半軌道位于210°E—330°E范圍內(nèi)． 這很可能受到南大西洋異常區(qū)的影響． 從全球地磁場(chǎng)強(qiáng)度分布圖（見Heirtzler（2002）一文中圖1）上可以看到， 在巴西南部海岸附近存在地磁場(chǎng)強(qiáng)度的極小值， 這一弱磁場(chǎng)區(qū)域被稱為南大西洋異常區(qū)（Heirtzler， 2002）． 該區(qū)域分布較廣， 緯度范圍為60°S—0°N， 經(jīng)度范圍為120°W—30°E． 而地磁場(chǎng)強(qiáng)度極小值集中在約20°S—40°S、 40°W—60°W（300°E —320°E）的中心區(qū)域內(nèi)． 這與上述智利地震文獻(xiàn)中報(bào)道的觀測(cè)到ULF電場(chǎng)擾動(dòng)的區(qū)域一致． 本文的研究結(jié)果發(fā)現(xiàn)ULF電場(chǎng)擾動(dòng)出現(xiàn)在比震中經(jīng)度更廣的范圍內(nèi)． 例如， 2010年2月17日ULF電場(chǎng)擾動(dòng)的范圍為210°E—330°E， 且都出現(xiàn)在20°S—40°S緯度帶內(nèi)． 這正好位于南大西洋異常區(qū)． 因此， 大范圍觀測(cè)到的ULF電場(chǎng)擾動(dòng)很可能受南大西洋異常區(qū)的影響． 智利地震剛好發(fā)生在南大西洋異常的中心區(qū)域， 很難判斷在此區(qū)域觀測(cè)到的ULF電場(chǎng)擾動(dòng)是否受該地震影響． 同時(shí)， 本文針對(duì)智利地震震前2—9天（該期間未見文獻(xiàn)報(bào)道ULF異常擾動(dòng)）， 利用本文提出的ULF波形預(yù)處理方法得到全球上行半軌道的ULF電場(chǎng)擾動(dòng)波形， 發(fā)現(xiàn)這8天中同樣觀測(cè)到了大范圍的ULF電場(chǎng)擾動(dòng)． 該擾動(dòng)出現(xiàn)在20°S—40°S及其共軛緯度帶內(nèi)， 擾動(dòng)的半軌道集中在210°E—330°E的經(jīng)度范圍內(nèi)． 這一結(jié)果與P?a等（2011）震前10天和Zhang等（2014）震前1天的結(jié)果相同， 即位于南大西洋異常區(qū)內(nèi)的半軌道均觀測(cè)到了ULF電場(chǎng)擾動(dòng)． 連續(xù)多天全球半軌道觀測(cè)的ULF電場(chǎng)擾動(dòng)分析表明， ULF電場(chǎng)擾動(dòng)在時(shí)間上連續(xù)出現(xiàn)， 擾動(dòng)的空間范圍位于20°S—40°S、 210°E—330°E． 這一擾動(dòng)特征表明觀測(cè)到的ULF電場(chǎng)擾動(dòng)在20°S—40°S、 210°E—330°E的區(qū)域內(nèi)具有穩(wěn)定的擾動(dòng)來源． 該擾動(dòng)空間分布范圍與南大西洋異常區(qū)重合， 結(jié)合目前對(duì)南大西洋異常區(qū)的認(rèn)識(shí)， 該擾動(dòng)的來源應(yīng)該是南大西洋異常區(qū)長(zhǎng)期存在的地磁場(chǎng)強(qiáng)度極小值．

此外， 本文發(fā)現(xiàn)了另一個(gè)有趣的現(xiàn)象． 在2010年2月17日觀測(cè)到ULF電場(chǎng)擾動(dòng)的半軌中， 301231， 301241和301251在20°S—40°S出現(xiàn)了ULF電場(chǎng)擾動(dòng)， 在其共軛區(qū)（約20°N—40°N）出現(xiàn)了稍弱的ULF電場(chǎng)擾動(dòng)． 其具體變化形態(tài)可參見圖1c， 圖中給出了301241半軌Ex，Ey和Ez分量的擾動(dòng)波形． 而301261和301271則沒有這種共軛擾動(dòng)現(xiàn)象， 其ULF電場(chǎng)擾動(dòng)只出現(xiàn)在20°S—40°S的區(qū)域． 同樣， 在分析2010年2月26日全球半軌的擾動(dòng)中， 發(fā)現(xiàn)302561和302571半軌出現(xiàn)了共軛擾動(dòng)現(xiàn)象， 302581半軌只在20°S—40°S區(qū)域內(nèi)出現(xiàn)ULF電場(chǎng)擾動(dòng)現(xiàn)象． 2010年2月17日和2010年2月26日擾動(dòng)半軌道的共同特點(diǎn)為： 位于南大西洋異常區(qū)的半軌道均在20°S—40°S緯度帶內(nèi)觀測(cè)到ULF電場(chǎng)擾動(dòng)； 位于南大西洋異常中心區(qū)域（地磁場(chǎng)強(qiáng)度極小的小區(qū)域）的半軌道觀測(cè)到了ULF共軛擾動(dòng)， 南大西洋異常中心區(qū)域外未見共軛擾動(dòng)現(xiàn)象． 這種特點(diǎn)可能是由于在南大西洋異常中心區(qū)域ULF電場(chǎng)擾動(dòng)幅度較大， 可以沿磁力線傳播至其共軛區(qū)； 而在南大西洋異常中心區(qū)域外， 由于ULF電場(chǎng)擾動(dòng)幅度較小， 在電離層中由于碰撞效應(yīng)使得波動(dòng)更易被衰減， 因此無法傳播至共軛區(qū)．

綜上所述， 本文得到如下結(jié)論：

1） 本文提出了ULF電場(chǎng)擾動(dòng)的預(yù)處理方法， 與P?a等（2011）提出的方法相比， 更簡(jiǎn)單易行． 通過預(yù)處理后， 可以高效地查看沿半軌道的ULF電場(chǎng)擾動(dòng)情況．

2） ULF電場(chǎng)擾動(dòng)并非只在最靠近智利地震震中的半軌觀測(cè)到， 在較大經(jīng)度范圍內(nèi)多條半軌均出現(xiàn)了ULF電場(chǎng)擾動(dòng)． 通過進(jìn)一步分析， 發(fā)現(xiàn)這一較大范圍的ULF電場(chǎng)擾動(dòng)很可能受到南大西洋異常的影響． 同時(shí)， 南大西洋異常中心區(qū)域表現(xiàn)為ULF共軛擾動(dòng)， 中心區(qū)域外ULF電場(chǎng)擾動(dòng)僅出現(xiàn)在20°S—40°S緯度帶內(nèi)．

本文研究表明， 空間中的ULF電場(chǎng)擾動(dòng)可能受到多種因素的影響， 并非在地磁平靜期震中上空或附近區(qū)域出現(xiàn)ULF電場(chǎng)擾動(dòng)即認(rèn)為可能與地震相關(guān)， 還需要對(duì)ULF電場(chǎng)擾動(dòng)的全球分布特征有更全面的認(rèn)識(shí)． 本文結(jié)果可為后續(xù)開展衛(wèi)星觀測(cè)的震前ULF電場(chǎng)擾動(dòng)研究提供參考．

何宇飛， 楊冬梅， 陳化然， 錢家棟， 朱榮， Parrot M. 2009. DEMETER衛(wèi)星探測(cè)到可能與汶川地震有關(guān)的地面VLF發(fā)射站信號(hào)的信噪比變化[J]. 中國(guó)科學(xué): D輯， 39（4）: 403--412.

He Y F， Yang D M， Chen H R， Qian J D， Zhu R， Parrot M. 2009. SNR changes of VLF radio signals detected onboard the DEMETER satellite and their possible relationship to the Wenchuan earthquake[J].ScienceinChina:SeriesD， 52（6）: 754--763.

劉靜， 萬(wàn)衛(wèi)星， 黃建平， 張學(xué)民， 趙庶凡， 歐陽(yáng)新艷， 澤仁志瑪. 2011. 智利8.8級(jí)地震的震前電子濃度擾動(dòng)[J]. 地球物理學(xué)報(bào)， 54（11）: 2717--2725.

Liu J， Wan W X， Huang J P， Zhang X M， Zhao S F， Ouyang X Y， Zeren Z M. 2011. Electron density perturbation before ChileM8.8 earthquake[J].ChineseJournalofGeophysics， 54（11）: 2717--2725 （in Chinese）.

歐陽(yáng)新艷， 張學(xué)民， 申旭輝， 劉靜， 錢家棟， 蔡晉安， 趙庶凡. 2008. 普洱地震前電離層電子密度擾動(dòng)變化研究[J]. 地震學(xué)報(bào)， 30（4）: 424--436.

Ouyang X Y， Zhang X M， Shen X H， Liu J， Qian J D， Cai J A， Zhao S F. 2008. Study on ionosphericNedisturbances before 2007 Pu’er， Yunnan of China， earthquake[J].ActaSeismologicaSinica， 30（4）: 424--436 （in Chinese）.

歐陽(yáng)新艷， 張學(xué)民， 申旭輝， 黃建平， 劉靜， 澤仁志瑪， 趙庶凡. 2011. DEMETER衛(wèi)星探測(cè)到的強(qiáng)震前O+濃度變化[J]. 空間科學(xué)學(xué)報(bào)， 31（5）: 607--617.

Ouyang X Y， Zhang X M， Shen X H， Huang J P， Liu J， Zeren Z M， Zhao S F. 2011. Disturbance of O+density before major earthquake detected by DEMETER satellite[J].ChineseJournalofSpaceScience， 31（5）: 607--617 （in Chinese）.

顏蕊， 王蘭煒， 胡哲， 劉大鵬， 張興國(guó)， 張宇. 2013. 利用DEMETER衛(wèi)星數(shù)據(jù)分析強(qiáng)震前后的電離層異常[J]. 地震學(xué)報(bào)， 35（4）: 498--511.

Yan R， Wang L W， Hu Z， Liu D P， Zhang X G， Zhang Y. 2013. Ionospheric disturbances before and after strong earthquakes based on DEMETER data[J].ActaSeismologicaSinica， 35（4）: 498--511 （in Chinese）.

澤仁志瑪， 申旭輝， 曹晉濱， 張學(xué)民， 黃建平， 劉靜， 歐陽(yáng)新艷， 趙庶凡. 2012. 強(qiáng)震前 ELF/VLF磁場(chǎng)的擾動(dòng)特征統(tǒng)計(jì)研究[J]. 地球物理學(xué)報(bào)， 55（11）: 3699--3708.

Zenren Z M， Shen X H， Cao J B， Zhang X M， Huang J P， Liu J， Ouyang X Y， Zhao S F. 2012. Statistical analysis of ELF/VLF magnetic field disturbances before major earthquakes[J].ChineseJournalofGeophysics， 55（11）: 3699--3708 （in Chinese）.

曾中超， 張蓓， 方廣有， 王東峰， 陰和俊. 2009. 利用DEMETER衛(wèi)星數(shù)據(jù)分析汶川地震前的電離層異常[J]. 地球物理學(xué)報(bào)， 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].ChineseJournalofGeophysics， 52（1）: 11--19 （in Chinese）.

張學(xué)民， 申旭輝， 歐陽(yáng)新艷， 蔡晉安， 黃建平， 劉靜， 趙庶凡. 2009a. 汶川8級(jí)地震前空間電離層VLF電場(chǎng)異?，F(xiàn)象[J]. 電波科學(xué)學(xué)報(bào)， 24（6）: 1024--1032.

Zhang X M， Shen X H， Ouyang X Y， Cai J A， Huang J P， Liu J， Zhao S F. 2009a. Ionosphere VLF electric field anomalies before WenchuanM8 earthquake[J].ChineseJournalofRadioScience， 24（6）: 1024--1032 （in Chinese）.

張學(xué)民， 錢家棟， 歐陽(yáng)新艷， 蔡晉安， 劉靜， 申旭輝， 趙庶凡. 2009b. 新疆于田7.2級(jí)地震前的電離層電磁擾動(dòng)[J]. 空間科學(xué)學(xué)報(bào)， 29（2）: 213--221.

Zhang X M， Qian J D， Ouyang X Y， Cai J A， Liu J， Shen X H， Zhao S F. 2009b. Ionospheric electro-magnetic distur-bances prior to Yutian 7.2 earthquake in Xinjiang[J].ChineseJournalofSpaceScience， 29（2）: 213--221 （in Chinese）.

中國(guó)地震臺(tái)網(wǎng)中心. 2010. 中國(guó)地震臺(tái)網(wǎng)地震目錄[EB/OL]. [2015-04-01]. http:∥www.csndmc.ac.cn/newweb/cgi-bin/csndmc/csn_catalog_p003.pl?mode=catalog&ot=20100227_0634140.

China Earthquake Networks Center. 2010. Earthquake catalogues from China Seismogical Networks[EB/OL]. [2015-04-01]. http://www.csndmc.ac.cn/newweb/cgi-bin/csndmc/csn_catalog_p003.pl?mode=catalog&ot=20100227_0634140 （in Chinese）.

朱濤， 王蘭煒. 2011. DEMETER衛(wèi)星觀測(cè)到的與汶川地震有關(guān)的LF電場(chǎng)異常[J]. 地球物理學(xué)報(bào)， 54（3）: 717--727.

Zhu T， Wang L W. 2011. LF electric field anomalies related to Wenchuan earthquake observed by DEMETER satellite[J].ChineseJournalofGeophysics， 54（3）: 717--727 （in Chinese）.

Akhoondzadeh M， Parrot M， Saradjian M R. 2010. Investigation of VLF and HF waves showing seismo-ionospheric anomalies induced by the 29 September 2009 Samoa earthquake （MW=8.1）[J].NatHazardsEarthSystSci， 10（5）: 1061--1067.

Athanasiou M A， Anagnostopoulos G C， Iliopoulos A C， Pavlos G P， David C N. 2011. Enhanced ULF radiation observed by DEMETER two months around the strong 2010 Haiti earthquake[J].NatHazardsEarthSystSci， 11（4）: 1091--1098.

Berthelier J J， Godefroy M， Leblanc F， Malingre M， Menvielle M， Lagoutte D， Brochot J Y， Colin F， Elie F， Legendre C， Zamora P， Benoist D， Chapuis Y， Artru J， Pfaff R. 2006. ICE， the electric field experiment on DEMETER[J].PlanetSpaceSci， 54（5）: 456--471.

Bhattacharya S， Sarkar S， Gwal A， Parrot M. 2007. Observations of ULF/ELF anomalies detected by DEMETER satellite prior to earthquakes[J].IndianJRadioSpacePhys， 36（2）: 103--113.

Chmyrev V M， Isaev N V， Bilichenko S V， Stanev G. 1989. Observation by space-borne detectors of electric fields and hydromagnetic waves in the ionosphere over an earthquake centre[J].PhysEarthPlanetInt， 57（1/2）: 110--114.

Cussac T， Clair M A， Ultré-Guerard P， Buisson F， Lassalle-Balier G， Ledu M， Elisabelar C， Passot X， Rey N. 2006. The DEMETER microsatellite and ground segment[J].PlanetSpaceSci， 54（5）: 413--427.

Fraser-Smith A C， Bernardi A， McGill P R， Ladd M E， Helliwell R A， Villard O G. 1990. Low-frequency magnetic field measurements near the epicenter of theMS7.1 Loma Prieta earthquake[J].GeophysResLett， 17（9）: 1465--1468.

Gokhberg M B， Morgounov V A， Yoshino T， Tomizawa I. 1982. Experimental measurement of electromagnetic emissions possibly related to earthquakes in Japan[J].JGeophysRes， 87（B9）: 7824--7828.

Hayakawa M， Kawate R， Molchanov O A， Yumoto K. 1996. Results of ultra-low-frequency magnetic field measurements during the Guam earthquake of 8 August 1993[J].GeophysResLett， 23（3）: 241--244.

Hayakawa M， Hattori K， Ohta K. 2007. Monitoring of ULF （ultra-low-frequency） geomagnetic variations associated with earthquakes[J].Sensors， 7（7）: 1108--1122.

Heirtzler J R. 2002. The future of the South Atlantic anomaly and implications for radiation damage in space[J].JAtmosSolar-TerrPhys， 64（16）: 1701--1708.

Molchanov O A， Kopytenko Y A， Voronov P M， Kopytenko E A， Matiashvili T G， Fraser-Smith A C， Bernardi A. 1992. Results of ULF magnetic field measurements near the epicenters of the Spitak （MS=6.9） and Loma Prieta （MS=7.1） earthquakes: Comparative analysis[J].GeophysResLett， 19（14）: 1495--1498.

Molchanov O A， Mazhaeva O A， Golyavin A N， Hayakawa M. 1993. Observation by the Intercosmos-24 satellite of ELF-VLF electromagnetic emissions associated with earthquakes[J].AnnGeophys， 11（5）: 431--440.

Němec F， Santolík O， Parrot M， Berthelier J J. 2008. Spacecraft observations of electromagnetic perturbations connected with seismic activity[J].GeophysResLett， 35（5）: L05109.

Onishi T， Berthelier J J, Kamogawa M. 2011. Critical analysis of the electrostatic turbulence enhancements observed by DEMETER over the Sichuan region during the earthquake preparation[J].NatHazardsEarthSystSci， 11（2）: 561--570.

Parrot M. 1994. Statistical study of ELF/VLF emissions recorded by a low-altitude satellite during seismic events[J].JGeophysRes:SpacePhys， 99（A12）: 23339--23347.

Parrot M. 1995. Electromagnetic noise due to earthquakes[G]∥Volland H eds.HandbookofAtmosphericElectro-dynamics,PartⅡ. Boca Raton: CRC Press: 95--116.

Parrot M， Benoist D， Berthelier J J， Blecki J， Chapuis Y， Colin F， Elie F， Fergeau P， Lagoutte D， Lefeuvre F， Legen-dre C， Lévêque M， Pin?on J L， Poirier B， Seran H C， Zamora P. 2006. The magnetic field experiment IMSC and its data processing onboard DEMETER: Scientific objectives， description and first results[J].PlanetSpaceSci， 54（5）: 441--455.

Walker S N， Kadirkamanathan V， Pokhotelov O A. 2013. Changes in the ultra-low frequency wave field during the precursor phase to the Sichuan earthquake: DEMETER observations[J].AnnGeophys， 31（9）: 1597--1603.

World Data Center for Geomagnetism， Kyoto. 2010. Geomagnetic Data Service，Kpindex[EB/OL]. [2015-04-01]. http:∥wdc.kugi.kyoto-u.ac.jp/kp/index.html.

Zhang X M， Zeren Z M， Parrot M， Battiston R， Qian J D， Shen X H. 2011. ULF/ELF ionospheric electric field and plasma perturbations related to Chile earthquakes[J].AdvSpaceRes， 47（6）: 991--1000.

Zhang X M， Shen X H， Parrot M， Zeren Z M， Ouyang X Y， Liu J， Qian J D， Zhao S F， Miao Y Q. 2012a. Phenomena of electrostatic perturbations before strong earthquakes （2005—2010） observed on DEMETER[J].NatHazardsEarthSystSci， 12（1）: 75--83.

Zhang X M， Chen H R， Liu J， Shen X H， Miao Y Q， Du X B， Qian J D. 2012b. Ground-based and satellite DC-ULF electric field anomalies around WenchuanM8.0 earthquake[J].AdvSpaceRes， 50（1）: 85--95.

Zhang X M， Shen X H， Zhao S F， Yao L， Ouyang X Y， Qian J D. 2014. The characteristics of quasistatic electric field perturbations observed by DEMETER satellite before large earthquakes[J].JAsianEarthSci， 79: 42--52.

A method for pre-processing ULF electric field disturbances observed by DEMETER and its case application analysis

1）InstituteofEarthquakeScience，ChinaEarthquakeAdministration，Beijing100036，China

2）SchoolofEarthandSpaceSciences，PekingUniversity，Beijing100871，China

3）InstituteofCrustalDynamics，ChinaEarthquakeAdministration，Beijing100085，China

This paper presents a method for pre-processing ULF electric field disturbances based on ULF waveforms observed by French DEMETER satellite. This method is simpler and easier to carry out. We analyze comprehensively ULF disturbances before the ChileMS8.8 earthquake on 27 February 2010， which have been reported in the literature， by using our method. The results show that reported ULF disturbances in 20°S—40°S geomagnetic latitude bands do not just appear on the half-orbit that is the closest to the epicenter， but ULF disturbances are observed in a wider longitude range. The large range with ULF disturbances is located in the region of South Atlantic anomaly （SAA）， and ULF disturbances may be affected by the geomagnetic field anomaly in this region. ULF disturbances in the center of SAA region show a conjugate distribution，i.e.， electric field fluctuates with the amplitude of 5 mV/m at 20°—40° north and south geomagnetic latitudes. The amplitude of ULF disturbances at northern latitudes is weaker than those at southern latitudes. However， ULF disturbances outside the center of SAA region are only observed at 20°S—40°S geomagnetic latitudes， and there are no ULF fluctuations in the conjugate region. As ULF fluctuations outside the center of SAA region with the weaker amplitude are easier to be attenuated because of collision effects in the ionosphere， ULF disturbance cannot propagate to the conjugate region. The method and results in this paper can provide references for subsequently thorough analyses on ULF fluctuations recorded by satellites in the ionosphere.

DEMETER satellite; ULF electric field disturbance; South Atlantic anomaly （SAA） region

10.11939/jass.2015.05.010.

中國(guó)地震局地震預(yù)測(cè)研究所基本科研業(yè)務(wù)費(fèi)專項(xiàng)（2015IES0405， 2014IES010203, 2012IES0401）和國(guó)家自然科學(xué)基金（41404126）共同資助.

2015-04-14收到初稿， 2015-07-02決定采用修改稿.

e-mail: oyxy@cea-ies.ac.cn

10.11939/jass.2015.05.010

P319.3+2

A

歐陽(yáng)新艷, 申旭輝. 2015. DEMETER衛(wèi)星觀測(cè)的超低頻電場(chǎng)擾動(dòng)預(yù)處理方法及震例應(yīng)用分析. 地震學(xué)報(bào), 37（5）: 820--829.

Ouyang X Y, Shen X H. 2015. A method for pre-processing ULF electric field disturbances observed by DEMETER and its case application analysis.ActaSeismologicaSinica, 37（5）: 820--829. doi:10.11939/jass.2015.05.010.