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

    同步軌道相對論電子通量長期倒空事件的統(tǒng)計(jì)研究

    2021-06-02 10:42:26吳晗陳濤KalegaevPanasyuk段素平何兆海李仁康
    地球物理學(xué)報(bào) 2021年6期

    吳晗, 陳濤, V. V. Kalegaev, M. I. Panasyuk,段素平, 何兆海, 李仁康

    1 中國科學(xué)院國家空間科學(xué)中心, 空間天氣學(xué)國家重點(diǎn)實(shí)驗(yàn)室, 北京 100190 2 中國科學(xué)院大學(xué), 北京 100049 3 莫斯科羅蒙諾索夫國立大學(xué), 核物理研究所, 俄羅斯 莫斯科 119234 4 云南師范大學(xué)物理與電子信息學(xué)院, 昆明 650500

    0 引言

    地球輻射帶是圍繞在地球周圍的充滿著大量被捕獲的高能帶電粒子區(qū)域,是許多衛(wèi)星飛行軌道所在區(qū)域,也是探索系外行星航天器飛離地球的必經(jīng)區(qū)域.其中,外輻射帶位于L=3~7之間,其內(nèi)部被捕獲的高能電子易受太陽風(fēng)、宇宙線等影響,處于不斷的動態(tài)變化中(Baker et al., 2018).MeV能量的高能電子也被稱為“殺手電子”,當(dāng)其通量增強(qiáng)到一定程度時(shí),會對在軌衛(wèi)星的正常工作和安全造成巨大威脅.外輻射帶中存在著一類特殊的物理現(xiàn)象,就是電子通量倒空現(xiàn)象(dropout),其特征是在一個(gè)廣闊的能段、投擲角和L范圍內(nèi),被捕獲的電子群通量顯著下降的現(xiàn)象(Turner et al., 2012a).這類現(xiàn)象發(fā)生時(shí)高能電子通量會降至極低水平,當(dāng)持續(xù)時(shí)間較長時(shí),可以為衛(wèi)星的安全和長時(shí)間艙外作業(yè)提供有效的保障(Wrenn, 1995).揭示這類現(xiàn)象的發(fā)生規(guī)律既是空間物理學(xué)研究的熱點(diǎn),也可以加深對電子在輻射帶中加速、損失機(jī)制的認(rèn)識和理解,為預(yù)報(bào)空間天氣提供強(qiáng)有力的科學(xué)支持.

    輻射帶內(nèi)相對論電子通量具有高度的可變性,可以在從幾分鐘到幾年的時(shí)間尺度上發(fā)生數(shù)量級的變化.長期的變化一般與太陽活動周不同階段太陽活動水平高低的影響有關(guān),短期的變化則由太陽風(fēng)和行星際條件決定(Hajra et al., 2014; Hajra and Tsurutani, 2018).目前已有許多工作對輻射帶電子長時(shí)間通量倒空現(xiàn)象進(jìn)行了研究.Lee等(2013)通過分析2009年外輻射帶電子通量變化發(fā)現(xiàn):外輻射帶在這一年當(dāng)中整體比較平靜,發(fā)生通量(幾百keV至MeV電子)倒空現(xiàn)象數(shù)次,最長持續(xù)時(shí)間高達(dá)2個(gè)月.這類現(xiàn)象與極弱的太陽風(fēng)條件有關(guān),在太陽活動低年,長時(shí)間平靜的太陽風(fēng)條件大大地減少磁層對流和注入過程,導(dǎo)致等離子層膨脹,甚至可以擴(kuò)大到地球同步軌道高度.等離子層內(nèi)的波動,如等離子體層嘶聲可以引起電子投擲角散射,使電子最終沉降到地球大氣從輻射帶中永久損失(Breneman et al., 2015).Pinto等(2018)通過時(shí)間序列疊加法對發(fā)生在1996—2006年間的21例相對論電子通量長期倒空事件分析得出:低速太陽風(fēng)、突然增加的質(zhì)子密度以及事件開始后緊隨的IMFBz分量平均值由南向轉(zhuǎn)為北向的過程是這一類事件的普遍特征.

    地磁暴發(fā)生期間,外輻射帶內(nèi)高能電子通量通常會發(fā)生變化.典型的變化是:相對論電子通量會隨著地磁活動的增強(qiáng),在磁暴主相期間下降,隨后在磁暴恢復(fù)相階段增長,在伴隨著高強(qiáng)度持續(xù)性亞暴活動的磁暴恢復(fù)相期間相對論電子通量甚至可能遠(yuǎn)超暴前水平(李柳元等,2006).Onsager等(2002)通過分析一次中等磁暴期間電子通量的變化發(fā)現(xiàn),在磁暴結(jié)束后,只有能量較低的高能電子(E<300 keV)可以恢復(fù)到暴前水平,而較高能量的相對論電子(>2 MeV)并沒有恢復(fù)到暴前水平.在外輻射帶相對論電子通量對磁暴的響應(yīng)這一研究領(lǐng)域中,比較有代表性的是Reeves等(2003)關(guān)于1989—2000年間對276個(gè)中、大磁暴前后,外輻射帶相對論電子通量變化的研究工作.其研究結(jié)果表明:有53%的磁暴會使同步軌道相對論電子通量增加,19%的磁暴會減少同步軌道相對論電子通量,而28%的磁暴對同步軌道相對論電子通量并沒有產(chǎn)生顯著的變化.近年來,Anderson等(2015)采用和Reeves等(2003)類似的方法,同樣研究了1989—2000年間磁暴前后相對論電子通量的變化情況.不同的是,他們的研究對象是342個(gè)小磁暴事件,但也得出了相似的結(jié)論,即小磁暴和中、大磁暴一樣可以使相對論電子通量增加、減少或者不變.由此可見,磁暴發(fā)生與否不能完全決定暴后相對論電子通量的走向,應(yīng)結(jié)合更多的參數(shù)綜合分析.

    以往的研究工作更多地集中在相對論電子通量產(chǎn)生變化的動態(tài)過程,即關(guān)注點(diǎn)在相對論電子通量顯著增強(qiáng)或者損失的階段,一些研究工作也會采用時(shí)間序列疊加法來進(jìn)行統(tǒng)計(jì)研究,或者對單例事件分析其相關(guān)的物理機(jī)制等(Green and Kivelson, 2004; Borovsky and Denton, 2009; Ni et al., 2016; Pinto et al., 2018).與前人工作不同的是,本文對相對論電子通量長期保持倒空狀態(tài)這一通量穩(wěn)定不變的階段研究其統(tǒng)計(jì)特征.首先通過對約1.5個(gè)太陽活動周內(nèi)共62例相對論電子通量長期倒空事件的持續(xù)時(shí)間和分布情況,研究其對太陽活動水平的響應(yīng).在此基礎(chǔ)上,對這62例事件進(jìn)行時(shí)間序列疊加法處理,分析開始和結(jié)束階段的普遍特征.同時(shí),也對這些相對論電子通量長時(shí)間維持在背景通量水平期間,太陽風(fēng)條件和地磁活動水平以及一些其他參數(shù)的特征進(jìn)行不同類別事件的對比研究.最后,本文結(jié)合相關(guān)的物理機(jī)制對這類事件特征的統(tǒng)計(jì)結(jié)果做出可能的、合理的物理解釋.

    1 數(shù)據(jù)來源與事件定義

    本文采用來自O(shè)mni數(shù)據(jù)庫5 min分辨率數(shù)據(jù)來研究太陽風(fēng)參數(shù)的變化情況(包括IMFBz分量、太陽風(fēng)密度、速度和動壓),和來自WDC(World Data Center)1 h分辨率數(shù)據(jù)來判斷地磁干擾強(qiáng)度(亞暴AE指數(shù)和磁暴Dst指數(shù)).文中使用來自SILSO(Sunspot Index and Long-term Solar Observations)的太陽黑子數(shù)據(jù)來衡量太陽活動水平的高低.

    ULF波是一種全球持續(xù)性地磁脈動,在磁層能量、質(zhì)量和動量的輸運(yùn)過程中起著重要的作用(Zong et al., 2017).ULF波動尤其是與輻射帶高能電子漂移運(yùn)動周期相近的Pc5頻段波動和粒子的相互作用一直被認(rèn)為是輻射帶電子加速的主要機(jī)制之一(Rostoker et al., 1998;李柳元等,2005;宗秋剛等,2011,2013).Regi(2016)通過分析2006—2010年間ULF波在Pc5頻段波動功率的漲落發(fā)現(xiàn),ULF波動功率在約2天左右的時(shí)間延遲上與地球同步軌道高能電子通量變化有強(qiáng)烈的相關(guān)性.Friedel等(2002)也展示相關(guān)的統(tǒng)計(jì)結(jié)果,表明磁暴期間至少持續(xù)12小時(shí)足夠強(qiáng)的ULF波功率,是隨后地球同步軌道相對論電子增強(qiáng)的最有效指標(biāo).此外,與分布范圍有限并且不長期存在的波動如合聲波、Electromagnetic Ion Cyclotron(EMIC)波等相比,ULF波動有全球廣泛分布的特性和可持續(xù)觀測的優(yōu)點(diǎn),利于統(tǒng)計(jì).因此,本文引入Kozyreva等(2007)和Pilipenko等(2017)介紹的ULF波指數(shù)(Pc5頻段),來研究ULF波強(qiáng)度可能對相對論電子通量產(chǎn)生的影響.為了更好地對應(yīng)全球ULF波動強(qiáng)度,避免受衛(wèi)星磁地方時(shí)等因素的影響,本文僅采用地面ULF波指數(shù),即從北半球幾十個(gè)地磁臺站數(shù)據(jù)重新計(jì)算的代表全球ULF波活動強(qiáng)度的指數(shù).此指數(shù)相比其他如同步軌道ULF波指數(shù)和行星際ULF波指數(shù)能更好地反映電子通量的變化(Pilipenko et al., 2017).

    受不同衛(wèi)星發(fā)射、運(yùn)行時(shí)間的限制,也為了得到更多有效的長期倒空事件,本文主要研究相對論電子通量在地球同步軌道上的變化,并采用地球同步軌道衛(wèi)星GOES (Geostationary Operational Environmental Satellite)對>2 MeV能段電子通量的觀測結(jié)果.其中,我們使用GOES-10,GOES-11和GOES-13衛(wèi)星數(shù)據(jù),分別統(tǒng)計(jì)2000—2005年,2006—2010年和2011—2016年相對論電子通量長期倒空事件.受Pinto等(2018)工作的啟發(fā),我們定義相對論電子通量長期倒空事件為日均通量下降至30 cm-2·sr-2·s-1以下,且持續(xù)時(shí)間達(dá)3天以上,通量無明顯隨磁地方時(shí)變化的現(xiàn)象.與Pinto等(2018)定義日均通量下降至10 cm-2·sr-2·s-1不同,我們放寬這個(gè)閾值至30 cm-2·sr-2·s-1,是因?yàn)镻into等(2018)采用5 min原始數(shù)據(jù)統(tǒng)計(jì)2007年前的相對論電子通量長期倒空事件,而我們?yōu)榱耸褂肎OES-13優(yōu)化過的1 min科學(xué)數(shù)據(jù),在文中采用1 min原始數(shù)據(jù)進(jìn)行統(tǒng)計(jì).受儀器精度、計(jì)數(shù)等的影響,2000—2010年間,5 min原始數(shù)據(jù)計(jì)算得到日均通量明顯較1 min原始數(shù)據(jù)計(jì)算得到日均通量低(約低10~20 cm-2·sr-2·s-1);而2011—2016年間,在使用GOES-13衛(wèi)星數(shù)據(jù)統(tǒng)計(jì)時(shí),二者相差不大,并且在此期間無論是5 min原始數(shù)據(jù)還是1 min原始數(shù)據(jù)計(jì)算得到的日均通量都沒有小于10 cm-2·sr-2·s-1的情況.因此,我們在Pinto等(2018)工作的基礎(chǔ)上做了優(yōu)化,并統(tǒng)計(jì)了跨度約1.5個(gè)太陽活動周更長時(shí)間的事例.在這種篩選條件下,相對論電子通量降低至衛(wèi)星背景觀測通量水平,即相對論電子近乎完全從地球同步軌道上損失.

    圖1 相對論電子通量長期倒空事件的示意圖Fig.1 The identification of a long-term relativistic electron flux dropout

    圖1為相對論電子通量長期倒空事件示意圖.圖中是2005年10月19日03 ∶ 00至10月28日16 ∶ 00期間>2 MeV電子通量的變化,紅色部分則為我們所選取的相對論電子通量長期倒空事件.我們以通量開始下降后達(dá)到通量最小值的時(shí)間點(diǎn)為事件開始的時(shí)間點(diǎn),以通量開始顯著增加并在隨后的2日內(nèi)有明顯日變化現(xiàn)象為事件結(jié)束的時(shí)間點(diǎn).在排除衛(wèi)星變軌等原因?qū)е码娮油繑?shù)據(jù)缺失和某些缺乏連續(xù)的、有效的太陽風(fēng)參數(shù)觀測數(shù)據(jù)的事件后,我們在2000—2016年間共挑選出62例相對論電子通量長期倒空事件,其具體的起止時(shí)間和持續(xù)時(shí)間見表1.

    表1 2000—2016年間地球同步軌道相對論電子通量長期倒空事件列表Table 1 List of long-term dropouts of relativistic electron flux in geosynchronous orbit from 2000 to 2016

    續(xù)表1

    2 相對論電子通量長期倒空事件在太陽活動周內(nèi)的年分布特征

    圖2所示為表1中列出的62例相對論電子通量長期倒空事件在不同年份的分布情況.圖中紅色柱狀圖為每年發(fā)生的長期倒空事件持續(xù)時(shí)間總和.這里對于起止時(shí)間分布在不同年份的事件,則按照其在不同年份實(shí)際持續(xù)的時(shí)間分別統(tǒng)計(jì).圖中柱狀圖上方的黑色數(shù)字表示每年發(fā)生事件的頻次,以事件開始的時(shí)間做統(tǒng)計(jì).藍(lán)色折線是以太陽黑子數(shù)年平均值表示的太陽活動年平均水平的變化曲線,其中2000—2007年屬于第23太陽活動周,2008—2016年屬于第24太陽活動周.本文共統(tǒng)計(jì)從第23太陽活動周峰年至第24太陽活動周下降期約1.5個(gè)太陽活動周期共17年內(nèi)的相對論電子通量長期倒空事件.

    從圖中可以看出,第23太陽活動周下降期的2003—2007年與第24太陽活動周上升期的2010—2014年相比,總體上在同等太陽活動水平條件下,第23太陽活動周下降期的相對論電子通量長期倒空事件持續(xù)時(shí)間較短,發(fā)生次數(shù)也較少.這可能是因?yàn)樵谝粋€(gè)太陽活動周的下降期,具有較多的HILDCAA(high-intensity long-duration continuous auroral activity)事件發(fā)生,這種持續(xù)性的高強(qiáng)度亞暴活動有助于相對論電子通量增強(qiáng)(Hajra et al., 2014, 2015; Hajra and Tsurutani, 2018),因此我們統(tǒng)計(jì)的長期倒空事件就會減少甚至沒有發(fā)生.

    2009年是我們觀測到長期倒空事件總持續(xù)時(shí)間最長的一年,累計(jì)持續(xù)時(shí)間共約5000小時(shí).這表示在2009年這一整年內(nèi),地球同步軌道上有一半以上的時(shí)間都處于相對論電子通量倒空狀態(tài),具體的原因Lee等(2013)已經(jīng)做出了詳細(xì)的物理解釋,并且我們在引言部分也介紹了他們的工作,故在此不做過多的解釋說明.

    總持續(xù)時(shí)間緊隨2009年,并且發(fā)生頻次最高的是2014年.這一年雖然是第24太陽活動周峰年,但我們?nèi)匀挥^測到許多長期倒空事件,這說明太陽活動水平較高不代表沒有長期倒空事件發(fā)生;與之相對的,在太陽活動水平較低的2008年,我們并沒有觀測到類似2009年長時(shí)間的通量倒空現(xiàn)象.綜上所述,在太陽活動周的下降期,相對論電子通量倒空現(xiàn)象較少發(fā)生;但在太陽活動周峰年和谷年,太陽活動水平的高低并沒有完全決定是否容易出現(xiàn)相對論電子通量倒空現(xiàn)象.

    3 相對論電子通量長期倒空事件產(chǎn)生條件和結(jié)束條件的統(tǒng)計(jì)特征

    這一節(jié),我們簡要對前述的62例事件開始和結(jié)束階段的IMFBz分量、太陽風(fēng)速度、密度、動壓、向陽面磁層頂高度MP(magnetopause standoff distance)、等離子體層頂高度LPP(plasmapause position)、亞暴AE指數(shù)、磁暴環(huán)電流Dst指數(shù)和ULF波指數(shù)進(jìn)行統(tǒng)計(jì)研究,發(fā)現(xiàn)并討論其開始、結(jié)束階段的普遍特征以及可能導(dǎo)致該特征產(chǎn)生的相關(guān)物理機(jī)制.其中,采用Shue等(1997)的方法計(jì)算向陽面磁層頂高度;采用Liu等(2015)的方法計(jì)算等離子體層頂高度.由于Liu等(2015)的方法無法重現(xiàn)等離子體層羽狀結(jié)構(gòu),在昏側(cè)的誤差較大,我們使用24個(gè)MLT區(qū)域的平均值來反映等離子體層頂?shù)恼w變化.經(jīng)驗(yàn)證,在MLT=6、12和24處等離子體層頂?shù)淖兓厔菖c所有MLT區(qū)域等離子體層頂平均值的變化趨勢保持一致,區(qū)別在于:當(dāng)MLT=6和12時(shí),等離子體層頂?shù)母叨缺绕骄德缘?;?dāng)MLT=24時(shí),等離子體層頂?shù)母叨缺绕骄德愿?

    圖2 相對論電子通量長期倒空事件在2000—2016年間的年分布特征Fig.2 The annual distribution of long-term dropouts of relativistic electron flux from 2000 to 2016

    我們采用相對論電子通量下降到最低值時(shí)為事件開始階段時(shí)間序列疊加的時(shí)間0點(diǎn),以通量從背景水平開始增長時(shí)為事件結(jié)束階段時(shí)間序列疊加的時(shí)間0點(diǎn),并分別統(tǒng)計(jì)兩個(gè)階段各種參數(shù)的變化情況,具體如圖3所示.圖中藍(lán)色曲線為上下四分位點(diǎn),紅色曲線為中位數(shù).此外,我們也檢驗(yàn)了平均值的變化,但受某些極高值的影響,部分參數(shù)的平均值會高于上四分位點(diǎn),并不能很好地反映整體的變化水平,因此本文僅展示中位數(shù)的變化.

    從圖3a的相對論電子通量倒空事件開始階段的特征圖可以看出,太陽風(fēng)速度一直維持在300~400 km·s-1之間,僅在事件開始前幾小時(shí)略微增加,但整體變化比較平穩(wěn).太陽風(fēng)密度和動壓在事件發(fā)生前呈現(xiàn)顯著增加趨勢,中位數(shù)最大值出現(xiàn)在事件開始前3~4 h附近并且約是事件發(fā)生前3天數(shù)值的兩倍.這樣高值的太陽風(fēng)密度和動壓一直持續(xù)到事件發(fā)生后半日才恢復(fù)到平穩(wěn)的狀態(tài).與之相對應(yīng)的向陽面磁層頂高度在太陽風(fēng)動壓顯著增加時(shí)向靠近地球的方向運(yùn)動,其中位數(shù)在事件開始時(shí)達(dá)到最小值,由事件開始前3天約11Re高度收縮至10Re以下.

    GSM坐標(biāo)系下IMFBz分量的中位數(shù)在事件開始前數(shù)小時(shí)有明顯北向轉(zhuǎn)南向的過程,在事件開始后數(shù)小時(shí)內(nèi)Bz分量逐漸由南向轉(zhuǎn)到北向并在此后幾乎一直維持在北向.與之相對應(yīng)的,Dst中位數(shù)在IMFBz分量達(dá)到南向最大值時(shí)呈現(xiàn)明顯的下降趨勢,Dst極小值出現(xiàn)在事件開始后約10 h左右.此外,在事件開始這一階段,AE指數(shù)和ULF波指數(shù)僅在事件開始前后數(shù)小時(shí)內(nèi)有明顯的增加,但整體增長的幅度并不高.

    在以上太陽風(fēng)活動水平和地磁活動水平條件下,等離子體層頂高度在事件發(fā)生前幾乎一直維持在高L區(qū)域,其中位數(shù)高于6Re,部分時(shí)刻上四分位點(diǎn)等離子體層頂位置甚至超過地球同步軌道高度.這種狀態(tài)一直持續(xù)到事件發(fā)生前數(shù)小時(shí),等離子體層頂才逐漸向內(nèi)收縮,并在事件開始后幾小時(shí)又開始緩慢向外擴(kuò)展.雖然本文相對論電子通量長期倒空事件與Pinto等(2018)的選取標(biāo)準(zhǔn)略有不同并且統(tǒng)計(jì)的事件數(shù)也更多,但得出了與Pinto等(2018)相似的結(jié)論,即事件開始時(shí)顯著增加的太陽風(fēng)動壓和密度,較弱的太陽風(fēng)速度和IMFBz分量中位數(shù)在事件發(fā)生前轉(zhuǎn)南向隨后一直維持在北向.我們在前人工作的基礎(chǔ)上增加了更多的參數(shù),這些參數(shù)在事件開始、結(jié)束階段呈現(xiàn)不同的變化,為研究相對論電子通量長期倒空現(xiàn)象提供了更多相關(guān)的統(tǒng)計(jì)理論支持,也可以為未來復(fù)雜的輻射帶建模等提供一些參考依據(jù).

    圖3 基于時(shí)間序列疊加法的62例相對論電子通量長期倒空事件開始(a)和結(jié)束階段(b)的不同特征Fig.3 Superposed epoch analysis results of the start time (a) and end time (b) of all 62 long-term dropouts of relativistic electron fluxes

    “磁層頂陰影”(magnetopause shadowing)是輻射帶電子真實(shí)損失機(jī)制的一種.當(dāng)太陽風(fēng)動壓增加時(shí),向內(nèi)壓縮的磁層頂會導(dǎo)致一些原來處于閉合的漂移軌道上的粒子運(yùn)動軌跡發(fā)生改變,沿著開放的磁力線逃逸并從輻射帶中大量流失(Shprits et al., 2006; Turner et al., 2012b; Kang et al., 2018).本文統(tǒng)計(jì)的相對論電子通量長期倒空事件開始前幾小時(shí),明顯增加的太陽風(fēng)動壓和密度引起向陽面磁層頂向靠近地球方向運(yùn)動,說明了某些事例相對論電子通量損失的原因與電子向外的徑向擴(kuò)散運(yùn)動并最終在壓縮的磁層頂處永久損失有關(guān).

    IMFBz分量南向時(shí),磁層頂重聯(lián),大量粒子注入使環(huán)電流強(qiáng)度增強(qiáng),當(dāng)增強(qiáng)的環(huán)電流粒子產(chǎn)生的磁場與地磁場疊加時(shí),地磁水平分量減小,磁暴發(fā)生.圖3a中IMFBz分量南向達(dá)到最大值后,Dst指數(shù)的下降則反映了這一點(diǎn).在磁暴主相期間,環(huán)電流的增強(qiáng)導(dǎo)致地磁場減小,電子為了保持第三絕熱不變量守恒會向著遠(yuǎn)離地心的方向徑向運(yùn)動.當(dāng)這些電子運(yùn)動到較弱的磁場時(shí),為了保持第一絕熱不變量守恒,電子的能量會降低.因此,指定能段電子在固定的位置會觀測到通量減小的現(xiàn)象(Kim and Chan, 1997).Kim 和 Chan(1997)通過計(jì)算確定磁暴期間三種絕熱不變量都保持守恒,證明了這個(gè)過程是一種完全的絕熱作用.本文的研究結(jié)果顯示,事件發(fā)生前Dst指數(shù)的下降表明某些事例與磁暴有關(guān),說明這些事例部分電子通量的降低可能是由“Dst效應(yīng)”造成的.

    但是,“Dst效應(yīng)”是一種完全的絕熱作用,相對論電子通量的降低不代表其真實(shí)地從輻射帶中損失.而磁暴發(fā)生期間,溫度各向異性的環(huán)電流離子的注入利于EMIC波的激發(fā),EMIC波與高能電子相互作用使其沉降到地球大氣進(jìn)而從輻射帶中永久損失是輻射帶快速損失相對論電子的主要機(jī)制(Millan and Thorne, 2007; Summers et al., 2007; Shprits et al., 2008).EMIC波主要分布在等離子體層內(nèi)昏側(cè)層頂附近和等離子體羽狀結(jié)構(gòu)中(Summers et al., 2007).地磁水平活躍時(shí),EMIC波在12—18 MLT區(qū)域最為常見(Meredith et al., 2014).本文對相對論電子通量長期倒空事件開始和結(jié)束時(shí)間0點(diǎn)的MLT分布進(jìn)行了統(tǒng)計(jì),如圖4所示.從圖4a中可以看出,事件開始時(shí)刻主要集中在午后至子夜MLT=14—22區(qū)域,這一區(qū)域剛經(jīng)過正午向陽面磁層頂被顯著壓縮的位置,也覆蓋了EMIC波的主要分布區(qū)域.Green等(2004)統(tǒng)計(jì)了同步軌道上相對論電子通量快速下降的區(qū)域,也呈現(xiàn)出在昏側(cè)分布最為顯著的特征.Xiang等(2017)列舉了三類輻射帶電子通量倒空事件,其主要的損失機(jī)制可以分為三種:磁層頂陰影占主導(dǎo)地位、EMIC波散射占主導(dǎo)地位以及二者綜合作用的損失機(jī)制.隨后他們利用Van Allen Probes上的相空間密度PSD(phase space density)數(shù)據(jù)統(tǒng)計(jì)了4年間輻射帶電子通量倒空事件,結(jié)果顯示EMIC波散射電子占主導(dǎo)地位的電子通量損失機(jī)制主要發(fā)生在低L*區(qū)域,并通常伴有活躍的地磁活動水平和太陽風(fēng)條件;在略高的L*區(qū)域,輻射帶電子通量損失更有可能是EMIC波散射和向外的徑向擴(kuò)散綜合作用導(dǎo)致的(Xiang et al., 2018).圖3a中,等離子體層頂在事件開始前幾乎一直維持在高L區(qū)域,部分時(shí)刻上四分位點(diǎn)位置可以超過地球同步軌道高度,這表明在某些事件,EMIC波對相對論電子的散射損失作用可能發(fā)生在高L區(qū)域甚至地球同步軌道高度.結(jié)合磁層頂高度下四分位點(diǎn)位置并沒有低至地球同步軌道高度,以及主要分布在午后至子夜MLT區(qū)域的事件開始0點(diǎn)(圖4a),表明地球同步軌道高度上,相對論電子的損失更有可能是向外的徑向擴(kuò)散作用和EMIC波在高L區(qū)域的投擲角散射引起的,但具體哪個(gè)機(jī)制占主導(dǎo)地位,需要具體事例具體分析.

    圖3b所示為相對論電子通量長期倒空現(xiàn)象結(jié)束階段的特征圖.與開始階段不同的是,在事件結(jié)束前后三天,并沒有觀察到明顯的IMFBz分量反轉(zhuǎn),也沒有觀察到Dst指數(shù)大幅下降的現(xiàn)象.這說明相對論電子通量結(jié)束長期倒空狀態(tài)時(shí),磁暴的發(fā)生不是必要條件.此外,盡管太陽風(fēng)動壓和密度的上四分位點(diǎn)在事件結(jié)束前1天有略微的增強(qiáng),但其中位數(shù)并沒有顯著改變;磁層頂高度的下四分位點(diǎn)在太陽風(fēng)動壓和密度略微增強(qiáng)時(shí),向靠近地球方向移動,但移動的幅度并不大,并且很快就恢復(fù)到原來的水平.

    事件結(jié)束階段的一個(gè)顯著特點(diǎn)就是太陽風(fēng)速度的增加,其中位數(shù)由事件結(jié)束前3天約350 km·s-1增加到事件結(jié)束后3天400 km·s-1以上,速度的上四分位點(diǎn)在事件結(jié)束前幾小時(shí)有著明顯的上升趨勢,并一直保持較高水平直到事件結(jié)束后2天才有輕微的下降.事件結(jié)束階段的另一個(gè)特點(diǎn)就是,長時(shí)間的亞暴活動和ULF波活動.AE指數(shù)和ULF波指數(shù)在事件結(jié)束前1天就開始增強(qiáng),并在事件結(jié)束后半日內(nèi)仍然維持在較高水平.盡管AE指數(shù)的中位數(shù)在事件開始階段的最大值(圖3a,約200 nT)和事件結(jié)束階段的最大值(圖3b,約200 nT)近乎相等,但在事件結(jié)束階段AE指數(shù)有著近1天半的時(shí)間都處于較高水平.而ULF波指數(shù)中位數(shù)在事件結(jié)束階段的最大值(圖3b,約20 nT),比事件開始階段的最大值(圖3a,約15 nT)明顯高.等離子體層頂高度在事件結(jié)束前1天就開始向內(nèi)收縮,隨后一直維持在較低高度上,直到事件結(jié)束后3天,其上四分位點(diǎn)也幾乎沒有超過6Re高度.

    Iles等(2002)通過統(tǒng)計(jì)相對論電子通量對磁暴響應(yīng)不同的三類事件,得出在磁暴恢復(fù)相期間,高速的太陽風(fēng)速度和IMFBz分量在0附近波動或主要是南向分布時(shí),是相對論電子顯著增強(qiáng)的重要條件.類似的,O′Brien等(2001)也提到了磁暴后相對論電子通量的增加,其實(shí)更取決于高速的太陽風(fēng)和恢復(fù)相期間長時(shí)間的ULF波活動.本文的研究結(jié)果也印證了高速太陽風(fēng)會打破相對論電子通量倒空現(xiàn)象,即相對論電子通量的增加,正如圖3b所示.

    圖4 相對論電子通量長期倒空事件開始(a)和結(jié)束(b)時(shí)刻MLT分布特征Fig.4 MLT distributions of the start time (a) and end time (b) of relativistic electron flux long-term dropouts

    亞暴活動也是相對論電子通量增加的一個(gè)重要原因.亞暴注入的源粒子群(數(shù)十keV)有助于合聲波的增長,當(dāng)種子粒子(數(shù)百keV)與合聲波相互作用時(shí),可以被加速到更高能量甚至相對論電子能量水平(Jaynes et al., 2015).一些強(qiáng)亞暴過程甚至可以直接向輻射帶注入MeV能量的電子(Dai et al., 2014, 2015).合聲波分布在等離子體層外側(cè)午夜至下午磁地方時(shí)區(qū)域(Summers et al., 2007).從圖4b中可以看出,事件結(jié)束時(shí)間0點(diǎn)主要分布在MLT=0—9午夜后至晨側(cè)區(qū)域,當(dāng)?shù)入x子體層向內(nèi)收縮時(shí),等離子體層外側(cè)此區(qū)域的合聲波可以在廣闊的L范圍內(nèi)加速由磁尾注入的粒子至更高能量.

    與亞暴發(fā)生期間增強(qiáng)的合聲波對電子的局地加速機(jī)制不同的是,ULF波動可以通過漂移-共振加速增強(qiáng)全球的相對論電子通量,因?yàn)閁LF波動中Pc5能段的波動頻率與相對論電子的漂移頻率接近(Green and Kivelson, 2004; 宗秋剛等, 2013; Su et al., 2015; Regi, 2016; Li et al., 2018).由磁尾注入的電子可以在午夜至晨側(cè)區(qū)域通過合聲波加速到較高能量,也可以通過與ULF波漂移共振繼續(xù)加速,并逐漸擴(kuò)散至所有MLT區(qū)域,因此,相對論電子通量長期倒空現(xiàn)象得以恢復(fù).本文所展示的在事件結(jié)束階段長達(dá)一天半時(shí)間內(nèi)高強(qiáng)度亞暴活動和ULF波活動與事件開始階段只有幾小時(shí)的活躍水平相比,可以為相對論電子通量的增長提供必要的源和加速機(jī)制,故而相對論電子通量從長期倒空狀態(tài)恢復(fù)時(shí)常伴隨著較強(qiáng)的亞暴活動和ULF波活動.

    4 維持相對論電子通量長期倒空現(xiàn)象的條件及其與磁暴強(qiáng)度的關(guān)系

    前文描述了相對論電子通量長期倒空事件開始和結(jié)束階段的不同特征.在開始階段,磁層頂向內(nèi)收縮,等離子體層頂一直維持在較高的高度上,向外的徑向擴(kuò)散結(jié)合磁層頂損失以及磁暴過程中波粒相互作用散射電子等損失機(jī)制占主導(dǎo)地位,相對論電子通量下降至背景通量水平.在結(jié)束階段,長時(shí)間的亞暴活動和ULF波活動等能為相對論電子的產(chǎn)生提供足夠的粒子源和加速過程,因此相對論電子通量可以從背景通量水平顯著增長.這一節(jié)則主要研究相對論電子可以持續(xù)3天以上時(shí)間維持在背景通量水平這一階段的特征.

    4.1 有無磁暴發(fā)生事件的選取標(biāo)準(zhǔn)

    由于磁暴發(fā)生期間的各種物理機(jī)制比較復(fù)雜,也為了研究有磁暴發(fā)生和無磁暴發(fā)生事件是否存在顯著的區(qū)別,本文把事件發(fā)生期間的過程分為兩類,即有磁暴發(fā)生事件和無磁暴發(fā)生事件.通常,Dst<-200 nT是大磁暴,-200≤Dst<-100 nT是強(qiáng)磁暴,-100 ≤Dst<-50 nT是中等磁暴,-50 ≤Dst<-30 nT是小磁暴(Gonzalez et al., 1999).我們以磁暴重要的兩個(gè)階段,即主相和恢復(fù)相(磁暴不同相位的劃分依照Partamies等(2013)的標(biāo)準(zhǔn))都出現(xiàn)在相對論電子通量長期倒空事件期間定義為有磁暴發(fā)生事件,如圖5a所示(圖5中所示參數(shù)與圖3中所示參數(shù)一致).圖中陰影區(qū)域表示相對論電子通量長期倒空現(xiàn)象.而如圖5b所示,盡管在9月12日事件開始前有中等磁暴發(fā)生,但當(dāng)相對論電子通量開始降低至背景通量水平時(shí)處于磁暴恢復(fù)相階段,其余陰影部分并沒有磁暴發(fā)生,我們把這類事件定義為無磁暴發(fā)生事件.關(guān)于圖5b所示的2014年9月相對論電子通量長期倒空事件,Alves等(2016)已進(jìn)行了相關(guān)研究,這個(gè)事件的產(chǎn)生是由在磁層深處L≤5區(qū)域電子通過共振波驅(qū)動的散射損失和在L>5區(qū)域絕熱和非絕熱作用導(dǎo)致的遠(yuǎn)離地心的徑向輸送過程將電子傳輸?shù)饺諅?cè)被壓縮的磁層區(qū)域,最終由磁層頂損失導(dǎo)致的.

    4.2 有無磁暴發(fā)生事件期間不同參數(shù)的統(tǒng)計(jì)對比

    我們把表1中的62例事件按4.1節(jié)的標(biāo)準(zhǔn)分類后,共得到18個(gè)有磁暴發(fā)生事件和44個(gè)無磁暴發(fā)生事件.為了更好的歸一化比較,并且清晰地陳列不同參數(shù)在不同數(shù)值區(qū)間的分布情況,我們以參數(shù)的具體數(shù)值區(qū)間為橫軸,以在不同區(qū)間范圍內(nèi)的事件個(gè)數(shù)占這類事件總數(shù)的百分比為縱軸,并對事件發(fā)生期間的太陽風(fēng)速度、密度、動壓、磁層頂高度、等離子體層頂高度、AE指數(shù)和ULF波指數(shù)的平均值和最值(除磁層頂高度為最小值外,其余均為最大值)分別統(tǒng)計(jì)比較,具體如圖6所示.圖中紅色柱狀圖表示有磁暴發(fā)生事件,藍(lán)色柱狀圖為無磁暴發(fā)生事件,每張小圖不同顏色的垂直線則表示這一類事件在不同參數(shù)統(tǒng)計(jì)下的平均水平(即單類事件在各個(gè)參數(shù)下的平均值).

    圖5 2002年2月有磁暴相對論電子通量長期倒空事件(a)和2014年9月無磁暴相對論電子通量長期倒空事件(b)Fig.5 Example of long-term dropout of relativistic electron flux with storm in February 2002 and without storm in September 2014

    圖6 有無磁暴發(fā)生事件期間各類參數(shù)的統(tǒng)計(jì)對比圖Fig.6 Comparison of various parameters during events with or without magnetic storm

    如圖6a1所示,有磁暴發(fā)生事件期間太陽風(fēng)速度平均值的平均水平(垂直紅線)為377 km·s-1,與無磁暴發(fā)生事件期間太陽風(fēng)速度平均值的平均水平(垂直藍(lán)線)376 km·s-1幾乎一致;這兩類事件期間太陽風(fēng)速最大值的平均水平相差也并不大,分別為536 km·s-1和522 km·s-1.這說明這兩類事件在太陽風(fēng)速度上并沒有較大差距,盡管大部分事件都有最高風(fēng)速在400~600 km·s-1范圍內(nèi),但整體上太陽風(fēng)速度比較低,事件發(fā)生期間太陽風(fēng)速度平均值幾乎都低于450 km·s-1以下.從圖6b1—b2、c1—c2可以看出,在太陽風(fēng)密度和動壓的平均值和最大值方面,有、無磁暴發(fā)生事件期間的平均水平相近,兩類事件在不同參數(shù)的不同區(qū)間分布情況也類似,只是有磁暴發(fā)生事件的平均水平比無磁暴發(fā)生事件的平均水平略高.圖6f1—f2、g1—g2關(guān)于AE指數(shù)和ULF波指數(shù)的統(tǒng)計(jì)結(jié)果也類似,即有磁暴發(fā)生事件的平均水平略高于無磁暴發(fā)生事件的平均水平.只有如圖6d1—d2、e1—e2所示的磁層頂高度與等離子體層頂高度的統(tǒng)計(jì)結(jié)果中,無磁暴發(fā)生事件期間的平均水平才比有磁暴發(fā)生事件期間的平均水平略高.這是因?yàn)闆]有磁暴發(fā)生時(shí),磁層相比地磁水平活躍時(shí)向外擴(kuò)展,因此磁層頂高度與等離子體層頂高度會較地磁水平活躍時(shí)偏高.

    綜上所述,相對論電子通量長期倒空現(xiàn)象發(fā)生時(shí),主要有兩種情況:第一種,通量倒空時(shí)有磁暴發(fā)生,太陽風(fēng)速度、密度、動壓都相比平靜時(shí)略高,與亞暴活動有關(guān)的AE指數(shù)偏高,ULF波指數(shù)也偏高.這類事件可能是由于磁暴發(fā)生期間,電子的損失機(jī)制(例如:動壓增強(qiáng)導(dǎo)致的磁層頂損失,磁暴期間波粒相互作用投擲角散射損失等)與亞暴注入、ULF波加速粒子機(jī)制等保持了一個(gè)動態(tài)平衡,所以相對論電子通量仍然可以維持在背景通量水平.第二種,與平靜的太陽風(fēng)條件和較弱的地磁活動水平、較少的ULF波活動水平相關(guān),這期間磁層頂高度與等離子體層頂高度都相比地磁水平活躍時(shí)偏高.這類事件可能是因?yàn)樘栵L(fēng)擾動較少,地球輻射帶一直保持在相對穩(wěn)定的狀態(tài),并且缺乏可以使相對論電子通量增長的源,故而這種情況下,相對論電子通量也可以一直保持在較低水平.盡管在不同種類的相對論電子通量長期倒空事件期間,不同參數(shù)不同數(shù)值區(qū)間的分布有些許差異,但這兩類事件在各個(gè)參數(shù)的平均水平十分相近,沒有較大差距,磁暴發(fā)生與否究竟與相對論電子通量長期倒空現(xiàn)象是否有關(guān),我們將在4.3節(jié)繼續(xù)探討.

    4.3 磁暴強(qiáng)度與維持相對論電子通量長期倒空現(xiàn)象的關(guān)系

    在引言部分,我們簡述了磁暴會對相對論電子通量產(chǎn)生三種影響,即增長、損失和近乎保持不變(Reeves et al., 2003; Anderson et al., 2015),并且在Kim等(2015)文中也提到了同步軌道電子通量的增強(qiáng)似乎與磁暴無關(guān),許多研究都表明磁暴并不是引起相對論電子通量變化的一個(gè)有效指標(biāo).但是,受Anderson等(2015)研究的小磁暴對相對論電子通量的損失能力要比中、高磁暴略強(qiáng)這一結(jié)果的啟發(fā),我們也對相對論電子通量能維持長期倒空狀態(tài)期間磁暴的強(qiáng)度進(jìn)行相關(guān)研究,來分析磁暴強(qiáng)度與維持這種長期倒空現(xiàn)象可能的聯(lián)系.

    為了對比事件開始、結(jié)束階段與事件發(fā)生期間出現(xiàn)磁暴的異同,我們統(tǒng)計(jì)了18例有磁暴發(fā)生事件和44例無磁暴發(fā)生事件,對事件開始、結(jié)束和相對論電子通量持續(xù)倒空期間三個(gè)不同時(shí)間階段出現(xiàn)的磁暴強(qiáng)度(S:小磁暴;M:中等磁暴;I:強(qiáng)磁暴和大磁暴)和相位(M:磁暴主相;R:磁暴恢復(fù)相)進(jìn)行統(tǒng)計(jì),如表2所示.從表中可以看出,無論是有磁暴發(fā)生事件還是無磁暴發(fā)生事件,在事件開始和結(jié)束階段,都有部分事例與磁暴過程相關(guān).并且,這兩類事件在事件開始時(shí),大部分事例處于磁暴主相階段;在事件結(jié)束時(shí),大部分事例處于磁暴恢復(fù)相階段.這是因?yàn)榇疟┲飨嗥陂g,與絕熱不變量有關(guān)的絕熱損失、沉降至地球大氣和磁層頂損失的非絕熱損失使相對論電子通量降低;在磁暴恢復(fù)相階段,通過亞暴活動和大規(guī)模太陽風(fēng)活動激發(fā)的感應(yīng)電場引起電子注入,伴隨著絕熱和非絕熱加速作用使相對論電子通量增長(Lazutin et al., 2018).

    表2 兩類相對論電子通量長期倒空事件開始、結(jié)束及事件期間磁暴強(qiáng)度、相位的統(tǒng)計(jì)Table 2 The intensity and phase of geomagnetic storms at the start time, end time and during the event occurred between two types of relativistic electron long-term dropouts

    縱觀這三個(gè)階段發(fā)生的磁暴強(qiáng)度可以看出,幾乎所有的磁暴都是小、中強(qiáng)度,并且小磁暴發(fā)生的次數(shù)也遠(yuǎn)比中等磁暴高.特別是在相對論電子通量維持長期倒空狀態(tài)時(shí),沒有一例大磁暴事件發(fā)生.雖然磁暴會引起相對論電子通量增加、減少或者保持不變,無法單從磁暴發(fā)生與否的角度判斷相對論電子通量的走向,但是我們的研究結(jié)果表明,相對論電子通量的變化在一定的程度上與磁暴的相位和磁暴的強(qiáng)度有關(guān).此外,相對論電子通量維持長期倒空時(shí),磁暴的強(qiáng)度并不高,這也是為什么我們4.2節(jié)統(tǒng)計(jì)的有磁暴發(fā)生事件和無磁暴發(fā)生事件在不同參數(shù)的平均水平相差不大的原因.

    磁暴效應(yīng)可以由兩種日地關(guān)系過程引發(fā),一種是日冕物質(zhì)拋射(CME)驅(qū)動的磁暴,另一種是共轉(zhuǎn)相互作用區(qū)(CIR)驅(qū)動的磁暴.大磁暴一般都是由CME驅(qū)動的,其磁暴主相的形成是由南向的IMF導(dǎo)致的,恢復(fù)相通常比較短暫;CIR驅(qū)動的磁暴主相由于IMFBz分量不斷波動所以呈現(xiàn)不規(guī)律性,其強(qiáng)度也比CME驅(qū)動的磁暴小很多,但恢復(fù)相可以從幾天持續(xù)到數(shù)周(Gonzalez et al., 1999; Borovsky and Denton, 2006).袁憧憬和宗秋剛(2011)通過研究磁暴發(fā)生期間外輻射帶1.5~6.0 MeV電子通量的變化發(fā)現(xiàn),對于CME驅(qū)動的磁暴無論是在主相還是恢復(fù)相期間電子通量在6≤L≤7區(qū)域都出現(xiàn)了增強(qiáng),而CIR驅(qū)動的磁暴,在相同的區(qū)域并沒有明顯地觀察到增強(qiáng)現(xiàn)象.Yuan和Zong(2012,2019)通過統(tǒng)計(jì)大量不同太陽風(fēng)/行星際結(jié)構(gòu)觸發(fā)的磁暴過程期間輻射帶電子總含量指數(shù)(Radiation Belt Content Index, RBC)的變化發(fā)現(xiàn),CME驅(qū)動的磁暴比CIR驅(qū)動的磁暴能夠更有效地加速整個(gè)外輻射帶電子,但在地球同步軌道上CIR驅(qū)動的磁暴比CME驅(qū)動的磁暴可以產(chǎn)生更多的相對論電子,這可能是因?yàn)镃ME驅(qū)動的磁暴在磁暴恢復(fù)相階段外輻射帶向內(nèi)移動,而CIR驅(qū)動的磁暴在磁暴恢復(fù)相階段外輻射帶向外移動.本文的研究結(jié)果表明,在維持同步軌道相對論電子通量長期倒空現(xiàn)象時(shí),沒有強(qiáng)磁暴發(fā)生,可能是因?yàn)榛謴?fù)相短暫的CME驅(qū)動的大磁暴在磁暴主相和恢復(fù)相間相對論電子通量都會出現(xiàn)增強(qiáng),不利于相對論電子通量長期維持在背景通量水平.雖然CIR磁暴能夠更有效地產(chǎn)生同步軌道相對論電子,但是并不是所有的磁暴都能加速電子,而我們的研究結(jié)果也僅代表地球同步軌道高度相對論電子通量的變化,不能代表整個(gè)外輻射帶相對論電子總含量的變化,也無從對比磁暴對整個(gè)外輻射帶產(chǎn)生的影響.當(dāng)整個(gè)外輻射帶電子通量倒空時(shí),同步軌道上相對論電子通量會出現(xiàn)倒空現(xiàn)象;當(dāng)輻射帶中心向內(nèi)移動,整個(gè)外輻射帶被壓縮時(shí),同步軌道相對論電子通量同樣會出現(xiàn)倒空現(xiàn)象.同步軌道相對論電子通量的變化僅是外輻射帶變化的一部分,整個(gè)外輻射帶的形態(tài)以及受磁暴的響應(yīng)仍需深入研究和探討.

    5 總結(jié)

    本文通過對2000—2016年GOES衛(wèi)星觀測的地球同步軌道上相對論電子通量長期倒空現(xiàn)象進(jìn)行統(tǒng)計(jì)研究,得出主要結(jié)論如下:

    (1)相對論電子通量長期倒空事件在太陽活動周不同階段的大尺度時(shí)間分布特征表明,在太陽活動周下降期,相對論電子通量長期倒空事件較少甚至沒有發(fā)生;在太陽活動周峰年和谷年,相對論電子通量倒空現(xiàn)象的發(fā)生率與太陽活動水平的高低并沒有直接聯(lián)系.

    (2)這類事件開始時(shí),主要發(fā)生在MLT=14—22區(qū)域.這類事件的產(chǎn)生條件,一方面與IMFBz分量在事件發(fā)生前北向轉(zhuǎn)南向,事件發(fā)生后由南向轉(zhuǎn)回北向并一直維持北向,Dst指數(shù)顯著下降發(fā)生磁暴過程,以及長時(shí)間擴(kuò)張的等離子體層有關(guān);另一方面與事件發(fā)生前后顯著增強(qiáng)太陽風(fēng)動壓、密度,磁層頂向內(nèi)壓縮有關(guān).

    (3)當(dāng)電子通量降低至背景通量水平時(shí),為了維持這種現(xiàn)象,無論長期倒空事件發(fā)生期間有無磁暴發(fā)生,都是在較弱的太陽風(fēng)條件、ULF波活動水平條件下.此外,長期倒空現(xiàn)象發(fā)生時(shí),不會出現(xiàn)強(qiáng)磁暴過程,并且小磁暴事件也比中等強(qiáng)度磁暴事件發(fā)生頻率高.

    (4)相對論電子通量打破倒空狀態(tài)時(shí),主要發(fā)生在MLT=0—9區(qū)域,通常會伴有長時(shí)間的高速太陽風(fēng)條件、長時(shí)間的高強(qiáng)度亞暴活動和ULF波活動以及向內(nèi)壓縮的等離子體層.

    本文研究的關(guān)注點(diǎn)在相對論電子通量可以長期保持倒空狀態(tài)這一類特殊現(xiàn)象,從統(tǒng)計(jì)方面加深了引起相對論電子通量變化的有關(guān)物理機(jī)制的認(rèn)識.本文的研究結(jié)果表明,等離子體層頂和磁層頂高度都有助于判斷相對論電子通量的變化,ULF波指數(shù)可以作為判斷相對論電子通量增長情況的一個(gè)有效指標(biāo),如果后續(xù)的預(yù)報(bào)模型可以增加更多的參數(shù),對加速機(jī)制和損失機(jī)制強(qiáng)弱的競爭關(guān)系充分建模,會對相對論電子通量的準(zhǔn)確預(yù)報(bào)起到一定的促進(jìn)作用.此外,本文的研究結(jié)果也充分說明了磁暴的強(qiáng)度可能會對相對論電子通量的變化產(chǎn)生影響,但有關(guān)磁暴過程中各種物理機(jī)制的綜合分析仍需深入的研究和探討.

    致謝特別感謝俄羅斯基礎(chǔ)研究基金會(19-05-00960)、國際空間科學(xué)研究所(ISSI-Beijing No.439)和國家留學(xué)基金的資助.感謝來自NOAA GOES提供的同步軌道相對論電子通量數(shù)據(jù),OMNI database提供的太陽風(fēng)參數(shù)數(shù)據(jù), WDC提供的地磁指數(shù)和ulf.gcras.ru/網(wǎng)站提供的ULF波指數(shù)數(shù)據(jù).特別感謝加利福尼亞大學(xué)伯克利分校Ilan Roth教授的建議和討論.

    91麻豆精品激情在线观看国产| 亚洲av免费在线观看| 中文字幕熟女人妻在线| 操出白浆在线播放| av专区在线播放| 99视频精品全部免费 在线| 18禁黄网站禁片免费观看直播| or卡值多少钱| 黄色片一级片一级黄色片| 黑人欧美特级aaaaaa片| av在线蜜桃| 久久久久久久久大av| 69人妻影院| 日本黄色片子视频| 国产伦在线观看视频一区| 亚洲精品亚洲一区二区| 天天添夜夜摸| 国产欧美日韩一区二区精品| 欧美一级a爱片免费观看看| 三级男女做爰猛烈吃奶摸视频| 亚洲精品成人久久久久久| 乱人视频在线观看| 少妇的丰满在线观看| 91在线观看av| 亚洲欧美日韩无卡精品| 69av精品久久久久久| 久久精品国产99精品国产亚洲性色| 少妇的逼水好多| 久久久久久人人人人人| 一个人免费在线观看的高清视频| 国产亚洲欧美98| 看黄色毛片网站| 国内少妇人妻偷人精品xxx网站| 亚洲国产欧洲综合997久久,| 国产真实乱freesex| 国产精品 欧美亚洲| 欧美日韩瑟瑟在线播放| 精品国产三级普通话版| 午夜老司机福利剧场| 亚洲av电影不卡..在线观看| 国产精品99久久久久久久久| 久久人人精品亚洲av| 亚洲狠狠婷婷综合久久图片| 首页视频小说图片口味搜索| 手机成人av网站| 亚洲精品日韩av片在线观看 | 在线播放国产精品三级| 动漫黄色视频在线观看| 国产精品一区二区三区四区久久| 俺也久久电影网| 天堂影院成人在线观看| 午夜福利在线观看免费完整高清在 | 在线天堂最新版资源| 久久精品国产亚洲av香蕉五月| 国产精品久久久久久人妻精品电影| 老司机午夜福利在线观看视频| 国产亚洲精品一区二区www| 国产av在哪里看| 嫩草影院精品99| 精品无人区乱码1区二区| 麻豆一二三区av精品| ponron亚洲| 亚洲内射少妇av| 亚洲欧美日韩高清专用| 一级毛片高清免费大全| 午夜福利18| 观看美女的网站| 首页视频小说图片口味搜索| 亚洲av二区三区四区| 久久精品国产亚洲av香蕉五月| 90打野战视频偷拍视频| av视频在线观看入口| 一个人免费在线观看电影| 听说在线观看完整版免费高清| 女人十人毛片免费观看3o分钟| 国产v大片淫在线免费观看| 国产亚洲精品一区二区www| 午夜福利18| 国产伦人伦偷精品视频| svipshipincom国产片| av天堂中文字幕网| 十八禁网站免费在线| av女优亚洲男人天堂| 国产精品电影一区二区三区| 亚洲美女黄片视频| 日本 欧美在线| 法律面前人人平等表现在哪些方面| 手机成人av网站| 色综合亚洲欧美另类图片| www.www免费av| 国产高清三级在线| 亚洲av成人不卡在线观看播放网| 亚洲欧美一区二区三区黑人| 精品久久久久久久久久久久久| 日本成人三级电影网站| 亚洲精品国产精品久久久不卡| 99精品久久久久人妻精品| 国产伦精品一区二区三区四那| 久久精品人妻少妇| 亚洲av第一区精品v没综合| 色av中文字幕| 久久精品91蜜桃| www.999成人在线观看| 亚洲国产色片| netflix在线观看网站| 村上凉子中文字幕在线| 免费人成视频x8x8入口观看| 制服人妻中文乱码| 久久久久亚洲av毛片大全| 欧美一区二区亚洲| 黄色丝袜av网址大全| 一区二区三区激情视频| 久久久久久大精品| 国产成人欧美在线观看| 亚洲成人免费电影在线观看| 欧美性猛交黑人性爽| 色视频www国产| 久久久久国产精品人妻aⅴ院| 一区福利在线观看| 国产精华一区二区三区| 在线播放无遮挡| 在线视频色国产色| 黄片大片在线免费观看| 日韩欧美在线乱码| 热99在线观看视频| 99精品久久久久人妻精品| 欧美高清成人免费视频www| 亚洲美女黄片视频| 色综合婷婷激情| 午夜影院日韩av| 日本a在线网址| svipshipincom国产片| 精品一区二区三区视频在线观看免费| 丰满乱子伦码专区| 亚洲不卡免费看| 日本黄大片高清| 综合色av麻豆| www国产在线视频色| 亚洲精品一卡2卡三卡4卡5卡| 宅男免费午夜| 97超级碰碰碰精品色视频在线观看| 一个人看视频在线观看www免费 | 91九色精品人成在线观看| 国产真人三级小视频在线观看| 免费看十八禁软件| 少妇的逼水好多| eeuss影院久久| 亚洲五月婷婷丁香| 国产亚洲精品综合一区在线观看| 久久精品国产亚洲av香蕉五月| 久久九九热精品免费| 成人精品一区二区免费| 搡女人真爽免费视频火全软件 | 岛国在线观看网站| 哪里可以看免费的av片| 成年免费大片在线观看| 亚洲av中文字字幕乱码综合| 成人国产一区最新在线观看| 91久久精品国产一区二区成人 | 美女免费视频网站| 给我免费播放毛片高清在线观看| 99久久99久久久精品蜜桃| 一级a爱片免费观看的视频| 国产精品 国内视频| 亚洲在线观看片| 国产三级在线视频| 中国美女看黄片| 美女被艹到高潮喷水动态| 欧美av亚洲av综合av国产av| 国产精华一区二区三区| 伊人久久精品亚洲午夜| 88av欧美| 欧美另类亚洲清纯唯美| 成人国产一区最新在线观看| 91久久精品国产一区二区成人 | 久久婷婷人人爽人人干人人爱| 最近在线观看免费完整版| 亚洲成av人片在线播放无| 国产真人三级小视频在线观看| 男人的好看免费观看在线视频| a级一级毛片免费在线观看| 一级a爱片免费观看的视频| 久久久久九九精品影院| 日本熟妇午夜| 国产黄a三级三级三级人| 欧美3d第一页| 91久久精品电影网| 在线观看午夜福利视频| 亚洲国产日韩欧美精品在线观看 | 99精品在免费线老司机午夜| 色在线成人网| av女优亚洲男人天堂| 国产精品乱码一区二三区的特点| www.999成人在线观看| 性欧美人与动物交配| 午夜两性在线视频| 69人妻影院| 人妻丰满熟妇av一区二区三区| 亚洲国产精品999在线| 99riav亚洲国产免费| 亚洲中文日韩欧美视频| 99热只有精品国产| 欧美黄色淫秽网站| 日本免费一区二区三区高清不卡| 在线观看午夜福利视频| 成人高潮视频无遮挡免费网站| 亚洲av美国av| 夜夜夜夜夜久久久久| 又紧又爽又黄一区二区| 日本a在线网址| 国产色婷婷99| 日本 av在线| 精品熟女少妇八av免费久了| 婷婷亚洲欧美| 叶爱在线成人免费视频播放| 丁香六月欧美| 免费电影在线观看免费观看| 91麻豆av在线| 国产爱豆传媒在线观看| 国产精品 国内视频| 中文字幕av在线有码专区| 搡女人真爽免费视频火全软件 | 一个人看的www免费观看视频| 亚洲熟妇熟女久久| 国内久久婷婷六月综合欲色啪| 中文字幕av成人在线电影| 欧美成人性av电影在线观看| 欧美高清成人免费视频www| 亚洲激情在线av| x7x7x7水蜜桃| 国产精品久久久久久久电影 | 变态另类成人亚洲欧美熟女| 69人妻影院| 成人欧美大片| 九九久久精品国产亚洲av麻豆| 一进一出抽搐动态| 久久精品91蜜桃| 亚洲av第一区精品v没综合| 99热精品在线国产| 18禁国产床啪视频网站| 国产精品久久电影中文字幕| 午夜影院日韩av| 久久久久九九精品影院| 午夜视频国产福利| 女人高潮潮喷娇喘18禁视频| 91久久精品电影网| 12—13女人毛片做爰片一| 噜噜噜噜噜久久久久久91| 99在线视频只有这里精品首页| 特大巨黑吊av在线直播| 免费av不卡在线播放| 日韩欧美三级三区| 久9热在线精品视频| 亚洲国产日韩欧美精品在线观看 | 少妇丰满av| 精品免费久久久久久久清纯| 国产单亲对白刺激| 香蕉久久夜色| 91av网一区二区| 一夜夜www| 校园春色视频在线观看| av福利片在线观看| 99热只有精品国产| av女优亚洲男人天堂| 午夜免费男女啪啪视频观看 | 成年女人永久免费观看视频| 18禁黄网站禁片午夜丰满| 99久久精品一区二区三区| 此物有八面人人有两片| 男人和女人高潮做爰伦理| 成人国产综合亚洲| 免费大片18禁| 久久久久久久久中文| 美女高潮喷水抽搐中文字幕| 国产在视频线在精品| 亚洲欧美日韩无卡精品| 18禁黄网站禁片午夜丰满| 亚洲人成网站在线播放欧美日韩| 色精品久久人妻99蜜桃| 日韩精品中文字幕看吧| 少妇裸体淫交视频免费看高清| av专区在线播放| 久久中文看片网| 亚洲avbb在线观看| 国产一区二区在线观看日韩 | 亚洲人成网站在线播放欧美日韩| 亚洲欧美日韩卡通动漫| 两性午夜刺激爽爽歪歪视频在线观看| 深爱激情五月婷婷| 中文资源天堂在线| 噜噜噜噜噜久久久久久91| 免费观看精品视频网站| 99精品欧美一区二区三区四区| 欧美另类亚洲清纯唯美| 国产不卡一卡二| 99久久精品国产亚洲精品| 成人av一区二区三区在线看| 欧美乱妇无乱码| 黄色成人免费大全| 日本a在线网址| 国产熟女xx| 久久久久久国产a免费观看| 精品久久久久久久久久免费视频| 日韩欧美三级三区| xxx96com| 国产精品免费一区二区三区在线| 日本在线视频免费播放| 国产一区二区在线观看日韩 | 午夜福利高清视频| 久9热在线精品视频| 全区人妻精品视频| 在线观看午夜福利视频| 国产精品亚洲一级av第二区| 中亚洲国语对白在线视频| 在线十欧美十亚洲十日本专区| 国产精品自产拍在线观看55亚洲| 91在线观看av| 俄罗斯特黄特色一大片| 国产伦在线观看视频一区| 色综合站精品国产| 午夜视频国产福利| 欧美性猛交黑人性爽| 99精品在免费线老司机午夜| 有码 亚洲区| 最近最新中文字幕大全电影3| 男女做爰动态图高潮gif福利片| 蜜桃亚洲精品一区二区三区| 亚洲精品成人久久久久久| 日本黄色片子视频| 在线视频色国产色| 老司机午夜福利在线观看视频| 最好的美女福利视频网| 97超级碰碰碰精品色视频在线观看| 91字幕亚洲| 国产精品嫩草影院av在线观看 | 国产精品一区二区免费欧美| 欧美黄色片欧美黄色片| 淫妇啪啪啪对白视频| 国产伦精品一区二区三区视频9 | 欧美在线黄色| 婷婷丁香在线五月| 午夜a级毛片| 成人亚洲精品av一区二区| 无遮挡黄片免费观看| 99热6这里只有精品| 国产99白浆流出| 国产精品爽爽va在线观看网站| 国产伦精品一区二区三区四那| 长腿黑丝高跟| 脱女人内裤的视频| 有码 亚洲区| 欧美最黄视频在线播放免费| 久久香蕉国产精品| 人妻丰满熟妇av一区二区三区| 深夜精品福利| 亚洲欧美日韩高清在线视频| 亚洲精品国产精品久久久不卡| 欧美最黄视频在线播放免费| 成人特级av手机在线观看| 国产精品野战在线观看| 国产精品爽爽va在线观看网站| 精品久久久久久久人妻蜜臀av| 国产国拍精品亚洲av在线观看 | 黄色女人牲交| av黄色大香蕉| 五月玫瑰六月丁香| 看免费av毛片| 国产精品日韩av在线免费观看| av专区在线播放| 国产一区二区在线av高清观看| 午夜精品一区二区三区免费看| 亚洲18禁久久av| 精品一区二区三区av网在线观看| 国产精品一区二区三区四区久久| 美女被艹到高潮喷水动态| 亚洲精品国产精品久久久不卡| 丁香六月欧美| 乱人视频在线观看| 亚洲欧美日韩东京热| av视频在线观看入口| 色老头精品视频在线观看| 亚洲aⅴ乱码一区二区在线播放| 国产精品一区二区三区四区免费观看 | 最近最新中文字幕大全电影3| 亚洲,欧美精品.| 观看免费一级毛片| 高清日韩中文字幕在线| 人妻丰满熟妇av一区二区三区| АⅤ资源中文在线天堂| 亚洲av中文字字幕乱码综合| 18禁在线播放成人免费| 成人高潮视频无遮挡免费网站| 亚洲aⅴ乱码一区二区在线播放| 欧美不卡视频在线免费观看| 中亚洲国语对白在线视频| 精品不卡国产一区二区三区| 精品欧美国产一区二区三| 一进一出抽搐gif免费好疼| 亚洲成人久久爱视频| 男人舔奶头视频| 99热6这里只有精品| 最好的美女福利视频网| 欧美极品一区二区三区四区| 高清日韩中文字幕在线| 啦啦啦韩国在线观看视频| 女生性感内裤真人,穿戴方法视频| 成人精品一区二区免费| 在线观看美女被高潮喷水网站 | 国产精品久久久久久久电影 | 又黄又爽又免费观看的视频| 母亲3免费完整高清在线观看| 日本成人三级电影网站| 久9热在线精品视频| 亚洲 国产 在线| 99久久精品国产亚洲精品| 亚洲欧美日韩无卡精品| 欧美国产日韩亚洲一区| 88av欧美| 在线国产一区二区在线| 午夜亚洲福利在线播放| av天堂中文字幕网| 国产欧美日韩一区二区精品| 老熟妇仑乱视频hdxx| 高清毛片免费观看视频网站| www日本在线高清视频| 国产国拍精品亚洲av在线观看 | 国产高清videossex| 精品一区二区三区人妻视频| 久久6这里有精品| 国产免费男女视频| 免费高清视频大片| 亚洲欧美精品综合久久99| 91久久精品电影网| 免费观看的影片在线观看| 欧美三级亚洲精品| 中文字幕精品亚洲无线码一区| 最好的美女福利视频网| 法律面前人人平等表现在哪些方面| 91麻豆精品激情在线观看国产| 国内精品久久久久精免费| 亚洲av日韩精品久久久久久密| 国产伦一二天堂av在线观看| 男女视频在线观看网站免费| 亚洲电影在线观看av| 可以在线观看的亚洲视频| 国产av在哪里看| 成熟少妇高潮喷水视频| 亚洲av免费在线观看| 精品久久久久久久人妻蜜臀av| 免费观看人在逋| 国产黄片美女视频| 757午夜福利合集在线观看| 99国产精品一区二区蜜桃av| 熟女电影av网| 亚洲成人久久爱视频| 亚洲成人精品中文字幕电影| 法律面前人人平等表现在哪些方面| 美女黄网站色视频| 国产伦精品一区二区三区四那| 国产精品女同一区二区软件 | 亚洲无线观看免费| 99久久无色码亚洲精品果冻| 老汉色∧v一级毛片| 精品无人区乱码1区二区| 精华霜和精华液先用哪个| 久久伊人香网站| 免费高清视频大片| 国产成+人综合+亚洲专区| 欧美在线一区亚洲| 亚洲 欧美 日韩 在线 免费| 一级黄片播放器| 搡老岳熟女国产| 欧美大码av| 久久精品综合一区二区三区| 久久精品国产综合久久久| 狠狠狠狠99中文字幕| 在线天堂最新版资源| 亚洲人与动物交配视频| 少妇丰满av| 淫妇啪啪啪对白视频| 舔av片在线| 国产毛片a区久久久久| 十八禁人妻一区二区| 国内少妇人妻偷人精品xxx网站| 成人鲁丝片一二三区免费| eeuss影院久久| 国产精品一区二区免费欧美| 18禁在线播放成人免费| 女同久久另类99精品国产91| 国产成人aa在线观看| www日本黄色视频网| 国产综合懂色| 乱人视频在线观看| 一级黄色大片毛片| 国产成人av教育| 啦啦啦免费观看视频1| 国产色爽女视频免费观看| 久久精品综合一区二区三区| 欧美不卡视频在线免费观看| 最近在线观看免费完整版| 久久国产乱子伦精品免费另类| 国产久久久一区二区三区| 床上黄色一级片| 国产伦一二天堂av在线观看| 桃红色精品国产亚洲av| 亚洲精品一区av在线观看| 日本在线视频免费播放| 成人特级av手机在线观看| 老熟妇乱子伦视频在线观看| 成人av在线播放网站| 亚洲精品日韩av片在线观看 | 草草在线视频免费看| 啦啦啦韩国在线观看视频| 51午夜福利影视在线观看| 中文亚洲av片在线观看爽| 狂野欧美白嫩少妇大欣赏| 啦啦啦免费观看视频1| 国产成人av激情在线播放| 香蕉丝袜av| 在线天堂最新版资源| 欧美色欧美亚洲另类二区| 国产精品久久久人人做人人爽| 波多野结衣巨乳人妻| 欧美zozozo另类| 欧美日韩黄片免| 女人被狂操c到高潮| 久久中文看片网| 欧美乱码精品一区二区三区| 嫁个100分男人电影在线观看| 此物有八面人人有两片| 丰满人妻熟妇乱又伦精品不卡| 成人特级av手机在线观看| 天天一区二区日本电影三级| 母亲3免费完整高清在线观看| 99久久99久久久精品蜜桃| 又黄又爽又免费观看的视频| 天美传媒精品一区二区| 国产探花极品一区二区| 亚洲性夜色夜夜综合| 内射极品少妇av片p| 99久久无色码亚洲精品果冻| 又紧又爽又黄一区二区| 日本一二三区视频观看| 中文亚洲av片在线观看爽| 国产精品嫩草影院av在线观看 | 婷婷六月久久综合丁香| 亚洲中文日韩欧美视频| 亚洲精品久久国产高清桃花| 精品人妻1区二区| 久久国产精品人妻蜜桃| 一夜夜www| 12—13女人毛片做爰片一| 校园春色视频在线观看| 丰满乱子伦码专区| 精品不卡国产一区二区三区| 两人在一起打扑克的视频| 欧美黄色片欧美黄色片| 亚洲欧美日韩东京热| 成人av在线播放网站| 国产主播在线观看一区二区| 18+在线观看网站| 国模一区二区三区四区视频| 亚洲乱码一区二区免费版| 久久久久久久亚洲中文字幕 | 亚洲久久久久久中文字幕| 色综合站精品国产| 国产精品亚洲一级av第二区| 国产成人欧美在线观看| 国产成人福利小说| 精品福利观看| 亚洲欧美激情综合另类| 久久久久久久久久黄片| www日本在线高清视频| 熟女人妻精品中文字幕| 免费av毛片视频| 91九色精品人成在线观看| 欧美成狂野欧美在线观看| 美女免费视频网站| 国产激情偷乱视频一区二区| 中国美女看黄片| 亚洲自拍偷在线| 夜夜躁狠狠躁天天躁| 最近最新中文字幕大全免费视频| 91字幕亚洲| 国产 一区 欧美 日韩| а√天堂www在线а√下载| 老熟妇仑乱视频hdxx| 99久久成人亚洲精品观看| 俄罗斯特黄特色一大片| 我的老师免费观看完整版| 欧洲精品卡2卡3卡4卡5卡区| or卡值多少钱| 欧美日韩黄片免| 国产成人啪精品午夜网站| 国产精品精品国产色婷婷| 麻豆久久精品国产亚洲av| 校园春色视频在线观看| 嫩草影院精品99| 国内精品美女久久久久久| 亚洲不卡免费看| 国内毛片毛片毛片毛片毛片| 免费在线观看影片大全网站| 内地一区二区视频在线| 久久久精品欧美日韩精品| 亚洲 欧美 日韩 在线 免费| 一卡2卡三卡四卡精品乱码亚洲| 九九热线精品视视频播放| 日韩中文字幕欧美一区二区| 国产v大片淫在线免费观看| 色综合站精品国产| 亚洲电影在线观看av| 精品久久久久久成人av| 俄罗斯特黄特色一大片| 91久久精品电影网| 欧美大码av|