馮婷 劉默然 王翔 周晨

（武漢大學(xué)電子信息學(xué)院，武漢 430072）

引 言

非相干散射雷達(dá)（incoherent scatter radar，ISR）是進(jìn)行電離層探測(cè)的重要地基無(wú)線電設(shè)備，可獲取電離層電子密度、電子溫度、離子溫度、等離子體漂移速度等多種參量[1-3]. 典型的ISR散射信號(hào)功率譜由離子譜線和等離子體譜線構(gòu)成. 標(biāo)準(zhǔn)的ISR電離層參數(shù)是通過(guò)測(cè)量離子聲波散射的離子譜線來(lái)獲取的，但其測(cè)量精度受到雷達(dá)系統(tǒng)的狀態(tài)、測(cè)量距離、積累時(shí)間、編碼方式等多種因素的影響，因此ISR測(cè)量的某些參量，如電子密度，通常需要校正和定標(biāo). 另外一種可靠且精度更高的電子密度測(cè)量可通過(guò)分析朗繆爾波散射的等離子體譜線來(lái)實(shí)現(xiàn)[4-5]，但等離子體譜線的信號(hào)強(qiáng)度通常低于雷達(dá)的噪聲水平，很難直接觀測(cè). 然而，某些情況下，ISR等離子體譜線的能量可以增強(qiáng)到噪聲水平之上，使得利用等離子體譜線提取電離層參量成為可能[6]. 例如，Djuth等人利用Arecibo站的日間觀測(cè)數(shù)據(jù)[7]，Carlson等人利用Arecibo站的夜間高頻（high frequency, HF）加熱實(shí)驗(yàn)數(shù)據(jù)[8]，Wang和Zhou利用EISCAT站的電離層加熱實(shí)驗(yàn)數(shù)據(jù)[9]，均從增強(qiáng)等離子體譜線中提取出了電子密度參數(shù).

Akbari等人在2017年的綜述文章[6]中指出，超熱電子有四種來(lái)源，分別是太陽(yáng)極紫外線輻射[6,10-11]、電子沉降[12-14]、電子束[6,15]和大功率高頻無(wú)線電波加熱電離層實(shí)驗(yàn)[8,16]. 源于太陽(yáng)遠(yuǎn)紫外線輻射的超熱電子被稱為光電子[17]，而源于極光粒子沉降的電子被稱為次級(jí)電子（secondary electrons）[18]. 此外，還有另外一種超熱電子是沿著磁力線從磁層運(yùn)動(dòng)至極區(qū)F層的粒子群，被稱為電子束[19-20]. 在不同緯度上，電離層等離子體譜線增強(qiáng)具有不同的機(jī)制. 在中低緯地區(qū)，由太陽(yáng)紫外線輻射所產(chǎn)生的光電子會(huì)導(dǎo)致當(dāng)?shù)鼗蛘叽殴曹椀貐^(qū)等離子體譜線增強(qiáng)，從而被Arecibo和Millstone Hill的ISR[21]觀測(cè)到. 在高緯極光區(qū)的電離層E層，Wickwar[22]和Vierinen等人[12]利用Sondrestorm ISR觀測(cè)到了強(qiáng)極光活動(dòng)期間，由強(qiáng)極光電子沉降引起的等離子體譜線增強(qiáng). 在高緯F區(qū)，電子束可以激發(fā)朗繆爾波[15]從而導(dǎo)致朗繆爾湍流產(chǎn)生，朗繆爾湍流通過(guò)等離子體不穩(wěn)定性導(dǎo)致等離子體譜線增強(qiáng)[6]. 此外，人工的電離層加熱實(shí)驗(yàn)也可以通過(guò)誘發(fā)和增強(qiáng)朗繆爾湍流來(lái)增強(qiáng)等離子體譜線[8].

增強(qiáng)的等離子體譜線為ISR研究電離層特性提供了一種更加準(zhǔn)確的方法，但其測(cè)量比較困難，觀測(cè)數(shù)據(jù)較為稀少. 本文針對(duì)歐洲非相干散射科學(xué)聯(lián)合會(huì)(European Incoherent Scatter Scientific Association，EISCAT)的超高頻(ultra high frequency, UHF)雷達(dá)觀測(cè)得到的增強(qiáng)等離子體譜線展開(kāi)討論，首先利用曲線擬合的方法將增強(qiáng)的譜線提取出來(lái)，然后通過(guò)朗繆爾色散關(guān)系反演得到電子密度剖面，最后由提取出的增強(qiáng)等離子體譜線反演得到隨時(shí)間變化的高精度電子密度剖面，希望可以通過(guò)進(jìn)一步的研究將這一方法應(yīng)用到電離層的參數(shù)估計(jì)中.

1 擬合方法

本文使用EISCAT的UHF雷達(dá)[2]（位于Troms?，69°N，19°E）在2015年10月31日測(cè)量得到的電離層數(shù)據(jù). 該雷達(dá)有Beata、Bella和Manda三種模式，文中使用Beata模式下的觀測(cè)數(shù)據(jù). 該模式的編碼方式為32位交替碼，碼元寬20 μs，采樣間隔為0.4 μs，脈沖重復(fù)周期為5.58 ms. Beata模式下的UHF雷達(dá)可以探測(cè)電離層50～700 km的等離子體參數(shù)，探測(cè)的時(shí)間和空間分辨率分別為5 s和3 km. 然而，由于設(shè)備的限制，雷達(dá)只能觀測(cè)到下移的高頻增強(qiáng)等離子體譜線，無(wú)法觀測(cè)到上移的等離子體譜線[23]，所以本文主要討論下移的增強(qiáng)等離子體譜線.

圖1為EISCAT的UHF雷達(dá)在2015年10月31日14:00—14:09UT觀測(cè)到的增強(qiáng)等離子體譜線，該時(shí)間段內(nèi)加熱機(jī)為關(guān)機(jī)狀態(tài). 從圖1中可以清晰地觀測(cè)到隨時(shí)間變化的下移增強(qiáng)等離子體譜線，對(duì)應(yīng)的頻率變化在?6 MHz左右. 由于雷達(dá)空間分辨率為3 km，為得到空間分辨率更高的等離子體譜線，本文將提取出不同高度對(duì)應(yīng)的等離子體譜線頻率值并進(jìn)行擬合處理. 以14:00:30UT時(shí)刻為例，圖2為EISCAT雷達(dá)在此時(shí)刻觀測(cè)到的增強(qiáng)等離子體譜線隨高度的變化. 在227.9～260.2 km高度上每3 km的頻譜中選取譜峰作為該高度等離子體譜線的頻率值，然后對(duì)選取的頻率值進(jìn)行擬合. 由圖2可以看到，該剖面的峰值對(duì)應(yīng)高度約為245.5 km，220～270 km的增強(qiáng)譜線的形狀類似于二次函數(shù)，因此我們使用二次曲線來(lái)擬合增強(qiáng)譜線隨高度的變化，擬合結(jié)果如圖2中白色虛線所示. 可以看到，白色虛線與增強(qiáng)的等離子體譜線在220～270 km基本重合，在235～260 km的擬合效果最好.

圖1 2015年10月31日14:00—14:09UT EISCAT ISR觀測(cè)到的209～383 km高度上的等離子體譜線Fig. 1 Undecoded plasma lines observed by the EISCAT ISR in altitude range 209?383 km at 14:00?14:09UT on October 31, 2015

圖2 2015年10月31日14:00:30UT高度200～300 km處的增強(qiáng)等離子體譜線Fig. 2 Plasma lines intensity profile obtained by the EISCAT ISR with altitude from 200 km to 300 km at 14:00:30UT on 31 October 2015, the white dashedline is the curve fitting with quadratic function

2 反饋結(jié)果

在磁化等離子體中，朗繆爾波的線性色散關(guān)系[24-25]為

式中：ωL是朗繆爾波頻率；是電子等離子體頻率，e、ne、me分別是電子電量、電子密度以及電子質(zhì)量，ε0是自由空間的電導(dǎo)率；k是朗繆爾波的波數(shù)，此處，它等于EISCAT UHF雷達(dá)的波數(shù)的兩倍，即是電子熱速度，kb是 玻爾茲曼常數(shù)，Te是 電子溫度；ωce=eB0/me是電子回旋頻率，B0是 地磁場(chǎng)強(qiáng)度；θ是雷達(dá)視線和地磁場(chǎng)B之間的夾角. 從方程(1)推導(dǎo)出電子密度為

在此次實(shí)驗(yàn)中，雷達(dá)的天線仰角等于Troms?的磁傾角78°，即雷達(dá)波束平行于地磁場(chǎng)，θ為0，因此方程(2)的第二項(xiàng)為0. 由于電子的熱速度vTe在計(jì)算電子密度時(shí)只會(huì)引起相對(duì)很小的誤差，例如，電子溫度為1000 K時(shí)，電子密度為6×1011m?3的誤差只有3.5%. 因此，方程(2)的第三項(xiàng)在電子密度的測(cè)量中可以忽略，于是電子密度的公式簡(jiǎn)化為

式中，朗繆爾波頻率 ωL由擬合出的等離子體譜線給出.

為檢驗(yàn)等離子體譜線反演結(jié)果的準(zhǔn)確性，我們對(duì)比了等離子體譜線和離子譜線的反演結(jié)果，如圖3所示. 由于目前已經(jīng)廣泛使用的離子譜線反演方法需要將離子譜線進(jìn)行一定時(shí)間的累積，才能得到較為準(zhǔn)確的電子密度，因此這里的離子譜線是經(jīng)過(guò)9 min（14:00—14:09UT）累積的離子譜線，而等離子體譜線為14:09UT時(shí)刻5 s的等離子體譜線. 可以看到，離子譜線與等離子體譜線反演的電子密度之間存在較小的波動(dòng)（在220～270 km，波動(dòng)小于5%），因此，可以認(rèn)為等離子體譜線反演的結(jié)果是準(zhǔn)確的.

圖3 14:09:00UT時(shí)刻的等離子體譜線和經(jīng)過(guò)9 min（14:00 —14:09UT）累積的離子譜線分析反演的電子密度剖面Fig. 3 Electron density profiles estimated by the plasma lines(solid blue line) at 14:09UT and the ion lines(solid red line) integrated over 9 minutes

圖3所示的反演結(jié)果也證明等離子體譜線反演方法具有更高的時(shí)間分辨率，即等離子體譜線只需要某個(gè)時(shí)刻5 s的數(shù)據(jù)即可反演出較為準(zhǔn)確的電子密度剖面，而離子譜線反演電子密度剖面則需要更長(zhǎng)的時(shí)間累積. 如果離子譜線用5 s的數(shù)據(jù)來(lái)反演，等離子體譜線反演方法的高時(shí)間分辨率優(yōu)勢(shì)則轉(zhuǎn)換為高精度優(yōu)勢(shì). 圖4給出了14:00:30UT時(shí)刻5 s數(shù)據(jù)下等離子體譜線和離子譜線反演的帶有誤差條（error bars）的電子密度剖面. 紅色的誤差條長(zhǎng)度代表GUISDAP的離子譜線分析得到的擬合誤差，而藍(lán)色的誤差條長(zhǎng)度是由二次曲線擬合方法得到的該時(shí)刻處關(guān)于頻率的預(yù)測(cè)標(biāo)準(zhǔn)誤差.可以看到離子譜線反演的電子密度剖面有較大的起伏，而等離子體譜線得到的電子密度結(jié)構(gòu)非常清晰和光滑，并且后者的誤差遠(yuǎn)小于前者. 兩種反演方法的效果不同主要是由于其反演原理上的差異，離子譜線給出的電子密度取決于離子譜線的強(qiáng)度，然而在實(shí)際測(cè)量中，由于非相干積累時(shí)間很短，導(dǎo)致信噪比低，故反演出的電子密度剖面效果不佳，如果想反演更好的電子密度剖面，必須增大非相干積累時(shí)間. 相比之下，等離子體譜線反演電子密度取決于方程(3)中朗繆爾波的頻率，誤差產(chǎn)生于頻點(diǎn)的選取過(guò)程. 此外，等離子體譜線具有較高的頻率分辨率（1.56 kHz），可以得到更好的結(jié)果.因此，增強(qiáng)的等離子體譜線反演電子密度剖面具有高精度或者高時(shí)間分辨率的優(yōu)點(diǎn).

圖4 14:00:30UT時(shí)刻等離子體譜線和離子譜線分析反演的帶有誤差條的電子密度剖面Fig. 4 Electron density profiles estimated by the plasma lines(solid blue line) and the ion lines (solid red line)with error bars at 14:00:30UT

圖5給出了從2015年10月31日14:00—14:09UT的譜線分析中得到的205～295 km的電子密度剖面隨時(shí)間的變化. 圖中時(shí)空分辨率分別為5 s和3 km.如圖5所示，離子譜線反演的電子密度噪點(diǎn)很多，而等離子體譜線反演出的電子密度清晰平滑. 從圖5下圖可以看出，電子密度的強(qiáng)度在240～260 km高度上隨時(shí)間減弱. 例如，在250 km處，14:00:05UT的電子密度為4.773×1011m?3，在14:09UT時(shí)減少為4.18×1011m?3. 在這9 min內(nèi)，電子密度強(qiáng)度減少了12.42%.

圖5 2015年10月31日14:00—14:09UT利用離子譜線（上）和等離子體譜線（下）分析反演出的電子密度剖面Fig. 5 Electron density estimated by ion lines (top) and plasma lines (bottom) with temporal resolution of 5 s and altitude resolution of 3 km from 205 km to 295 km on 31 October 2015 (14:00?14:09UT)

3 結(jié)果分析

等離子體譜線是朗繆爾波對(duì)ISR的散射回波[12]，在其測(cè)量中，起決定性作用的是散射效率以及朗繆爾波的激勵(lì)和阻尼[6]. 在一定條件下朗繆爾波的激發(fā)和阻尼會(huì)對(duì)等離子體譜線產(chǎn)生足夠的影響，使其增強(qiáng)到可以被雷達(dá)觀測(cè)到. 在不同的緯度上，相應(yīng)的超熱電子會(huì)通過(guò)朗道阻尼激發(fā)朗繆爾波從而增強(qiáng)等離子體譜線[14]. 在中低緯電離層，主要考慮光電子的作用[26-27]，而在高緯電離層，超熱電子主要由次級(jí)電子和來(lái)自磁層的電子束產(chǎn)生[6,15]，這兩種源產(chǎn)生的超熱電子分別在E層和F層高度上增強(qiáng)等離子體譜線. 本文中選取的增強(qiáng)等離子體譜線產(chǎn)生高度位于F層（如圖2所示），且不在電離層加熱實(shí)驗(yàn)期間，所以主要考慮電子束的作用. 圖6是2015年10月31日13:00—16:00UT的AE指數(shù)和SYM-H指數(shù)，其中，虛線框內(nèi)是14:00—14:09UT的指數(shù)變化. 在13:00—14:00 UT期間，AE指數(shù)的值為150 ～250 nT，有一個(gè)較小的增強(qiáng)；而SYM-H指數(shù)由10 nT減少至5 nT，變化緩慢，但不足以說(shuō)明此時(shí)沒(méi)有電子束的作用. 圖7是2015年10月31日國(guó)防氣象衛(wèi)星計(jì)劃（defense meteorological satellite program, DMSP）F16衛(wèi)星探測(cè)的粒子能譜，從上到下分別為：電子與離子的總能通量（eV·(cm2·s·sr)?1）、平均能量（eV）的變化以及電子和離子的能譜，該衛(wèi)星的運(yùn)行高度距地面約為848 km.在13:43UT時(shí)刻，該衛(wèi)星經(jīng)過(guò)Troms?，可以看到此時(shí)圖中相應(yīng)區(qū)域具有較大的粒子能量（電子總能通量大于1010eV·(cm2·s·sr)?1），與AE指數(shù)有相同的趨勢(shì)，說(shuō)明此時(shí)有來(lái)自磁層的粒子沉降至F層. 雖然在14:00UT時(shí)刻沒(méi)有衛(wèi)星觀測(cè)數(shù)據(jù)，但此時(shí)的增強(qiáng)等離子體譜線仍可能是電子束的作用，這一推論還需更加精確的數(shù)據(jù)來(lái)證明.

圖6 2015年10月31日13:00—16:00UT的AE指數(shù)（上）和SYM-H指數(shù)（下）Fig. 6 AE index (top) and SYM-H index (bottom) during 13:00?16:00UT on 31 October 2015, where the dot lines indicate the time range of 14:00?14:09UT

圖7 2015年10月31日13:40—14:01UT DMSP探測(cè)的粒子能譜Fig. 7 Spectrum obtained by DMSP during 13:40?14:01UT on 31 October 2015

4 結(jié) 論

利用等離子體譜線的頻偏獲得電離層電子密度是目前最精確的地面無(wú)線電遙測(cè)方法. 本文通過(guò)曲線擬合的方法，提取出了位于Troms?的 EISCAT UHF雷達(dá)所觀測(cè)到的持續(xù)時(shí)長(zhǎng)為9 min的增強(qiáng)等離子體譜線，然后利用朗繆爾波色散關(guān)系反演得到電子密度剖面. 與離子譜線計(jì)算的電子密度剖面相比，這種方法得到的電子密度剖面具有高時(shí)空分辨率的特點(diǎn). 對(duì)于增強(qiáng)等離子體譜線具體機(jī)制的討論，還需要大量的實(shí)驗(yàn)數(shù)據(jù)分析和理論研究，可以作為下一步的研究工作. 未來(lái)還可以進(jìn)一步開(kāi)發(fā)自主的等離子體譜線測(cè)量及檢測(cè)技術(shù)，并應(yīng)用于國(guó)內(nèi)新的ISR.

致謝：感謝EISCAT觀測(cè)站提供的數(shù)據(jù). 感謝中國(guó)國(guó)家自然科學(xué)基金會(huì)（ 42074187，41774162和41704155）及中國(guó)博士后科研基金（2017M622504）對(duì)本工作的支持. 感謝中國(guó)極地研究中心胡澤駿研究員提供的DMSP數(shù)據(jù). 感謝評(píng)審專家的建議及編輯的辛勤勞動(dòng).