馬?？? 常紅芳, 孫 琳, 郭欣興

(1. 西安工程大學(xué) 理學(xué)院， 陜西 西安 710048; 2. 西安電子科技大學(xué) 物理與光電工程學(xué)院, 陜西 西安 710071)

極區(qū)電離層GPS信號閃爍特征統(tǒng)計分析

馬?？?,2, 常紅芳1, 孫 琳1, 郭欣興1

(1. 西安工程大學(xué) 理學(xué)院， 陜西 西安 710048; 2. 西安電子科技大學(xué) 物理與光電工程學(xué)院, 陜西 西安 710071)

對意大利高緯電離層閃爍監(jiān)測網(wǎng)2004年全年的電離層閃爍數(shù)據(jù)進行分析，通過計算GPS 信號振幅和相位閃爍在給定強度變化范圍內(nèi)的發(fā)生率，研究太陽活動低年GPS 信號L1頻標(biāo)電離層振幅、相位閃爍的發(fā)生率及其時空分布隨電離層總電子含量、不同季節(jié)的氣候?qū)W變化特征。結(jié)果表明:在高緯地區(qū)，電離層相位閃爍比振幅閃爍更容易發(fā)生；在極蓋區(qū)，相位閃爍的季節(jié)性特征較為明顯；不同的電離層總電子含量變化率可導(dǎo)致產(chǎn)生不同類型的電離層閃爍。

電離層閃爍；全球定位系統(tǒng)；統(tǒng)計特性

電離層及其電子密度不均勻體的存在，使得穿越于其中的無線電波的幅度、相位、到達(dá)角以及極化狀態(tài)等發(fā)生快速起伏，即為電離層閃爍。電離層閃爍主要集中于磁赤道和高緯地區(qū)[1-2]，高緯地區(qū)的電離層閃爍，其并發(fā)的物理機制異常復(fù)雜，如來自太陽的極紫外輻射電離、X射線輻射、太陽耀斑以及日冕物質(zhì)拋射的高速帶電粒子和太陽風(fēng)等，致使整個高緯、極蓋及極隙區(qū)的電離層電子密度結(jié)構(gòu)發(fā)生變化，產(chǎn)生一系列復(fù)雜的時空分布的電子密度不均勻體，進而形成高緯地區(qū)的電離層閃爍[2-3]。

全球定位系統(tǒng)(Global Positioning System, GPS)的出現(xiàn)為電離層閃爍研究提供了新途徑。利用GPS衛(wèi)星監(jiān)測臺網(wǎng)，能夠監(jiān)測局域或全球范圍內(nèi)的電離層閃爍及電離層電子濃度總量(Total Electron Content， TEC)的變化。

文[4]利用在澳大利亞以北至東南亞地區(qū)建立的GPS電離層閃爍監(jiān)測網(wǎng)，對該區(qū)域的電離層閃爍進行了分析。文[5]將GPS電離層閃爍的觀測結(jié)果同電離層閃爍寬帶模型(Wide Band Model, WBMOD)進行了對比，文[6]總結(jié)了多年的觀測結(jié)果，分析了太陽活動對電離層閃爍的影響。文[7-8]利用海南站的電離層閃爍數(shù)據(jù)，分析了太陽活動高到低年，擴展F和電離層閃爍的關(guān)系。文[9]利用海南站的電離層數(shù)據(jù)分析了太陽活動低年，電離層GPS信號的閃爍特征和譜特征。

受地理位置以及數(shù)據(jù)獲取途徑的制約，高緯電離層閃爍的研究在我國相對甚少，而高緯電離層閃爍的研究，是對電離層閃爍全球形態(tài)及其并發(fā)物理機制研究的有力補充，具有重要意義。

本文對意大利電離層閃爍監(jiān)測網(wǎng)(Istituto Nazionale di Geofisica e vulcanologia，INGV)[10]的電離層閃爍數(shù)據(jù)進行統(tǒng)計分析，以2004年為例，對比研究太陽活動低年，高緯電離層GPS信號閃爍的極區(qū)特征。

1 電離層閃爍數(shù)據(jù)獲取

位于挪威Svalbard群島的電離層觀測站與我國南極中山站(69.37°S,76.38°E)在磁正午附近分別位于南北半球的極隙區(qū)，在夜側(cè)則分別位于南北半球的極蓋區(qū)，這樣的地理位置分布為共軛觀測高緯電離層閃爍提供了絕佳的條件。INGV的4個電離層閃爍監(jiān)測陣列中[10]，Nya0和Nya1分別位于挪威Svalbard的NyAlesund(78.9°N,11.9°E)，Lyb0位于Svalbard的Longyearbyen(78.2°N,16.0°E)，而另一個觀測陣列Btn0則位于南極洲的Mario Zucchelli觀測站(74.7°S,164.1°E)。Nya0、Nya1和Lyb0三站的觀測陣列通常可對亞極光、極光、極隙區(qū)和極蓋區(qū)的電離層閃爍活動進行觀測，而Btn0和我國中山站通常可觀測極隙區(qū)及極蓋區(qū)電離層活動。另外，INGV的監(jiān)測設(shè)備為GSV4004雙頻GPS閃爍和TEC監(jiān)測儀，該設(shè)備以50 Hz的速率采樣GPS L1(1 575.42 MHz)信標(biāo)的幅度和相位閃爍。

2 數(shù)據(jù)處理

為了避免人為因素引起的誤差，將數(shù)據(jù)的間斷段以及大的調(diào)變值作丟棄處理，且選用前后漲落相差不是太大的數(shù)據(jù)段。為了減少斜入射、多徑效應(yīng)及低仰角效應(yīng)對統(tǒng)計結(jié)果的影響，選取GPS衛(wèi)星仰角大于25°的數(shù)據(jù)。定義閃爍事件的發(fā)生率為一段時間內(nèi)接收到振幅、相位閃爍滿足某強度要求的閃爍發(fā)生個數(shù)與該時間段內(nèi)接收到閃爍總發(fā)生個數(shù)的比率。將振幅閃爍劃分為3個等級,即0.20.6為強閃爍。將相位閃爍劃分為3個等級,即0.1 rad<σφ≤0.5 rad為弱閃爍, 0.5 rad<σφ≤1 rad為中等閃爍, σφ>1 rad為強閃爍。

振幅閃爍S4定義為減去系統(tǒng)噪聲每分鐘信號強度的歸一化標(biāo)準(zhǔn)差[3-8]，即

(1)

其中〈〉表示一定時間間隔內(nèi)的均值，I′為信號強度I經(jīng)過濾波和消趨勢后的值，C/N0為系統(tǒng)輸出的信噪比，一般取均值。

相位閃爍σφ定義為消趨后信號載波相位的標(biāo)準(zhǔn)偏差，不妨使用濾波器函數(shù)對原始相位數(shù)據(jù)進行消趨處理。

INGV還可提供15 s、30 s、45 s和60 s四種不同間隔的總電子含量TEC或相對TEC值ETEC。為便于分析，將斜路徑數(shù)據(jù)投影到了垂直路徑上，即

(2)

(3)

ETEC，Vert=ETEC,Slantsinαelev,

(4)

其中傾斜因子

(5)

Re為地球半徑，HIPP代表電離層穿刺點高度，αelev為衛(wèi)星仰角。式(2)中的指數(shù)a依賴于相位閃爍的譜指數(shù)p，對于各向異性電離層，根據(jù)譜指數(shù)的變化，a通常可取0.9。式(3)中的指數(shù)b一般取為0.5[11]。

3 統(tǒng)計分析

(1) 太陽低年北極GPS信號電離層閃爍隨月份的變化。

2004年為太陽低年，太陽黑子數(shù)的月均值在20～50之間，太陽活動在2004年基本表現(xiàn)為平靜。

2004年北極NyAlesund地區(qū)的Nya0和Nya1兩個GPS電離層閃爍和TEC觀測站觀測到的GPS L1頻標(biāo)信號振幅、相位閃爍的統(tǒng)計結(jié)果如圖1所示。由圖1(a)可見，在對GPS 11、21和5號衛(wèi)星連續(xù)一年的數(shù)據(jù)觀測中，三顆衛(wèi)星的L1頻標(biāo)在北極全年主要表現(xiàn)為弱閃爍，且春季振幅閃爍的發(fā)生率相對較小，同時，閃爍指數(shù)S4>0.4的中等強度閃爍主要集中在夏季和秋季。由圖1(b)可見，相比于振幅閃爍，相位閃爍0.1 rad<σφ<0.5 rad在北極全年均表現(xiàn)為高發(fā),相位閃爍σφ>0.5的中等閃爍主要發(fā)生在春季、夏季和秋季，其較高的發(fā)生率可能與這一時期電離層不均勻體的形成及其空間結(jié)構(gòu)的分布有關(guān)。

(a) 振幅閃爍S4的發(fā)生率隨月份的變化

(b) 相位閃爍σφ的發(fā)生率隨月份的變化

(2) 太陽低年北極GPS信號振幅、相位閃爍隨地理經(jīng)緯度的變化。

2004年北極GPS L1頻標(biāo)電離層閃爍的發(fā)生率在地理坐標(biāo)下的分布如圖2所示。從圖2(a)可見，2004年，振幅閃爍S4的發(fā)生率具有明顯的極區(qū)分布特征，北半球振幅閃爍的高發(fā)區(qū)在北緯70°以上。從圖2(b)可見，北極GPS信號的相位閃爍更加頻繁，在北緯70°以上，相位閃爍σφ極為頻繁，且受觀測站點位置的限制，北半球的相位閃爍發(fā)生率在地理坐標(biāo)下則呈三角狀分布。

(a) 北極振幅閃爍 S4

(b) 北極相位閃爍 σφ

(3) 太陽低年北極GPS 信號相位閃爍的季節(jié)性變化。

北極相位閃爍在冬季和夏季的分布變化如圖3所示。在北極的冬季和夏季，北緯70°以上，相位閃爍的發(fā)生率很高，地理位置從西經(jīng)20°一直持續(xù)到東經(jīng)53°左右。相位閃爍與該區(qū)域總電子含量的變化率(Rate Of TEC, ROT)密切相關(guān)，北極冬季相位閃爍的高發(fā)率說明，冬季這一區(qū)域電子密度的漲落很大，理論上可以存在各種尺度的電子密度不均勻體，而高的電子濃度變化率及其均方根是導(dǎo)致GPS L波段相位閃爍的主要原因。這與采用無線電層析成像技術(shù)得到的結(jié)論一致[12-13]。

(a) 北極冬季相位閃爍σφ

(b) 北極夏季相位閃爍 σφ

(4) 太陽低年北極電離層電子濃度總量變化率的分布。

為進一步說明電離層閃爍與電離層電子密度不均勻體及其空間分布之間的聯(lián)系，圖4給出了北極總電子含量TEC變化率(ROT)的空間分布。由圖4(a)可見，北半球總電子濃度變化率的分布不管在量值上，還是地理位置上均較大。北半球總電子濃度變化率較大的區(qū)域呈徑向分布，北緯60°左右，總電子濃度變化率最大，這也預(yù)示著，由于在該區(qū)域各種并發(fā)物理機制的出現(xiàn)，如高緯槽，以及高密度的太陽輻射電離，等離子體向著或背離極蓋區(qū)輸運，從而使得該區(qū)域電離層總電子濃度變化率快速起伏，導(dǎo)致GPS信號振幅和相位發(fā)生閃爍。由圖4(b)可見，在冬季，北半球總電子濃度變化率的最大值主要分布在北緯57°以上，且北緯60°左右電子濃度總量的變化率很大。由圖4(c)可見，而到了夏季，北半球總電子濃度變化率的最大值主要分布在北緯75°左右以及東經(jīng)25°～50°的一個呈徑向分布的帶上。

夏季出現(xiàn)的E層電導(dǎo)，使得增強的橫過磁場傳播的離子產(chǎn)生的電場發(fā)生短路，進而導(dǎo)致這一區(qū)域不均勻體強度減弱，從而使得閃爍的強度譜變得更加陡峭，譜指數(shù)增大[14]。電子濃度總量變化率的這種分布，也預(yù)示著該區(qū)域可能存在尺度上從幾米到數(shù)十千米的電子密度不均勻體，正是因為它們空間上的這種分布，導(dǎo)致了穿越其中的GPS信號發(fā)生閃爍。

(a) 北極ROT的分布

(b) 北極冬季ROT的分布

(c) 北極夏季ROT的分布

4 結(jié)語

從太陽活動低年對高緯電離層閃爍特征的統(tǒng)計分析來看，在高緯地區(qū)，GPS L波段電離層閃爍的相位閃爍比振幅閃爍更容易發(fā)生。而就電離層閃爍發(fā)生的季節(jié)性特征來看，北極冬季和夏季，在北緯60°左右，相位閃爍與季節(jié)基本無關(guān)，在這一區(qū)域，電離層閃爍主要表現(xiàn)為相位閃爍。而在極蓋區(qū)，相位閃爍的季節(jié)性特征比較明顯，夏季，相位閃爍的發(fā)生率相對偏高。冬季，小尺度的電離層電子密度不均勻體在北極幾乎不出現(xiàn)，到了夏季，在觀測到振幅閃爍的地方，小尺度的電子密度不均勻體才趨于出現(xiàn)。統(tǒng)計結(jié)果表明，總電子含量變化率的不同，將可能導(dǎo)致產(chǎn)生不同類型的電離層閃爍。如在北緯60°以上，總電子含量的變化率很大，則在太陽活動低年，該區(qū)域主要表現(xiàn)為相位閃爍，振幅閃爍的發(fā)生率相對較小。

[1] 熊年祿,唐存琛,李行建.電離層物理概論[M]. 武漢：武漢大學(xué)出版社, 1999: 32-40.

[2] Li Guozhu, Ning Baiqi, Ren Zhipeng, et al.. Statistics of GPS ionospheric scintillation and irregularities over polar regions at solar minimum[J/OL].GPS Solution,2010,14:331-341. (2014-08-24) [2014-09-01]. http://www.polar.gov.cn/paperDetail/?id=7114.

[3] 劉鈍,馮健,鄧忠新,等. 電離層閃爍對GPS系統(tǒng)定位性能的影響研究[J]. 電波科學(xué)學(xué)報,2010,25(4):702-710.

[4] Thomas R M ,Cervera M A, Eftaxiadis K, et al.. A regional GPS receiver network for monitoring equatorial scintillation and total electron content [J]. Radio Science, 2001, 36(6):1545-1557.

[5] Cervera M A, Thomas R M . Groves K M, et al.. Validation of WBMOD in the southeast asian region [J]. Radio Science, 2001,36(6):1557-1572.

[6] Babayev E S. Influence of the solar and geomagnetic activity on trans-ionospheric radio waves scintillation[C]//Regional Meeting on Solar Physics Solar Researches in the South-Eastern European Countries. Bucharest: Proceedings of the Regional Meeting on Solar Physics, 2001: 29-34.

[7] 王霄,史建魁,肖佐.中國低緯(海南)電離層漂移特性[J].電波科學(xué)學(xué)報,2005,20(4):482-486.

[8] 尚社平,史建魁,郭善兼.海南地區(qū)電離層閃爍監(jiān)測及初步統(tǒng)計分析[J].空間科學(xué)學(xué)報,2005,25(1):23-28.

[9] 甄衛(wèi)民,龍其利,馬健敏.磁赤道異常區(qū)電離層F區(qū)不均勻體發(fā)展過程中的閃爍譜研究[J].空間科學(xué)學(xué)報, 1995,15(2):143-147.

[10] Istituto Nazionale di Geofisica e vulcanologia-Rome. Electrinc space weather upper atmosphere[EB/OL].(2001-01-10) [2013-10-15] http://www.eswua.ingv.it/ingv/home.php?res=1024.

[11] Spogli L, Alfonsi L, Franceschi G De, et al.. Climatology of GPS ionospheric scintillation over high and mid latitude European region[J]. Annals of Geophysics, 2009,27(3):3429-3437.

[12] Pryse S, Kersley L, Malan D, et al.. Parameterization of the main ionospheric trough in the European sector[J]. Radio Science, 2006,41(5):1-9.

[13] Pryse S E, Dewis K L, Balthazor R L, et al.. The dayside high‐latitude trough under quiet geomagnetic conditions: Radio tomography and the CTIP model [J]. Annals of Geophysics, 2005, 23(4): 1199-1206.

[14] Kivanc O, Heelis R A. Structures in ionospheric number density and velocity associated with polar cap ionization patches [J]. Journal Geophys Research, 1997,102(A1): 307-318.

[責(zé)任編輯:瑞金]

Statistical analyses on characteristics of GPS signal scintillation in high latitude ionosphere

MA Baoke1,2, CHANG Hongfang1, SUN Lin1, GUO Xinxing1

(1. School of Science, XI’an Polytechnic University, Xi’an 710048, China; 2. School of Physics and Optoelectronic Engineering, Xidian University, Xi’an 710071, China)

The ionospheric scintillation data which were obtained from the monitoring network of INGV (Istituto Nazionale di Geofisica e vulcanologia) of Italy in 2004 are analysed in this paper. The change characteristics of spatial and temporal distribution of amplitude and phase scintillation of L1-band frequency of GPS with different variation of ionospheric total electron content (TEC), different seasons as well as other related space parameters at the low solar activity are compared by calculating the incidence of amplitude and phase scintillation of signals of Global Position System (GPS) within a given intensity range. Results show that the phase scintillation is more likely to take place than the amplitude scintillation at high-latitude, the seasonal change feature for the phase scintillation is more obviously in polar cap, and meanwhile different rate of change of the total electron content may lead to different types of ionospheric scintillation.

ionospheric scintillation, global position system (GPS), statistical characteristic

2014-10-15

陜西省自然科學(xué)基金資助項目(2013JM8011)；陜西省教育廳專項科研計劃資助項目(2014JK1295)；西安工程大學(xué)博士科研啟動基金資助項目(BS1418)；國家大學(xué)生創(chuàng)新創(chuàng)業(yè)基金資助項目(201310709007)

馬?？?1972-)，男，博士，教授，從事復(fù)雜媒質(zhì)中的電磁波傳播與散射特性研究。E-mail: Baokema2006@126.com 常紅芳(1977-)，女，碩士，講師，從事電磁波的傳播特性研究。E-mail: chfang1018@163.com

10.13682/j.issn.2095-6533.2015.03.016

P352

A

2095-6533(2015)03-0092-05