中國大陸14次強震前電離層異常統計分析

徐 彤 胡艷莉 吳 健 李春斌吳振森 索玉成 馮 健

（1.中國電波傳播研究所 電波環境特性及?；夹g重點實驗室，山東 青島266107；2.西安電子科技大學，陜西 西安710071）

引 言

2008年5月12日汶川發生的8.0級大地震給人類的生命和財產造成了重大損失，觀測事實表明：汶川地震前電離層發生了異常明顯的擾動，很多學者利用不同類型的電離層測量數據對汶川地震前電離層異常進行了報道［1－5］。

實際上地震電離層異常的研究始于20世紀60年代中期，Barnes等［6］在1964年阿拉斯加大地震時，發現電離層出現異常擾動現象，而此次的分析似乎是第一次發現電離層的異常擾動與地震兩者之間存在著某種關聯。之后全球大量的研究者開始對地震電離層擾動進行統計研究，試圖找出震前電離層擾動的規律。Furumoto等［7］發現1969年的Kurile島地震前電離層也出現類似的異常擾動現象。臺灣中央大學以劉正彥教授為首的研究團隊［8－10］所做的地震電離層擾動統計工作最具代表性。他們對臺灣地震面波震級Ms≥5.0級上百次地震進行了統計分析，得出了許多重要的統計規律。Dabas等［11］對低緯地區數次大地震進行了電離層擾動分析。Singh等人［12］對2006年9月到2007年11月印度43次5.0級以上地震進行了電離層變化統計分析。然而這些結論不盡一致。由于樣本貧乏，地震電離層擾動的統計規律仍然沒有定論。

目前俄羅斯、印度、日本、美國、墨西哥、保加利亞、法國、中國臺灣等多個國家和地區開展了地震電離層異常擾動研究［13－16］。法國和俄羅斯分別于2004年、2006年發射了DEMETER和COMPASS－2地震衛星，進行地震電離層異常擾動監測。

2005年，中國以國家863項目為標志，正式開始地震電磁衛星的研究工作，并制訂了相應發展計劃，目前該工作正在進行中，地震電離層擾動研究越來越受到重視。

我國從20世紀40年代開始建立電離層垂測站，目前中國電波傳播研究所已經建立了遍及我國中低緯地區13個電離層垂測站（各站分布如圖1所示），積累了數個太陽周的電離層觀測數據。

我們對汶川地震前電離層垂測參數進行了詳細分析，發現電離層F2層臨界頻率foF2震前響應最為明顯［17－19］。本文利用我國電離層觀測數據，對1960年以來，中國大陸地區14次Ms≥7.0級大地震前電離層foF2進行了統計分析，希望借此形成7級以上大地震電離層擾動的一般規律，并初步得到了一些有意義的結論。

圖1 我國電離層垂測站分布

1.震例及分析方法

本文選擇了我國大陸地區發生的部分7級以上地震，共14次。表1列出了這些地震的發生時間、地域、震級等信息。

表1 60年代至今發生在中國的Ms≥7.0級地震

分析方法為四分位數法［10，17，20］，即 利用震 中附近電離層觀測站foF2數據，確定震前15天以及地震當天電離層foF2的上下邊界，然后將觀測值與其進行比較以確定是否發生電離層異常擾動。由于電離層foF2的變化很復雜，與很多因素有關，在此，我們考慮了太陽和地磁活動的影響，濾除了地磁指數Dst＜－15nT對應的foF2數據。

2.結果與分析

2.1 震前電離層時間變化

地震電離層前兆研究，必須統計震前電離層的時間變化。Liu等［10］采用中位數的統計方法對1994－1999年間發生在臺灣的5級以上的地震震前的電離層數據進行研究，發現在震前5天內foF2會有顯著的變化。Dabas［11］對2003－2005年發生在低緯地區的11次6級以上地震分析發現，地震電離層擾動發生在震前2～25天之間，時間尺度較大。震前電離層的時間變化目前沒有定論。我們對選用的14次地震電離層擾動時間進行了分析。由于地震事件較少，不能用統計學上的方法來確定地震擾動時間。設X為foF2異常矩陣，且滿足

式中：i為震前天數，i∈［－15，0］；j為地震事例，j∈［1，14］（與表1中地震序號對應）。i＝0為地震當天。矩陣X生成的圖像如圖2所示。

電離層異常距離地震當日的時間尺度較小，比較集中在震前幾天，只有昆侖山地震在前14天發生了異常。其中有6次在震前第3天發生了異常，3次在震前第7天和第1天發生了異常。初步可以看出，時間尺度大概在數天時間內，71.4%的異常發生在震前7天內。Liu等［10］統計發現地震電離層擾動大多發生在前5天內，和我們統計的結果較為接近。超過10天時間尺度的，只有一例，與Dabas［11］的分析結果存在較大差異，還需進一步證實。

此外，我們分析了地震電離層發生擾動的地方時間分布。設y為各個擾動序列，有

式中：i為地方時；j為地震事例。

圖3給出了地震電離層擾動地方時分布?？梢郧逦乜闯觯?0∶00LT～23∶00LT整個時間段內，都有擾動發生，但是擾動時間大部分集中在11∶00LT～17∶00LT.其中16∶00LT和17∶00 LT發生的次數最多。汶川地震也是這個時間段內發生了異常明顯的擾動［3－4］。Liu［10］等分析得出地震電離層擾動大部分發生在12∶00LT～17∶00LT之間的結論，和我們的分析結果較為一致。

2.2 foF2異常擾動與地震震級的關系

既然地震的規模（震級）代表能量的大小，而電離層foF2異常擾動的偏移量可能與外力影響的強度有關。因此，可探討電離層foF2擾動偏移量與地震震級之間的關聯性。

利用兩種方法來考察電離層foF2擾動與地震震級之間的關系。一是利用發生擾動時間內相對最大百分比%Dev來分析，二是利用最大絕對擾動量ΔfoF2來分析。對于一次地震，有多個電離層觀測站發生異常情況，我們選取離震中相對較近的站進行分析。

圖4給出了電離層foF2絕對擾動量ΔfoF2、相對擾動量%Dev與地震震級之間的關系。對昆侖山8.1級地震的分析采用了距震中約1 000km處的烏魯木齊站的數據，距離較遠可能是造成擾動幅度偏小的主要原因，因此，分別擬合。其中二次擬合的紅色曲線和線性擬合的藍色直線為不包括昆侖山8.1級地震數據的結果?？梢钥闯?，隨著地震震級的增加，電離層的擾動量開始增加，呈現明顯的正相關關系。其中汶川地震擾動最為明顯，ΔfoF2達到4.5MHz.然而昆侖山地震，震級雖達到8.1級，絕對擾動ΔfoF2卻不是特別明顯，為2.15MHz.由于電離層的擾動量除了可能和地震的震級有關外，可能還與震源的深度、觀測站距震中的距離有關，這些因素可能共同決定著電離層擾動。從圖4（b）也大體能夠看出地震震級增大，電離層擾動越發增強，但是這種關系沒有圖4（a）明顯，數據比較散亂。雖然電離層絕對擾動很大，但是由于上邊界的參考值隨觀測站位置，特別是當地時間有顯著的變化，可能造成了百分比擾動%Dev與絕對擾動ΔfoF2存在很大的差異。這里似乎說明在討論電離層foF2擾動量時，分析絕對偏差ΔfoF2較相對偏差%Dev更合理，這還需大量的統計事例分析證實。

2.3 foF2異常擾動與震中距離的關系

由于地震震級對電離層的擾動幅度有著重要的影響，所以分析電離層響應與震中距離關系時，應該剔除地震震級的影響。另外，震源深度可能也是電離層響應的一個重要因素。這里僅以汶川地震這一次地震事例，初步說明震中距離與電離層擾動之間的關系。

地震孕育期間，電離層擾動局限在一定的空間尺度內［1，3，20］。在對汶川地震分析時，發現距離震中較遠的拉薩站電離層沒有異常擾動發生，而距震中最近的重慶站電離層擾動最為明顯。

圖5為我們利用全國13個站生成的電離層擾動二維分布圖。可見離震中較遠的區域，電離層沒有擾動。卓裕榮［20］對臺灣地區地震事件進行統計分析發現，電離層擾動幅度與震中距離成反比。圖6給出了汶川地震發生前2008年5月9日重慶、蘭州、昆明和拉薩站電離層foF2最大擾動量與震中距離之間的關系。從圖6可以看出，隨著與震中距離的增加，電離層foF2擾動幅度顯著減小，呈現出非線性的關系。

2.4 foF2異常擾動與震源深度的關系

對同一震級的地震，震源深度不同，釋放到地球表面的能量會有所不同，可能對電離層產生不同程度的影響。卓裕榮［20］通過對臺灣地區184個5級以上地震統計分析發現，深度80km以上的地震，其foF2異常擾動量隨著深度的增加有減小的趨勢。由于地震事件樣本有限，這14次地震大多是淺表地震，所以很難進行統計分析。但是我們發現在吉林琿春7.0級地震前，離震中僅有400km的長春站電離層沒有出現異常擾動，而吉林琿春地震的震源深度為578km，震源深度太深，可能是這次地震電離層沒有異常擾動現象的一個重要原因。這似乎說明震源深度可能與電離層foF2異常擾動幅度之間存在著某種復雜的負相關關系。

3.總結與討論

地震電離層異常擾動現象報道眾多，我們對我國14次大地震的分析進一步證實了地震電離層前兆的存在，然而地震電離層前兆規律目前仍然沒有定論。

就這14次大地震而言，有85.7%的震例出現了電離層擾動現象，而且基本為正相擾動；電離層擾動集中發生在震前7天以內的11∶00LT～17∶00LT；電離層foF2最大絕對擾動量基本大于1.0MHz，與即將發生的地震震級存在明顯正相關關系，與震中距離、震源深度呈現負相關關系。

眾所周知，電離層受太陽、地磁、氣象等多種因素的影響，變化復雜［21］。在磁寧靜無震條件下，電離層亦有明顯的逐日變化。Mikhailov等［22］對磁靜日電離層正負相擾動特征進行了深入的研究，發現中低緯地區電離層磁靜日正相擾動主要發生在夜間，擾動空間尺度大。這與我們統計的震前電離層擾動特征（白天，局部性）有明顯不同。同時，局部擾動特征也是震前電離層擾動與磁暴期間電離層擾動的一個重要區別。

然而，地震激發電離層異常擾動的機制目前尚不明確。主要有兩種假說：一是地震前地殼運動產生的重力波（AGW）影響；二是地震孕育區上方的異常電場透入電離層產生電離層擾動。前者存在著諸多局限［23］。Pulinets等［23］指出，聲重波向上傳輸過程中受到中性風、熱傳導及其他擴散過程的影響，很難傳到電離層高度。由于地震孕育區上方出現異常電場己經多次被觀測到［24－25］，異常電場的假說逐漸被越來越多的研究者接受。但是目前電動力學模型得到的電離層高度穿透電場強度偏小，不足以解釋震前電離層異常明顯的擾動，異常電場假說仍需研究。

地震電離層異常的研究對形成地震預測預報方法有重要的意義。但是由于地震和電離層耦合的復雜性和電離層異常狀態的多樣性，若將電離層異常作為地震前兆用于實際的地震預報中，在機制和前兆信息識別和提取方面有許多工作要做。本文所選的地震樣本較少，統計規律還有待進一步完善。地震電離層前兆規律沒有定論，還需要進行大量的統計，另外識別地震電離層前兆的方法也需深入研究。應用空間技術開展大區域范圍的電離層的實時監測，如GPS臺網、地震電離層衛星，并與電離層垂測、斜測儀觀測相結合，建立天—地基一體化的地震電離層前兆監測系統，無疑會加速這項研究工作的進程。

［1］LIU J Y，CHEN Y I，CHEN C H，et al.Seismoionospheric GPS total electron content anomalies observed before the 12May 2008Mw7.9Wenchuan earthquake［J］.J Geophys Res，2009，114 （A04320）.doi：10.1029/2008JA013698.

［2］ZHAO B Q，WANG M，YU T，et al.Is an unusual large enhancement of ionospheric electron density linked with the 2008great Wenchuan earthquake？［J］.J Geophys Res，2008，113（A11304）.doi：10.1029/2008JA013613.

［3］PULINETS S A，BONDUR V G.，TSIDILINA M N，et al.Verification of the concept of seismo－ionospheric coupling under quiet helio－geomagnetic conditions using the Wenchuan （China）earthquake of May 12，2008，as an example［J］.Geomagn Aeron，2010，50（2）：231－242.

［4］ZHOU Y Y，Wu Y，QIAO X J.Ionospheric anomalies detected by ground－based GPS before the Mw 7.9 Wenchuan earthquake of May 12，2008，China［J］.J Atmos Solar－Terr Phys，2009，71（8/9）：959－966.

［5］張學民，申旭輝，歐陽新艷，等.汶川8級地震前空間電離層VLF電場異常現象［J］.電波科學學報，2009，24（6）：1024－1032.ZHANG Xuemin，SHEN Xuhui，OUYANG Xiayan，et al.Ionosphere VLF electric field anomalies before Wenchuan M8earthquake［J］.Chinese Journal of Radio Science，2009，24（6）：1024－1032.（in Chinese）

［6］BARNES R A，LEONARD R S.Observations of ionospheric disturbances following the Alaska earthquake［J］.J Geophys Res，1965，70（9）：1250－1253.

［7］FURUMOTO A S，PROLSS G W，WEAVER P F，et al.Acoustic coupling into the ionosphere from seismic waves of the earthquake at Kuile islands on August 11，1969［J］.Nature，1971，226（5252）：1239－1241.doi：10.1038/2261239a0.

［8］HSIA C C，LIU J Y，OYAMA K I，et al.Ionospheric electron density anomaly prior to the December 26，2006M7.0Pingtung earthquake doublet observed by FORMOSAT－3/COSMIC ［J］.Phys Chem Earth，2009，34（6/7）：474－478.

［9］LIU J Y，TSAI Y B，CHEN S W，et al.Giant ionospheric disturbances excited by the M9.3Sumatra earthquake of 26December 2004［J］.Geophys Res Let，33：L02103.doi：10.1029/2005GL023963.

［10］LIU J Y，CHEN Y I，CHUO Y J.A Statistical investigation of pre－earthquake ionospheric anomaly［J］.J Geophys Res，2006，111：A05304.doi：10.1029/2005JA 011333.

［11］DABAS R S，RUPESH M D，SHARMA K，et al.I－onospheric precursors observed over low latitudes during some of the recent major earthquakes［J］.J Atmos Solar－Terr Phys，2007，69（15）：1813－1824.

［12］SINGH O P，CHAUHAN V，SINGH V，et al.A－nomalous variation in total electron content（TEC）associated with earthquakes in India during September 2006－November 2007［J］.J Phys Chem Earth，2009，34（6/7）：479－484.

［13］HEGAI V V，KIM V P，LIU J Y.The Ionospheric Effect of atmospheric gravity waves excited prior to strong earthquakes［J］.Adv Space Res，2006，37（4）：653－659.

［14］TSOLIS G S，XENOS T D.A qualitative study of the seismo－ionospheric precursors prior to the 6April 2009earthquake in L＇Aquila.Italy［J］.Nat Hazards Earth Syst Sci，2010（10）：133－137.doi：10.5194/nhess－10－133－2010.

［15］SHALIMOV S L，GOKHBERG M B.Lithosphereionosphere coupling mechanism and its application to the earthquake in Iran on June 20，1990：interpretation of its ionosphere effects［J］.J Earthquake Pred Res，1998，7：98－111.

［16］OMORI Y.Preseismic alteration of atmospheric electrical conditions due to anomalous radon emanation［J］.J Phys Chen Earth，2008，34（6/7）：435－440.

［17］XU T，HU Y L，WU J，et al.Abnormal perturbations in the ionospheric F2Region before Wenchuan earthquake on 12May 2008 ［J］.Sci China （Earth Sci），2010，53（11）：1671－1674.

［18］XU T，HU Y L，WU J，et al.Giant disturbance in the ionospheric F2region prior to the M8.0Wenchuan earthquake on 12May 2008［J］.Ann Geophysicae，2010，28（8）：1533－1538.

［19］丁宗華，吳 健，孫樹計，等.汶川大地震前電離層參量的變化特征與分析［J］.地球物理學報，2010，53（1）：30－38.DING Zonghua，WU Jian，SUN Shuji，et al.The variation of ionosphere on some days before the Wenchuan earthquake［J］.Chinese Journal of Geophysics 2010，53（1）：30－38.（in Chinese）

［20］卓裕榮.電離層地震前兆之研究［D］.臺灣中壢：國立中央大學，2002年.CHUO Yurong.Study on Seismo－ionospheric Precursor［D］.Chungli Taiwan：National Central University，2002.（in Chinese）

［21］XU T，WU Z S，WU J，et al.Solar cycle variation of the monthly medianfoF2at Chongqing station ［J］.Adv Space Res，2008，42（1）：213－218.

［22］MIKHAILOV A V，DEPUEVA A K，LESCHINSKAYA T Y.Morphology of quiet time F2－layer disturbances：high to lower latitudes［J］.Int J Geomagn Aeron，2003（5）：GI1006.doi：10.1029/2003GI 000058.

［23］PULINETS S.Ionospheric precursors of earthquakes：recent advances in theory and practical applications［J］.TAO，2004，15（3）：413－435.

［24］MIKHAILOV Y M，MIKHAILOV G A，KAPUSTINA O V，et al.Variations of power spectra of electric field of the near－ground atmosphere on Kamchatka［J］.Geomagn Aeron，2003，4（3）：422－428.

［25］RULENKO O P.Operative precursors of earthquakes in the electricity of near－earth atmosphere［J］.Vokanol Seismol，2000（4）：57－68.