姚宜斌，翟長治，孔　建，劉　磊

1. 武漢大學(xué)測繪學(xué)院，湖北 武漢 430079； 2. 武漢大學(xué)地球空間環(huán)境與大地測量教育部重點實驗室，湖北 武漢 430079； 3. 地球空間信息技術(shù)協(xié)同創(chuàng)新中心，湖北 武漢 430079； 4. 武漢大學(xué)中國南極測繪研究中心，湖北 武漢 430079

?

2015年尼泊爾地震的震前電離層異常探測

姚宜斌1,2,3，翟長治1，孔建4，劉磊1

1. 武漢大學(xué)測繪學(xué)院，湖北 武漢 430079； 2. 武漢大學(xué)地球空間環(huán)境與大地測量教育部重點實驗室，湖北 武漢 430079； 3. 地球空間信息技術(shù)協(xié)同創(chuàng)新中心，湖北 武漢 430079； 4. 武漢大學(xué)中國南極測繪研究中心，湖北 武漢 430079

摘要：提出了一種基于奇異譜分析的電離層異常探測的方法。通過對尼泊爾地震震中周圍GIM格網(wǎng)點TEC時間序列的探測，發(fā)現(xiàn)在2015年4月23日震中東部區(qū)域出現(xiàn)電離層正異常。進(jìn)一步利用二維電離層地圖分析異?？臻g分布，發(fā)現(xiàn)出現(xiàn)電離層正異常的區(qū)域為25°N—37.5°N，90°E—110°E，時間為2015年4月23日UT9:00—15:00。利用中國陸態(tài)網(wǎng)數(shù)據(jù)計算異常區(qū)域衛(wèi)星穿刺點軌跡STEC變化情況，發(fā)現(xiàn)2015年4月23日穿刺點軌跡進(jìn)入異常區(qū)域后STEC值比前后幾天明顯增大，而離開異常區(qū)域后又恢復(fù)正常。采用CIT(computerized ionosphere tomography)方法詳細(xì)地呈現(xiàn)了電離層異常的三維形態(tài)，發(fā)現(xiàn)4月23日UT9:00—15:00在震中東部區(qū)域出現(xiàn)電離層正異常，峰值位于約30°N，115°E，高度范圍為100～500 km，且異常峰值隨高度變化與電離層本身垂直密度分布規(guī)律相一致。

關(guān)鍵詞：地震-電離層異常; 奇異譜分析; STEC; CIT

2015年4月25日，尼泊爾發(fā)生了MS8.1級大地震，震中位于28.147°N，84.708°E，震源深度為15 km。地震是對人類造成傷害最大的自然災(zāi)害之一，但由于地震產(chǎn)生原因的復(fù)雜性，地震預(yù)報一直處于探索階段[1]。許多學(xué)者致力于研究震前地質(zhì)構(gòu)造變化以及與地震有關(guān)的其他震前異?，F(xiàn)象，其中震前電離層異常是研究熱點之一。

1964年美國阿拉斯加大地震時，文獻(xiàn)[2]發(fā)現(xiàn)在震中區(qū)上空出現(xiàn)了電離層異常擾動，第一次將電離層異常與地震聯(lián)系到一起。文獻(xiàn)[3]發(fā)現(xiàn)在1966年Tashkent地震時電離層TEC有增加的現(xiàn)象。文獻(xiàn)[4]發(fā)現(xiàn)在1969年Kurile島地震期間電離層出現(xiàn)擾動異常。文獻(xiàn)[5—6]發(fā)現(xiàn)1979—1981年間發(fā)生的大地震之前電離層參數(shù)有擾動現(xiàn)象，隨后統(tǒng)計了50個M>5.0的地震期間電離層變化情況，結(jié)果表明電離層F2臨界頻率(f0F2)會明顯降低，并且電離層異常區(qū)域與震區(qū)地理位置相對應(yīng)，但并不重合。文獻(xiàn)[7]發(fā)現(xiàn)1999年Chi Chi大地震引起的同震電離層擾動現(xiàn)象，并利用波束形成和光線追蹤技術(shù)確定了CID(coseismic ionospheric disturbances)傳播速度以及發(fā)源地位置。文獻(xiàn)[8]利用電離層臺站數(shù)據(jù)和地磁臺網(wǎng)數(shù)據(jù)對比分析了1997年瑪尼7.5級地震和2001年昆侖山口西8.1級地震前電磁異常和電離層異常情況。分析結(jié)果表明，兩次地震前電磁異常和電離層異?？臻g分布均具有較好的一致性，震中周圍出現(xiàn)明顯的電離層f0F2異常。

傳統(tǒng)電離層觀測手段主要有電離層探測儀、非相干散射雷達(dá)等。但是這些技術(shù)成本高且只能獲得站點上空電離層信息，不能對電離層進(jìn)行長時間、大范圍監(jiān)測。GNSS技術(shù)作為一種全天候、全自動的高精度觀測手段，能夠準(zhǔn)確獲得信號傳播路徑上的TEC值，并且在全球范圍內(nèi)共有數(shù)千個GNSS站連續(xù)觀測，相比于其他觀測技術(shù)，其時空分辨率大大提高。文獻(xiàn)[9]最早利用GPS技術(shù)對1994年Northridge地震期間電離層進(jìn)行探測，發(fā)現(xiàn)震后幾分鐘出現(xiàn)了電離層TEC擾動。文獻(xiàn)[10]利用地基GPS以及COSMIC數(shù)據(jù)對汶川地震期間TEC以及電離層F2層峰值電子濃度(NmF2)進(jìn)行功率譜分析，發(fā)現(xiàn)TEC和NmF2時空分布變化存在電離層震前擾動和同震響應(yīng)，并且TEC震后的異常變化相比震前更加顯著。文獻(xiàn)[11]分別利用Demeter衛(wèi)星和GPS對四個大地震電離層異常擾動進(jìn)行對比分析，兩種觀測技術(shù)獲得的電離層正負(fù)異常變化能夠很好地吻合。

基于GNSS技術(shù)有多種電離層異常分析方法。文獻(xiàn)[12]等利用四分位距法和滑動時窗法分析了汶川地震電離層擾動。文獻(xiàn)[13]利用GIM分析了2011年日本Tohoku Oki地震電離層異常，發(fā)現(xiàn)地震前第3 d震中偏向赤道方向電離層有明顯正異常，在赤道共軛區(qū)域也存在正異常。文獻(xiàn)[14]利用計算GPS衛(wèi)星信號電離層穿刺點軌跡STEC的方法，分析了2007年Bengkulu地震期間電離層變化情況，文獻(xiàn)[15]分析了2011年日本Tohoku-Oki地震電離層變化情況，均發(fā)現(xiàn)了震前TEC正異常。文獻(xiàn)[16]利用GIM數(shù)據(jù)對2002—2010年全球范圍內(nèi)的736次M≥6.0地震前電離層電子含量做了統(tǒng)計分析。結(jié)果表明，震前1～21 d內(nèi)電離層異常出現(xiàn)頻率與地震震級大小以及震源深度有關(guān)。近年來，隨著全球GNSS觀測站的增多，電離層層析技術(shù)(CIT)逐漸得到發(fā)展和完善，成為一種能夠重建三維甚至四維電離層形態(tài)的新方法[17-21]。文獻(xiàn)[22]利用電離層層析技術(shù)分析了中緯度地區(qū)夏季夜間電離層異?，F(xiàn)象。文獻(xiàn)[23]利用GPS和GLONASS觀測數(shù)據(jù)建立了電離層模型，分析了俄羅斯地區(qū)電離層情況。文獻(xiàn)[24]基于日本高密度的GNSS監(jiān)測網(wǎng)絡(luò)GEONET建立了電離層三維模型。

本文提出利用奇異譜分析進(jìn)行電離層異常探測的新方法。對尼泊爾地震震中周圍格網(wǎng)點TEC時間序列進(jìn)行異常探測，并結(jié)合太陽活動和地磁活動情況，初步認(rèn)定2015年4月23日在震中東部區(qū)域出現(xiàn)電離層正異常。進(jìn)而利用二維電離層地圖以及電離層3DCIT技術(shù)對此次異常的時空變化進(jìn)行研究。

1震前電離層異常探測原理

1.1基于奇異譜分析的TEC時間序列異常探測

奇異譜分析(SSA)是一種廣義的功率譜分析，不受正弦波假定的約束，對信號的識別和描述采用時域性的頻率域分析方式，能夠穩(wěn)定識別和強(qiáng)化周期信號[25]。SSA的分析對象是一維的時間序列。將時間序列x1、x2、x3、…、xN中心化處理后，按照相同的序列延遲排列得到

(1)

式中，N為時間序列長度；M為嵌入維數(shù)。經(jīng)過大量試驗證明，M取N/3左右時，分析效果比較理想。式(1)的變量間協(xié)方差是原序列x不同滯后的自協(xié)方差。構(gòu)造滯后自協(xié)方差陣Tx

(2)

Tx是Toeplitz矩陣，其主對角線元素是時間序列x的方差(或時遲為0的自協(xié)方差)。C(j)為時間序列x遲后為j的自協(xié)方差，0≤j≤M-1，用Yule-Walke估計法得到

(3)

然后根據(jù)公式

TxEk=λkEk,k=1,2,3,…,M

(4)

求得Tx的特征值λk和特征向量Ek。Ek就是M個分量構(gòu)成的一個時間序列，它反映時間序列x中的時間演變型。定義狀態(tài)向量在第M個特征向量上的投影

(5)

可以由其中一部分特征向量和時間系數(shù)來重建x的成分

(6)

與主成分分析類似，將Tx的特征值λk從大到小排列，λ1≥λ2≥λ3≥…≥λM≥0。截取前p個較大的特征值，由其所對應(yīng)的xk之和重建出可充分反映原序列的整體特征，即

(7)

由于電離層變化受太陽活動、地球自轉(zhuǎn)和公轉(zhuǎn)、地磁條件等多種因素的影響，因此電離層變化具有很強(qiáng)的日周期性。而奇異譜分析能夠識別和強(qiáng)化周期信號，利用奇異譜分析的方法可以提取出電離層TEC時間序列中除去異常擾動以及觀測噪聲部分的日周期部分，作為主要成分TECmain。與時間序列ARMA等異常探測方法相比，奇異譜方法提取的TECmain中包含了電離層隨著地球公轉(zhuǎn)的季節(jié)變化、太陽27 d自轉(zhuǎn)周期引起的電離層9 d、13.5 d、27 d周期變化等的影響，排除了背景場時間段與異常探測時間段地球外界環(huán)境不同造成的干擾。選取震前一段時間內(nèi)不受太陽活動和地磁異常擾動的TEC數(shù)據(jù)計算背景噪聲ε，取兩倍ε作為上下限值[26-29]，由此得到電離層TEC的正常變化范圍

(8)

以up和low作為TEC序列變化的上下界，超過該范圍就認(rèn)為是異常值。

1.2基于衛(wèi)星穿刺點軌跡STEC的區(qū)域電離層異常探測

GNSS信號的電離層延遲與信號的頻率的平方有反比的關(guān)系，根據(jù)此關(guān)系利用GNSS雙頻信號可以計算得到信號傳播路徑上的總電子含量TEC[30-32]。利用GNSS觀測數(shù)據(jù)計算電離層TEC的方法主要有偽距觀測值法、相位觀測值法、相位平滑偽距觀測值法等，目前國際上通常采用載波相位平滑的方法求解電離層TEC。文獻(xiàn)[33—34]提出非組合PPP方法求解TEC可以合理利用IGS發(fā)布的衛(wèi)星軌道、鐘差產(chǎn)品，減少待估參數(shù)，提高電離層估計精度。但是該方法依賴于外部高精度的衛(wèi)星軌道、鐘差等信息[35]，計算結(jié)果的可靠性、穩(wěn)定性有待驗證。本文采用載波相位平滑的方法求解電離層TEC，在實際計算中將小于30min的弧段舍去，保證計算結(jié)果具有足夠的精度。載波相位平滑的方法求解電離層TEC具體表達(dá)式為

(9)

1.3基于CIT技術(shù)的三維電離層異常分析

電離層三維層析技術(shù)(CIT)是利用GNSS無線電波投影反演電離層電子密度空間分布的新技術(shù)。

GNSS信號傳播路徑上的總電子含量可表示為

TEC=∫lNe(r,t)ds

(10)

式中，TEC為總電子含量；Ne為時間t測站與衛(wèi)星之間向徑r處的電子密度；l為信號的傳播路徑。由于測站的稀疏性以及信號角度的有限性，CIT技術(shù)一般采用離散化的方法求解，利用合適的函數(shù)基來表示反演區(qū)域電子密度

TEC=Ax+ε

(11)

式中，A為構(gòu)造矩陣，元素A(i,j)表示第i條射線在第j個格網(wǎng)中的截距；x為所有反演區(qū)域格網(wǎng)電子密度組成的列向量；ε為觀測誤差向量。圖2為電離層空間離散化格網(wǎng)示意圖。在實際反演過程中，由于構(gòu)造矩陣A為巨型稀疏矩陣，一般為秩虧矩陣，無法對矩陣求逆。

圖1　衛(wèi)星穿刺點軌跡示意圖Fig.1　IPP trajectories and SIP trajectories

圖2　電離層空間離散化格網(wǎng)示意圖Fig.2　Schematic diagram of ionospheric space discretization grid

本文采用國際參考電離層(IRI2012)作為迭代初始值，采用ART迭代算法對反演區(qū)域進(jìn)行迭代重構(gòu)。ART迭代算法的迭代公式為

(12)

式中，ai為矩陣構(gòu)造矩陣A的第i行；λ為迭代松弛因子，λ∈(0,1)，對于含有誤差的觀測數(shù)據(jù)，選擇合理的松弛因子至關(guān)重要[36]。

2分析與討論

2.1TEC時間序列的奇異譜分析

選取震前30 d內(nèi)不受太陽活動和地磁異常擾動的TEC數(shù)據(jù)計算背景噪聲ε，由此得到電離層TEC變化的上下界范圍。本文采用GIM提供的全球電離層格網(wǎng)數(shù)據(jù)進(jìn)行電離層異常探測。圖3依次給出了2015年3月26日至4月27日33 d的F10.7指數(shù)、Dst指數(shù)、Kp指數(shù)變化圖，從圖3中可以看出，3月26日至4月10日15天內(nèi)，F(xiàn)10.7均超過100但是小于150，太陽輻射較為強(qiáng)烈。Dst指數(shù)均大于-40，Kp指數(shù)除了4月3日短時間內(nèi)超過了4.0，其余均小于4.0，說明該段時間太陽活動較少，地磁場平靜，屬于電離層平靜時期，可以用于計算背景場噪聲。

圖3　2015年3月26日至4月27日F10.7、Dst、Kp指數(shù)變化Fig.3　The fluctuation of F10.7, Dst, Kp indices from 26, March to 27, April

利用2015年3月26日至4月10日的TEC數(shù)據(jù)計算得到電離層背景噪聲后，利用奇異譜方法對震中周圍電離層格網(wǎng)TEC序列進(jìn)行異常探測。發(fā)現(xiàn)4月23日在震中東南部，緯度約22.5°N—32.5°N，經(jīng)度約85°E—110°E范圍內(nèi)出現(xiàn)電離層正異常。圖4給出了電離層異常分布圖。圖中紅五角星為震中位置，紅底橢圓形區(qū)域為電離層異常區(qū)域，藍(lán)色正三角形為異常區(qū)域和非異常區(qū)域格網(wǎng)點。圖5和圖6只給出了2個非異常區(qū)域格網(wǎng)點和3個異常區(qū)域格網(wǎng)點的TEC時間序列圖。

圖5中給出了兩個非異常區(qū)域格網(wǎng)點和3個異常區(qū)域格網(wǎng)點2015年4月18日至4月27日的TEC時間序列圖。圖5(a)、(b)兩圖分別對應(yīng)圖4中a、b兩格網(wǎng)點。紅色曲線為實際TEC值，綠色為上界曲線，藍(lán)色為下界曲線，紅色豎線為地震發(fā)生時間。從圖中可以看出，這兩個格網(wǎng)點TEC時間序列均處于上下界范圍內(nèi)，沒有出現(xiàn)異常。圖5(c)、(d)、(e)分別對應(yīng)圖4中c、d、e格網(wǎng)點，紅色箭頭指出了異常出現(xiàn)的時間。從圖5可以看出，2015年4月23日前后幾天c、d、e格網(wǎng)點實際TEC序列均處于正常的上下界范圍之內(nèi)，而在4月23日當(dāng)天，實際TEC序列超出了上界約15個TECu，e點超出約10個TECu，說明當(dāng)天出現(xiàn)了明顯的正異常。

從圖3可以看出，在2015年4月23日，太陽活動指數(shù)F10.7低于150，且與前后幾天持平，沒有明顯變化；Dst指數(shù)接近0值，Kp指數(shù)小于3.5，基本可以排除太陽活動、地磁擾動引起電離層異常的可能性，表明此次電離層異常與4月25日發(fā)生的尼泊爾地震有較大的相關(guān)性。

2.2二維電離層地圖異常探測

為了進(jìn)一步驗證上文電離層異常分析結(jié)果的可靠性，對緯度15°N—40°N，經(jīng)度70°N—115°N范圍內(nèi)GIM電離層地圖進(jìn)行異常探測。同樣利用2015年3月26日至4月10日電離層地圖作為背景值，2倍標(biāo)準(zhǔn)差作為上下界，對4月23日電離層地圖進(jìn)行異常探測。如果TEC值與背景值之差絕對值小于2ε，表示沒有明顯異常，如果大于2ε，ΔTEC>0表示有正異常，ΔTEC<0表示有負(fù)異常。圖6給出了2015年4月23日電離層異常的探測結(jié)果。從圖中可以看出，UT06:00時電離層沒有出現(xiàn)異?，F(xiàn)象，電離層非常平靜；UT09:00時在(25°N，110°E)開始出現(xiàn)電離層異常。隨著時間推移，異常區(qū)域開始向西移動，異常峰值逐漸增大。在UT12:00時異常峰值達(dá)到最大，ΔTEC≥20TECu，異常范圍約為25°N—37.5°N，90°E—110°E，異常區(qū)域整體呈橢圓形。隨后異常減小并擴(kuò)散，UT15:00時異常峰值減小為15TECu左右，UT18:00后異常逐漸消失，UT21:00時電離層完全恢復(fù)平靜。異常整體持續(xù)時間為9h。文獻(xiàn)[37—41]均發(fā)現(xiàn)了在大型地震前有大范圍的、與震中距離較遠(yuǎn)的電離層異常，說明此次探測到的電離層異常是完全有可能的。利用二維電離層地圖的異常探測結(jié)果與上文中奇異譜時間序列探測結(jié)果吻合較好，進(jìn)一步確認(rèn)了異常的存在性。

圖4　電離層異常區(qū)域示意圖Fig.4　Schematic diagram of ionospheric anomalise region

圖5　電離層異常和非異常區(qū)域格網(wǎng)點TEC時間序列圖Fig.5　TEC time series from 18 to 28, April of detection grid points in ionospheric abnormal and normal area

圖6　電離層地圖異常分布Fig.6　Ionospheric abnormal area disribution on ionosphere map

2.3衛(wèi)星穿刺點軌跡TEC異常分析

本文利用中國大陸構(gòu)造環(huán)境監(jiān)測網(wǎng)絡(luò)數(shù)據(jù)，選取了靠近震區(qū)的3個測站(XZAR、XZCD、XZCY)的觀測數(shù)據(jù)對上文中電離層異常區(qū)域進(jìn)行研究。圖7給出了6條穿刺點軌跡的位置以及3個陸態(tài)網(wǎng)測站的分布圖。圖中紅色五角星為震中位置，藍(lán)色圓點為陸態(tài)網(wǎng)測站，帶箭頭曲線為穿刺點軌跡，箭頭方向指示穿刺點時間順序，不同曲線代表不同測站或不同衛(wèi)星形成的穿刺點軌跡，灰色橢圓為電離層異常區(qū)域。從圖7中可以看出，6條穿刺點軌跡均穿過了上文中探測到的異常區(qū)域，能夠利用這些曲線對該區(qū)域電離層狀況進(jìn)行研究。

圖7　測站和穿刺點軌跡分布Fig.7　The distribution of crustal movement observation network stations and IPP trajectories

圖8給出了圖7中6條穿刺點軌跡2015年4月21日至27日STEC的對比圖。圖8中紅色曲線均為4月23日STEC變化曲線，黑色為其余

天數(shù)的STEC。以圖7中深紅色曲線XZAR-SAT1和圖8中對應(yīng)的XZAR-SAT1子圖為例，隨著穿刺點進(jìn)入電離層異常區(qū)域，23日STEC逐漸高于其他天穿刺點STEC，經(jīng)過異常區(qū)域中心時，差值達(dá)到最大。然后隨著穿刺點遠(yuǎn)離異常峰值中心并離開異常區(qū)域，差值逐漸減小直至消失。其余穿刺點軌跡均也出現(xiàn)類似情況。6條穿刺點軌跡23日STEC值均高于其他日期10～20個TECu左右，且時間為UT9:00—13:00左右，與上文電離層地圖探測到的異常出現(xiàn)時間十分吻合。

2.4基于CIT的4月23日電離層異常時空分布分析

本文采用了中國大陸構(gòu)造環(huán)境監(jiān)測網(wǎng)絡(luò)青海、西藏、四川、云南地區(qū)共73個測站數(shù)據(jù)，表1給出了CIT反演范圍，圖9給出了青藏川滇地區(qū)測站的分布圖。利用2015年3月26日至4月9日15d的每小時的反演結(jié)果作為背景場，對4月23日電離層時空分布情況進(jìn)行異常探測。

表1　CIT反演范圍

圖8　4月21日至4月27日穿刺點軌跡STEC對比圖Fig.8　STEC of IPP trajectories from 21 to 27 April

圖9　中國陸態(tài)網(wǎng)青藏川滇地區(qū)站點分布圖Fig.9　The distribution of crustal movement observation network stations in Tibet,Qinghai,Sichuan and Yunnan

圖10給出了2015年4月23日UT03：00、UT06：00、UT09：00、UT12：00、UT15:00、UT18:00，100～500 km高度面的探測結(jié)果。從圖10中可以看出， UT3:00時，電離層各個高度面都非常平靜，沒有出現(xiàn)異常現(xiàn)象；在UT6:00時，在100 km高度面出現(xiàn)少量異常，但是比較分散，而到了UT9:00時，在30°N，105°E附近開始出現(xiàn)明顯的電離層異常。其中100～300 km隨著高度的增加，電離層異常峰值不斷增大，異常峰值最大約為15 TECu左右；300～500 km，隨著高度增加峰值逐漸減小，符合電離層密度垂直分布規(guī)律；UT12:00時異常范圍擴(kuò)大并稍微向西移動，峰值減小。UT15:00時，峰值進(jìn)一步減小，異常開始消散；UT18:00時異常消失，電離層恢復(fù)平靜。電離層三維層析結(jié)果與2.1節(jié)奇異譜分析結(jié)果、2.2節(jié)電離層地圖異常分布以及2.3節(jié)衛(wèi)星穿刺點軌跡異常探測結(jié)果的區(qū)域位置和時間都十分吻合，并且更加清晰、立體地展現(xiàn)出電離層異常的空間水平、垂直分布，以及隨時間的變化情況。

圖10　2015年4月23日15°N—40°N，75°E—115°E區(qū)域不同時間段、不同高度面的異常分布Fig.10　The ionospheric anomaly distribution in the inversion range 15°N—40°N，75°E—115°E in 23 April,2015 of different time and height

3結(jié)論

本文提出利用奇異譜方法對震前TEC時間序列進(jìn)行異常探測。利用GIM提供的全球電離層地圖數(shù)據(jù)，對震中附近格網(wǎng)點探測結(jié)果表明，2015年4月23日在震中正東方向出現(xiàn)大范圍電離層正異常。同時利用電離層地圖對電離層異常區(qū)域進(jìn)行探測，確定了異常出現(xiàn)的時間為2015年4月23日UT9:00—15:00，主要分布在25°N—37.5°N，90°E—110°E，異常區(qū)域整體呈橢圓形。為了進(jìn)一步驗證該結(jié)論的正確性，本文利用中國大陸構(gòu)造環(huán)境監(jiān)測網(wǎng)絡(luò)西藏地區(qū)測站觀測數(shù)據(jù)計算4月21日至4月27日異常區(qū)域穿刺點軌跡STEC隨時間變化情況，同樣在上述區(qū)域UT9:00—15:00出現(xiàn)了電離層正異常。

為了更加精細(xì)、立體地研究電離層的時空分布，本文利用CIT技術(shù)對15°N—40°N，75°E—115°E，100～500 km范圍內(nèi)電離層進(jìn)行反演，以2015年3月26日至4月9日反演結(jié)果作為背景場，對4月23日反演結(jié)果進(jìn)行異常探測，探測結(jié)果與上文中奇異譜分析以及衛(wèi)星穿刺點軌跡STEC方法的探測結(jié)果基本吻合。三維層析結(jié)果進(jìn)一步反映出電離層異常隨高度的變化情況：電離層異常出現(xiàn)在4月23日UT9:00—15:00，峰值位置約為30°N，105°E，主要分布在200～400 km高度范圍內(nèi)；100～300 km隨著高度增加，異常峰值逐漸增大，異常范圍擴(kuò)大；300～500 km，隨著高度增加，異常峰值逐漸減小，異常范圍縮小，異常隨高度變化規(guī)律與電離層電子密度垂直分布規(guī)律相一致。CIT反演結(jié)果從三維角度提供了電離層異常的空間分布情況，為分析此次電離層異常與地震之間的關(guān)系提供了進(jìn)一步的參考。

參考文獻(xiàn)：

[1]祝芙英,吳云,林劍,等. 汶川地震前電離層VTEC的異常響應(yīng)[J]. 地震學(xué)報, 2009, 31(2): 180-187.

ZHU Fuying, WU Yun, LIN Jian, et al. Anomalous Response of Ionospheric VTEC before the Wenchuan Earthquake[J]. Acta Seismologica Sinica, 2009, 31(2): 180-187.

[2]LEONARD R S, BARNES JR R A. Observation of Ionospheric Disturbances Following the Alaska Earthquake[J]. Journal of Geophysical Research, 1965, 70(5): 1250-1253.

[3]ANTSELEVICH M G. The Influence of Tashkent Earthquake on the Earth’s Magnetic Field and the Ionosphere[C]∥Tashkent Earthquake 26 April 1966. Tashkent: FAN, 1971: 187-188.

[4]WEAVER P F, YUEN P C, PROLSS G W. Acoustic Coupling into the Ionosphere from Seismic Waves of the Earthquake at Kurile Islands on August 11, 1969[J]. Nature, 1970, 226(5252): 1239-1241.

[5]PULINETS S A. Seismic Activity as a Source of the Ionospheric Variability[J]. Advances in Space Research, 1998, 22(6): 903-906.

[6]PULINETS S A. Natural Radioactivity, Earthquakes, and the Ionosphere[J]. EOS Transactions American Geophysical Union, 2007, 88(20): 217-218.

[7]LIU J Y, CHEN Y I, CHUO Y J, et al. Variations of Ionospheric Total Electron Content during the Chi Chi Earthquake[J]. Geophysical Research Letters, 2001, 28(7): 1383-1386.

[8]丁鑒海, 索玉成, 余素榮. 地磁場與電離層異?，F(xiàn)象及其與地震的關(guān)系[J]. 空間科學(xué)學(xué)報, 2005, 25(6): 536-542.

DING Jianhai, SUO Yucheng, YU Surong. Phenomena of Geomagnetic and Ionospheric Anomalies and Their Relation to Earthquakes[J]. Chinese Journal of Space Science, 2005, 25(6): 536-542.

[9]CALAIS E, MINSTER J B. GPS Detection of Ionospheric Perturbations Following the January 17, 1994, Northridge Earthquake[J]. Geophysical Research Letters, 1995, 22(9): 1045-1048.

[10]馬新欣, 林湛, 陳化然, 等. 基于GPS和COSMIC數(shù)據(jù)分析汶川地震TEC和NmF2擾動[J]. 地球物理學(xué)報, 2014, 57(8): 2415-2422.

MA Xinxin, LIN Zhan, CHEN Huaran, et al. Analysis on Ionospheric Perturbation of TEC andNmF2Based on GPS and COSMIC Data before and after the Wenchuan Earthquake[J]. Chinese Journal of Geophysics, 2014, 57(8): 2415-2422.

[11]AKHOONDZADEH M, PARROT M, SARADJIAN M R. Electron and Ion Density Variations before Strong Earthquakes (M>6.0) Using DEMETER and GPS Data[J]. Natural Hazards and Earth System Sciences, 2010, 10(1): 7-18.

[12]祝芙英, 吳云, 林劍, 等. 汶川Ms8.0地震前電離層TEC異常分析[J]. 大地測量與地球動力學(xué), 2008, 28(6): 16-21.

ZHU Fuying, WU Yun, LIN Jian, et al. Study on Ionospheric TEC Anomaly Prior to WenchuanMs8.0 Earthquake[J]. Journal of Geodesy and Geodynamics, 2008, 28(6): 16-21.

[13]姚宜斌, 陳鵬, 吳寒, 等. 2011年3月11日日本地震震前電離層異常變化分析[J]. 科學(xué)通報, 2012, 57(5): 355-365.

YAO Yibin, CHEN Peng, WU Han, et al. Analysis of Ionospheric Anomalies before the 2011Mw9.0 Japan Earthquake[J]. Chinese Science Bulletin, 2012, 57(5): 500-510.

[14]CAHYADI M N, HEKI K. Ionospheric Disturbances of the 2007 Bengkulu and the 2005 Nias Earthquakes, Sumatra, Observed with a Regional GPS Network[J]. Journal of Geophysical Research, 2013, 118(4): 1777-1787.

[15]HEKI K. Ionospheric Electron Enhancement Preceding the 2011 Tohoku Oki Earthquake[J]. Geophysical Research Letters, 2011, 38(17): L17312.

[16]LE H, LIU J Y, LIU L. A Statistical Analysis of Ionospheric Anomalies before 736M6.0+ Earthquakes During 2002—2010[J]. Journal of Geophysical Research, 2011, 116(A2): A02303.

[17]聞德保. 基于GPS的電離層層析算法及其應(yīng)用研究[D]. 武漢: 中國科學(xué)院測量與地球物理研究所, 2007.

WEN Debao. Investigation of GPS-based Ionospheric Tomographic Algorithms and Their Applications[D]. Wuhan: Institute of Geodesy and Geophysics Chinese Academy of Sciences, 2007.

[18]楊劍, 吳云, 周義炎. 基于電離層層析成像技術(shù)探測汶川地震前電離層異常[J]. 大地測量與地球動力學(xué), 2011, 31(1): 9-14.

YANG Jian, WU Yun, ZHOU Yiyan. Probe into Seismo-ionospheric Anomaly of WenchuanMs8.0 Earthquake Based on Computerized Ionospheric Tomography[J]. Journal of Geodesy and Geodynamics, 2011, 31(1): 9-14.

[19]鄒玉華. GPS地面臺網(wǎng)和掩星觀測結(jié)合的時變?nèi)S電離層層析[D]. 武漢: 武漢大學(xué), 2004.

ZOU Yuhua. A Study of Time-dependent 3-D Ionospheric Tomography with Ground-based GPS Network and Occultation Observations[D]. Wuhan: Wuhan University, 2004.

[20]孟泱, 安家春, 王澤民, 等. 基于GPS的南極電離層電子總含量空間分布特征研究[J]. 測繪學(xué)報, 2011, 40(1): 37-40.

MENG Yang, AN Jiachun, WANG Zemin, et al. Spatial Distribution of Antarctic Ionosphere TEC Based on GPS[J]. Acta Geodaetica et Cartographica Sinica, 2011, 40(1): 37-40.

[21]RIUS A, RUFFINI G, CUCURULL L. Improving the Vertical Resolution of Ionospheric Tomography with GPS Occultations[J]. Geophysical Research Letters, 1997, 24(18): 2291-2294.

[22]THAMPI S V, LIN C, LIU H, et al. First Tomographic Observations of the Midlatitude Summer Nighttime Anomaly over Japan[J]. Journal of Geophysical Research, 2009, 114(A10): A10318.

[23]KUNITSYN V E, NESTEROV I A, PADOKHIN A M, et al. Ionospheric Radio Tomography Based on the GPS/GLONASS Navigation Systems[J]. Journal of Communications Technology and Electronics, 2011, 56(11): 1269-1281.

[24]SEEMALA G K, YAMAMOTO M, SAITO A, et al. Three-dimensional GPS Ionospheric Tomography over Japan Using Constrained Least Squares[J]. Journal of Geophysical Research, 2014, 119(4): 3044-3052.

[25]王解先, 連麗珍, 沈云中. 奇異譜分析在GPS站坐標(biāo)監(jiān)測序列分析中的應(yīng)用[J]. 同濟(jì)大學(xué)學(xué)報(自然科學(xué)版), 2013, 41(2): 282-288.

WANG Jiexian, LIAN Lizhen, SHEN Yunzhong. Application of Singular Spectral Analysis to GPS Station Coordinate Monitoring Series[J]. Journal of Tongji University (Natural Science), 2013, 41(2): 282-288.

[26]杜品仁, 蔣和榮, 郭兼善. 電離層異常能否作為地震前兆的研究[J]. 地震, 1998, 18(2): 119-126.

DU Pinren, JIANG Herong, GUO Jianshan. Research on Possibility of Ionospheric Anomalies as an Earthquake Precursor[J]. Earthquake, 1998, 18(2): 119-126.

[27]LIU J Y, CHUO Y J, SHAN S J, et al. Pre-earthquake Ionospheric Anomalies Registered by Continuous GPS TEC Measurements[J]. Annales Geophysicae, 2004, 22(5): 1585-1593.

[28]OYAMA K I, KAKINAMI Y, LIU J Y, et al. Reduction of Electron Temperature in Low-latitude Ionosphere at 600 km before and after Large Earthquakes[J]. Journal of Geophysical Research, 2008, 113(A11): A11317.

[29]KIM V P, HEGAI V V. A Possible Presage of Strong Earthquakes in the Night-time Mid-latitudeF2Region Ionosphere[C]∥Atmospheric and Ionospheric Electromagnetic Phenomena Associated with Earthquakes. Tokyo: Terra Scientific Publishing Company, 1999: 619-627.

[30]李征航,黃勁松.GPS測量與數(shù)據(jù)處理[M].武漢: 武漢大學(xué)出版社, 2005: 102-105.

LI Zhenghang, HUANG Jinsong. GPS Surveying and Data Processing[M]. Wuhan: Wuhan University Press, 2005: 102-105.

[31]劉經(jīng)南, 陳俊勇, 張燕平, 等. 廣域差分GPS原理和方法[M]. 北京: 測繪出版社, 1999.

LIU Jingnan, CHEN Junyong, ZHANG Yanping, et al. Theory and Technique of Wide Area Differential GPS[M]. Beijing: Surveying and Mapping Press, 1999.

[32]袁運斌. 基于GPS的電離層監(jiān)測及延遲改正理論與方法的研究[D]. 武漢: 中國科學(xué)院測量與地球物理研究所, 2002.

YUAN Yunbin. Study on Theories and Methods of Correcting Ionospheric Delay and Monitoring Ionosphere Based on GPS[D]. Wuhan: Institute of Geodesy and Geophysics of Chinese Academy of Sciences, 2002.

[33]張寶成, 歐吉坤, 袁運斌, 等. 利用非組合精密單點定位技術(shù)確定斜向電離層總電子含量和站星差分碼偏差[J]. 測繪學(xué)報, 2011, 40(4): 447-453.

ZHANG Baocheng, OU Jikun, YUAN Yunbin, et al. Calibration of Slant Total Electron Content and Satellite-receiver’s Differential Code Biases with Uncombined Precise Point Positioning Technique[J]. Acta Geodaetica et Cartographica Sinica, 2011, 40(4): 447-453.

[34]張寶成, 歐吉坤, 李子申, 等. 利用精密單點定位求解電離層延遲[J]. 地球物理學(xué)報, 2011, 54(4): 950-957.

ZHANG Baocheng, OU Jikun, LI Zishen, et al. Determination of Ionospheric Observables with Precise Point Positioning[J]. Chinese Journal of Geophysics, 2011, 54(4): 950-957.

[35]李子申. GNSS/Compass電離層時延修正及TEC監(jiān)測理論與方法研究[D]. 武漢: 中國科學(xué)院測量與地球物理研究所, 2012.

LI Zishen. Study on the Mitigation of Ionospheric Delay and the Monitoring of Global Ionospheric TEC Based on GNSS/Compass[D]. Wuhan: Institute of Geodesy and Geophysics Chinese Academy of Sciences, 2012.

[36]PRYSE S E, KERSLEY L, MITCHELL C N, et al. A Comparison of Reconstruction Techniques Used in Ionospheric Tomography[J]. Radio Science, 1998, 33(6): 1767-1779.

[37]CHOI B K, PARK J U, LEE S J. Ionospheric Anomalies Observed over South Korea Preceding the Great Tohoku Earthquake of 2011[J]. Advances in Space Research, 2012, 50(3): 311-317.

[38]HAO Y Q, XIAO Z, ZHANG D H. Teleseismic Magnetic Effects (TMDs) of 2011 Tohoku Earthquake[J]. Journal of Geophysical Research, 2013, 118(6): 3914-3923.

[39]張學(xué)民, 劉靜, 申旭輝, 等. 2005年3月28日蘇門答臘8.6級地震前的電離層擾動[J]. 地球物理學(xué)報, 2010, 53(3): 567-575.

ZHANG Xuemin, LIU Jing, SHEN Xuhui, et al. Ionospheric Perturbations Associated with theM8.6 Sumatra Earthquake on 28 March 2005[J]. Chinese Journal of Geophysics, 2010, 53(3): 567-575.

[40]KLIMENKO M V, KLIMENKO V V, ZAKHARENKOVA I E, et al. Formation Mechanism of Great Positive TEC Disturbances Prior to Wenchuan Earthquake on May 12, 2008[J]. Advances in Space Research, 2011, 48(3): 488-499.

[41]HEKI K, ENOMOTO Y. Preseismic Ionospheric Electron Enhancements Revisited[J]. Journal of Geophysical Research, 2013, 118(10): 6618-6626.

(責(zé)任編輯：叢樹平)

修回日期： 2015-12-16

First author： YAO Yibin(1976—), male,professor, majors in geodetic data processing,GNSS space environment science.

E-mail： ybyao@sgg.whu.edu.cn

E-mail： czzhai@whu.edu.cn

The Pre-earthquake Ionosphere Anomaly of the 2015 Nepal Earthquake

YAO Yibin1,2,3，ZHAI Changzhi1，KONG Jian4，LIU Lei1

1. School of Geodesy and Geomatics, Wuhan University, Wuhan 430079, China; 2. Key Laboratory of Geospace Environment and Geodesy, Ministry of Education, Wuhan University, Wuhan 430079, China; 3. Collaborative Innovation Center for Geospatial Technology, Wuhan 430079, China; 4. Chinese Antarctic Center of Surveying and Mapping, Wuhan 430079, China

Abstract：An ionospheric anomaly detection method based on singular spectrum analysis is introduced firstly.With the TEC time sequence of GIM grid points around the Nepal earthquake epicenter, the positive ionospheric anomalies are determined on 23 April around the eastern region of the epicenter. Using the two-dimensional ionospheric anomalies map, further analysis is taken to comfirm the ionospheric positive anomaly area distributing from 25°N to 37.5° N, 90°E to 110°E, and the anomaly time ranging UT9:00 to UT15:00 on 23 April, 2015. The data of crustal movement observation network of China issued to analyze the STEC fluctuation changes of IPP(ionosphere piece piont) trajectories in the abnormal area, and it is found that on 23 April the STEC values of IPP trajectories passing through the abnormal area are obviously larger than the ones a few days after and before the earthquake time, and returned to the normal level after leaving the abnormal area. Lastly, the CIT (computerized ionosphere tomography) method is used to present the 3D ionospheric anomalies distribution, and it is found that, at UT9: 00—UT15:00 on 23 April the peak of ionospheric positive anomaly on the eastern region of the epicenter is approximately at (30°N, 115°E), with the altitude ranging from 100 to 500 km. Also the anomalous peak varies with height is consistent with ionosphere density vertical distribution.

Key words：seismic-ionospheric disturbance; singular spectrum analysis; STEC; CIT

Corresponding author：ZHAI Changzhi

通信作者：翟長治

第一作者簡介：姚宜斌(1976—)，男，教授，研究方向為測量數(shù)據(jù)處理理論與方法、GNSS空間環(huán)境學(xué)。

收稿日期：2015-07-20

基金項目：國家自然科學(xué)基金(41274022;41574028);湖北省杰出青年科學(xué)基金(2015CFA036)

中圖分類號：P228

文獻(xiàn)標(biāo)識碼：A

文章編號：1001-1595(2016)04-0385-11

Foundation support： The National Natural Science Foundation of China(Nos.41274022; 41574028); Natural Science Foundation for Distinguished Young Scholars of Hubei Province of China(No.2015CFA036)

引文格式：姚宜斌，翟長治，孔建,等.2015年尼泊爾地震的震前電離層異常探測[J].測繪學(xué)報，2016,45(4)：385-395. DOI:10.11947/j.AGCS.2016.20150384.

YAO Yibin，ZHAI Changzhi，KONG Jian,et al.The Pre-earthquake Ionosphere Anomaly of the 2015 Nepal Earthquake[J]. Acta Geodaetica et Cartographica Sinica,2016,45(4):385-395. DOI:10.11947/j.AGCS.2016.20150384.