李涌濤,李建文,龐鵬,劉德智,李婧

(1.信息工程大學(xué) 地理空間信息學(xué)院,河南 鄭州 450001;2.西安測繪總站,陜西 西安 710054)

0 引 言

電離層是距離地球60～1000 km的高層大氣受太陽高能輻射以及宇宙射線的激勵而電離形成的,因此電離層總電子含量(TEC)的變化很大程度上受到太陽活動的影響[1]。電離層除受太陽活動(如太陽耀斑、太陽黑子、太陽質(zhì)子事件等)的影響外還受地磁場等條件的影響,并隨之表現(xiàn)出不同的變化特性。電離層是空間環(huán)境研究的重要對象之一,它與人類生產(chǎn)、生活等息息相關(guān)。TEC是描述電離層形態(tài)和結(jié)構(gòu)的重要參數(shù),電離層TEC的異常變化會對地空無線電通信、衛(wèi)星導(dǎo)航定位、雷達(dá)等無線電信息系統(tǒng)電波信號傳播產(chǎn)生重要影響[2-3]。地磁暴是全球性地球磁場劇烈擾動的現(xiàn)象,發(fā)生強(qiáng)地磁暴時,會導(dǎo)致地球低軌道大氣密度的急劇增大,對衛(wèi)星的拖曳效應(yīng)迅速上升,引起衛(wèi)星軌道的快速衰減;地磁暴及其后數(shù)天,可能引發(fā)高能電子暴,高軌衛(wèi)星的深層充電效應(yīng)概率增大;可能發(fā)生電離層暴,引起全球?qū)Ш蕉ㄎ痪认陆?。全球?qū)Ш叫l(wèi)星的發(fā)展和廣泛應(yīng)用為探測、研究電離層帶來了革命性的變化,利用GNSS信號經(jīng)過電離層時發(fā)生折射效應(yīng)從而導(dǎo)出電離層總電子含量TEC[4],并利用GNSS全球跟蹤站長期觀測數(shù)據(jù)探測電離層TEC變化具有全球覆蓋、數(shù)據(jù)連續(xù)、時空分辨率和測量精度高、運行方便等顯著優(yōu)勢[3]。GNSS信號在經(jīng)過電離層時發(fā)生延遲,從而產(chǎn)生測距誤差,因此研究分析太陽活動、磁暴對TEC的影響[1]及TEC的變化規(guī)律對提高導(dǎo)航系統(tǒng)定位和導(dǎo)航精度具有十分重要的意義[5-6]。

本文利用國際GNSS監(jiān)測評估系統(tǒng)(iGMAS)分析中心解算的全球電離層TEC格網(wǎng)數(shù)據(jù)(GIM),iGMAS網(wǎng)址為:http://www.igmas.org/,分析了2017年9月6日太陽爆發(fā)X9.3級特大耀斑后引發(fā)的磁暴對全球電離層TEC變化的影響,以及全球不同緯度帶TEC與磁暴Dst指數(shù)的相關(guān)關(guān)系和變化特性。

1 數(shù)據(jù)來源

選取全球均勻分布的177個GNSS跟蹤站的觀測數(shù)據(jù),采樣率為30 s,衛(wèi)星截止高度角為10°,觀測量為L4非差無幾何線性組合相位平滑偽距,采用15階15次的球諧函數(shù)建立全球性的VTEC模型,球諧函數(shù)模型具體形式如公式(1)所示:

(1)

利用以上公式生成450 km高度,采樣間隔1 h,經(jīng)度差為5°和緯度差為2.5°的每天24幅全球電離層TEC格網(wǎng)數(shù)據(jù)(GIM),精度為8 TECU(電離層TEC單位,1TECU=1016/m2),足以滿足實驗精度要求。

本文地磁數(shù)據(jù)采用中國科學(xué)院國家空間科學(xué)中心(NSSC)空間環(huán)境預(yù)報中心(SEPC)提供的Dst指數(shù)數(shù)據(jù)。

2 數(shù)據(jù)分析

太陽耀斑爆發(fā)經(jīng)常向行星際空間噴射等離子體云,其徑向速度達(dá)100～500 km/s,經(jīng)過1.5～3天到達(dá)地球附近,與地磁場作用引起磁暴[7],2017年9月6日11:53 UTC太陽爆發(fā)X9.3級特大耀斑,13 h后發(fā)生地磁暴,39 h后磁暴主相達(dá)到峰值,取耀斑爆發(fā)后到磁暴恢復(fù)相之間即9月6日(DOY 249)12:00 UTC至9月9日(DOY 252)24:00 UTC之間85 h數(shù)據(jù)的TEC數(shù)據(jù)(如圖1所示)和Dst指數(shù)(如圖2所示)進(jìn)行分析,其中TEC數(shù)據(jù)采用每天24幅GIM數(shù)據(jù)計算電離層TEC平均值,采樣間隔為1 h.

磁暴環(huán)電流指數(shù)(Dst指數(shù))是在地球赤道附近,按大致均勻的經(jīng)度間隔選取四個地磁臺站,取這四個臺站每小時地磁水平強(qiáng)度變化的平均值,單位為納特斯拉(nT)[8]。磁暴時地磁水平分量H變化最大,最能代表磁暴過程特點,磁暴的大部分形態(tài)學(xué)和統(tǒng)計學(xué)特征是依據(jù)中低緯度H分量的變化得到的。典型磁暴的發(fā)展過程也是按照H分量的變化來劃分的,通?？煞譃槿齻€階段:初相、主相和恢復(fù)相[8]。

磁暴開始發(fā)生的典型特征是水平分量H突然正向增大,稱為磁暴急始(SSC),磁暴發(fā)生后,在一段時間內(nèi)磁場沒有明顯起伏變化,稱為初相。在初相后,水平分量H大幅度下降,直到最低值,并伴隨著劇烈擾動,即為主相[10]。磁暴的主要特點反映在主相期間,主相的最低點幅度是衡量磁暴大小的標(biāo)準(zhǔn)[9],其值逐漸減小表示磁擾幅度逐漸增大。主相后磁場較弱,逐漸恢復(fù)至正常水平,稱為恢復(fù)相[10]。由圖2可知13-34 h為初相,34-39 h為主相,39-85 h為恢復(fù)相。

兩個變量的相關(guān)系數(shù)r,當(dāng)0.3<|r|≤0.5時,表示兩個變量低度相關(guān);當(dāng)0.5<|r|≤0.8時,表示兩個變量顯著相關(guān);當(dāng)0.8<|r|≤1時,表示兩個變量高度相關(guān);當(dāng)|r|=1時,表示兩個變量完全相關(guān)[11];

初相13-34 h期間Dst指數(shù)與TEC相關(guān)系數(shù)為-0.579 0,表現(xiàn)為顯著負(fù)相關(guān)性;主相39-85 h期間Dst指數(shù)與TEC相關(guān)系數(shù)為-0.960 5,Dst指數(shù)與延后1 h的TEC相關(guān)系數(shù)為-0.999 7,表現(xiàn)出非常強(qiáng)的負(fù)相關(guān)性,此過程中TEC迅速增大3.3 TECU,達(dá)到峰值后5 h內(nèi)迅速減小4.5 TECU,表明磁暴主相期間H分量的變化引起了TEC的劇烈變化,且變化幅度較大。在恢復(fù)相39-85 h期間,TEC變化趨于穩(wěn)定,Dst指數(shù)與TEC相關(guān)系數(shù)為-0.8561,50-85 h期間Dst指數(shù)與TEC相關(guān)系數(shù)為-0.9589,表現(xiàn)出了非常強(qiáng)的負(fù)相關(guān)性。從初相到恢復(fù)相整個磁暴過程中Dst指數(shù)與TEC變化總體呈現(xiàn)強(qiáng)相關(guān)性,兩者整體變化一致,因此可以認(rèn)為磁暴的變化對TEC的變化產(chǎn)生了重大的影響。從相關(guān)系數(shù)分析,磁暴對TEC變化的影響主要是在主相階段和恢復(fù)相階段,并且主相期間的影響比恢復(fù)相期間更大。電離層和磁層之間通過許多不同的途徑相互耦合在一起,幾乎在磁層中發(fā)生的每一個過程都會對電離層造成影響,而在電離層中發(fā)生的過程也都可以在磁層中有所反應(yīng)[12]。并且從圖1和圖2的圖型上可以看出極強(qiáng)的對稱性,因此可以認(rèn)為在磁暴期間TEC的變化也能間接反映出磁暴的發(fā)生和磁場的變化。

按緯度劃分低緯度帶(0°～30°),中緯度帶(30°～60°),高緯度帶(60°～90°)統(tǒng)計各緯度帶的1 h采樣間隔的TEC平均值,主相和恢復(fù)相34-85 h期間各緯度帶TEC如圖3所示。

從圖中可看出主相34-39 h期間Dst達(dá)到峰值時(第39 h, DOY 251 02:00 UTC),高緯度帶TEC也達(dá)到峰值,即高緯帶上TEC隨磁暴Dst指數(shù)同時發(fā)生變化,1 h后(第40 h,DOY 251 03:00 UTC)中緯度帶TEC達(dá)到峰值,再1 h后(第41h, DOY 251 04:00 UTC)低緯度帶TEC達(dá)到峰值,并在恢復(fù)相開始5 h后(第44 h,DOY 251 07:00 UTC)不同緯度帶TEC變化時延消失,變化也趨于一致。從而也解釋了不分緯度帶時全球1 h采樣間隔的TEC平均值與Dst指數(shù)存在1 h時延。說明磁暴對不同緯度帶的影響不是同步的,而是從高緯到低緯的順序,由不同緯度帶TEC隨磁暴整體變化一致,說明磁暴在不同緯度帶的影響機(jī)制相同。其中可能的原因是日冕物質(zhì)拋射(CME)引發(fā)磁暴,尤其是中等以上非重現(xiàn)性地磁暴的主要起因[9],當(dāng)磁暴或者亞暴時,在磁尾的能量會通過場向電流、粒子沉降和對流電場等方式向地球兩極傳輸,產(chǎn)生極光,引起電離層暴[13],從而導(dǎo)致TEC變化異常。

統(tǒng)計磁暴不同階段Dst指數(shù)與TEC相關(guān)系數(shù)如表1所示。

表1 Dst指數(shù)與TEC相關(guān)系數(shù)

由表1分析可知磁暴主相階段對高緯度帶TEC的影響最大,其次是中緯度帶,最后是低緯度帶,可能的原因是磁暴對電離層的影響由高緯至低緯的順序時延造成了微小的相關(guān)性差異。

如圖3顯示磁暴期間同一時刻不同緯度帶TEC不同,高緯4～9 TECU,中緯6～14 TECU,低緯16～27 TECU,低緯明顯高于中緯和高緯,并不能說明磁暴對低緯TEC影響最大,這是因為在低緯地區(qū)太陽輻射能量比中緯、高緯大很多,同時在低緯地區(qū)存在赤道異常現(xiàn)象,即高度電子濃度的分布和峰值電子濃度在地磁赤道上空有極小值,在地磁赤道南北緯度15°～20°區(qū)域的兩側(cè)上空出現(xiàn)極大值,形成雙峰分布[14-15]。即使不發(fā)生磁暴赤道異常現(xiàn)象也存在,圖4示出了主相開始之后的15 h全球電離層TEC圖(DOY 250 21:00UTC-DOY 251 11:00UTC,對應(yīng)圖1中34-48 h)。DOY 251 01:00-05:00UTC,當(dāng)Dst指數(shù)達(dá)到峰值后一小時(DOY 251 03:00UTC)全球TEC達(dá)到最大值,且Dst指數(shù)峰值前后表現(xiàn)出明顯的赤道異?，F(xiàn)象,如圖4所示。

恢復(fù)相從DOY 251 02:00UTC開始后一直延續(xù)幾天,在這段期間磁暴逐漸恢復(fù)平靜水平,由圖4所示,從DOY 251 02:00UTC后數(shù)小時內(nèi)(由于篇幅所限,僅展示至DOY251 08:00UTC的TEC圖)隨磁暴減弱,TEC逐漸減少并趨于穩(wěn)定,沒有較大的變化和起伏。

在恢復(fù)相后期DOY 251 20:00-23:00UTC四小時期間,對應(yīng)圖1和圖2中的57-60 h,圖1中TEC小時均值與圖2中的Dst指數(shù)變化平穩(wěn),說明期間沒有發(fā)生磁暴起伏擾動,但是在低緯度帶出現(xiàn)大面積的TEC增大現(xiàn)象,且TEC最大值比磁暴最劇烈時高5.5 TECU,如圖5所示,從圖例看即高出了一個色度,可能是因為磁暴恢復(fù)相期間的某種特性變化引起的這種變化,具體原因尚不明確,有待進(jìn)一步分析研究。

3 結(jié)束語

本文基于2017年9月發(fā)生的X9.3級特大耀斑并在之后引發(fā)的磁暴等活動,利用GIM的小時數(shù)據(jù)與磁暴期間初相、主相和恢復(fù)相三個階段的Dst指數(shù)進(jìn)行相關(guān)性分析,重點分析了磁暴期間不同緯度帶電離層TEC的變化情況,結(jié)果表明:1) 此次太陽耀斑爆發(fā)13 h后引發(fā)磁暴,磁暴對TEC變化影響顯著,主相和恢復(fù)相階段Dst指數(shù)與TEC具有很強(qiáng)的相關(guān)性,主相階段Dst指數(shù)與時延1 h的TEC相關(guān)系數(shù)為-0.999 7,表現(xiàn)出完全負(fù)相關(guān)性,且磁暴對電離層TEC的影響時延約為1 h;2) 磁暴強(qiáng)度決定了TEC的變化程度,在磁暴主相階段TEC迅速增加,影響也最大,恢復(fù)相階段TEC隨磁暴減弱逐漸減小,變化趨于穩(wěn)定,在主相和恢復(fù)相階段Dst指數(shù)與TEC變化總體呈強(qiáng)相關(guān)性,兩者整體變化一致,在此次磁暴中可用TEC間接描述磁暴的變化水平;3) 磁暴期間對全球不同緯度帶電離層變化趨于一致,磁暴對TEC變化影響機(jī)制在不同緯度帶相同,影響程度大小由高緯至低緯逐漸遞減;4) 磁暴對電離層的影響并非全球同步,而是從高緯到低緯逐漸影響,此次磁暴對TEC高緯到低緯的影響時延約為2 h,并在磁暴的恢復(fù)相階段時延逐漸消失。磁暴與電離層TEC變化關(guān)系密切,磁暴恢復(fù)相后期幾小時內(nèi),低緯度TEC的劇烈變化原因還有待進(jìn)一步分析研究。

致謝: 感謝信息工程大學(xué)iGMAS分析中心和中國科學(xué)院空間環(huán)境預(yù)報中心給予的幫助和支持。