• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    地震電離層探測(cè)技術(shù)及其應(yīng)用研究進(jìn)展*

    2016-11-28 10:56:15張學(xué)民申旭輝趙庶凡歐陽新艷婁文宇澤仁志瑪何建輝
    地震學(xué)報(bào) 2016年3期

    張學(xué)民 申旭輝 趙庶凡 劉 靜 歐陽新艷 婁文宇 澤仁志瑪 何建輝 錢 庚

    1) 中國北京100036中國地震局地震預(yù)測(cè)研究所2) 中國北京100085中國地震局地殼應(yīng)力研究所

    ?

    地震電離層探測(cè)技術(shù)及其應(yīng)用研究進(jìn)展*

    1) 中國北京100036中國地震局地震預(yù)測(cè)研究所2) 中國北京100085中國地震局地殼應(yīng)力研究所

    2017年即將發(fā)射的中國電磁監(jiān)測(cè)試驗(yàn)衛(wèi)星將填補(bǔ)地震電離層立體監(jiān)測(cè)體系中不可或缺的空白區(qū)域, 也將為天地一體化地震電磁對(duì)比校驗(yàn)及圈層耦合機(jī)理認(rèn)識(shí)提供重要的科學(xué)支撐. 針對(duì)近期地震電離層立體監(jiān)測(cè)體系發(fā)展的需求, 本文主要介紹了目前國內(nèi)用于地震研究的地基及空基電離層探測(cè)技術(shù), 包括電離層垂測(cè)/斜測(cè)、 甚低頻(VLF)電波觀測(cè)、 舒曼諧振觀測(cè)、 GPS及空間衛(wèi)星電磁等, 并總結(jié)了各種探測(cè)技術(shù)在國內(nèi)外地震應(yīng)用研究中的進(jìn)展; 最后結(jié)合不同探測(cè)手段的優(yōu)勢(shì), 探討了地震電磁立體探測(cè)系統(tǒng)的構(gòu)建, 并就未來的多手段綜合應(yīng)用發(fā)展提出了建議.

    電離層垂測(cè)/斜測(cè) VLF電波 舒曼諧振 GPS TEC 電磁衛(wèi)星

    引言

    電離層作為日地空間的重要組成部分, 一直是與人類社會(huì)發(fā)展和生活息息相關(guān)的重要環(huán)境區(qū)域. 電離層探測(cè)主要針對(duì)總電子含量、 電子密度、 電子溫度、 離子組成成分、 離子密度、 離子溫度、 離子漂移運(yùn)動(dòng)速度、 電場(chǎng)強(qiáng)度、 磁場(chǎng)強(qiáng)度、 不均勻結(jié)構(gòu)的尺度、 高度分布等特征參量及其變化進(jìn)行觀測(cè), 分為直接探測(cè)和間接探測(cè)兩種. 直接探測(cè)是利用火箭、 衛(wèi)星等空間飛行器, 將探測(cè)裝置攜帶至電離層中, 通過探測(cè)電離層等離子體或環(huán)境對(duì)裝置的直接作用獲得電離層特性參量; 間接探測(cè)則是依據(jù)天然輻射或人工發(fā)射機(jī)發(fā)射的電磁波通過電離層傳播時(shí)與等離子體相互作用所產(chǎn)生的電磁效應(yīng)或傳播特征, 推算電離層特性參量. 隨著科學(xué)技術(shù)的進(jìn)步, 電離層探測(cè)水平不斷提高, 特別是人造地球衛(wèi)星的應(yīng)用使得從地基觀測(cè)向天基觀測(cè)發(fā)展成為可能, 并取得了許多關(guān)于電離層結(jié)構(gòu)及其變化規(guī)律的研究成果(梁百先等, 1994; 肖佐, 1997), 其中大多數(shù)與太陽活動(dòng)對(duì)電離層的影響相關(guān). 近年來, 隨著地震預(yù)測(cè)的需求日益迫切, 地基電離層探測(cè)技術(shù)被廣泛應(yīng)用于地震前兆及同震效應(yīng)研究; 與此同時(shí), 從事空間科學(xué)研究的人員也開始關(guān)注到地震電離層效應(yīng). 值得一提的是, 2004年法國科學(xué)家在綜合世界各國衛(wèi)星觀測(cè)資料的基礎(chǔ)上, 成功發(fā)射了DEMETER (detection of electric-magnetic emissions transmitted from earthquake regions)衛(wèi)星, 專門用于地震電離層效應(yīng)的探測(cè), 開創(chuàng)了天基地震電磁觀測(cè)的新局面. 天基觀測(cè)的空間覆蓋率遠(yuǎn)遠(yuǎn)高于地基觀測(cè), 基于全球觀測(cè)的優(yōu)勢(shì)使其可以獲得更多的震例積累, 該方面地基觀測(cè)則難以企及. DEMETER衛(wèi)星運(yùn)行期間發(fā)生了2004年12月印尼MS9.0, 2008年中國汶川MS8.0, 2010年海地MS7.2和智利MS8.8等多次災(zāi)難性大地震, 該衛(wèi)星獲取了大量的震例觀測(cè)資料, 受到地球科學(xué)界的廣泛關(guān)注.

    我國在借鑒DEMETER衛(wèi)星科學(xué)載荷配置的基礎(chǔ)上, 提出了發(fā)射地震電磁監(jiān)測(cè)衛(wèi)星的計(jì)劃, 這將是繼法國之后世界上第二個(gè)以服務(wù)地震預(yù)測(cè)為特點(diǎn)的專項(xiàng)衛(wèi)星計(jì)劃, 其亮點(diǎn)在于天基觀測(cè)與地基觀測(cè)的相互結(jié)合.

    本文將圍繞目前國內(nèi)主要的幾種地基電離層探測(cè)技術(shù), 對(duì)其探測(cè)原理和應(yīng)用成果進(jìn)行總結(jié); 然后結(jié)合衛(wèi)星探測(cè)技術(shù)及其應(yīng)用進(jìn)展, 分析地基觀測(cè)和天基觀測(cè)各自的優(yōu)勢(shì)和缺陷; 最后闡述我國地震電磁衛(wèi)星發(fā)展計(jì)劃的新亮點(diǎn)----地震電磁效應(yīng)立體觀測(cè)體系, 并就如何綜合發(fā)揮各自效能提出相關(guān)建議.

    1 地基電離層探測(cè)技術(shù)

    地基電離層探測(cè)基本都屬于間接探測(cè)的范疇, 如目前常用的電離層測(cè)高、 GPS 總電子含量(total electric content, 簡寫為TEC)、 甚低頻(very low frequency, 簡寫為VLF)電波等均是利用人工發(fā)射源所發(fā)射的電磁波信號(hào)的電離層效應(yīng)來反映或反演電離層特征參量, 而舒曼諧振觀測(cè)則是利用閃電源所激發(fā)的地球電離層波導(dǎo)系統(tǒng)內(nèi)電磁波的諧振效應(yīng)來推算諧振腔的尺度變化. 下面將對(duì)上述幾種地基電離層探測(cè)技術(shù)逐一介紹.

    1.1 電離層垂測(cè)/斜測(cè)

    電離層垂測(cè)/斜測(cè)是利用電離層測(cè)高儀所發(fā)射的高頻無線電波在地面接收其反射回波實(shí)現(xiàn)對(duì)電離層進(jìn)行日常觀測(cè)的一種技術(shù). 垂測(cè)是垂直向上發(fā)射頻率隨時(shí)間變化的無線電脈沖, 在同一地點(diǎn)接收這些脈沖的電離層反射信號(hào); 斜測(cè)則是將垂直探測(cè)方法中的發(fā)射和接收設(shè)備分別置于地面上相隔一定距離的兩個(gè)點(diǎn)位, 使用高精度GPS時(shí)間標(biāo)定實(shí)現(xiàn)收、 發(fā)同步, 然后測(cè)量電離層反射回波時(shí)延隨頻率的變化. 斜測(cè)可以實(shí)現(xiàn)較大范圍內(nèi)的電離層探測(cè)和研究, 通過測(cè)量電波往返的傳遞時(shí)延, 獲得反射高度與頻率的關(guān)系曲線, 即頻高圖.

    電離層測(cè)高的基本原理是發(fā)射不同頻率(1—30 MHz)的高頻信號(hào), 對(duì)回波信號(hào)進(jìn)行分析計(jì)算以獲取電離層特征參數(shù). 電離層E, F1, F2, Es等層的臨界頻率、 虛高等參數(shù)可以從頻高圖或電離圖(中國科學(xué)院地質(zhì)與地球物理研究所, 2015)中直接讀出, 其中虛高h(yuǎn)′由h′=ct/2計(jì)算而得, 式中c為真空電波傳播速度,t為回波時(shí)間, 由于實(shí)際電波速度總是小于真空光速, 因此由該式計(jì)算出的h′不是反射點(diǎn)的真正高度h, 而是比h大, 通常稱h′為等效高度或虛高. 另一方面, 電離層測(cè)高也存在其局限性, 例如, 無法探測(cè) D層的電離狀態(tài), 難于獲得E層與F層之間的谷區(qū)(高度120—140 km)信息, 不能研究F層電子密度峰值以上的電離層, 這些問題可由其它探測(cè)手段予以解決.

    由于電離層測(cè)高站費(fèi)用較高, 全球僅有百余個(gè), 國內(nèi)十幾個(gè)測(cè)高站大多由中國電子科技集團(tuán)公司第二十二研究所建設(shè)管理. 2008年汶川地震后, 鑒于距離震中較近的重慶站和昆明站所觀測(cè)到的foF2大幅度擾動(dòng)(張學(xué)民等, 2009a), 為加強(qiáng)首都圈地震監(jiān)測(cè)工作, 中國地震局與中電集團(tuán)22所利用華北地區(qū)已有的5個(gè)測(cè)高站, 在河北、 山東、 山西、 遼寧、 江蘇、 陜西等多省市布設(shè)完成20個(gè)接收站, 形成100條鏈路的電離層斜測(cè)網(wǎng), 使大華北地區(qū)電離層監(jiān)測(cè)的空間分辨率達(dá)到100 km, 大大提高了電離層測(cè)高站的工作效率, 改善了單站監(jiān)測(cè)空間有限的狀況. 為進(jìn)一步加強(qiáng)對(duì)南北地震帶南部川滇地震危險(xiǎn)區(qū)的監(jiān)測(cè), 中國地震局地震預(yù)測(cè)研究所自2013年起在樂山、 普洱和大理先后布設(shè)了3個(gè)電離層測(cè)高站, 其中樂山站和普洱站的觀測(cè)已有兩年之久, 我們計(jì)劃在未來幾年內(nèi), 通過增設(shè)電離層斜測(cè)站進(jìn)一步提高川滇地區(qū)的空間監(jiān)測(cè)分辨率.

    利用電離層測(cè)高研究地震始于20世紀(jì). Pulinets和Boyarchuk(2004)最早綜合介紹了國際上在地基電離層和衛(wèi)星電離層測(cè)高中觀測(cè)到的強(qiáng)震前電離層擾動(dòng)現(xiàn)象, 并分析了特征參量foF2在強(qiáng)震前后的時(shí)空演化特征, 例如, 1964年3月27日阿拉斯加MS9.2大地震前30小時(shí)直至震后18小時(shí), Alouette-1衛(wèi)星和地基觀測(cè)均記錄到震中上空電離層F2層臨界頻率foF2異常持續(xù)增強(qiáng)的現(xiàn)象. 蔡軍濤等(2007)總結(jié)了國外文獻(xiàn)中多次地震前E層臨界頻率增強(qiáng)及Es擴(kuò)展現(xiàn)象. 我國關(guān)于電離層測(cè)高的地震應(yīng)用研究起步較晚. 丁鑒海等(2004)最先開展了地磁觀測(cè)資料與電離層測(cè)高數(shù)據(jù)的聯(lián)合分析研究, 在2001年昆侖山口MS8.1和1997年西藏瑪尼MS7.5(余素榮等, 2004)等地震前均觀測(cè)到臨近電離層臺(tái)站上空的foF2相對(duì)背景月中值異常增強(qiáng), 且與地磁低點(diǎn)位移異常出現(xiàn)時(shí)間有很好的對(duì)應(yīng)關(guān)系. 徐彤等 (2012) 利用電離層垂測(cè)資料分析了20世紀(jì)以來的14次中國大陸強(qiáng)震, 結(jié)果顯示85.7%的地震前出現(xiàn)foF2異常, 大多數(shù)異常出現(xiàn)在震前7天內(nèi), 并集中在當(dāng)?shù)貢r(shí)段11:00—17:00. 由此可見, 異常擾動(dòng)強(qiáng)度與震級(jí)、 震源深度及震中距均有一定的相關(guān)性.

    鑒于電離層測(cè)高自身的觀測(cè)技術(shù)特點(diǎn), 從頻高圖中直接測(cè)量得到的foF2精度較高, 可以作為天地對(duì)比觀測(cè)中校驗(yàn)其它數(shù)據(jù)的參考源; 但如果換算到電子密度等參量進(jìn)行對(duì)比分析, 則需要借助相關(guān)公式進(jìn)行反演計(jì)算, 同時(shí)會(huì)導(dǎo)致精度大大降低. 另外, 電離層測(cè)高技術(shù)還可以獲得電離層D層和E層的峰值信息, 不同層位的峰區(qū)信息對(duì)于約束電離層結(jié)構(gòu)以及研究地震擾動(dòng)電離層效應(yīng)耦合機(jī)理非常有意義, 即使僅有各層峰值信息, 其精度也是其它探測(cè)技術(shù)所無法達(dá)到的.

    1.2 基于GPS觀測(cè)的TEC

    當(dāng)GPS衛(wèi)星發(fā)射的高頻信號(hào)在電離層中傳播時(shí), 由于受到電離層介質(zhì)的折射影響, 產(chǎn)生了附加的信號(hào)傳播時(shí)延, 使所接收信號(hào)的傳播時(shí)間產(chǎn)生誤差. 在影響GPS定位精度的各種誤差因素中, 電離層是僅次于衛(wèi)星星歷誤差的誤差源, 電離層所引起的GPS誤差主要與沿衛(wèi)星至接收機(jī)收視方向上的TEC有關(guān). 基于電磁波傳播與TEC的相關(guān)性, 利用電離層對(duì)電磁波的折射效應(yīng)即可研究電離層本身, 即利用GPS衛(wèi)星信號(hào)獲得的偽距ρ(ρ=cΔt, 其中Δt為某顆衛(wèi)星發(fā)送信號(hào)至接收到該信號(hào)的時(shí)間差,c為真空中電波速度, 實(shí)際電波速度低于該值, 因此據(jù)此計(jì)算的距離ρ稱為偽距)和載波相位差分可計(jì)算得到信號(hào)傳播路徑上的TEC, 具體公式(夏淳亮, 2004)為

    (1)

    式中: Δρ為偽距差分,f1和f2為GPS衛(wèi)星工作頻率,f1=1.57542 GHz,f2=1.2276 GHz.

    同樣, 利用載波的相位測(cè)量, 考慮到整周模糊度(因?yàn)檩d波在空間傳輸?shù)恼芷跀?shù)無法通過觀測(cè)直接獲得, 接收機(jī)只能記錄這個(gè)差值的非整周部分, 其整周部分未知, 故稱為整周模糊度)可計(jì)算相對(duì)TEC, 即

    (2)

    式中, Ф1和Ф2分別為對(duì)應(yīng)GPS衛(wèi)星工作頻率f1和f2的載波相位, Δn為整周模糊度. 因式(2)中含有一定的未知參量, 只能得到TEC的相對(duì)變化, 因此稱其為相對(duì)TEC, 其測(cè)量精度為1013; 而通過偽距差分得到的TEC稱為絕對(duì)TEC, 其測(cè)量精度可達(dá)1016.

    GPS觀測(cè)由于具有觀測(cè)及時(shí)、 覆蓋范圍廣、 觀測(cè)精度高、 成本低、 效率高等優(yōu)勢(shì)而被廣泛應(yīng)用于電離層研究中. 美國噴氣動(dòng)力實(shí)驗(yàn)室(Jet Propulsion Laboratory, 簡寫為JPL)利用IGS(International GNSS Service)發(fā)布的全球GPS觀測(cè)站數(shù)據(jù), 建立了TEC現(xiàn)報(bào)系統(tǒng)GIM(global ionospheric maps), 其空間分辨率為2.5°×5°, 時(shí)間分辨率為2小時(shí), 是國際上認(rèn)可度較高的電離層現(xiàn)報(bào)系統(tǒng). 利用JPL發(fā)布的GIM數(shù)據(jù), 我國科研人員也開展了大量的地震電離層應(yīng)用研究, 例如: Zhao等(2008)分析了2008年汶川MS8.0地震前的TEC異常, 結(jié)果顯示震前3天(5月9日)震中東南方向出現(xiàn)局部TEC明顯增強(qiáng)現(xiàn)象, 且在南半球的磁共軛區(qū)觀測(cè)到同步擾動(dòng); Liu等(2013)研究了1998—2012年我國56次M>6地震的TEC擾動(dòng)特征, 結(jié)果顯示震前2—9天TEC在下午時(shí)段呈快速下降異常, 與位于低緯度的印尼地區(qū)和我國臺(tái)灣地區(qū)的地震TEC擾動(dòng)具有相似的特征.

    中國地殼運(yùn)動(dòng)觀測(cè)網(wǎng)絡(luò)是“九五”國家重大科學(xué)工程建設(shè)項(xiàng)目, 自1998年投入開工建設(shè), 目前已形成具有260個(gè)GPS基準(zhǔn)站、 千余個(gè)不定期復(fù)測(cè)GPS站的區(qū)域觀測(cè)網(wǎng)絡(luò)(甘衛(wèi)軍等, 2007), 這些臺(tái)站的建設(shè)為反演產(chǎn)出我國大陸地區(qū)GPS TEC數(shù)據(jù)奠定了良好的基礎(chǔ). 利用國內(nèi)GPS網(wǎng)絡(luò)觀測(cè)數(shù)據(jù), 張學(xué)民等(2009a)關(guān)于汶川地震前后電離層擾動(dòng)異常的研究顯示, 距離震中較近的瀘州站TEC在震前3天(5月9日)與重慶站測(cè)高站觀測(cè)到的foF2均呈快速增強(qiáng). 林劍等(2009)利用中國25個(gè)基準(zhǔn)站結(jié)合國際IGS網(wǎng)站發(fā)布的全球數(shù)據(jù), 研究了汶川地震前后的TEC擾動(dòng)空間異常, 結(jié)果顯示地震前后一周內(nèi)震區(qū)上空連續(xù)出現(xiàn)異常擾動(dòng)且具備共軛特性, 并向磁赤道漂移. 李建勇(2010)利用全國GPS基準(zhǔn)站信息對(duì)2000—2009年中國大陸33次M>6地震的統(tǒng)計(jì)分析結(jié)果表明: 中強(qiáng)地震電離層TEC異常多出現(xiàn)在震前數(shù)天, 為非持續(xù)性變化, 異常持續(xù)時(shí)間多為幾個(gè)小時(shí); TEC異常的變化形態(tài)有正有負(fù), 相對(duì)變化幅度一般為20%—30%; 中國大陸TEC異常區(qū)多偏震中東南方向, 可能與電離層擾動(dòng)傳播受地磁場(chǎng)控制向磁赤道偏移有關(guān). 另一方面, TEC數(shù)據(jù)還被用于同震及震后效應(yīng)研究. 例如, Song等(2015)利用中國GPS觀測(cè)網(wǎng)提供的華南地區(qū)46個(gè)GPS接收臺(tái)站的TEC觀測(cè)數(shù)據(jù), 對(duì)汶川地震所引起的電離層同震效應(yīng)進(jìn)行分析研究, 結(jié)果顯示震后數(shù)分鐘發(fā)生兩個(gè)擾動(dòng)事件, 擾動(dòng)振幅為0.5—1 TECU, 水平相速度為500—800 m/s, 周期為6—10分鐘, 且分別向南、 東兩個(gè)方向傳播, 分析認(rèn)為這是由于震后地面垂直運(yùn)動(dòng)而激發(fā)的聲重力波傳播到達(dá)電離層高度所致. 隨著陸態(tài)網(wǎng)絡(luò)的完善及數(shù)據(jù)的實(shí)時(shí)發(fā)布更新(蔡華等, 2014), 數(shù)據(jù)時(shí)空覆蓋率不斷提高, 我國區(qū)域GPS網(wǎng)的TEC必將在地震前后的局地電離層時(shí)空演化特征研究中發(fā)揮更重要的作用.

    隨著全球GPS網(wǎng)絡(luò)的快速發(fā)展, 全球及區(qū)域性TEC分布圖的時(shí)空分辨率會(huì)持續(xù)提高, 為地震前后的電離層擾動(dòng)變化研究提供更有力的數(shù)據(jù)支撐. GPS臺(tái)站的密集分布不僅能為地震監(jiān)測(cè)預(yù)測(cè)提供時(shí)間信息, 也能很好地研究異常的局地特性和空間傳播特性. 就觀測(cè)技術(shù)而言, TEC需要通過一系列數(shù)據(jù)預(yù)處理技術(shù)剔除儀器本身所攜帶的以及觀測(cè)中所帶來的干擾因素, 反演計(jì)算獲得沿整個(gè)傳播路徑的積分TEC, 因此測(cè)量精度會(huì)受到一定制約, 一些微弱擾動(dòng)(1—2 TECU)很容易落在正常的反演誤差內(nèi), 導(dǎo)致與背景噪聲難以區(qū)分. 這種情況在一些地震臨震前及震后的電離層效應(yīng)研究中尤為突出, 反過來這也對(duì)TEC的反演精度和算法提出了更高的要求. 雖然, 可以通過增設(shè)臺(tái)站改善TEC數(shù)據(jù)的水平空間分辨率, 但該參量無法反映不同高度的信息, 僅是各高度電子濃度的積分綜合效應(yīng).

    1.3 甚低頻(VLF)電波

    甚低頻(3—30 kHz)電磁波主要通過大地與低電離層之間形成的波導(dǎo)進(jìn)行傳播, 其傳播損耗近似于自由空間, 只存在因能量擴(kuò)散引起的傳播損耗, 距離可達(dá)數(shù)千千米乃至覆蓋全球, 而且該頻段電磁波穿透海水的能力較強(qiáng), 適用于水下艦艇的遠(yuǎn)距離通信等. 迄今為止, 全球仍保留20—30個(gè)VLF導(dǎo)航通信臺(tái)(Wikipedia, 2015), 服務(wù)于陸地和海上通信.

    地面VLF電波接收觀測(cè)包含兩種模式, 一種是該頻段的天然信號(hào)觀測(cè), 另一種是利用全球已有的VLF電波導(dǎo)航信號(hào)發(fā)射站配合地面接收站接收該固定站點(diǎn)的固定頻段信號(hào), 并分析信號(hào)傳播鏈路上所發(fā)生的擾動(dòng). 對(duì)于第一種觀測(cè)模式而言, 存在于頻帶范圍內(nèi)的各類噪聲源會(huì)損害觀測(cè)數(shù)據(jù), 只有當(dāng)觀測(cè)站距離震中非常近時(shí)才能獲得一些與地震相關(guān)的局地異常; 第二種觀測(cè)模式由于VLF發(fā)射站頻率固定、 發(fā)射功率已知, 因此可以在接收站觀測(cè)到發(fā)射—接收站傳播路徑上來自地震-大氣層和地震-電離層對(duì)電磁波傳播造成的影響, 從而使觀測(cè)區(qū)域覆蓋更大的范圍, 在地震監(jiān)測(cè)應(yīng)用中多被采用.

    日本在VLF電波觀測(cè)領(lǐng)域開展的地震應(yīng)用研究較早, 目前意大利、 南美等國家和地區(qū)也陸續(xù)建立了AWESOME和SAVNET等VLF電波觀測(cè)網(wǎng), 用于研究地震所引起的低電離層擾動(dòng)與電波傳播的相關(guān)性. Hayakawa等(2010)全面總結(jié)了VLF/LF (3—300 kHz)電波觀測(cè)和異常提取算法方面的研究發(fā)展, 其對(duì)日本2001—2007年長達(dá)7年的VLF電波觀測(cè)資料的分析結(jié)果顯示, 日本境內(nèi)淺部(<40 km)M>6.0地震的低電離層VLF/LF異常波動(dòng)幅度均超過背景數(shù)據(jù)的兩倍標(biāo)準(zhǔn)差.

    2010年前后, 中國地震局地震預(yù)測(cè)研究所在北京、 通海和雅安建立了3個(gè)VLF電波接收站, 分別接收位于俄羅斯東部、 中部和西部3個(gè)阿爾法導(dǎo)航站發(fā)射的3個(gè)頻率(11.9, 12.6, 14.9 kHz)的VLF電波信號(hào), 形成了9條鏈路覆蓋中國北部大部分區(qū)域的電波監(jiān)測(cè)網(wǎng), 可通過電波幅度和相位變化提取擾動(dòng)信息. 3個(gè)VLF電波接收站運(yùn)行期間, 覆蓋區(qū)域內(nèi)發(fā)生多次M>6.0地震, Zeren等(2014) 通過分析VLF電波傳播路徑上的擾動(dòng)信號(hào), 結(jié)果顯示2010年玉樹和2013年蘆山兩次M>7.0地震前在地磁和太陽活動(dòng)平靜的情況下均監(jiān)測(cè)到震中鄰近傳播路徑上的低電離層擾動(dòng). 由于我國目前所使用的VLF電波接收儀僅能接收來自俄羅斯3個(gè)站點(diǎn)3個(gè)頻率的信息, 受鏈路空間所限能獲得的震例很有限. 為加強(qiáng)我國西南部地區(qū)的地震監(jiān)測(cè), 計(jì)劃在未來1—2年內(nèi)配置寬頻VLF電波接收儀, 以獲取更多東部和南半球國家的VLF發(fā)射站信號(hào), 拓展國內(nèi)VLF電波的空間監(jiān)測(cè)范圍.

    VLF人工源電波觀測(cè)的發(fā)射站、 接收站位置固定, 發(fā)射頻率和功率穩(wěn)定, 因此在接收站較易鎖定信號(hào), 易于獲取觀測(cè)相對(duì)穩(wěn)定的信息記錄. 但是, 電波傳播受路徑上各類電離層背景狀態(tài)及不規(guī)則體擾動(dòng)影響, 其傳播距離越長, 擾動(dòng)源的位置越難判定, 故需提高空間覆蓋分辨率, 通過多點(diǎn)多源對(duì)未來孕震區(qū)進(jìn)行預(yù)測(cè).

    1.4 舒曼諧振

    地球和電離層可以構(gòu)成一個(gè)諧振腔體, 腔體中存在一組特殊的諧振頻率(Schumann, 1952)

    (3)

    式中,fn為第n階頻率,c為光速,a為地球半徑,n=1, 2, 3…. 諧振頻率主要由地球及諧振腔的尺寸決定, 并由全球的閃電放電激發(fā), 這個(gè)現(xiàn)象被稱為舒曼(Schumann)諧振. Balser和Wagner(1960)在觀測(cè)自然噪音時(shí), 真實(shí)記錄了位于7.8, 14.2, 19.6, 25.9和32 Hz的舒曼波譜峰, 證實(shí)了地球腔體的舒曼諧振效應(yīng). 采用常規(guī)的超低頻段(super-low frequency, 簡寫為SLF)(30—300 Hz)電磁儀器進(jìn)行觀測(cè), 均能在電磁頻譜分析圖上觀察到該信號(hào). 舒曼諧振多用于全球閃電活動(dòng)、 全球變暖、 低電離層屬性等方面的研究. 在地震應(yīng)用領(lǐng)域, Hayakawa等(2005)最先報(bào)道了1999年臺(tái)灣集集地震前出現(xiàn)高階舒曼諧振幅度增強(qiáng)以及諧振信號(hào)頻移等現(xiàn)象. Ohta等(2006)通過分析1999—2004年臺(tái)灣發(fā)生的33次M>5.0地震開展了舒曼諧振與地震的統(tǒng)計(jì)關(guān)系研究, 結(jié)果表明所有陸地地震前均有舒曼諧振異常發(fā)生. Hayakawa等(2008)通過改變電離層D層高度模擬地震引起的舒曼諧振異常, 并將其與全球三大閃電源(非洲中部的剛果盆地、 東南亞的島嶼地區(qū)、 南美洲的亞馬遜雨林地區(qū))疊加, 模擬震源上空電離層變化引起的高階舒曼諧振幅度增強(qiáng)現(xiàn)象, 該模擬結(jié)果與實(shí)測(cè)結(jié)果比較一致.

    我國的舒曼諧振觀測(cè)試驗(yàn)始于2010年. 中國地震局地震預(yù)測(cè)研究所先后在云南省架設(shè)了巧家、 永勝、 通海和芒市等4個(gè)臺(tái)站, 均采用SRG-M01舒曼諧振儀進(jìn)行觀測(cè), 數(shù)據(jù)采樣率為100 Hz, 主要通過感應(yīng)線圈觀測(cè)磁場(chǎng)三分量中所攜帶的舒曼諧振信號(hào). 苗園青等(2011)的實(shí)測(cè)結(jié)果顯示, 觀測(cè)到的前四階舒曼諧振波分布在8, 14, 20, 26 Hz頻率附近, 類似的舒曼諧振信號(hào)在中國地震局地質(zhì)研究所布設(shè)的超低頻電磁觀測(cè)網(wǎng)(頻率范圍含SLF/ELF/TLF, 0.1—300 Hz)中也被記錄到(范曄等, 2013), 反映了舒曼諧振廣泛存在的特征. 2011年3月日本MS9.0地震后, 周洪娟等(2013)通過分析云南地區(qū)的舒曼諧振磁場(chǎng)觀測(cè)認(rèn)為, 日本地震前3—4天低階舒曼諧振出現(xiàn)比較明顯的幅度增強(qiáng); Zhou等(2013)的模擬計(jì)算結(jié)果表明, 日本地震確實(shí)能引起云南地區(qū)1—3階舒曼異常, 進(jìn)一步確認(rèn)了地-電離層空腔中地震上空不均勻體與舒曼諧振異常之間的相關(guān)性. 歐陽新艷等(2015)利用永勝臺(tái)資料研究了中國大陸M>6.0地震的舒曼諧振異常特征, 結(jié)果表明震前諧振幅度增強(qiáng), 并伴隨諧振頻率和半寬的變化. 由此可見, 舒曼諧振確實(shí)可以作為分析地震前電磁和電離層信號(hào)擾動(dòng)的一種探測(cè)技術(shù).

    舒曼諧振的產(chǎn)生與地球三大閃電源密切相關(guān), 在全球范圍內(nèi)產(chǎn)生共振效應(yīng), 異常信號(hào)產(chǎn)生后疊加在諧振信號(hào)中, 因此相距較遠(yuǎn)的臺(tái)站也能接收到來自震源的異常信號(hào), 也正因如此, 信號(hào)的來源往往難以確定, 目前國內(nèi)的觀測(cè)結(jié)果也顯示震前一般出現(xiàn)低階諧振信號(hào)增強(qiáng)現(xiàn)象, 但遠(yuǎn)距離強(qiáng)震與近區(qū)中小地震激發(fā)效應(yīng)比較類似. 對(duì)于震源激發(fā)信號(hào)與舒曼諧振信號(hào)的疊加調(diào)制及其在傳播過程中隨空腔結(jié)構(gòu)變化的耦合機(jī)理尚待進(jìn)一步研究. 相對(duì)VLF電波的鏈路觀測(cè), 舒曼諧振的全球影響效應(yīng)更為顯著, 如何使之在地震監(jiān)測(cè)中發(fā)揮更有效的作用, 尚需更多資料積累和更深入的理論研究.

    2 衛(wèi)星探測(cè)

    衛(wèi)星的飛行高度一般在幾百千米以上, 直接位于電離層甚至磁層中. 與地基電離層間接探測(cè)模式不同, 衛(wèi)星觀測(cè)多屬于直接探測(cè)或就位探測(cè), 其電磁場(chǎng)觀測(cè)技術(shù)和原位等離子體參量觀測(cè)技術(shù)使用較多. 為了獲得更多電離層精細(xì)結(jié)構(gòu), 利用GPS及衛(wèi)星本體發(fā)射電波信號(hào)反演電離層結(jié)構(gòu)的探測(cè)技術(shù)也在快速發(fā)展, 如掩星探測(cè)、 信標(biāo)探測(cè)等. 下面將針對(duì)不同種類的衛(wèi)星探測(cè)技術(shù)作簡單介紹.

    2.1 衛(wèi)星電磁場(chǎng)探測(cè)數(shù)據(jù)

    空間電磁場(chǎng)觀測(cè)基于與地基觀測(cè)相同的觀測(cè)原理, 探測(cè)儀器也較接近. 一般采用磁強(qiáng)計(jì)、 磁通門磁力儀和感應(yīng)式磁力儀分別觀測(cè)地磁總場(chǎng)、 地磁場(chǎng)三分量、 變化磁場(chǎng)三分量, 利用一定長度的不共面電極對(duì)組合觀測(cè)電場(chǎng)矢量. 由于電磁波頻段不同, 其在電離層的傳播模式及電離層擾動(dòng)響應(yīng)都存在較大差別, 因此在地震監(jiān)測(cè)應(yīng)用中須針對(duì)不同頻段的電磁波分別進(jìn)行分析.

    對(duì)于ELF頻段(3—30 Hz)的地震電磁信號(hào), 由于其采樣率較低, 一般直接從電磁波形記錄中判識(shí)地震異常. Chmyrev等(1989)利用位于800—900 km高度的Intercosmos-Bulgaria-1300衛(wèi)星資料, 在一次MS4.8地震前15分鐘衛(wèi)星準(zhǔn)靜電場(chǎng)垂直分量上出現(xiàn)3—7 mV/m的異常信號(hào). Gousheva等(2008)基于該衛(wèi)星數(shù)據(jù)的統(tǒng)計(jì)分析結(jié)果顯示, 在低、 中、 高緯度地震前, 電離層電場(chǎng)的水平分量和垂直分量均會(huì)出現(xiàn)2—15 mV/m的異常增強(qiáng). Zhang等(2012b)通過研究汶川地震前大地電場(chǎng)和衛(wèi)星電場(chǎng)擾動(dòng), 認(rèn)為地基電場(chǎng)與DEMETER衛(wèi)星記錄的電場(chǎng)矢量在震前幾個(gè)月的異常信號(hào)呈同步增強(qiáng). Zhang等(2014)利用DEMETER衛(wèi)星記錄分析了2010年全球M>7.0地震, 并針對(duì)赤道地區(qū)的印尼和中緯度的智利地區(qū)27次強(qiáng)震進(jìn)行了統(tǒng)計(jì)分析, 結(jié)果顯示電場(chǎng)擾動(dòng)幅值為1.5—16 mV/m, 其中10次地震前1天有明顯的電場(chǎng)擾動(dòng)信號(hào)出現(xiàn).

    對(duì)于SLF(30—300 Hz)/ULF頻段(0.3—3 kHz)的地震電磁擾動(dòng), Gokhberg等(1983)報(bào)道1970年3月30日發(fā)生的一次地震前OGO-6衛(wèi)星在100, 216, 467 Hz均記錄到該頻段磁場(chǎng)強(qiáng)度增加的現(xiàn)象. 2004年法國DEMETER衛(wèi)星發(fā)射升空后, Parrot等 (2006) 利用波矢量分析技術(shù)對(duì)DEMETER衛(wèi)星記錄的電磁波形數(shù)據(jù)進(jìn)行分析, 結(jié)果顯示, 在2005年1月23日一次MS6.2地震前2.5天, 震中上空軌道記錄到130 Hz處出現(xiàn)異常電磁波信號(hào). Němec等(2009)關(guān)于DEMETER衛(wèi)星的全球地震統(tǒng)計(jì)分析結(jié)果顯示, 震前4小時(shí)1.7 kHz附近的電磁場(chǎng)功率譜強(qiáng)度快速降低. Pí?a等(2013) 利用DEMETER運(yùn)行期間2004年7月—2010年11月長達(dá)6.5年的資料, 再次確認(rèn)了震前幾小時(shí)1.7 kHz附近出現(xiàn)電磁場(chǎng)功率譜強(qiáng)度下降異常, 并提出這類異常在3月份和8月份的高緯度地區(qū)和海域更為突出. Zhang等(2014)總結(jié)了DEMETER衛(wèi)星運(yùn)行期間(2005—2010年)全球中低緯度區(qū)總共發(fā)生的69次M>7.0地震的電磁異常, 結(jié)果顯示19.5—250 Hz頻段內(nèi)的電場(chǎng)功率譜在32次地震前會(huì)增加1—2個(gè)數(shù)量級(jí)以上. 于海雁等(2010)利用該衛(wèi)星資料對(duì)汶川地震前SLF/ULF頻段電磁波矢特征的分析結(jié)果顯示, 震前6天存在窄帶電磁輻射波(約380 Hz), 該電磁輻射信號(hào)具有很強(qiáng)的平面波特性, 且為上行傳播, 與汶川地震孕育有一定的相關(guān)性. 澤仁志瑪?shù)?2012)統(tǒng)計(jì)分析了2005—2009年北半球M>7.0強(qiáng)震前后370—897 Hz頻段的磁場(chǎng)功率譜密度變化, 結(jié)果顯示26次強(qiáng)震中77%的地震震前磁場(chǎng)增強(qiáng), 在異常最高點(diǎn)或者其轉(zhuǎn)折下降時(shí)發(fā)生地震. 因此, 空間SLF頻段電磁觀測(cè)的地震異常既有相對(duì)背景場(chǎng)的頻譜能量增強(qiáng)現(xiàn)象, 也有單純的具有平面波傳播特性的電磁波信號(hào)出現(xiàn).

    VLF頻段(3—30 kHz)的電磁擾動(dòng), 由于其采樣率較高, 主要以電磁頻譜特性分析為主; 而且由于一些地面導(dǎo)航站的發(fā)射信號(hào)也在十幾到幾十kHz頻段內(nèi), 分別有針對(duì)天然源和人工源的電磁擾動(dòng)信號(hào)分析. Molchanov等(2006)分析了衛(wèi)星記錄的人工源信號(hào)的信噪比變化, 結(jié)果顯示幾次強(qiáng)震前震中上空信噪比大范圍大幅度降低. 張學(xué)民等(2009b)針對(duì)汶川地震前DEMETER衛(wèi)星記錄的VLF電場(chǎng)異常的研究結(jié)果表明, 震前1—5天2—6 kHz電場(chǎng)頻譜在震中上空相對(duì)5—10天的背景信息顯著增強(qiáng). 朱濤和王蘭煒(2011)的研究則顯示汶川地震上空當(dāng)?shù)貢r(shí)間白天10—14 kHz和夜間10—18 kHz的電場(chǎng)在震中周圍300 km范圍內(nèi), 2008年4—5月數(shù)據(jù)相對(duì)其它無震年份呈現(xiàn)“震前一個(gè)半月平均功率譜明顯增強(qiáng)—震前半個(gè)月減小—震時(shí)達(dá)到最小—震后恢復(fù)”的時(shí)空演化特征. 何宇飛等(2009)應(yīng)用信噪比方法得到了汶川地震前一個(gè)月震中上空也是澳大利亞??怂姑┧筃WC (North West Cape)人工源發(fā)射站(發(fā)射頻率為19.8 kHz)的共軛區(qū), 而且2008年4月的電場(chǎng)信噪比相比2006和2007年顯著降低.

    目前SLF頻段的波矢分析應(yīng)用相對(duì)較多, 但更高頻段的電磁波后向追蹤技術(shù)仍處于初步發(fā)展中, 故未能有效地應(yīng)用于實(shí)測(cè)VLF電磁波數(shù)據(jù)分析中. 震例研究結(jié)果已證明TLF—VLF頻段衛(wèi)星觀測(cè)在震前的響應(yīng)事實(shí), 但由于目前地基相對(duì)應(yīng)的高頻段電磁場(chǎng)觀測(cè)很少, 或運(yùn)行時(shí)段不一致, 因此未能有可相互驗(yàn)證的天地同步地震電磁擾動(dòng)信號(hào)被觀測(cè)到. 為了配合衛(wèi)星觀測(cè)和解釋地震電離層電磁場(chǎng)擾動(dòng)機(jī)理, 地基系統(tǒng)需要配備同頻段觀測(cè)、 同步監(jiān)測(cè)電磁信號(hào)從巖石層、 大氣層至電離層的激發(fā)和傳播過程, 以實(shí)現(xiàn)數(shù)據(jù)的相互校驗(yàn)和圈層耦合理論的發(fā)展完善.

    2.2 原位等離子體參量

    原位等離子體參量主要指衛(wèi)星高度的電子、 離子成分、 高能粒子通量等多種等離子體參量, 是衛(wèi)星觀測(cè)中維度最低的一些參量, 僅與觀測(cè)時(shí)間和衛(wèi)星軌道位置有關(guān). 其中電子密度是聯(lián)系天基觀測(cè)與地基觀測(cè)的重要參量, 是天地對(duì)比分析的橋梁, 電離層測(cè)高、 GPS反演TEC等最終均與電離層中電子含量密切相關(guān).

    Lebreton等(2006)詳細(xì)介紹了DEMETER衛(wèi)星上朗繆爾(Langmuir)探針的觀測(cè)原理, 使用朗繆爾探針獲得的伏安曲線進(jìn)行反演即可獲得衛(wèi)星原位電子密度和溫度信息; Bertherlier等(2006)介紹了等離子體分析儀觀測(cè)技術(shù), 利用阻滯儀和離子漂移計(jì)兩套設(shè)備可分別獲得離子成分、 離子溫度和離子漂移速度等信息.

    原位等離子體探測(cè)由于設(shè)備相對(duì)輕便, 很多空間探測(cè)衛(wèi)星上均搭載有此類載荷, 因此其地震應(yīng)用實(shí)例非常豐富. 基于Intercosmos-24衛(wèi)星(高度2300—2500 km), Bo?ková等(1994) 研究了1990年6月20日一次伊朗地震前的離子成分異常, 結(jié)果顯示震中上空軌道記錄的輕離子密度明顯增加. Pulinets等(2003) 利用AE-C衛(wèi)星上搭載的質(zhì)譜儀也發(fā)現(xiàn)震中區(qū)域內(nèi)平均離子質(zhì)量下降的現(xiàn)象. 2004年DEMETER衛(wèi)星發(fā)射升空后, Parrot 等(2006)利用該衛(wèi)星數(shù)據(jù)最先開展了多次強(qiáng)震前的研究工作, 結(jié)果顯示震前電子密度和離子密度均有增大. 之后, 世界各國研究人員利用DEMETER衛(wèi)星資料開展了大量的研究, 中國也通過國際合作開始介入衛(wèi)星地震電磁數(shù)據(jù)分析研究, 并取得了一定的成果. 張學(xué)民等(2008)、 歐陽新艷等(2008)和劉靜等(2011)先后開展了關(guān)于我國西藏改則MS6.9、 云南普洱MS6.4、 智利MS8.8地震的電離層異常識(shí)別的研究, 通過對(duì)比相鄰軌道與重訪軌道數(shù)據(jù), 獲得了離子溫度、 電子密度等多個(gè)參量的震前擾動(dòng)現(xiàn)象. 除了單個(gè)震例的分析研究, 多震例統(tǒng)計(jì)近年來也得到快速發(fā)展, 統(tǒng)計(jì)結(jié)果(Li, Parrot, 2013; Zhangetal, 2013; Liuetal, 2014)反映了電離層高能粒子沉降、 離子密度及電子密度等擾動(dòng)與地震在時(shí)間、 空間上的相關(guān)性.

    原位等離子體參量是地震研究中異常提取最為豐富的參量, 但由于電離層本身存在的赤道電離層異常以及逐日、 季節(jié)等多種擾動(dòng)特性, 衛(wèi)星觀測(cè)原位等離子體參量的異常閾值界定在不同空間區(qū)域內(nèi)會(huì)有較大波動(dòng), 也會(huì)出現(xiàn)較多的虛假異常信息, 因此一個(gè)合理的電離層背景場(chǎng)構(gòu)建及異常閾值的界定一直是這方面的研究重點(diǎn).

    2.3 電離層結(jié)構(gòu)探測(cè)

    由于地面電離層測(cè)高一般只能獲得觀測(cè)站所在固定地理位置上電離層臨界頻率隨時(shí)間變化的信息, 難以獲得大尺度電離層的空間結(jié)構(gòu)信息, 也難以精確地探測(cè)中性分子密度大而電子密度小的電離層D層、 E層和F層峰值過后的下降區(qū), 無法全面滿足電離層研究的需求, 因此更多的電離層結(jié)構(gòu)探測(cè)技術(shù)得以發(fā)展, 下面主要介紹掩星和三頻信標(biāo)(tri-band beacon, 簡寫為TBB)探測(cè)技術(shù).

    2.3.1 掩星探測(cè)技術(shù)與方法

    衛(wèi)星掩星探測(cè)技術(shù)是指GPS衛(wèi)星發(fā)射的電波信號(hào)被地球大氣層所遮掩, 經(jīng)過地球大氣層和電離層折射后到達(dá)近地軌道衛(wèi)星, 也即對(duì)GPS衛(wèi)星的臨邊觀測(cè). 衛(wèi)星接收的GPS信號(hào)中包含了地球大氣層和電離層的信息, 可通過相關(guān)的反演算法得到大氣層和電離層的剖面結(jié)構(gòu).

    阿貝爾(Abel)積分是目前運(yùn)用最廣泛的電離層掩星反演方法, 對(duì)于GPS載波頻率, 即使在電離層變化最劇烈的狀態(tài)下, 電離層所造成的彎曲角仍然相當(dāng)小, 因此可以把掩星事件發(fā)生時(shí)的傳播路徑近似為直線傳播, 反演公式為(Schreineretal, 1999):

    (4)

    式中:Ne為電子密度;p為掩星切點(diǎn)到地心的距離, 稱為碰撞高度; STEC為斜TEC, 其實(shí)是路徑TEC, 因?yàn)椴皇谴怪睖y(cè)點(diǎn)下方而被稱為斜TEC;r為剖面高度.

    需要說明的是, 阿貝爾積分反演算法包含了3個(gè)假設(shè): ① 假設(shè)信號(hào)傳播路徑為直線; ② 假設(shè)局部球?qū)ΨQ; ③ 假設(shè)低軌衛(wèi)星軌道高度以上的電子含量忽略不計(jì). 此外, 阿貝爾積分反演算法中存在積分上下限問題. 為了消除以上因素對(duì)反演精度的影響, 電子密度的反演算法仍在不斷發(fā)展中(劉經(jīng)南等, 2010), 如“洋蔥分層”反演算法(Leietal, 2007)等.

    美國GPS/MET計(jì)劃于1995年發(fā)射MicroLab-1低軌衛(wèi)星, 首次實(shí)現(xiàn)地球大氣層和電離層的無線電GPS掩星觀測(cè). 目前已發(fā)射的用于GNSS掩星技術(shù)研究的低軌衛(wèi)星計(jì)劃還有德國的CHAMP計(jì)劃、 阿根廷的SAC-C計(jì)劃、 丹麥的Orsted計(jì)劃、 美國和中國臺(tái)灣合作的COSMIC計(jì)劃、 澳大利亞的FedSat計(jì)劃等(王也英等, 2009), 國內(nèi)已有氣象衛(wèi)星如風(fēng)云3號(hào), 未來的電磁衛(wèi)星CSES上也將搭載GNSS無線電掩星設(shè)備. 在上述掩星計(jì)劃中, COSMIC計(jì)劃是迄今為止世界上最大的掩星計(jì)劃. COSMIC系統(tǒng)由6顆相同的微衛(wèi)星組成, 每顆衛(wèi)星裝載有GPS接收機(jī)、 小型電離層光度計(jì)和三頻信標(biāo)儀等科學(xué)探測(cè)儀. 這6顆衛(wèi)星分布在約800 km 軌道高度, 傾角為72°, 相鄰衛(wèi)星軌道平面夾角為30°. 掩星接收機(jī)通過接收GPS與衛(wèi)星之間的掩星事件, 計(jì)算電波在電離層傳播過程中產(chǎn)生的相位延時(shí)量, 利用相位差分反演得到電離層TEC及電子密度剖面.

    COSMIC衛(wèi)星已經(jīng)在全球電離層變化特征、 天氣和氣候預(yù)報(bào)以及地震電離層異常探測(cè)中取得許多成果. Hsiao等(2010)通過FORMOSAT-3/COSMIC掩星觀測(cè)對(duì)汶川地震前后電離層擾動(dòng)的詳細(xì)分析結(jié)果顯示, 地震前5天內(nèi)正午左右, 靠近震源區(qū)上空F2層峰值高度(hmF2)下降約25 km, 峰值密度(NmF2)下降2×105el/cm3, TEC在250—300 km高度間下降2 TECU. Kakinami等(2010)基于COSMIC掩星數(shù)據(jù)構(gòu)建了電離層電子含量(ionospheric electron content, 簡寫為IEC)經(jīng)驗(yàn)?zāi)P停?模型數(shù)據(jù)與實(shí)測(cè)數(shù)據(jù)的對(duì)比結(jié)果顯示, 汶川地震前6天內(nèi)除了第3天即5月9日震中附近IEC顯著增強(qiáng), 其余5天IEC均呈下降異常. 楊劍等(2008)利用COSMIC掩星觀測(cè)的電離層電子密度剖面研究了2006年11月15日本千島MS8.0地震和2007年2月25日青海海西MS5.3地震, 結(jié)果表明臨震前地震上空附近的電子密度存在異常變化. 耿光宇(2010)利用COSMIC掩星數(shù)據(jù)研究了2009年9月30日蘇門答臘MS7.7地震前的電離層異?,F(xiàn)象, 采用人工神經(jīng)網(wǎng)絡(luò)插值技術(shù)獲得了研究區(qū)域電離層電子密度峰值的空間分布, 結(jié)果顯示震前5天在地磁和太陽活動(dòng)平靜的情況下震中上空NmF2出現(xiàn)顯著下降, 可能與該地震的孕育有關(guān). 姚璐等(2014)結(jié)合JPL發(fā)布的GPS TEC和COSMIC掩星數(shù)據(jù)研究了2010年4月13日玉樹MS7.1地震前的電離層異常變化, 結(jié)果顯示地震當(dāng)天出現(xiàn)明顯的TEC增強(qiáng)現(xiàn)象, 而掩星探測(cè)電子密度剖面在4月9—10日的異常更為突出, 兩者不能完全吻合可能與掩星探測(cè)剖面出現(xiàn)的概率較低及其時(shí)空分布太過稀疏有關(guān).

    利用GPS衛(wèi)星和近地軌道衛(wèi)星的掩星觀測(cè)具有全球覆蓋、 全天候、 高垂直分辨率和高精度等特點(diǎn), 可提供大氣(近地面至60 km)參量和電離層電子密度剖面, 但觀測(cè)數(shù)據(jù)在空間上隨機(jī)分布, 無法保證研究區(qū)域的長時(shí)間連續(xù)覆蓋. 當(dāng)研究某一區(qū)域時(shí), 如經(jīng)常出現(xiàn)數(shù)據(jù)缺乏的情形, 則需借助插值技術(shù)、 現(xiàn)有經(jīng)驗(yàn)?zāi)P图捌渌^測(cè)數(shù)據(jù)進(jìn)行時(shí)空分布的補(bǔ)充. 受低電離層反演精度的制約, 地震應(yīng)用研究中一般選擇可信度和測(cè)量精度較高的TEC或電離層峰值高度F2層的特征參量如峰值密度、 峰值高度等進(jìn)行分析, 其它層位的貢獻(xiàn)則相對(duì)較少.

    2.3.2 三頻信標(biāo)探測(cè)技術(shù)與方法

    Austen 等(1988)最早將層析成像(computerized tomography, 簡寫為CT)方法應(yīng)用于電離層研究, 此后電離層層析成像(computerized ionospheric tomography, 簡寫為CIT)技術(shù)成為衛(wèi)星信標(biāo)和地面接收站配套使用的探測(cè)手段. 三頻信標(biāo)的工作機(jī)制是發(fā)射機(jī)在VHF(30—300 MHz)和UHF(300—3000 MHz )頻段上輸出頻率穩(wěn)定且相位相關(guān)的3個(gè)載頻信號(hào), 并經(jīng)搭載在衛(wèi)星上的全向天線向預(yù)定覆蓋區(qū)域輻射, 位于地面的三頻信標(biāo)接收機(jī)跟蹤鎖定星載三頻信標(biāo)信號(hào)后, 實(shí)時(shí)輸出3個(gè)頻段的同相/正交(in-phase/quadrature)兩路信號(hào). 對(duì)3個(gè)頻段數(shù)據(jù)分別進(jìn)行差分處理, 然后利用下式計(jì)算相對(duì)TEC:

    (5)

    其中,

    式中,f1和f2表示3個(gè)頻率中的任意兩個(gè)頻點(diǎn). 進(jìn)而利用不同頻率的差分相位結(jié)合素?cái)?shù)論原理即可計(jì)算絕對(duì)TEC, 即

    (6)

    式中: Δφ12和Δφ13分別為兩個(gè)頻率之間的差分相位;k2為一正整數(shù), 通常稱為三頻相位積分常數(shù), 可利用多站數(shù)據(jù)結(jié)合雙頻計(jì)算TEC粗差擬合得出;x12和x13可由3個(gè)頻率間的倍數(shù)素?cái)?shù)論原理獲得. 將相對(duì)和絕對(duì)TEC數(shù)據(jù)作為輸入源, 利用代數(shù)重建法、 代數(shù)迭代重建法和奇異值分解法等層析成像技術(shù)即可生成電離層電子密度圖像, 得到電子密度的二維或三維結(jié)構(gòu)分布.

    2006年COSMIC衛(wèi)星發(fā)射成功, 星上搭載了TBB發(fā)射機(jī), 發(fā)射頻率分別為150, 400, 1067 MHz的頻率穩(wěn)定且相位相關(guān)的3組載頻信號(hào), 經(jīng)全向天線向預(yù)定覆蓋區(qū)域輻射. 我國臺(tái)灣地區(qū)據(jù)此建立了低緯電離層CT探測(cè)臺(tái)鏈, 通過接收COSMIC和OSCAR等低軌道衛(wèi)星三頻發(fā)射信號(hào), 對(duì)低緯電離層的形態(tài)和擾動(dòng)進(jìn)行研究, 以便有效監(jiān)測(cè)地震危險(xiǎn)區(qū)的電離層結(jié)構(gòu)變化. 美國從2011年開始計(jì)劃用三頻信標(biāo)裝備24顆衛(wèi)星, 印度和我國臺(tái)灣地區(qū)也啟動(dòng)采用三頻信標(biāo)裝備衛(wèi)星的計(jì)劃, 我國2017年即將發(fā)射的電磁監(jiān)測(cè)衛(wèi)星CSES上也將搭載三頻信標(biāo)發(fā)射機(jī). 這些衛(wèi)星的發(fā)射將為全面開展電離層三維結(jié)構(gòu)探測(cè)提供支撐.

    三頻信標(biāo)技術(shù)研究目前集中于電離層層析成像技術(shù)的更新發(fā)展, 觀測(cè)數(shù)據(jù)被應(yīng)用于電離層數(shù)據(jù)同化等, 截至目前為止未見利用信標(biāo)觀測(cè)得到的地震震例研究結(jié)果, 其地震應(yīng)用效能尚待檢驗(yàn).

    3 地震電磁立體觀測(cè)體系

    隨著地基和空間電離層探測(cè)技術(shù)的快速發(fā)展, 地震前后電磁多參量對(duì)比分析研究陸續(xù)展開. 丁鑒海等(2005)分析了地磁場(chǎng)變化與電離層擾動(dòng)在多次地震前的同步響應(yīng)關(guān)系. 安張輝等(2011)利用由希爾伯特-黃(Hilbert-Huang)變換計(jì)算得到的邊際譜對(duì)天基和地基電場(chǎng)觀測(cè)數(shù)據(jù)予以分析, 結(jié)果顯示DEMETER衛(wèi)星記錄電場(chǎng)和震中周邊地區(qū)臺(tái)站觀測(cè)電場(chǎng)的邊際譜在汶川地震前半個(gè)月左右增強(qiáng)了兩個(gè)數(shù)量級(jí), 天、 地電場(chǎng)擾動(dòng)時(shí)間有較好的一致性. 張學(xué)民等(2014)基于地基測(cè)高、 GPS TEC和DEMETER衛(wèi)星探測(cè)原位等離子體參量, 對(duì)比分析了2010年4月13日玉樹地震前后不同高度電離層特征參量的時(shí)空分布, 確認(rèn)地震前一天自低電離層到衛(wèi)星高度(670 km)出現(xiàn)多參量同步擾動(dòng)現(xiàn)象. 基于汶川地震豐富的各類研究成果, 張學(xué)民等(2009a)從地磁低點(diǎn)位移、 地基GPS TEC、 電離層測(cè)高以及DEMETER衛(wèi)星觀測(cè)等離子體參量, 綜合分析了汶川地震前的電磁多參量異常演化特征, 并綜合應(yīng)力、 紅外、 高光譜等其它探測(cè)手段, 總結(jié)了各探測(cè)參量在地震孕育不同階段的敏感性, 由此可見多參量的協(xié)同監(jiān)測(cè)有助于理解地震孕育發(fā)展的全過程(Zhangetal, 2012a). 對(duì)于同震及震后電離層異常研究, Hao等(2012)針對(duì)2011年日本Tohoku大地震的研究表明, 利用地表高頻多普勒頻移和GPS TEC均觀測(cè)到傳播至電離層的地表振蕩所激發(fā)的地震次聲波, 而電離層中電子密度與電流的相互作用同樣被地磁場(chǎng)觀測(cè)結(jié)果所證實(shí).

    綜上所述, 地基和空基電離層探測(cè)技術(shù)各有其優(yōu)勢(shì)和劣勢(shì), 具體列于表1. 可以看到, 地基探測(cè)可以在固定區(qū)域連續(xù)觀測(cè), 但一般受臺(tái)站位置所限, 震例積累相對(duì)有限; 而衛(wèi)星觀測(cè)可以覆蓋全球獲取更多震例資料, 但是單顆衛(wèi)星受軌道所限對(duì)同一區(qū)域的時(shí)間間隔往往在幾天以上, 期間可能錯(cuò)失很多信息. 為了更好地發(fā)揮不同探測(cè)手段的應(yīng)用效能, 圖1給出了地震電離層立體監(jiān)測(cè)體系的基本框架. 可以看出, 綜合利用現(xiàn)有探測(cè)技術(shù)即可實(shí)現(xiàn)自地表、 電離層底界面D區(qū)、 E區(qū)、 F2區(qū)至衛(wèi)星高度的全面立體監(jiān)測(cè), 并可通過不同探測(cè)手段獲得相同參量的觀測(cè), 從而改善觀測(cè)參量的時(shí)空分辨率和三維成像技術(shù), 為研究電磁擾動(dòng)激發(fā)機(jī)制以及從地表到空間的傳播耦合機(jī)理提供不同層位的基礎(chǔ)觀測(cè)結(jié)果和實(shí)際約束條件.

    關(guān)于地基電離層探測(cè), 在國內(nèi)建立了測(cè)高GPS TEC, VLF電波和舒曼諧振等相關(guān)技術(shù)觀測(cè)臺(tái)站, 并開展了地震監(jiān)測(cè)應(yīng)用研究, 也取得了一定的成效, 但尚有很多問題值得我們深思.

    針對(duì)電離層測(cè)高技術(shù), 我國在華北地區(qū)建立了電離層斜測(cè)網(wǎng), 明顯提高了該區(qū)域的電離層空間分辨率, 但電離層斜測(cè)接收的是電離層垂測(cè)站的高頻信號(hào), 通過反射到達(dá)一定距離之外的接收站后其能量相對(duì)較弱, 導(dǎo)致頻高圖上反射信號(hào)有一定程度的損失, 對(duì)頻高圖

    表1 電離層探測(cè)技術(shù)各項(xiàng)指標(biāo)對(duì)比

    圖1 地震電離層立體監(jiān)測(cè)體系基本框架

    上關(guān)鍵數(shù)據(jù)的判讀也造成較大影響, 因此斜測(cè)數(shù)據(jù)的準(zhǔn)確性和可靠性在國際上的認(rèn)可度不高. 根據(jù)目前的觀測(cè)積累, 在每年的夏季和雷雨季節(jié), 頻高圖上信號(hào)中斷特別突出, 數(shù)據(jù)缺失嚴(yán)重, 因此斜測(cè)發(fā)射接收的觀測(cè)技術(shù)本身尚需改進(jìn). 此外, 自2009年電離層斜測(cè)網(wǎng)運(yùn)行以來, 華北地區(qū)未發(fā)生M>6.0強(qiáng)震, 整個(gè)監(jiān)測(cè)網(wǎng)未有典型震例分析的經(jīng)驗(yàn)積累, 只在磁暴等空間天氣事件中與GPS TEC有相關(guān)的對(duì)比分析結(jié)果, 故其在地震中的應(yīng)用效能還需更長期的數(shù)據(jù)積累.

    TEC在國際上應(yīng)用非常廣泛, 我國的GPS數(shù)據(jù)預(yù)處理和TEC反演技術(shù)基本保持與國際同步發(fā)展的水平, 但目前的地震監(jiān)測(cè)應(yīng)用多以單站時(shí)間序列分析為主. JPL GIM模型僅使用了我國約40個(gè)GPS站點(diǎn)的數(shù)據(jù), 在我國的空間分辨率很低, 而國內(nèi)目前尚無分辨率較高的區(qū)域TEC現(xiàn)報(bào)數(shù)據(jù)開放, 因此利用國內(nèi)現(xiàn)有基礎(chǔ)資源提供高精度區(qū)域TEC現(xiàn)報(bào)數(shù)據(jù), 進(jìn)而加強(qiáng)地震電離層擾動(dòng)時(shí)空分布研究是近期需要完成的主要任務(wù)之一.

    VLF電波在國內(nèi)處于試驗(yàn)研究階段, 目前使用設(shè)備僅能接收俄羅斯阿爾法3個(gè)導(dǎo)航站的信號(hào), 發(fā)射與接收之間以電磁波的南北向傳播為主, 地震異常分析結(jié)果與日本震例應(yīng)用中較多使用的日出日落時(shí)間(terminator time, 簡寫為TT)異常幾乎無相似之處, 且鏈路距離比較長, 因此在觀測(cè)技術(shù)、 鏈路設(shè)計(jì)以及理論電波傳播計(jì)算模擬等方向亟需加強(qiáng).

    舒曼諧振在國際上較多應(yīng)用于閃電研究, 地震應(yīng)用只是其很小的應(yīng)用研究方向. 在國內(nèi)云南地區(qū)的地震監(jiān)測(cè)應(yīng)用中, 雖然一些地震震前顯示了各階舒曼諧振信號(hào)的增強(qiáng), 但目前整體強(qiáng)震積累不足, 資料說服力不強(qiáng), 加之其全球影響的特性, 單純利用諧振信號(hào)來進(jìn)行地震監(jiān)測(cè)的應(yīng)用價(jià)值并不高. 考慮到其屬于ELF頻段的磁場(chǎng)觀測(cè), 未來與其它資料配合研究變化磁場(chǎng)或地下結(jié)構(gòu)探測(cè)等有一定價(jià)值(趙國澤等, 2015).

    伴隨空間電磁探測(cè)技術(shù)的發(fā)展和探測(cè)水平的不斷提高, 以及關(guān)于大量震前顯著電離層短臨擾動(dòng)特性的深入研究, 地震電離層探測(cè)受到世界各國的廣泛關(guān)注, 尤其是中國電磁衛(wèi)星計(jì)劃的確立為構(gòu)建地震電磁三維立體監(jiān)測(cè)體系提供了一條無可替代的空間探測(cè)技術(shù)途徑, 使我國的地震立體監(jiān)測(cè)計(jì)劃進(jìn)一步向電離層頂部拓展.

    張學(xué)民等(2009c)基于法國DEMETER衛(wèi)星開展的震例研究結(jié)果進(jìn)一步證實(shí)了地震電離層現(xiàn)象的存在及其顯著的短臨特性, 但不同衛(wèi)星、 不同載荷、 不同參量之間的相關(guān)性分析在國內(nèi)外的研究程度明顯不足, 其背后的理論支撐亟需加快發(fā)展. 相對(duì)DEMETER衛(wèi)星而言, 在保證多數(shù)載荷相同的基礎(chǔ)上, 中國電磁監(jiān)測(cè)試驗(yàn)衛(wèi)星專門增加了電離層結(jié)構(gòu)探測(cè)GNSS和三頻信標(biāo)兩個(gè)載荷, 對(duì)電離層的系統(tǒng)觀測(cè)在原位探測(cè)基礎(chǔ)上, 增加了對(duì)衛(wèi)星下方電離層結(jié)構(gòu)的探測(cè), 這樣整個(gè)觀測(cè)系統(tǒng)相對(duì)更為完善, 同時(shí)我國地基電磁、 電離層探測(cè)手段的不斷加強(qiáng)和完善也為天地一體化數(shù)據(jù)校驗(yàn)及地震電離層信息立體監(jiān)測(cè)提供了更多的途徑. 法國DEMETER衛(wèi)星運(yùn)行期間, 法國國內(nèi)地面也相應(yīng)建立了一些電磁觀測(cè)站, 主要使用ULF, ELF和VLF頻段電磁場(chǎng)觀測(cè), 在活動(dòng)斷層和火山活動(dòng)監(jiān)測(cè)中予以驗(yàn)證(Zlotnickietal, 2006). 我國不僅在電磁觀測(cè)手段方面比較完備, 同時(shí)多種地基電離層等離子體參量觀測(cè)手段也基本齊備, 因此基于中國地震電磁衛(wèi)星項(xiàng)目, 我國的地震電磁立體觀測(cè)體系相比其它國家更系統(tǒng), 也更完整.

    地震是發(fā)生在地下深部的構(gòu)造活動(dòng), 在其孕育過程中由于能量的大量積累激發(fā)產(chǎn)生相關(guān)的形變、 電磁輻射以及地下流體、 化學(xué)物質(zhì)的放射等, 這些信息在地表積累進(jìn)而傳播影響到電離層衛(wèi)星觀測(cè), 結(jié)合國內(nèi)外研究人員提出的各類模型綜合分析認(rèn)為, 其傳播耦合途徑主要有3種: 電磁波傳播、 聲重力波傳播以及地球化學(xué)等因素激發(fā)的直流電場(chǎng)耦合. 從理論驗(yàn)證的角度而言, 立體監(jiān)測(cè)將為不同的模型校驗(yàn)提供更多的約束資料, 從而推動(dòng)理論研究的深化. 在電磁波傳播的立體監(jiān)測(cè)體系中, 可在地面利用天然源場(chǎng)或者人工源場(chǎng)等探測(cè)地下電性結(jié)構(gòu)以及與地震可能相關(guān)的異常變化, 又可以作為波導(dǎo)傳播模型中地下介質(zhì)設(shè)定的重要參數(shù), 直接觀測(cè)地表人工源及可能來自孕震區(qū)的天然源電磁輻射; 加上衛(wèi)星觀測(cè)電磁場(chǎng)可以提供在電離層高度的電磁波傳播特性, 因此寬頻電磁場(chǎng)觀測(cè)可以實(shí)現(xiàn)從巖石層到電離層的有效覆蓋, 并通過地面與空間相同頻率范圍的觀測(cè)結(jié)果進(jìn)行對(duì)比研究.

    地基和空間觀測(cè)缺一不可, 其中人工源電波作為其中電磁波的特例可以為發(fā)展完善傳播模型提供最好的范本; 配合電離層等離子體參量觀測(cè)(測(cè)高、 TEC、 掩星、 衛(wèi)星原位觀測(cè)等), 可以研究電磁波傳播在電離層不均勻體等結(jié)構(gòu)下的擾動(dòng), 形成一套比較完善的監(jiān)測(cè)體系. 就目前而言, 由于甚低頻段電磁波的反射、 折射能量損失主要發(fā)生在電離層底界面附近, 而目前國內(nèi)使用的地震電離層探測(cè)技術(shù)在此界面嚴(yán)重不足, 這也是監(jiān)測(cè)體系中缺失的重要環(huán)節(jié). 在該問題上, 三頻信標(biāo)和掩星觀測(cè)可以發(fā)揮更大的作用, 同時(shí)可考慮利用雷達(dá)、 空間站等彌補(bǔ)低電離層觀測(cè)的缺乏.

    聲重力波傳播模型目前較多應(yīng)用于同震、 震后以及海嘯、 臺(tái)風(fēng)等大氣波動(dòng)事件引起的電離層擾動(dòng), Hao等(2012)的研究顯示了地震波、 地磁場(chǎng)、 GPS TEC、 多普勒頻移、 次聲波等觀測(cè)技術(shù)在日本2011年3月11日大地震后的聯(lián)動(dòng)效應(yīng). 聲重力波傳播偏重于大氣層的機(jī)械波傳播, 其到達(dá)電離層底界面后與電離層電子、 離子的響應(yīng)運(yùn)移機(jī)制目前尚未有較好的機(jī)理模型解釋; 國外研究人員在解釋人工源電波在震前電離層信號(hào)強(qiáng)度降低的現(xiàn)象時(shí)也認(rèn)為震前聲重力波擾動(dòng)電離層引起電波強(qiáng)度改變(Molchanovetal, 2006), 但截至目前未有同震或震后與GPS TEC等觀測(cè)同步的人工源電波擾動(dòng)現(xiàn)象發(fā)布, 兩者之間的關(guān)聯(lián)尚待在立體監(jiān)測(cè)體系中進(jìn)一步驗(yàn)證.

    直流電場(chǎng)耦合模型最直接的應(yīng)用就是GPS TEC的地震擾動(dòng)異常及其共軛效應(yīng), 不同高度的等離子體參量響應(yīng)也可以對(duì)直流電場(chǎng)激發(fā)的電離層不均勻體特性進(jìn)行約束, 但目前該模型僅利用電磁觀測(cè)在地表和大氣層的校驗(yàn), 所以其約束性還比較弱. 為了驗(yàn)證該模型, 需要地基、 大氣層、 電離層化學(xué)物質(zhì)的探測(cè), 地表電荷的積累, 大氣電場(chǎng)監(jiān)測(cè), 大氣電導(dǎo)率結(jié)構(gòu)探測(cè)等結(jié)果, 目前的探測(cè)技術(shù)中已有地表殼幔源氣體, 衛(wèi)星高光譜氣體探測(cè), 電磁探測(cè)手段中掩星觀測(cè)可以提供大氣層中氣壓、 溫度等剖面結(jié)構(gòu), 因此需要充分利用現(xiàn)有的地球化學(xué)探測(cè)資源, 配合大氣電場(chǎng)等電磁觀測(cè), 研究地表-大氣層多參量的耦合響應(yīng), 為直流電場(chǎng)模型提供更多支撐.

    4 討論與展望

    地震預(yù)測(cè)仍然是世界性難題, 針對(duì)地震監(jiān)測(cè)應(yīng)用的每一種探測(cè)技術(shù)除了其獨(dú)特的優(yōu)勢(shì)外, 均存在片面性和局限性. 鑒于電磁尤其是地震電離層監(jiān)測(cè)領(lǐng)域的顯著短臨擾動(dòng)特性, 國內(nèi)外科研人員提出了地震立體監(jiān)測(cè)體系的發(fā)展理念(Bleier, Freund, 2005; 趙國澤等, 2007, 2015; 張學(xué)民等, 2009a). 本文基于國內(nèi)外現(xiàn)有的電離層探測(cè)技術(shù)及其在國內(nèi)的主要應(yīng)用進(jìn)展, 分析了各類探測(cè)技術(shù)的優(yōu)勢(shì)和不足, 并與地震電離層耦合機(jī)理理論相結(jié)合, 希望通過立體監(jiān)測(cè)體系深化不同探測(cè)參量之間的關(guān)聯(lián), 驗(yàn)證完善地震電離層耦合理論模型, 以數(shù)據(jù)發(fā)展校驗(yàn)?zāi)P停?以理論帶動(dòng)技術(shù)的發(fā)展思路, 最終實(shí)現(xiàn)理論模型下的多參量有機(jī)串聯(lián)和組合, 在地震電離層立體監(jiān)測(cè)體系的框架下, 為地基和衛(wèi)星地震電離層立體監(jiān)測(cè)提供理論指導(dǎo).

    根據(jù)目前國內(nèi)電離層探測(cè)技術(shù)及其應(yīng)用研究進(jìn)展, 建議進(jìn)一步開展以下幾方面的研究:

    1) 完善觀測(cè)體系, 主要包括地基現(xiàn)有電離層探測(cè)技術(shù)的改進(jìn)、 反演算法的優(yōu)化、 多頻段電磁場(chǎng)觀測(cè)的加強(qiáng)以及新技術(shù)的引進(jìn)等.

    2) 校驗(yàn)圈層耦合機(jī)理, 通過在地震前后及其它災(zāi)害事件中的天地一體化電磁場(chǎng)以及不同高度等離子體參量的同步響應(yīng), 利用其時(shí)空強(qiáng)關(guān)聯(lián)特性校驗(yàn)圈層耦合模型.

    3) 改進(jìn)地球物理場(chǎng)精細(xì)建模技術(shù), 包括電離層現(xiàn)報(bào)系統(tǒng)、 電離層層析成像以及地磁場(chǎng)建模等, 為地震監(jiān)測(cè)應(yīng)用分析提供精細(xì)電離層結(jié)構(gòu)和基礎(chǔ)背景支撐.

    4) 發(fā)展地震電磁預(yù)測(cè)模型. 挖掘多源多參量數(shù)據(jù)中的相關(guān)性, 初步搭建地震電離層立體監(jiān)測(cè)體系下的地震電磁短臨預(yù)測(cè)模型, 充分發(fā)揮立體監(jiān)測(cè)效能.

    安張輝, 杜學(xué)彬, 范瑩瑩, 劉君, 譚大誠, 陳軍營, 解滔. 2011. 汶川MS8.0級(jí)大震前天基與陸基電場(chǎng)資料聯(lián)合應(yīng)用研究[J]. 地球物理學(xué)報(bào), 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 8.0 earthquake of 2008 using both space-based and ground-based observational data[J].ChineseJournalofGeophysics, 54(11): 2876--2884 (in Chinese).

    蔡華, 孫漢榮, 李子申. 2014. 利用陸態(tài)網(wǎng)絡(luò)實(shí)時(shí)監(jiān)測(cè)電離層TEC及在地震監(jiān)測(cè)中的應(yīng)用分析[J]. 地震, 34(1): 41--48.

    Cai H, Sun H R, Li Z S. 2014. Real-time ionospheric TEC monitoring based on CMONOC and its application[J].Earthquake, 34(1): 41--48 (in Chinese).

    蔡軍濤, 趙國澤, 詹艷, 湯吉,陳小斌. 2007. 地震期間電離層擾動(dòng)現(xiàn)象研究[J]. 地球物理學(xué)進(jìn)展, 22(3): 695--701.

    Cai J T, Zhao G Z, Zhan Y, Tang J, Chen X B. 2007. The study on ionospheric disturbances during earthquakes[J].ProgressinGeophysics, 22(3): 695--701 (in Chinese).

    丁鑒海, 索玉成, 余素榮, 肖武軍. 2004. 昆侖山口西8.1級(jí)地震前電離層與磁場(chǎng)短期異常對(duì)比研究[J]. 地震, 23(1): 104--111.

    Ding J H, Suo Y C, Yu S R, Xiao W J. 2004. Comparison study on short-term anomalies of ionosphere and geomagnetic field before KunlunshanMS8.1 earthquake[J].Earthquake, 23(1): 104--111 (in Chinese).

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

    Ding J H, Suo Y C, Yu S R. 2005. Phenomena of geomagnetic and ionospheric anomalies and their relation to earthquakes[J].ChineseJournalofSpaceScience, 25(6): 536--542 (in Chinese).

    范曄, 湯吉, 趙國澤, 王立鳳, 吳江星, 李細(xì)順, 黃天保, 劉廣寬. 2013. 電磁監(jiān)測(cè)臺(tái)站觀測(cè)的舒曼諧振背景變化[J]. 地球物理學(xué)報(bào), 56(7): 2369--2377.

    Fan Y, Tang J, Zhao G Z, Wang L F, Wu J X, Li X S, Huang T B, Liu G K. 2013. Schumann resonances variation observed from electromagnetic monitoring stations[J].ChineseJournalofGeophysics, 56(7): 2369--2377 (in Chinese).

    甘衛(wèi)軍, 張銳, 張勇, 唐方頭. 2007. 中國地殼運(yùn)動(dòng)觀測(cè)網(wǎng)絡(luò)的建設(shè)及應(yīng)用[J]. 國際地震動(dòng)態(tài), (7): 32--37.

    Gan W J, Zhang R, Zhang Y, Tang F T. 2007. Development of the crustal movement observation network of China and its application[J].RecentDevelopmentsinWorldSeismology, (7): 32--37 (in Chinese).

    耿光宇. 2010. 利用COSMIC掩星數(shù)據(jù)研究震前電離層異常[J]. 桂林電子科技大學(xué)學(xué)報(bào), 30(2): 154--157.

    Geng G Y. 2010. Using COSMIC occultation data to evaluate seismo-ionospheric anomalies[J].JournalofGuilinUniversityofElectronicTechnology, 30(2): 154--157 (in Chinese).

    何宇飛, 楊冬梅, 陳化然, 錢家棟, 朱榮, Parrot M. 2009. DEMETER衛(wèi)星探測(cè)到可能與汶川地震有關(guān)的地面VLF發(fā)射站信號(hào)的信噪比變化[J]. 中國科學(xué): 地球科學(xué), 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].ScienceChinaEarthSciences, 52(6): 754--763.

    李建勇. 2010. 利用GPS數(shù)據(jù)研究地震電離層TEC異常[D]. 北京: 中國地震局地震預(yù)測(cè)研究所: 1--107.

    Li J Y. 2010.ResearchonSeismicIonosphericAnomaliesBasedonGPSData[D]. Beijing: Institute of Earthquake Science, China Earthquake Administration: 1--107 (in Chinese).

    梁百先, 李鈞, 馬淑英. 1994. 我國的電離層研究[J]. 地球物理學(xué)報(bào), 37(S1): 51--73.

    Liang B X, Li J, Ma S Y. 1994. Progress of ionospheric research in China[J].ChineseJournalofGeophysics, 37(S1): 51--73 (in Chinese).

    林劍, 吳云, 祝芙英, 喬學(xué)軍, 周義炎. 2009. 基于GPS探測(cè)汶川地震電離層TEC的異常[J]. 地球物理學(xué)報(bào), 52(1): 297--300.

    Lin J, Wu Y, Zhu F Y, Qiao X J, Zhou Y Y. 2009. Wenchuan earthquake ionosphere TEC anomaly detected by GPS[J].ChineseJournalofGeophysics, 52(1): 297--300 (in Chinese).

    劉靜, 萬衛(wèi)星, 黃建平, 張學(xué)民, 趙庶凡, 歐陽新艷, 澤仁志瑪. 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).

    劉經(jīng)南, 趙瑩, 張小紅. 2010. GNSS無線電掩星電離層反演技術(shù)現(xiàn)狀與展望[J]. 武漢大學(xué)學(xué)報(bào): 信息科學(xué)版, 35(6): 631--635.

    Liu J N, Zhao Y, Zhang X H. 2010. Current situation and expectation of inversion of ionospheric GNSS occultation[J].GeomaticsandInformationScienceofWuhanUniversity, 35(6): 631--635 (in Chinese).

    苗園青, 張學(xué)民, 申旭輝, 歐陽新艷. 2011. 云南地區(qū)震前舒曼諧振異常初步分析[J]. 中國地震, 27(3): 290--299.

    Miao Y Q, Zhang X M, Shen X H, Ouyang X Y. 2011. Initial results on anomalous Schumann resonance before earthquakes in Yunnan[J].EarthquakeResearchinChina, 27(3): 290--299 (in Chinese).

    歐陽新艷, 張學(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).

    歐陽新艷, 婁文宇, 申旭輝, 張學(xué)民. 2015. 中國大陸區(qū)域6及以上地震舒曼諧振異常特征分析[J]. 地震, 35(1): 123--130.

    Ouyang X Y, Lou W Y, Shen X H, Zhang X M. 2015. Extraction and analysis of Schumann resonance anomaly features based on earthquakes larger thanM6.0 in mainland China[J].Earthquake, 35(1): 123--130 (in Chinese).

    王也英, 符養(yǎng), 杜曉勇, 薛震剛. 2009. 全球GNSS掩星計(jì)劃進(jìn)展[J]. 氣象科技, 37(1): 74--78.

    Wang Y Y, Fu Y, Du X Y, Xue Z G. 2009. Advances in global GNSS occultation projects[J].MeteorologicalScienceandTechnology, 37(1): 74--78 (in Chinese).

    夏淳亮. 2004. GPS臺(tái)網(wǎng)觀測(cè)中電離層TEC的解算方法及TEC現(xiàn)報(bào)系統(tǒng)的初步研制[D]. 武漢: 中國科學(xué)院武漢物理與數(shù)學(xué)研究所: 1--61.

    Xia C L. 2004.TECIonosphericComputingTechniqueandInitialDevelopmentofTECNowcastingSysteminGPSObservationalNetwork[D]. Wuhan: Wuhan Institute of Physics and Mathematics, Chinese Academy of Sciences: 1--61 (in Chinese).

    肖佐. 1997. 近年來中國電離層物理研究進(jìn)展[J]. 地球物理學(xué)報(bào), 40(S1): 21--28.

    Xiao Z. 1997. Research on ionospheric physics in recent years in China[J].ChineseJournalofGeophysics, 40(S1): 21--28 (in Chinese).

    徐彤, 胡艷莉, 吳健, 李春斌, 吳振森, 索玉成, 馮健. 2012. 中國大陸14次強(qiáng)震前電離層異常統(tǒng)計(jì)分析[J]. 電波科學(xué)學(xué)報(bào), 27(3): 507--512.

    Xu T, Hu Y L, Wu J, Li C B, Wu Z S, Suo Y C, Feng J. 2012. Statistical analysis of seismo-ionospheric perturbation before 14MS≥7.0 strong earthquakes in Chinese subcontinent[J].ChineseJournalofRadioScience, 27(3): 507--512 (in Chinese).

    楊劍, 吳云, 周義炎. 2008. 利用GPS無線電掩星數(shù)據(jù)研究震前電離層異常[J]. 大地測(cè)量與地球動(dòng)力學(xué), 28(1): 16--22.

    Yang J, Wu Y, Zhou Y Y. 2008. Research on seismo-ionospheric anomalies using GPS radio occultation data[J].JournalofGeodesyandGeodynamics, 28(1): 16--22 (in Chinese).

    姚璐, 申旭輝, 張學(xué)民. 2014. 玉樹MS7.1地震前電離層異常擾動(dòng)分析[J]. 地震, 34(3): 74--85.

    Yao L, Shen X H, Zhang X M. 2014. Analysis of ionospheric anomalies preceding the 2010 YushuMS7.1 earthquake[J].Earthquake, 34(3): 74--85 (in Chinese).

    于海雁, 周洪娟, 喬曉林. 2010. 汶川MS8.0地震前ELF異常電磁輻射傳播分析[J]. 地震學(xué)報(bào), 32(6): 641--648.

    Yu H Y, Zhou H J, Qiao X L. 2010. Study on wave propagation of ELF emission anomaly beforeMS8.0 Wenchuan earthquake[J].ActaSeismologicaSinica, 32(6): 641--648 (in Chinese).

    余素榮, 丁鑒海, 索玉成. 2004. 瑪尼7.5級(jí)地震前磁場(chǎng)與電離層異常特征研究[J]. 地震, 24(增刊): 77--83.

    Yu S R, Ding J H, Suo Y C. 2004. Research on anomaly characteristics of geomagnetic field and ionosphere before ManiMS7.5 earthquake[J].Earthquake, 24(Suppl): 77--83 (in Chinese).

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

    Zeren 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).

    張學(xué)民, 劉靜, 錢家棟, 申旭輝, 蔡晉安, 歐陽新艷, 趙庶凡. 2008. 西藏改則6.9級(jí)地震前的電離層電磁擾動(dòng)[J]. 地震, 28(3): 14--22.

    Zhang X M, Liu J, Qian J D, Shen X H, Cai J A, Ouyang X Y, Zhao S F. 2008. Ionospheric electromagnetic disturbance before Gaize earthquake withMS6.9, Tibet[J].Earthquake, 28(3): 14--22 (in Chinese).

    張學(xué)民, 丁鑒海, 申旭輝, 王敏, 劉靜, 余素榮, 王亞麗, 歐陽新艷. 2009a. 汶川8級(jí)地震前電磁擾動(dòng)與電磁立體監(jiān)測(cè)體系[J]. 電波科學(xué)學(xué)報(bào), 24(1): 1--8.

    Zhang X M, Ding J H, Shen X H, Wang M, Liu J, Yu S R, Wang Y L, Ouyang X Y. 2009a. Electromagnetic perturbations before WenchuanM8 earthquake and stereo electromagnetic observation system[J].ChineseJournalofRadioScience, 24(1): 1--8 (in Chinese).

    張學(xué)民, 申旭輝, 歐陽新艷, 蔡晉安, 黃建平, 劉靜, 趙庶凡. 2009b. 汶川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. 2009b. Ionospheric VLF electric field anomalies before WenchuanM8 earthquake[J].ChineseJournalofRadioScience, 24(6): 1024--1032 (in Chinese).

    張學(xué)民, 申旭輝, 錢家棟, 路立, 陳化然, 歐陽新艷. 2009c. 我國地震電磁衛(wèi)星數(shù)據(jù)分析及應(yīng)用研究進(jìn)展[J]. 地震, 29(Z1): 34--45.

    Zhang X M, Shen X H, Qian J D, Lu L, Chen H R, Ouyang X Y. 2009c. Advances in the analysis and application of seismo-electromagnetic satellite data in China[J].Earthquake, 29(Z1): 34--45 (in Chinese).

    張學(xué)民, 劉靜, 趙必強(qiáng), 徐彤, 申旭輝, 姚璐. 2014. 玉樹地震前的電離層異常現(xiàn)象分析[J]. 空間科學(xué)學(xué)報(bào), 34(6): 822--829.

    Zhang X M, Liu J, Zhao B Q, Xu T, Shen X H, Yao L. 2014. Analysis on ionospheric perturbations before Yushu earthquake[J].ChineseJournalofSpaceScience, 34(6): 822--829 (in Chinese).

    趙國澤, 陳小斌, 蔡軍濤. 2007. 電磁衛(wèi)星和地震預(yù)測(cè)[J]. 地球物理學(xué)進(jìn)展, 22(3): 667--673.

    Zhao G Z, Chen X B, Cai J T. 2007. Electromagnetic observation by satellite and earthquake prediction[J].ProgressinGeophysics, 22(3): 667--673 (in Chinese).

    趙國澤, Bi Y X, 王立鳳, 韓冰, 汪曉, 肖騎彬, 蔡軍濤, 詹艷, 陳小斌, 湯吉, 王繼軍. 2015. 中國地震交變電磁場(chǎng)觀測(cè)數(shù)據(jù)處理技術(shù)新進(jìn)展[J]. 中國科學(xué): 地球科學(xué), 45(1): 22--33.

    Zhao G Z, Bi Y X, Wang L F, Han B, Wang X, Xiao Q B, Cai J T, Zhan Y, Chen X B, Tang J, Wang J J. 2015. Advances in alternating electromagnetic field data processing for earthquake monitoring in China[J].ScienceChinaEarthSciences, 58(2): 172--182.

    中國科學(xué)院地質(zhì)與地球物理研究所. 2015. DPS4D digisonde[EB/OL].[2015-11-25]. http:∥space.iggcas.ac.cn/SYIonogram.html.

    Institute of Geology and Geophysics, Chinese Academy of Sciences. 2015. DPS4D digisonde[EB/OL].[2015-11-25]. http:∥space.iggcas.ac.cn/SYIonogram.html (in Chinese).

    周洪娟, 于海雁, 曹丙霞, 喬曉林. 2013. 日本MW9.0級(jí)地震前的舒曼諧振異常分析[J]. 地震學(xué)報(bào), 35(3): 400--409.

    Zhou H J, Yu H Y, Cao B X, Qiao X L. 2013. Schumann resonance anomaly before theMW9.0 Japan earthquake[J].ActaSeismologicaSinica, 35(3): 400--409 (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).

    Austen J R, Franke S J, Liu C H. 1988. Ionospheric imaging using computerized tomography[J].RadioSci, 23(3): 299--307.

    Balser M, Wagner C. 1960. Observations of earth-ionosphere cavity resonances[J].Nature, 188(4751): 638--641.

    Berthelier J J, Godefroy M, Leblanc F, Seran E, Peschard D, Gilbert P, Artru J. 2006. IAP, the thermal plasma analyzer on DEMETER[J].PlanetSpaceSci, 54(5): 487--501.

    Bleier T, Freund F. 2005. Impending earthquakes have been sending US warning signals and people are starting to listen[J].IEEESpectrumINT, 3: 3--7.

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

    Gousheva M, Danov D, Hristov P, Matova M. 2008. Quasi-static electric fields phenomena in the ionosphere associated with pre- and post earthquake effects[J].NatHazardsEarthSystSci, 8(1): 101--107.

    Gokhberg M B, Pilipenko V A, Pokhotelov O A. 1983. Observation from a satellite of electromagnetic radiation above the region of an earthquake in preparation[J].DokladyoftheUSSRAcademyofSciences:SciencesSections, 268(1): 5--7 (in Russian).

    Hao Y Q, Xiao Z, Zhang D H. 2012. Multi-instrument observation on co-seismic ionospheric effects after great Tohoku earthquake[J].JGeophysRes, 117(A2): A02305. doi:10.1029/2011JA017036.

    Hayakawa M, Ohta K, Nickolaenko A P, Ando Y. 2005. Anomalous effect in Schumann resonance phenomena observed in Japan, possibly associated with the Chi-Chi earthquake in Taiwan[J].AnnGeophys, 23(4): 1335--1346.

    Hayakawa M, Nickolaenko A P, Sekiguchi M, Yamshita K, Ida Y, Yano M. 2008. Anomalous ELF phenomena in the Schumann resonance band as observed at Moshiri (Japan) in possible association with an earthquake in Taiwan[J].NatHazardsEarthSystSci, 8(6): 1309--1316.

    Hayakawa M, Kasahara Y, Nakamura T, Muto F, Horie T, Maekawa S, Hobara Y, Rozhnoi A A, Solovieva M, Molchanov O A. 2010. A statistical study on the correlation between lower ionospheric perturbations as seen by subionospheric VLF/LF propagation and earthquakes[J].JGeophysRes, 115(A9): A09305.

    Hsiao C, Liu J Y, Oyama K I, Yen N L, Liou Y A, Chen S S, Miau J J. 2010. Seismo-ionospheric precursor of the 2008MW7.9 Wenchuan earthquake observed by FORMOSAT-3/COSMIC[J].GPSSolut, 14(1): 83--89.

    Kakinami Y, Liu J Y, Tsai L C, Oyama K I. 2010. Ionospheric electron content anomalies detected by a FORMOSAT-3/COSMIC empirical model before and after the Wenchuan earthquake[J].IntJRemoteSens, 31(13): 3571--3578.

    Lebreton J P, Stverak S, Travnicek P, Mksimovic M, Klinge D, Merikallio S, Lagoutte D, Poirier B, Blelly P L, Kozacek Z, Salaquarda M. 2006. The ISL Langmuir probe experiment and its data procession onboard DEMETER: Scientific objectives, description and first results[J].PlanetSpaceSci, 54(5): 472--486.

    Lei J H, Syndergaard S, Burns A G, Solomon S C, Wang W B, Zeng J, Roble R G, Wu Q, Kuo Y H, Holt J M, Zhang S R, Hysell D L, Rodrigues F S, Lin C H. 2007. Comparison of COSMIC ionospheric measurements with ground-based observations and model predictions: Preliminary results[J].JGeophysRes, 112(A7): A07308.

    Li M, Parrot M. 2013. Statistical analysis of an ionospheric parameter as a base for earthquake prediction[J].JGeophysRes, 118(6): 3731--3739.

    Liu J, Huang J P, Zhang X M. 2014. Ionospheric perturbations in plasma parameters before global strong earthquakes[J].AdvSpaceRes, 53(5): 776--787.

    Liu J Y, Chen C H, Tsai H F, Le H. 2013. A statistical study on seismo-ionospheric anomalies of the total electron content for the period of 56M≥6.0 earthquakes occurring in China during 1998--2012[J].ChinJSpaceSci, 33(3): 258--269.

    Molchanov O, Rozhnoi A, Solovieva M, Akentieva O, Berthelier J J, Parrot M, Lefeuvre F, Biagi P F, Castellana L, Hayakawa M. 2006. Global diagnostics of the ionospheric perturbations related to the seismic activity using the VLF radio signals collected on the DEMETER satellite[J].NatHazardsEarthSystSci, 6(5): 745--753.

    Němec F, Santolík O, Parrot M. 2009. Decrease of intensity of ELF/VLF waves observed in the upper ionosphere close to earthquakes: A statistical study[J].JGeophysRes, 114(A4): A04303. doi:10.1029/2008JA013972.

    Ohta K, Watanabe N, Hayakawa M. 2006. Survey of anomalous Schumann resonance phenomena observed in Japan, in possible association with earthquakes in Taiwan[J].PhysChemEarth, 31(4/5/6/7/8/9): 397--402.

    Parrot M, Berthelier J J, Lebreton J P, Sauvaud J A, Santolik O, Blecki J. 2006. Examples of unusual ionospheric observations made by the DEMETER satellite over seismic regions[J].PhysChemEarth, 31(4/5/6/7/8/9): 486--495.

    Pí?a D, Němec F, Santolík O, Parrot M, Rycroft M. 2013. Additional attenuation of natural VLF electromagnetic waves observed by the DEMETER spacecraft resulting from preseismic activity[J].JGeophysRes, 118(8): 5286--5295.

    Pulinets S A, Boyarchuk K. 2004.IonosphericPrecursorsofEarthquakes[M]. Berlin: Springer: 75--169.

    Pulinets S A, Legen′ka A D, Gaivoronskaya T V, Depuev V K. 2003. Main phenomenological features of ionospheric precursors of strong earthquakes[J].JAtmosSol-TerrPhy, 65(16/17/18): 1337--1347.

    Schreiner W S, Sokolovskiy S V, Rocken C, Hunt D C. 1999. Analysis and validation of GPS/MET radio occultation data in the ionosphere[J].RadioSci, 34(4): 949--966.

    Schumann W O. 1952. über die strahlungslosen Eigenschwingungen einer leitenden Kugel, die von einer Luftschicht und einer Ionosph?renhülle umgeben ist[J].ZeitschriftundNaturfirschung, 7a: 149--154.

    Song Q, Ding F, Yu T, Wan W X, Ning B Q, Liu L B, Zhao B Q. 2015. GPS detection of the coseismic ionospheric disturbances following the 12 May 2008M7.9 Wenchuan earthquake in China[J].ScienceChinaEarthSciences, 58(1): 151--158. doi: 10.1007/s11430-014-5000-7.

    Wikipedia. 2015. List_of_VLF-transmitters[EB/OL]. [2015-11-25]. https:∥en.wikipedia.org/wiki/List_of_VLF-transmitters.

    Zeren Z M, Zhang X M, Shen X H, Sun W H, Ning D M, Ruzhin Y. 2014. VLF radio signal anomalies associated with strong earthquakes[C]∥Proceedingsofthe2014XXXIthURSIIEEEGeneralAssemblyandScientificSymposium. Beijing: IEEE: 1--4.

    Zhang X M, Shen X H, Miao Y Q. 2012a. Electromagnetic anomalies around Wenchuan earthquake and their relationship with earthquake preparation[J].ProcediaEnvironSci, 12: 693--701.

    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, Fidani C, Huang J P, Shen X H, Zeren Z M, Qian J D. 2013. Burst increases of precipitating electrons recorded by the DEMETER satellite before strong earthquakes[J].NatHazardsEarthSystSci, 13(1): 197--209.

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

    Zhao B Q, Yu T, Wang M, Wan W X, Lei J H, Liu L B, Ning B Q. 2008. Is an unusual large enhancement of ionospheric electron density linked with the 2008 great Wenchuan earthquake?[J].JGeophysRes, 113(A11): A11304.

    Zhou H J, Zhou Z Q, Qiao X L, Yu H Y. 2013. Anomalous phenomena in Schumann resonance band observed in China before the 2011 magnitude 9.0 Tohoku-Oki earthquake in Japan[J].JGeophysRes, 118(23): 13338--13345.

    Zlotnicki J, Le Mou?l J L, Kanwar R, Yvetot P, Vargemezis G, Menny P, Fauquet F. 2006. Ground-based electro-magnetic studies combined with remote sensing based on Demeter mission: A way to monitor active faults and volcanoes[J].PlanetSpaceSci, 54(5): 541--557.

    The seismo-ionospheric monitoring technologies and their application research development

    1)InstituteofEarthquakeScience,ChinaEarthquakeAdministration,Beijing100036,China2)InstituteofCrustalDynamics,ChinaEarthquakeAdministration,Beijing100085,China

    China Seismo-Electromagnetic Satellite (CSES) will be launched in 2017, which will fill in the gap in the stereo seismo-ionospheric monitoring system, and provide the scientific support for integrated space-ground comparison and validation in seismo-electromagnetic field, and for further understanding geosphere coupling mechanism. In order to meet the development requirements in stereo seismo-ionospheric monitoring system and fully understand each tool, this paper introduces the major ground-based and space-borne ionospheric monitoring technologies in China, including vertical/slant ionosounding, VLF (very low frequency) electromagnetic wave observation, Shumann resonance detection, GPS TEC (total electron content), and space electromagnetic satellite. Meanwhile, the application achievements in earthquake research are summarized, especially the developing status in China. Finally based on the integration of advantages for each technology, the basic framework of stereo electromagnetic monitoring system is discussed, and future comprehensive research fields are also suggested based on multi-source observation data from this system.

    vertical/slant ionosounde; VLF electromagnetic wave; Shumann resonance; GPS TEC; electromagnetic satellite

    綜 述

    中國地震局地震預(yù)測(cè)研究所基本科研業(yè)務(wù)經(jīng)費(fèi)專項(xiàng)(2015IES0101)和亞太空間合作組織國際科技合作項(xiàng)目“亞太地區(qū)地基電離層監(jiān)測(cè)及地震擾動(dòng)特征研究”(APSCO-SP/PM-EARTHQUAKE)共同資助.

    2015-11-17收到初稿, 2016-02-23決定采用修改稿.

    e-mail: zhangxm96@126.com

    10.11939/jass.2016.03.004

    P352.7

    A

    張學(xué)民, 申旭輝, 趙庶凡, 劉靜, 歐陽新艷, 婁文宇, 澤仁志瑪, 何建輝, 錢庚. 2016. 地震電離層探測(cè)技術(shù)及其應(yīng)用研究進(jìn)展. 地震學(xué)報(bào), 38(3): 356--375. doi:10.11939/jass.2016.03.004.

    Zhang X M, Shen X H, Zhao S F, Liu J, Ouyang X Y, Lou W Y, Zeren Z M, He J H, Qian G. 2016. The seismo-ionospheric monitoring technologies and their application research development.ActaSeismologicaSinica, 38(3): 356--375. doi:10.11939/jass.2016.03.004.

    啦啦啦 在线观看视频| 在线观看一区二区三区激情| 亚洲免费av在线视频| 19禁男女啪啪无遮挡网站| 久久国产精品影院| 欧美日韩福利视频一区二区| 亚洲国产看品久久| 欧美成人午夜精品| 9色porny在线观看| 成人手机av| 成年人免费黄色播放视频| 极品少妇高潮喷水抽搐| 一本大道久久a久久精品| 电影成人av| www.自偷自拍.com| 性色av一级| 欧美 亚洲 国产 日韩一| 后天国语完整版免费观看| 亚洲成色77777| 18禁观看日本| 美女国产高潮福利片在线看| 国产精品 欧美亚洲| 十八禁人妻一区二区| 看免费成人av毛片| 亚洲欧美色中文字幕在线| 亚洲一区中文字幕在线| 精品久久蜜臀av无| 18在线观看网站| 女人久久www免费人成看片| 欧美精品亚洲一区二区| 99国产精品99久久久久| 亚洲天堂av无毛| 欧美精品一区二区大全| 国产亚洲av片在线观看秒播厂| 老司机深夜福利视频在线观看 | h视频一区二区三区| 美女扒开内裤让男人捅视频| 中国美女看黄片| 午夜免费观看性视频| 精品高清国产在线一区| 久久精品亚洲av国产电影网| 国产在线视频一区二区| 亚洲精品一区蜜桃| 亚洲精品国产av成人精品| 高清欧美精品videossex| 欧美少妇被猛烈插入视频| 交换朋友夫妻互换小说| av视频免费观看在线观看| 丁香六月天网| 欧美精品亚洲一区二区| 国产伦理片在线播放av一区| 热re99久久精品国产66热6| 精品少妇久久久久久888优播| 黄色片一级片一级黄色片| 国产精品久久久久久人妻精品电影 | 欧美+亚洲+日韩+国产| 欧美 亚洲 国产 日韩一| 手机成人av网站| 精品福利观看| 久久人人97超碰香蕉20202| 国产精品麻豆人妻色哟哟久久| 男人添女人高潮全过程视频| 黄片小视频在线播放| 九草在线视频观看| 另类亚洲欧美激情| 男男h啪啪无遮挡| 一区二区三区乱码不卡18| 亚洲三区欧美一区| 乱人伦中国视频| 欧美日韩av久久| 午夜两性在线视频| 欧美 日韩 精品 国产| 欧美+亚洲+日韩+国产| 成年动漫av网址| 亚洲,欧美,日韩| 一级毛片黄色毛片免费观看视频| 欧美激情高清一区二区三区| 久久性视频一级片| 亚洲欧美成人综合另类久久久| 国产精品三级大全| 女人久久www免费人成看片| 老汉色av国产亚洲站长工具| 成年动漫av网址| 国产伦理片在线播放av一区| 美女脱内裤让男人舔精品视频| 亚洲欧美精品自产自拍| 51午夜福利影视在线观看| 精品国产乱码久久久久久男人| 成人亚洲欧美一区二区av| 90打野战视频偷拍视频| 手机成人av网站| 成年人午夜在线观看视频| 成年人免费黄色播放视频| 久久精品国产a三级三级三级| 电影成人av| 亚洲国产精品一区二区三区在线| 19禁男女啪啪无遮挡网站| 国产精品二区激情视频| 狠狠精品人妻久久久久久综合| 久久青草综合色| 精品一区在线观看国产| 又大又黄又爽视频免费| 亚洲成色77777| 亚洲精品美女久久久久99蜜臀 | 日本av免费视频播放| 精品人妻一区二区三区麻豆| 亚洲三区欧美一区| 咕卡用的链子| 国产免费一区二区三区四区乱码| 精品欧美一区二区三区在线| 成年美女黄网站色视频大全免费| 一区二区av电影网| 国产欧美日韩一区二区三区在线| 男人爽女人下面视频在线观看| 高清av免费在线| 狂野欧美激情性xxxx| √禁漫天堂资源中文www| 中国国产av一级| 高潮久久久久久久久久久不卡| 精品少妇久久久久久888优播| 国产麻豆69| 精品少妇黑人巨大在线播放| 999精品在线视频| 久久精品久久精品一区二区三区| 亚洲欧美成人综合另类久久久| 一区在线观看完整版| 观看av在线不卡| 在线观看免费日韩欧美大片| 欧美变态另类bdsm刘玥| 国产免费一区二区三区四区乱码| 亚洲熟女精品中文字幕| 国产有黄有色有爽视频| 80岁老熟妇乱子伦牲交| 欧美另类一区| 国产日韩欧美亚洲二区| 国产亚洲欧美在线一区二区| 久久精品熟女亚洲av麻豆精品| 精品人妻1区二区| 麻豆乱淫一区二区| 亚洲精品国产一区二区精华液| 国产亚洲精品久久久久5区| 黄网站色视频无遮挡免费观看| 一区二区三区四区激情视频| 一级毛片电影观看| 国产精品久久久人人做人人爽| 精品高清国产在线一区| 亚洲精品自拍成人| 精品少妇内射三级| 捣出白浆h1v1| kizo精华| 久久久久久久大尺度免费视频| 天堂8中文在线网| 别揉我奶头~嗯~啊~动态视频 | 亚洲中文日韩欧美视频| 亚洲欧美日韩高清在线视频 | 三上悠亚av全集在线观看| 黑人猛操日本美女一级片| 国产成人一区二区三区免费视频网站 | 9色porny在线观看| 高清视频免费观看一区二区| 男人爽女人下面视频在线观看| 国产高清不卡午夜福利| 日韩中文字幕欧美一区二区 | 美女高潮到喷水免费观看| 老司机深夜福利视频在线观看 | 黄色 视频免费看| 免费在线观看视频国产中文字幕亚洲 | 亚洲精品成人av观看孕妇| 久久国产精品人妻蜜桃| 电影成人av| 女人爽到高潮嗷嗷叫在线视频| 黄色视频在线播放观看不卡| 少妇人妻久久综合中文| √禁漫天堂资源中文www| 亚洲精品av麻豆狂野| 午夜日韩欧美国产| 9热在线视频观看99| 中文字幕制服av| 男女午夜视频在线观看| 亚洲精品一二三| 中文字幕另类日韩欧美亚洲嫩草| 免费看十八禁软件| 一二三四社区在线视频社区8| 丝袜喷水一区| 亚洲国产日韩一区二区| 精品福利永久在线观看| 国产成人av教育| 狂野欧美激情性xxxx| 多毛熟女@视频| 丝袜脚勾引网站| 欧美在线黄色| 欧美激情高清一区二区三区| 老司机午夜十八禁免费视频| 天天躁日日躁夜夜躁夜夜| 亚洲av在线观看美女高潮| 老司机影院毛片| 久久这里只有精品19| 亚洲国产精品一区三区| 久久人人爽人人片av| 亚洲国产看品久久| 亚洲国产欧美一区二区综合| 亚洲精品国产一区二区精华液| 在现免费观看毛片| 丝袜喷水一区| 久久国产精品大桥未久av| 狠狠婷婷综合久久久久久88av| 人成视频在线观看免费观看| 少妇人妻久久综合中文| 成人18禁高潮啪啪吃奶动态图| 免费看不卡的av| 国产av一区二区精品久久| 国产高清不卡午夜福利| 欧美成狂野欧美在线观看| 欧美中文综合在线视频| 丰满迷人的少妇在线观看| 欧美精品亚洲一区二区| 国产免费又黄又爽又色| 女性生殖器流出的白浆| 精品福利永久在线观看| 男女之事视频高清在线观看 | 亚洲成人免费电影在线观看 | 母亲3免费完整高清在线观看| 久久精品久久久久久噜噜老黄| 18禁裸乳无遮挡动漫免费视频| 亚洲 欧美一区二区三区| www日本在线高清视频| 狂野欧美激情性bbbbbb| videos熟女内射| 亚洲国产精品国产精品| 人人妻人人添人人爽欧美一区卜| 免费日韩欧美在线观看| 国产在线观看jvid| 国产免费现黄频在线看| 久久久久久久久久久久大奶| www.熟女人妻精品国产| 亚洲精品一卡2卡三卡4卡5卡 | 美女扒开内裤让男人捅视频| 亚洲成人手机| 亚洲欧美精品综合一区二区三区| 老熟女久久久| 午夜福利乱码中文字幕| 国产精品av久久久久免费| 国产成人一区二区三区免费视频网站 | 久久国产精品人妻蜜桃| 亚洲国产欧美日韩在线播放| 九草在线视频观看| 亚洲一区二区三区欧美精品| 99国产精品一区二区蜜桃av | 精品亚洲成a人片在线观看| 久久久精品免费免费高清| 老司机影院毛片| 熟女少妇亚洲综合色aaa.| 午夜两性在线视频| 欧美黑人欧美精品刺激| 久久精品亚洲av国产电影网| 国产99久久九九免费精品| av天堂在线播放| 国产成人精品久久二区二区91| 国产精品三级大全| 99热全是精品| 国产午夜精品一二区理论片| 狂野欧美激情性xxxx| 国产欧美亚洲国产| 丝袜人妻中文字幕| 国产激情久久老熟女| a 毛片基地| 人人妻人人爽人人添夜夜欢视频| 中文欧美无线码| 91精品伊人久久大香线蕉| 国产在线一区二区三区精| 欧美xxⅹ黑人| 国产成人av教育| 精品国产一区二区三区久久久樱花| 肉色欧美久久久久久久蜜桃| 日韩中文字幕视频在线看片| 亚洲国产精品成人久久小说| 在线观看免费视频网站a站| 国产成人精品在线电影| 一级片'在线观看视频| 亚洲欧美精品综合一区二区三区| 国产又色又爽无遮挡免| 免费人妻精品一区二区三区视频| 免费在线观看日本一区| 亚洲 国产 在线| 国产熟女欧美一区二区| 亚洲精品乱久久久久久| 一本久久精品| 国产淫语在线视频| 欧美成狂野欧美在线观看| 看免费av毛片| 在线亚洲精品国产二区图片欧美| 看十八女毛片水多多多| 国产亚洲欧美精品永久| 国产精品国产三级国产专区5o| 亚洲精品美女久久久久99蜜臀 | 欧美性长视频在线观看| 18禁国产床啪视频网站| bbb黄色大片| 成年人午夜在线观看视频| 国产欧美日韩精品亚洲av| 国产一区二区在线观看av| 精品久久久精品久久久| 亚洲国产看品久久| 后天国语完整版免费观看| 丰满少妇做爰视频| 午夜福利一区二区在线看| 免费观看a级毛片全部| 久久99热这里只频精品6学生| 欧美 日韩 精品 国产| 精品人妻熟女毛片av久久网站| 七月丁香在线播放| 一二三四在线观看免费中文在| 大香蕉久久网| 亚洲天堂av无毛| 大话2 男鬼变身卡| 色网站视频免费| e午夜精品久久久久久久| 老司机深夜福利视频在线观看 | 国语对白做爰xxxⅹ性视频网站| 精品视频人人做人人爽| 视频区图区小说| 亚洲精品成人av观看孕妇| 一级毛片黄色毛片免费观看视频| 9色porny在线观看| 亚洲第一av免费看| 亚洲天堂av无毛| 欧美日韩av久久| 国产亚洲欧美在线一区二区| e午夜精品久久久久久久| 不卡av一区二区三区| av视频免费观看在线观看| 国产精品秋霞免费鲁丝片| 国产片特级美女逼逼视频| 亚洲精品美女久久久久99蜜臀 | 性色av一级| 欧美激情极品国产一区二区三区| 亚洲人成电影免费在线| 真人做人爱边吃奶动态| 精品一区二区三卡| 嫁个100分男人电影在线观看 | 黄色视频不卡| 一级毛片 在线播放| 美女脱内裤让男人舔精品视频| 国产男人的电影天堂91| 国产在视频线精品| 爱豆传媒免费全集在线观看| 国产亚洲av高清不卡| 在线精品无人区一区二区三| 久久人人97超碰香蕉20202| 久久精品国产亚洲av涩爱| 日本a在线网址| 久久av网站| 欧美精品一区二区大全| 青春草亚洲视频在线观看| 纯流量卡能插随身wifi吗| 免费观看av网站的网址| 男女国产视频网站| 三上悠亚av全集在线观看| 国产成人精品久久久久久| 成人影院久久| 日本五十路高清| 人成视频在线观看免费观看| 亚洲 欧美一区二区三区| 男女午夜视频在线观看| 亚洲国产欧美一区二区综合| 丁香六月欧美| 亚洲国产精品成人久久小说| 王馨瑶露胸无遮挡在线观看| 精品人妻1区二区| 国产福利在线免费观看视频| 久久热在线av| 成人国产一区最新在线观看 | 18禁裸乳无遮挡动漫免费视频| 亚洲欧美中文字幕日韩二区| 国产欧美日韩一区二区三 | 国产亚洲欧美在线一区二区| 国产精品亚洲av一区麻豆| 欧美精品亚洲一区二区| 欧美激情极品国产一区二区三区| 久久久国产一区二区| 国产高清视频在线播放一区 | 18在线观看网站| 午夜福利影视在线免费观看| 日本黄色日本黄色录像| 国产1区2区3区精品| 一二三四社区在线视频社区8| 色视频在线一区二区三区| 在线精品无人区一区二区三| 久久国产亚洲av麻豆专区| 五月开心婷婷网| 丁香六月欧美| 日韩熟女老妇一区二区性免费视频| 久9热在线精品视频| 麻豆乱淫一区二区| 午夜av观看不卡| 老司机靠b影院| av福利片在线| 日韩制服丝袜自拍偷拍| 女警被强在线播放| 欧美亚洲日本最大视频资源| 人成视频在线观看免费观看| 亚洲国产欧美网| 国产午夜精品一二区理论片| 午夜av观看不卡| 亚洲精品日本国产第一区| 午夜91福利影院| 不卡av一区二区三区| 91精品三级在线观看| 国产成人91sexporn| 一区二区三区激情视频| 午夜福利在线免费观看网站| 日本猛色少妇xxxxx猛交久久| 十八禁高潮呻吟视频| 久久精品国产亚洲av高清一级| 久久青草综合色| 国产又爽黄色视频| 又大又黄又爽视频免费| 国产精品久久久久久精品电影小说| 久久女婷五月综合色啪小说| 成人18禁高潮啪啪吃奶动态图| 精品人妻一区二区三区麻豆| 免费高清在线观看视频在线观看| www日本在线高清视频| 国产一区二区激情短视频 | 国产精品成人在线| 黄色a级毛片大全视频| 9191精品国产免费久久| 高清视频免费观看一区二区| 日本猛色少妇xxxxx猛交久久| 99久久人妻综合| 日本色播在线视频| 侵犯人妻中文字幕一二三四区| 成人免费观看视频高清| 我的亚洲天堂| 国产精品香港三级国产av潘金莲 | 国产亚洲av高清不卡| 90打野战视频偷拍视频| 国产片内射在线| 亚洲第一av免费看| 又大又黄又爽视频免费| 蜜桃国产av成人99| 日韩人妻精品一区2区三区| 热re99久久国产66热| 国产精品亚洲av一区麻豆| 国产在线一区二区三区精| 狂野欧美激情性xxxx| h视频一区二区三区| 男女高潮啪啪啪动态图| 精品人妻熟女毛片av久久网站| 男男h啪啪无遮挡| 中文字幕人妻丝袜制服| 精品福利观看| 成人手机av| 欧美 日韩 精品 国产| 精品久久久久久电影网| 久久精品熟女亚洲av麻豆精品| 激情五月婷婷亚洲| 国产欧美日韩精品亚洲av| 色播在线永久视频| 国产成人精品无人区| 一级a爱视频在线免费观看| 成年人黄色毛片网站| 欧美+亚洲+日韩+国产| 爱豆传媒免费全集在线观看| 免费看不卡的av| av在线播放精品| 精品国产一区二区三区久久久樱花| 五月天丁香电影| 久久精品成人免费网站| 在线观看www视频免费| 视频区图区小说| 午夜影院在线不卡| 欧美另类一区| 久热这里只有精品99| 免费在线观看黄色视频的| √禁漫天堂资源中文www| 在线观看人妻少妇| 国产成人一区二区三区免费视频网站 | 2018国产大陆天天弄谢| 人人妻,人人澡人人爽秒播 | 成人国语在线视频| 久久这里只有精品19| 亚洲欧美一区二区三区国产| 99精国产麻豆久久婷婷| av在线播放精品| 国产免费又黄又爽又色| 两人在一起打扑克的视频| 丝袜在线中文字幕| 电影成人av| 人妻人人澡人人爽人人| 亚洲精品久久成人aⅴ小说| 一级毛片 在线播放| 日韩制服丝袜自拍偷拍| 一区二区av电影网| 国产精品一二三区在线看| 中国美女看黄片| 欧美精品亚洲一区二区| 欧美激情高清一区二区三区| 美女福利国产在线| 黄片小视频在线播放| 亚洲国产中文字幕在线视频| 熟女少妇亚洲综合色aaa.| 又紧又爽又黄一区二区| 国产免费福利视频在线观看| 王馨瑶露胸无遮挡在线观看| 无遮挡黄片免费观看| 久久影院123| 激情视频va一区二区三区| 国产精品 国内视频| 不卡av一区二区三区| 大陆偷拍与自拍| 天天躁日日躁夜夜躁夜夜| 91精品三级在线观看| 日本av手机在线免费观看| 在线精品无人区一区二区三| 91老司机精品| 1024香蕉在线观看| av天堂久久9| 久久亚洲国产成人精品v| 精品一区在线观看国产| 无遮挡黄片免费观看| 精品久久久久久电影网| 亚洲专区中文字幕在线| 日本av免费视频播放| 国产精品久久久久久精品古装| 嫁个100分男人电影在线观看 | 国产精品.久久久| 国产男女内射视频| 人人澡人人妻人| kizo精华| 欧美少妇被猛烈插入视频| 国产又爽黄色视频| 亚洲综合色网址| h视频一区二区三区| √禁漫天堂资源中文www| 亚洲人成网站在线观看播放| 国产精品欧美亚洲77777| 国产又色又爽无遮挡免| 欧美久久黑人一区二区| 亚洲,欧美,日韩| 免费一级毛片在线播放高清视频 | 免费看av在线观看网站| 国产一区二区三区综合在线观看| 日韩制服骚丝袜av| 久久久久久久久免费视频了| 久久毛片免费看一区二区三区| 青青草视频在线视频观看| 久久天堂一区二区三区四区| 欧美久久黑人一区二区| 19禁男女啪啪无遮挡网站| netflix在线观看网站| 国产精品av久久久久免费| 一区二区三区乱码不卡18| 你懂的网址亚洲精品在线观看| 热99国产精品久久久久久7| 亚洲一码二码三码区别大吗| 欧美日韩精品网址| a级毛片在线看网站| 国产精品免费视频内射| 视频在线观看一区二区三区| 久久久久精品国产欧美久久久 | 日本av免费视频播放| 国产精品免费视频内射| 亚洲精品日韩在线中文字幕| 国产精品一区二区精品视频观看| 视频在线观看一区二区三区| 秋霞在线观看毛片| 后天国语完整版免费观看| 精品人妻一区二区三区麻豆| 男女边吃奶边做爰视频| 十八禁网站网址无遮挡| 日本av手机在线免费观看| 成人亚洲欧美一区二区av| 五月天丁香电影| 免费看十八禁软件| 中文字幕亚洲精品专区| 国产精品亚洲av一区麻豆| 考比视频在线观看| 久久青草综合色| 色网站视频免费| 欧美国产精品va在线观看不卡| 交换朋友夫妻互换小说| 麻豆乱淫一区二区| 欧美精品一区二区免费开放| 久久久久久久国产电影| 一本色道久久久久久精品综合| 国产不卡av网站在线观看| 国产99久久九九免费精品| 欧美激情高清一区二区三区| 99热全是精品| 中文字幕人妻熟女乱码| 亚洲av在线观看美女高潮| 日韩大片免费观看网站| 成人午夜精彩视频在线观看| 天堂8中文在线网| videosex国产| 我要看黄色一级片免费的| 午夜福利免费观看在线| 一本久久精品| 高清黄色对白视频在线免费看| 国产老妇伦熟女老妇高清| 天天躁夜夜躁狠狠久久av| 国产日韩欧美视频二区| 精品熟女少妇八av免费久了| 国产亚洲精品第一综合不卡| 亚洲精品国产一区二区精华液| 久久av网站| 亚洲av片天天在线观看| 久久精品国产亚洲av高清一级| 亚洲精品成人av观看孕妇| 日韩伦理黄色片| 国产91精品成人一区二区三区 | 久久青草综合色| 亚洲一区二区三区欧美精品| 日本91视频免费播放|