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

    高能電子爆發(fā)與繞月衛(wèi)星表面電位大幅下降的聯(lián)動(dòng)效應(yīng)

    2016-11-08 03:04:55王馨悅張愛兵荊濤Reme孔令高張珅毅李春來(lái)
    地球物理學(xué)報(bào) 2016年10期
    關(guān)鍵詞:太陽(yáng)風(fēng)高能等離子體

    王馨悅, 張愛兵, 荊濤, H.Reme, 孔令高, 張珅毅, 李春來(lái)

    1 中國(guó)科學(xué)院國(guó)家空間科學(xué)中心, 北京 100190 2 天基空間環(huán)境探測(cè)北京市重點(diǎn)實(shí)驗(yàn)室, 北京 100190 3 法國(guó)國(guó)家科研中心, 法國(guó) 圖盧茲 31028 4 中國(guó)科學(xué)院國(guó)家天文臺(tái), 北京 100012

    ?

    高能電子爆發(fā)與繞月衛(wèi)星表面電位大幅下降的聯(lián)動(dòng)效應(yīng)

    王馨悅1,2, 張愛兵1,2, 荊濤1,2, H.Reme3, 孔令高1,2, 張珅毅1,2, 李春來(lái)4

    1 中國(guó)科學(xué)院國(guó)家空間科學(xué)中心, 北京100190 2 天基空間環(huán)境探測(cè)北京市重點(diǎn)實(shí)驗(yàn)室, 北京100190 3 法國(guó)國(guó)家科研中心, 法國(guó) 圖盧茲31028 4 中國(guó)科學(xué)院國(guó)家天文臺(tái), 北京100012

    “嫦娥”一號(hào)、二號(hào)繞月飛行經(jīng)歷地球磁尾邊界層區(qū)域時(shí),分別在2007年11月26日—2008年2月5日和2010年10月3日—2011年2月28日,發(fā)現(xiàn)了15次月球軌道0.1~2 MeV電子急劇增加(Bursts of 0.1~2 MeV Energetic Electrons, BEE),衛(wèi)星周圍等離子體離子加速的現(xiàn)象.統(tǒng)計(jì)研究表明,這類現(xiàn)象發(fā)生在穩(wěn)定太陽(yáng)風(fēng)和弱行星際磁場(chǎng)條件下,且無(wú)顯著空間環(huán)境擾動(dòng)事件發(fā)生時(shí),離子的加速滯后于高能電子爆發(fā),離子能量的變化與高能電子通量的時(shí)間演化正相關(guān),地球磁鞘內(nèi)側(cè)或邊界層過(guò)渡區(qū)域是該類現(xiàn)象的高發(fā)區(qū),離子能量增加時(shí)衛(wèi)星表面電位大幅下降可達(dá)負(fù)幾千伏.為了研究高能電子爆發(fā)與繞月衛(wèi)星表面電位變化的關(guān)系及其對(duì)月球表面電位的影響,本文用電流平衡法建立繞月衛(wèi)星和月球表面充電模型,并假設(shè)能量電子(>2 eV~2 MeV)滿足冪律譜的分布,模擬急劇增加的能量電子對(duì)衛(wèi)星和月球表面電位的影響.模擬結(jié)果表明,能量電子急劇增加使得繞月衛(wèi)星和月球表面電位大幅下降;能量電子總流量>1011cm-2時(shí),繞月衛(wèi)星和月球表面充電電位可達(dá)負(fù)上千伏;月球充電到大的負(fù)電位的時(shí)間僅為衛(wèi)星充電時(shí)間的1/10.鑒于高能電子急劇增加事件的高發(fā)生率(~125次/年),能量電子急劇增加使得繞月衛(wèi)星表面電位大幅下降的發(fā)生率應(yīng)大于實(shí)測(cè)等離子體離子加速現(xiàn)象的發(fā)生率(~25次/年).

    嫦娥一號(hào); 嫦娥二號(hào);月球;表面充電;高能電子急劇增加事件;等離子體

    1 引言

    月球沒(méi)有內(nèi)稟磁場(chǎng)和稠密的大氣層,每月約4天時(shí)間,月球穿越地球磁尾,其余時(shí)間則位于太陽(yáng)風(fēng)中.帶電粒子和太陽(yáng)光輻射直接轟擊月球,使得月球表面帶電.等離子體、太陽(yáng)光、二次電子和高能帶電粒子的撞擊使得繞月衛(wèi)星表面帶電(Stubbs et al., 2007).Lunar Prospector Electron Reflectometer (LP ER)的觀測(cè)表明(Halekas et al., 2005, 2008),通常情況下,月球表面向陽(yáng)側(cè)電位約為+10~+20 V,背陽(yáng)側(cè)電位為-100 V左右;在地球等離子體片區(qū)域,月球背陽(yáng)側(cè)電位可達(dá)負(fù)上千伏.Halekas等(2008,2009)利用LP ER實(shí)測(cè)的電子能譜模擬了太陽(yáng)質(zhì)子事件(Solar Energetic Proton Events, SEP)發(fā)生后繞月衛(wèi)星的表面充電電位,結(jié)果表明,SEP時(shí)高能粒子累積在衛(wèi)星表面可使繞月衛(wèi)星表面充電電位大幅下降.

    “嫦娥一號(hào)”衛(wèi)星發(fā)射入軌后,在月球軌道200 km高度上,首次發(fā)現(xiàn)了0.1~2 MeV電子急劇增加(BEE),隨后衛(wèi)星周圍等離子體離子加速的現(xiàn)象,這類現(xiàn)象通常發(fā)生在穩(wěn)定太陽(yáng)風(fēng)、弱行星際磁場(chǎng)(Interplanetary Magnetic Field, IMF)、空間環(huán)境弱擾動(dòng)條件下,伴隨BEE發(fā)生的等離子體離子加速現(xiàn)象可能是由于衛(wèi)星表面電位下降到大的負(fù)電位引起的,BEE可能是除了等離子體片和SEP外第三種使得繞月衛(wèi)星表面充電電位達(dá)負(fù)上千伏的誘因(Wang et al.,2012;王馨悅等,2012).“嫦娥二號(hào)”衛(wèi)星入軌后,在月球軌道100 km高度上也發(fā)現(xiàn)了同類現(xiàn)象.這類現(xiàn)象發(fā)生時(shí),0.1~2 MeV電子峰值通量比背景電子至少增長(zhǎng)了一個(gè)數(shù)量級(jí),且持續(xù)時(shí)間超過(guò)1 min.已有觀測(cè)和研究表明,能量電子急劇增長(zhǎng)事件可能出現(xiàn)在太陽(yáng)風(fēng)、地球磁尾、磁鞘及邊界層過(guò)渡區(qū)等空間環(huán)境區(qū)域.太陽(yáng)風(fēng)中的高能電子可能來(lái)源于弓激波上行粒子或“超級(jí)暈狀(super-halo)”電子(Lin, 1998).磁尾中的高能電子則可能由磁重聯(lián)引起(Egedal et al., 2010; Lu et al., 2010; Wang et al., 2010a, 2010b, 2014; Wu et al., 2015).在磁鞘內(nèi)側(cè)或邊界層過(guò)渡區(qū)域,持續(xù)幾分鐘并具有45°傾角特征的高能電子可能來(lái)自于弓激波上行粒子,持續(xù)時(shí)間較長(zhǎng)的高能電子可能起源于磁層或內(nèi)磁鞘(Anderson et al., 1979; Baker and Stone,1978; Bieber and Stone,1982; Formisano, 1979; Klassen et al., 2009).上述研究指出,急劇增加的能量電子能譜通常符合冪律譜的分布特征.

    為了研究“嫦娥”衛(wèi)星發(fā)現(xiàn)的高能電子急劇增加、隨后衛(wèi)星周圍等離子體離子加速、繞月衛(wèi)星表面電位下降到大的負(fù)電位現(xiàn)象的機(jī)制,本文統(tǒng)計(jì)研究了“嫦娥”一號(hào)、二號(hào)衛(wèi)星繞月飛行時(shí),在空間環(huán)境相對(duì)寧?kù)o的219天內(nèi),觀測(cè)到的75次0.1~2 MeV電子急劇增加現(xiàn)象與15次等離子體離子加速、衛(wèi)星表面電位大幅下降現(xiàn)象的特征與相關(guān)性;用電流平衡法建立繞月衛(wèi)星和月球表面充電模型,并將“嫦娥”實(shí)測(cè)高能電子數(shù)據(jù)擬合滿足冪律譜分布的能量電子能譜與實(shí)測(cè)等離子體能譜作為輸入條件,模擬BEE時(shí)繞月衛(wèi)星和月球表面充電電位的變化,研究高能電子急劇增加與繞月衛(wèi)星表面充電電位大幅下降之間的關(guān)系及其對(duì)月球表面充電電位的影響,為深入了解月球帶電粒子環(huán)境、進(jìn)一步開展月球科學(xué)探測(cè)提供參考.

    2 “嫦娥”衛(wèi)星空間環(huán)境探測(cè)儀

    “嫦娥一號(hào)”、“二號(hào)”是兩顆三軸穩(wěn)定的繞月極軌衛(wèi)星,飛行高度分別為200 km和100 km,軌道周期分別為127 min和118 min.空間環(huán)境探測(cè)儀由兩臺(tái)太陽(yáng)風(fēng)離子探測(cè)器(Solar Wind Ion Detectors,SWIDs)和一臺(tái)太陽(yáng)高能粒子探測(cè)器(High Energetic Particles Detector,HPD)組成,用于探測(cè)月球軌道帶電粒子環(huán)境(歐陽(yáng)自遠(yuǎn), 2010).“嫦娥一號(hào)”和“二號(hào)”空間環(huán)境探測(cè)儀具有相同的技術(shù)指標(biāo)和安裝位置.SWIDs用于觀測(cè)等離子體環(huán)境,主要觀測(cè)對(duì)象是太陽(yáng)風(fēng)離子.HPD用于監(jiān)測(cè)太陽(yáng)質(zhì)子事件,主要觀測(cè)對(duì)象是高能電子、離子.SWIDs和HPD由中國(guó)科學(xué)院國(guó)家空間科學(xué)中心(前空間科學(xué)與應(yīng)用研究中心)空間環(huán)境探測(cè)研究室研制.

    SWIDs是一對(duì)(SWID-A/B)采用相同設(shè)計(jì)的探測(cè)儀器.探測(cè)器采用半球形靜電分析器原理,并根據(jù)能量-電荷比(E/Q)的關(guān)系,使得具有一定能量的正離子入射到微通道板上,激發(fā)出電信號(hào).SWIDs各有12個(gè)具備相同探測(cè)能力的極角,稱為極角1~12,每個(gè)極角有48個(gè)能量通道,標(biāo)號(hào)為1~48(Kong et al.,2011; Wang et al., 2012).SWIDs的技術(shù)指標(biāo)如表1.衛(wèi)星發(fā)射前,SWIDs利用法國(guó)IRAP(former CESR of Toulouse in France, Institut de Recherche en Astrophysique et Planétologie)的5~800 keV離子源完成了地面定標(biāo)試驗(yàn).SWIDs安裝在衛(wèi)星頂角處指向天空,SWID-A和B的視場(chǎng)幾乎相互垂直.衛(wèi)星繞月飛行時(shí),SWID-A和B的探測(cè)視場(chǎng)分別垂直和平行于衛(wèi)星前進(jìn)方向.

    表1 SWIDs的技術(shù)指標(biāo)Table 1 Performance of SWIDs

    HPD的傳感器由3片半導(dǎo)體硅探測(cè)器組成,探測(cè)器視場(chǎng)角為60°,粒子進(jìn)入傳感器后,以電離方式損失能量,根據(jù)傳感器D1、D2和D3輸出能量的不同,區(qū)分粒子成分與能量(王馨悅等,2012).HPD安裝在衛(wèi)星的頂部,視場(chǎng)方向?yàn)槌煜?P1-P6為質(zhì)子(4~400 MeV),E1-E2為電子(0.1~2 MeV, ≥2.0 MeV).He(13~105 MeV)表示氦離子,C(117~590 MeV)能道包含的主要離子種類為C、N、O離子,Li(34~210 MeV)能道包含的主要離子種類為L(zhǎng)i、B、Be離子.加速器定標(biāo)試驗(yàn)分別在中國(guó)原子能科學(xué)研究院核物理所及中國(guó)科學(xué)院近代物理所完成.放射源試驗(yàn)和模擬信號(hào)源試驗(yàn)在中國(guó)科學(xué)院國(guó)家空間科學(xué)中心完成.

    3 探測(cè)結(jié)果

    “嫦娥一號(hào)”和“嫦娥二號(hào)”繞月飛行的2007年11月—2008年2月與2010年11月—2011年2月,在無(wú)顯著空間環(huán)境擾動(dòng)事件發(fā)生時(shí),月球軌道附近0.1~2 MeV的高能電子平均流量約40 (cm-2·sr-1·s-1·MeV-1),≥2 MeV高能質(zhì)子平均流量約0.2~3.0(cm-2·sr-1·s-1·MeV-1),遠(yuǎn)低于周圍等離子體的流量.

    3.10.1~2 MeV高能電子急劇增加(BEE)

    我們對(duì)峰值流量大于500 (cm-2·sr-1·s-1·MeV-1) 的0.1~2 MeV高能電子急劇增加事件(BEE)開展了統(tǒng)計(jì)研究.結(jié)果表明,2007年11月26日—2008年2月5日 ,“嫦娥一號(hào)”繞月飛行的70天中有15天發(fā)現(xiàn)了0.1~2 MeV高能電子急劇增加事件,電子峰值流量最高可達(dá) 2.4×104(cm-2·sr-1·s-1·MeV-1).其中有4天月球位于太陽(yáng)風(fēng),共發(fā)現(xiàn)5次持續(xù)時(shí)間僅數(shù)分鐘的高能電子通量增加.其他11天月球位于地球磁尾、磁鞘或邊界層過(guò)渡區(qū)域,觀測(cè)到的電子通量增加事件次數(shù)達(dá)36次,持續(xù)時(shí)間從幾分鐘到近10小時(shí).2010年10月3日—2011年2月28日,“嫦娥二號(hào)”衛(wèi)星繞月飛行的149 天中有17天發(fā)現(xiàn)了高能電子通量增加事件,電子峰值流量最高可達(dá)0.9×104(cm-2·sr-1·s-1·MeV-1).其中2天(1次/天)月球位于太陽(yáng)風(fēng),電子通量持續(xù)增加時(shí)間超過(guò)12小時(shí).其余15天月球位于地球磁尾、磁鞘內(nèi)側(cè)或磁層頂邊界層過(guò)渡區(qū)域,觀測(cè)到的電子通量增加事件次數(shù)約32次,持續(xù)時(shí)間從幾分鐘到近1天.綜上所述,BEE的年發(fā)生率為125次/年;月球每次穿越磁層時(shí),都會(huì)在磁鞘內(nèi)側(cè)或磁層頂邊界層過(guò)渡區(qū)域附近發(fā)現(xiàn)高能電子通量增加事件,約占全部BEE事件的70%,圖1是不同空間區(qū)域BEE事件的年發(fā)生率.我們統(tǒng)計(jì)了BEE期間衛(wèi)星在月心坐標(biāo)的位置,發(fā)現(xiàn)約70%的BEE事件發(fā)生在月球向陽(yáng)側(cè).SWIDs對(duì)等離子體離子的觀測(cè)顯示了當(dāng)?shù)氐入x子體區(qū)域的變化,當(dāng)衛(wèi)星從太陽(yáng)風(fēng)穿越到磁尾時(shí),可以觀測(cè)到離子的密度和速度下降、溫度上升的現(xiàn)象.圖3c和圖4c是兩次典型的BEE事件.

    圖1 太陽(yáng)風(fēng)、磁鞘內(nèi)側(cè)與磁層頂邊界層過(guò)渡區(qū)域和磁尾區(qū)域高能電子急劇增長(zhǎng)事件的年發(fā)生率Fig.1 Annual occurrence ratios of BEE in solar wind, inner terrestrial magnetosheath or the transition region of the boundaries and geomagnetotail

    3.2離子加速與衛(wèi)星表面充電電位大幅下降

    衛(wèi)星繞月飛行時(shí)發(fā)現(xiàn)的0.1~2 MeV電子通量急劇增加后,衛(wèi)星周圍等離子體離子加速、整體能量顯著抬升的現(xiàn)象通常發(fā)生在日側(cè)、月球穿越磁尾、經(jīng)歷磁鞘內(nèi)側(cè)區(qū)域附近;它們發(fā)生的時(shí)段則處于太陽(yáng)活動(dòng)水平較低的2007年11月到2011年2月之間,穩(wěn)定太陽(yáng)風(fēng)和弱行星際磁場(chǎng)(<10 nT)條件下,且空間環(huán)境較為寧?kù)o,無(wú)地磁暴發(fā)生.SWIDs與HPD的探測(cè)視場(chǎng)如圖2.BEE發(fā)生后,大量電子累積在衛(wèi)星表面,可使衛(wèi)星表面帶負(fù)電;根據(jù)圖2中探測(cè)器視場(chǎng)的幾何關(guān)系,進(jìn)入SWIDs視場(chǎng)的離子運(yùn)動(dòng)方向指向衛(wèi)星,當(dāng)衛(wèi)星表面電位下降到大的負(fù)電位時(shí),如圖7所示,向衛(wèi)星運(yùn)動(dòng)的離子將被星表電場(chǎng)加速,離子能量整體抬升,離子能量的增加正比于衛(wèi)星電位變化,由于通常光照區(qū)域衛(wèi)星表面帶電僅為正幾伏(Halekas et al., 2005, 2008, 2009), 當(dāng)大的充負(fù)電事件發(fā)生時(shí),衛(wèi)星表面充電電位可估算為U≈-|ΔE/q|(對(duì)H+,q=1e),其中ΔE為衛(wèi)星周圍等離子體離子能量的增量(eV),U是衛(wèi)星表面電位的增量(V).我們研究了引言提到的75次BEE事件,“嫦娥”一號(hào)、二號(hào)分別發(fā)現(xiàn)8次和7次衛(wèi)星周圍等離子體離子的顯著加速現(xiàn)象,占總BEE事件的約20%,該現(xiàn)象的年發(fā)生率約為25次/年,等離子體離子能量的增加通常比高能電子通量增加滯后幾分鐘,離子能量的變化趨勢(shì)與高能電子通量的急劇增加正相關(guān).

    圖2 高能電子急劇增加和等離子體離子加速聯(lián)動(dòng)現(xiàn)象發(fā)生、繞月衛(wèi)星位于月球向陽(yáng)面時(shí), SWIDs和HPD的探測(cè)視場(chǎng)示意圖Fig.2 Schematic diagram showing observed geometry of SWIDs and HPD during the BEE events and satellite charging in selenocentric solar ecliptic (SSE) coordinate system

    上述15次等離子體離子加速事件都發(fā)生在磁鞘內(nèi)側(cè)或磁層頂邊界層附近,離子速度方向指向SWIDs視場(chǎng)時(shí),且電子通量增加的持續(xù)時(shí)間超過(guò)了10 min.在未發(fā)現(xiàn)顯著等離子體離子加速的BEE事件中,約5%的BEE發(fā)生時(shí)SWIDs沒(méi)有數(shù)據(jù)下傳;由于太陽(yáng)風(fēng)、磁鞘附近離子分布具有很強(qiáng)的方向性,約20%的BEE發(fā)生時(shí),離子運(yùn)動(dòng)方向恰好位于SWIDs視場(chǎng)外;約18%的BEE中,衛(wèi)星經(jīng)歷了離子密度低于SWIDs計(jì)數(shù)下限(103(cm-2·sr-1·s-1·keV-1))的磁尾區(qū)域;約9%的BEE中,高能電子急劇增加在離子進(jìn)入SWIDs視場(chǎng)前發(fā)生, 無(wú)法判斷是否存在等離子體離子加速的過(guò)程;另外約28%的BEE中,高能電子通量增加的持續(xù)時(shí)間很短,僅幾分鐘電子通量即下降到背景值.

    圖3—4是“嫦娥一號(hào)”于2007年12月22日、“嫦娥二號(hào)”于2010年10月20日在月球向陽(yáng)側(cè)觀測(cè)到的兩次典型0.1~2.0 MeV 高能電子急劇增加與等離子體離子加速聯(lián)動(dòng)現(xiàn)象.F1,F(xiàn)2分別對(duì)應(yīng)了“嫦娥”一號(hào)、二號(hào)觀測(cè)的高能電子通量達(dá)到峰值的時(shí)間點(diǎn).F1時(shí)刻月球坐標(biāo)(-50RE, 25RE, 4.7RE)GSE,太陽(yáng)入射角(太陽(yáng)光和HPD探測(cè)視場(chǎng)中軸線夾角)約45°; F2時(shí)刻月球坐標(biāo)(-44RE, 39RE, -2.8RE)GSE,太陽(yáng)入射角約50°.SWIDs對(duì)衛(wèi)星周圍等離子體的觀測(cè)表明F1,F(xiàn)2時(shí)刻衛(wèi)星位于磁鞘內(nèi)側(cè)磁層頂邊界層附近,月球從太陽(yáng)風(fēng)向磁尾運(yùn)動(dòng)的過(guò)程中.F1前后,行星際磁場(chǎng)IMF |B|為2.5~5 nT,太陽(yáng)風(fēng)速度570~660 km·s-1,溫度7.4×104~1.9×105K.F2前后,行星際磁場(chǎng)IMF |B|為 5.5~6.2 nT,太陽(yáng)風(fēng)速度425~433 km·s-1,溫度 4.1×104~1.4×105K.SWIDs和ACE衛(wèi)星的觀測(cè)結(jié)果表明,兩次事件均具有相對(duì)穩(wěn)定的太陽(yáng)風(fēng)條件和弱的IMF(|B|<10 nT),無(wú)明顯空間環(huán)境擾動(dòng).

    圖3 “嫦娥一號(hào)”觀測(cè)的高能電子急劇增加和等離子體離子加速的聯(lián)動(dòng)現(xiàn)象(a) SWID-A觀測(cè)的離子能譜;(b) SWID-B觀測(cè)的離子能譜;(c) HPD觀測(cè)的0.1~2 MeV高能電子急劇增加現(xiàn)象,其中F1是電子通量峰值,T1對(duì)應(yīng)了加速前的離子能量,T2對(duì)應(yīng)了加速后的離子能量.Fig.3 Bursts of energetic electrons and the acceleration of the ambient plasma ions synchronously observed by Chang′E-1(a) Energy-time spectrogram of ions measured by SWID-A; (b) Energy-time spectrogram of ions measured by SWID-B; (c) 0.1~2 MeV electrons bursting observed by HPD. F1 is the peak flux of the bursting electrons. T1 is the ions energy before charging. T2 is the peak energy of accelerated ions after electron bursting of F1.

    圖4 “嫦娥二號(hào)”觀測(cè)的高能電子急劇增加和等離子體離子加速的聯(lián)動(dòng)現(xiàn)象(a) SWID-A觀測(cè)的離子能譜; (b) SWID-B觀測(cè)的離子能譜; (c) HPD觀測(cè)的0.1~2 MeV高能電子急劇增加現(xiàn)象,其中F2是電子通量峰值,T3對(duì)應(yīng)了加速前的離子能量,T4對(duì)應(yīng)了加速后的離子能量.Fig.4 Bursts of energetic electrons and the acceleration of the ambient plasma ions synchronously observed by Chang′E-2(a) Energy-time spectrogram of ions measured by SWID-A; (b) Energy-time spectrogram of ions measured by SWID-B; (c) 0.1~2 MeV electrons bursting observed by HPD. F2 is the peak flux of the bursting electrons. T3 is the ions energy before charging. T4 is the peak energy of accelerated ions after electron bursting of F1.

    圖3中,F(xiàn)1時(shí)刻 (UT 05∶09) 電子通量達(dá)約103(cm-2·sr-1·s-1·MeV-1),10 min后,離子能量上升到最大值,電子通量高于背景值的持續(xù)時(shí)間約27 min.圖4中,F(xiàn)2時(shí)刻(UT 14∶38) 電子通量達(dá)約103(cm-2·sr-1·s-1·MeV-1),7 min后,離子能量上升到最大值,電子通量高于背景值的持續(xù)時(shí)間約29 min.圖3—4中,T1,T3是等離子體離子加速前離子能量最大值的分布范圍,T2和T4是離子加速后離子能量最大值的分布范圍,白色圓點(diǎn)表示最大能量平均值,結(jié)果表明,T2時(shí)刻衛(wèi)星表面電位最大降幅為-5.4 kV,T4時(shí)刻衛(wèi)星表面電位最大降幅為-2 kV.圖5是等離子加速事件發(fā)生前和發(fā)生期間的離子能譜,結(jié)果表明,衛(wèi)星表面的電位增加提高了衛(wèi)星周圍等離子體離子的能量.

    月球穿越磁尾,將遭遇磁鞘、磁層頂邊界層、磁尾瓣區(qū)域,磁尾瓣區(qū)域離子流量遠(yuǎn)低于太陽(yáng)風(fēng)離子流量,使得觀測(cè)計(jì)數(shù)率低于SWIDs觀測(cè)下限,各區(qū)域的典型離子能譜分布的觀測(cè)特征如圖6.與太陽(yáng)風(fēng)相比,磁鞘離子溫度顯著上升,能譜分布觀測(cè)特征則與太陽(yáng)風(fēng)相似,各向異性特征顯著;磁層頂邊界層過(guò)渡區(qū)域附近離子的觀測(cè)特征與太陽(yáng)風(fēng)、磁鞘區(qū)域離子存在較大差異,它的密度顯著下降,能譜分布各向同性特征較強(qiáng).根據(jù)SWIDs實(shí)測(cè)等離子體分布特征隨時(shí)間的變化,圖3—4中兩次離子加速事件前后,月球均處于穿越磁鞘向磁尾方向的運(yùn)動(dòng)過(guò)程中;F1與F2對(duì)應(yīng)時(shí)刻太陽(yáng)入射角45°~50°,偏離太陽(yáng)風(fēng)入射方向;加速事件前后離子能譜分布具較強(qiáng)的各向同性特征,上述結(jié)果表明,離子加速期間衛(wèi)星始終處于磁層頂邊界層過(guò)渡區(qū)域附近.

    圖5 等離子加速事件發(fā)生前(2007-12-22 03∶06,2010-10-20 14∶02)和等離子加速事件發(fā)生期間(2007-12-22 05∶12,2010-10-20 14∶46)的離子能譜分布Fig.5 Ions distribution of preceding orbit (2007-12-22 03∶06,2010-10-20 14∶02) and in charging events (2007-12-22 05∶12,2010-10-20 14∶46) when the spacecraft passed through identical location. The ions spectrums of two samples show same characteristics.

    4 繞月衛(wèi)星表面充電建模與分析

    4.1建模方法

    在月球向陽(yáng)面,空間環(huán)境相對(duì)較為平靜、BEE發(fā)生時(shí),繞月衛(wèi)星表面收集的電流,主要包括環(huán)境等離子體電子和離子電流、光照作用下產(chǎn)生的光電子電流、BEE注入的高能電子電流、離子和電子撞擊表面產(chǎn)生的二次電子電流等等效電流,如圖7是衛(wèi)星表面電流成分示意圖,在等離子體環(huán)境中,表面電位將在衛(wèi)星周圍產(chǎn)生等離子體鞘,影響周圍等離子體分布,進(jìn)而影響環(huán)境產(chǎn)生的入射電流.衛(wèi)星表面的總電流代數(shù)和為零,即達(dá)到了電流平衡的狀態(tài),可以表示為ΣJk=0,其中,Jk表示衛(wèi)星表面收集的電流,包含光電流Jv、離子電流JI、電子電流JE、高能電子電流JEI和二次電子電流JSEC,每一項(xiàng)都是衛(wèi)星表面電位U的函數(shù).根據(jù)電流平衡方程建立繞月衛(wèi)星表面充電模型,模式的計(jì)算流程為假定衛(wèi)星表面電勢(shì)為U,通過(guò)改變U值,計(jì)算表面凈電流,凈電流約0時(shí)的U就是表面平衡電位.

    圖6 2010-10-16—2010-10-21月球穿越太陽(yáng)風(fēng)、遭遇磁鞘、磁層頂邊界層過(guò)渡區(qū)域,進(jìn)入磁尾過(guò)程中,“嫦娥二號(hào)”SWID-B觀測(cè)的離子能譜(a) 太陽(yáng)風(fēng); (b) 磁鞘; (c) 磁層頂邊界層過(guò)渡區(qū)域.Fig.6 Ion spectrum observed by SWID-B of Chang′E-2 with the Lunar crossing from the solar wind, the magnetosheath and the transition region of the boundaries to the magnetotail(a) Solar wind; (b) Magnetosheath; (c) Transition region of the boundaries.

    圖7 衛(wèi)星表面電流成分示意圖Fig.7 Sketch of surface charging of spacecraft

    磁鞘內(nèi)側(cè)、邊界層附近等離子體的德拜半徑約幾十米到幾百米,大于繞月衛(wèi)星的特征尺寸(1 m),與衛(wèi)星表面相互作用的等離子體可以近似為厚鞘層.對(duì)于厚鞘限制,表面收集的電子電流JE、離子電流JI與表面電壓關(guān)系如下(Whipple,1981;師立勤,2011):

    (1)

    JI=JI0(1-eU/E),

    (2)

    e是單位電荷電量,E是帶電粒子能量,JE0和JI0分別是能量為E的入射電子和入射離子的自由電流密度.月球特征尺寸遠(yuǎn)大于德拜半徑,建模時(shí)考慮薄鞘限制.對(duì)于薄鞘限制,表面收集的電子電流JE、離子電流JI與表面電壓關(guān)系如下(Whipple,1981;師立勤,2011):

    (3)

    JE=JE0, U≥0,

    (4)

    JI=JI0, U≤0,

    (5)

    JI=JI0(1-eU/E), U>0

    (6)

    將繞月衛(wèi)星看作是一個(gè)直徑1 m的球,表面覆蓋材料是鋁,則光照產(chǎn)生的光電子電流約50 μA·m-2.對(duì)于月球,光照產(chǎn)生的光電子電流約5 μA·m-2(Halekas et al., 2008, 2009).入射電子產(chǎn)生的二次電子電流為JSECE=YSEJE;入射質(zhì)子產(chǎn)生的二次電子電流表達(dá)為JSECI=YSIJI,YSI和YSE分別是入射質(zhì)子和電子的二次電子產(chǎn)生率(Whipple,1981;師立勤,2011).

    (7)

    YSE=2.228δm(QE-1+e-QE)(Em/E)0.35/QE.

    (8)

    (7)式中QI=1/E-0.1, Y1為能量1 keV入射質(zhì)子的產(chǎn)生率,Em是產(chǎn)生最大入射率的入射質(zhì)子能量;(8)式中,δm是入射電子產(chǎn)生的最大二次電子的產(chǎn)生率,QE=2.28(Em/E)1.35,Em是產(chǎn)生最大入射率的入射電子能量.衛(wèi)星或月球周圍等離子體的電流密度來(lái)自于SWIDs的實(shí)測(cè)結(jié)果.已有研究認(rèn)為太陽(yáng)風(fēng)或磁層中偶發(fā)的高能電子流通常符合冪律譜(J=AEγ,其中E表示能量,J表示電子通量,A和γ是方程的系數(shù))的分布(Wu et al., 2011),本文假設(shè)BEE時(shí)高能電子滿足冪律譜分布,用HPD的觀測(cè)結(jié)果,反演獲得BEE時(shí),>2 eV~2 MeV高能電子的時(shí)間積分譜.

    表面充電的總時(shí)間可以累加每一步的充電時(shí)間估算得到,即Δt=CAΔV/Jmean,其中,CA是單位面積上的電容,Jmean是每一步的平均電流,ΔV是每一步的電壓變化值.衛(wèi)星采用直徑1 m球形電容的假設(shè),月球看作一個(gè)無(wú)限大的平面,采用平板電容的假設(shè).

    4.2模擬與分析

    衛(wèi)星經(jīng)過(guò)月球向陽(yáng)面時(shí),SWIDs和HPD指向天空,它們的探測(cè)平面均位于衛(wèi)星向陽(yáng)側(cè),模擬表明,通常在圖2情況下,電流平衡時(shí),月球向陽(yáng)側(cè)光照區(qū)表面充電電位約+10 V,嫦娥衛(wèi)星的光照面表面充電電位約+14.5 V,與風(fēng)云衛(wèi)星和LP ER的觀測(cè)結(jié)果基本一致(Halekas et al., 2008, 2009; Wang et al., 2012).

    模擬高能電子通量急劇增加與繞月衛(wèi)星、月球表面電位變化之間的關(guān)系時(shí),為了簡(jiǎn)化計(jì)算過(guò)程,定量地給出高能電子急劇增加后衛(wèi)星和月球表面電位的變化,我們忽略了高能電子隨時(shí)間演化的復(fù)雜過(guò)程,僅反演了>2 eV~2 MeV高能電子的時(shí)間積分譜(從高能電子通量增加開始到等離子體能量增加達(dá)到峰值為止)作為輸入?yún)?shù),并假設(shè)充負(fù)電過(guò)程的發(fā)生從高能電子注入結(jié)束、衛(wèi)星或月球表面累積大量負(fù)電荷后開始.模擬選擇衛(wèi)星位于月球向陽(yáng)側(cè)和月球光照區(qū),初始表面電位均設(shè)置為+15 V.

    模擬表明,太陽(yáng)風(fēng)速度約400 km·s-1,密度約1 cm-3,溫度約1×105K條件下,0.1~2 MeV電子通量急劇增加期間,>2 eV~2 MeV高能電子總流量的時(shí)間積分>1011cm-2時(shí),繞月衛(wèi)星和月球表面電位可能達(dá)-1 kV.對(duì)F1, 電子通量急劇增加期間0.1~2 MeV電子平均流量約180(cm-2·sr-1·s-1·MeV-1),充電時(shí)間約13 min,2 MeV電子平均流量保持在0.2(cm-2·sr-1·s-1·MeV-1),高能電子總流量約2×1014cm-2,衛(wèi)星表面充電電位最高-5500 V,月球表面充電電位最高-5788 V,繞月衛(wèi)星表面電流平衡時(shí)間約0.002 s,月球表面電流平衡時(shí)間約0.0002 s.對(duì)F2, 電子通量急劇增加期間0.1~2 MeV電子平均流量約550(cm-2·sr-1·s-1·MeV-1),充電時(shí)間約14 min,2 MeV電子平均流量保持在約1.5(cm-2·sr-1·s-1·MeV-1),高能電子總流量約3.0×1011cm-2,衛(wèi)星表面充電電位最高約-1995 V,月球表面充電電位最高-1960 V,繞月衛(wèi)星表面電流平衡時(shí)間約0.009 s,月球表面電流平衡時(shí)間約0.0006 s.繞月衛(wèi)星表面充電電位的模擬計(jì)算結(jié)果與圖3—4的實(shí)測(cè)結(jié)果基本一致.圖8—9模擬了“嫦娥一號(hào)”2007年12月22日F1時(shí)刻0.1~2 MeV電子急劇增加通量達(dá)峰值后,高能電子總積分流量達(dá)2.0×1014cm-2時(shí),繞月衛(wèi)星和月球表面電流平衡狀態(tài)下的表面充電電位和充電時(shí)間,圖8—9中,表面收集電流約0 μA時(shí),衛(wèi)星和月球表面處于電流平衡狀態(tài),此時(shí)的表面電勢(shì)即為最大表面充電電位.上述結(jié)果表明,空間環(huán)境相對(duì)寧?kù)o時(shí)期,月球軌道高能電子急劇增加是繞月衛(wèi)星表面充電電位大幅下降,導(dǎo)致環(huán)境等離子體離子加速的主要誘因;月球表面在高能電子急劇增加時(shí)期也可能充電到大的負(fù)電位;同一條件下,月球表面充電電位達(dá)平衡時(shí)間僅為衛(wèi)星表面充電電位平衡時(shí)間的1/10左右,這是由于衛(wèi)星采用了厚鞘層限制以及球形電容、組成材料為鋁的假設(shè),月球采用了薄鞘層限制以及平板電容、組成材料為SiO2的假設(shè),衛(wèi)星的單位電容比月球大約10倍,使得月球表面充電時(shí)間遠(yuǎn)遠(yuǎn)小于衛(wèi)星充電時(shí)間.

    圖8 衛(wèi)星表面充電模擬Fig.8 Simulation of satellite surface charging

    圖9 月球表面充電模擬Fig.9 Simulation of Lunar surface charging

    5 結(jié)論

    本文利用“嫦娥”一號(hào)、二號(hào)衛(wèi)星對(duì)月球附近等離子體和高能電子的觀測(cè)結(jié)果,統(tǒng)計(jì)研究了0.1~2 MeV電子急劇增加時(shí),衛(wèi)星周圍等離子體離子加速,繞月衛(wèi)星表面充電電位大幅下降達(dá)負(fù)上千伏的現(xiàn)象.采用電流平衡法,建立月球向陽(yáng)側(cè)光照區(qū)表面、衛(wèi)星處于該區(qū)域時(shí)光照面的表面充電模型,結(jié)合“嫦娥”一號(hào)、二號(hào)對(duì)月球附近帶電粒子環(huán)境的實(shí)測(cè)結(jié)果,模擬了空間環(huán)境相對(duì)寧?kù)o時(shí),月球軌道高能電子急劇增加對(duì)繞月衛(wèi)星和月球表面電位的影響,結(jié)論如下:

    (1) 0.1~2 MeV電子急劇增加時(shí),衛(wèi)星周圍等離子體離子加速、衛(wèi)星表面電位從正幾伏快速下降到負(fù)上千伏的現(xiàn)象通常發(fā)生在穩(wěn)定太陽(yáng)風(fēng)和弱行星際磁場(chǎng)條件下,且無(wú)顯著空間環(huán)境擾動(dòng)發(fā)生時(shí),地球磁鞘內(nèi)側(cè)及邊界層附近區(qū)域是該類現(xiàn)象的高發(fā)區(qū).離子加速滯后于高能電子增加,離子能量的變化與高能電子通量的時(shí)間演化正相關(guān),加速前后離子能譜分布不變.

    (2) 月球軌道高能電子急劇增加是繞月衛(wèi)星表面電位大幅下降、加速周圍等離子體離子的主要誘因.

    (3) 鑒于高能電子急劇增加事件的高發(fā)生率(約125次/年)與等離子體觀測(cè)的局限性,能量電子急劇增加使得衛(wèi)星表面電位大幅下降的發(fā)生率應(yīng)大于實(shí)測(cè)的等離子體離子加速現(xiàn)象的發(fā)生率(25次/年).

    (4) 與繞月衛(wèi)星表面充電現(xiàn)象類似,月球附近高能電子的急劇增加可使月球向陽(yáng)側(cè)表面電位大幅下降.

    (5)高能電子總流量的時(shí)間積分>1011cm-2時(shí),繞月衛(wèi)星和月球表面充電電位可達(dá)-1 kV.

    (6) 同樣條件下,相比繞月衛(wèi)星,月球表面充電電位達(dá)平衡時(shí)間更短,是繞月衛(wèi)星表面充電電位達(dá)平衡時(shí)間的約1/10.

    高能電子急劇增加事件總是發(fā)生在穩(wěn)定太陽(yáng)風(fēng)和弱行星際磁場(chǎng)條件下,且與地磁暴、地磁亞暴沒(méi)有顯著相關(guān)性;月球每次經(jīng)歷磁鞘內(nèi)側(cè)、磁層頂邊界層附近區(qū)域時(shí),都能夠觀測(cè)到高能電子急劇增加事件的發(fā)生;磁層頂邊界層的厚度超過(guò)1000 km,月球的運(yùn)動(dòng)速度約1 km·s-1,則月球穿越磁層頂邊界層的時(shí)間可能超過(guò)100 min;此外,磁層頂邊界層的磁場(chǎng)強(qiáng)度比等離子體片大.上述事實(shí)表明,可能存在著月球和磁層頂邊界層相互作用的特殊加速區(qū)域,使得磁能轉(zhuǎn)化為粒子動(dòng)能,加速了來(lái)源于太陽(yáng)風(fēng)或磁層內(nèi)部的電子,引起了地球磁鞘內(nèi)側(cè)和邊界層區(qū)域高能電子的急劇增加.

    地球磁鞘內(nèi)側(cè)和邊界層區(qū)域是月球和繞月衛(wèi)星表面充電到負(fù)上千伏事件的高發(fā)區(qū),對(duì)探月活動(dòng)可能產(chǎn)生災(zāi)害性影響(Bedingfield et al., 1996; Koons et al., 1999).月球表面電位大幅下降,可造成當(dāng)?shù)貕m埃等離子體環(huán)境劇烈擾動(dòng),加劇月表環(huán)境中的靜電塵埃污染(Stubbs et al., 2006, 2007).當(dāng)高能電子急劇增加在月球背陽(yáng)面發(fā)生時(shí),由于沒(méi)有光電流,月表可能達(dá)到更高的負(fù)電位.

    致謝感謝SWIDs 和 HPD工作團(tuán)隊(duì)設(shè)計(jì)生產(chǎn)了空間環(huán)境探測(cè)儀.感謝探月工程地面應(yīng)用系統(tǒng)提供了“嫦娥”衛(wèi)星科學(xué)數(shù)據(jù).感謝ACE 工作團(tuán)隊(duì)和NASA的 CDAWEB提供的空間環(huán)境數(shù)據(jù).

    Anderson K A. 1979. A review of upstream and bow shock energetic-particle measurements.IlNuovoCimentoC, 2(6): 747-771.Baker D N, Stone E C. 1978. The magnetopause energetic electron layer, 1. Observations along the distant magnetotail.J.Geophys.Res., 83(A9): 4327-4338. Bedingfield K L, Leach R D, Alexander M B. 1996. Spacecraft system failures and anomalies attributed to the natural space environment.NASAReferencePublication, 1390: 1-5. Bieber J W, Stone E C. 1982. Energetic electrons in the magnetosheath and upstream of the bow shock.J.Geophys.Res., 87(A1): 85-94. Egedal J, Lê A, Zhu Y, et al. 2010. Cause of super-thermal electron heating during magnetotail reconnection.Geophys.Res.Lett., 37(10): 253-265.Formisano V. 1979. Observations of energetic electrons of magnetospheric origin in the magnetosheath and in the solar wind.IlNuovoCimentoC, 2(6): 781-788. Halekas J S, Lin R P, Mitchell D L. 2005. Large negative lunar surface potentials in sunlight and shadow.Geophys.Res.Lett., 32(9): 302-317.Halekas J S, Delory G T, Lin R P, et al. 2008. Lunar Prospector observations of the electrostatic potential of the lunar surface and its response to incident currents.J.Geophys.Res., 113(A9), doi: 10.1029/2008JA013194.

    Halekas J S, Delory G T, Lin R P, et al. 2009. Lunar surface charging during solar energetic particle events: Measurement and prediction.J.Geophys.Res., 114(A5), doi:10.1029/2009JA014113.Klassen A, Gómez-Herrero R, Müller-Mellin R, et al. 2009. STEREO/SEPT observations of upstream particle events: almost monoenergetic ion beams.Ann.Geophys., 27(5): 2077-2085.Kong L G, Wang S J, Wang X Y, et al. 2011. In-flight performance and preliminary observational results of Solar Wind Ion Detectors (SWIDs) on Chang′E-1.PlanetaryandSpaceScience, 2012, 62(1): 23-30. Koons H C, Mazur J E, Selesnick R S, et al. 1999. The impact of the space environment on space systems. Aerospace Report No. TR-99(1670)-1.

    Lin R P. 1998. WIND observations of suprathermal electrons in the interplanetary medium.SpaceScienceReviews, 86(1-4): 61-78.Lu Q M, Zhou L H, Wang S. 2010. Particle-in-cell simulations of whistler waves excited by an electronκdistribution in space plasma.J.Geophys.Res., 115(A2), doi: 10.1029/2009JA014580.Ouyang Z Y. 2010. Science results of Chang′E-1 lunar orbiter and mission goals of Chang′E-2.SpacecraftEngineering(in Chinese), 19(5): 1-6, doi: 10.3969/j.issn.1673-8748.2010.05.001.Shi L Q. 2011. Research on radiation environment and surface charging effect for LEO spacecraft [Ph. D. thesis] (in Chinese). Hefei: University of Science and Technology of China, doi: 10.7666/d.y2029992. Stubbs T J, Vondrak R R, Farrell W M. 2006. A dynamic fountain model for Lunar dust.AdvancesinSpaceResearch, 37(1): 59-66. Stubbs T J, Halekas J S, Farrell W M, et al. 2007. Lunar surface charging: A global perspective using Lunar Prospector data. ∥Krueger H, Graps A. Workshop on Dust in Planetary Systems (ESA SP-643). Kauai, Hawaii: ADS, 181-184.

    Wang R S, Lu Q M, Du A M, et al. 2010a. In situ observations of a Secondary Magnetic Island in an ion diffusion region and associated energetic electrons.Phys.Rev.Lett., 104(17): 39-87.Wang R S, Lu Q M, Li X, et al. 2010b. Observations of energetic electrons up to 200 keV associated with a secondary island near the center of an ion diffusion region: A Cluster case study.J.Geophys.Res., 115(A11), doi:10.1029/2010JA015473.

    Wang R S, Lu Q M, Khotyaintsev Y V, et al. 2014. Observation of double layer in the separatrix region during magnetic reconnection.Geophys.Res.Lett., 41(14): 4851-4858.Wang X Y, Jing T, Zhang S Y, et al. 2012. The first results of Chang′E-1 high energetic particles detector.ProgressinGeophysics(in Chinese), 27(6): 2289-2295, doi: 10.6038/j.issn.1004-2903.2012.06.003. Wang X Y, Zhang A B, Zhang X G, et al. 2012. Bursts of energetic electron induced large surface charging observed by Chang′E-1.PlanetaryandSpaceScience, 71(1): 1-8.

    Whipple E C. 1981. Potentials of surfaces in space.ReportsonProgressinPhysics, 44(11): 1197-1250.

    Wu M Y, Lu Q M, Du A M, et al. 2011. The evolution of the magnetic structures in electron phase-space holes: Two-dimensional particle-in-cell simulations.J.Geophys.Res., 116(A10): 1441-1458.

    Wu M Y, Huang C, Lu Q M, et al. 2015. In situ observations of multistage electron acceleration driven by magnetic reconnection.J.Geophys.Res., 120(8): 6320-6331.

    附中文參考文獻(xiàn)

    歐陽(yáng)自遠(yuǎn). 2010. 嫦娥一號(hào)衛(wèi)星的初步科學(xué)成果與嫦娥二號(hào)衛(wèi)星的使命. 航天器工程, 19(5): 1-6, doi: 10.3969/j.issn.1673-8748.2010.05.001.

    師立勤. 2011. 低軌道航空器輻射環(huán)境和表面充電效應(yīng)研究[博士論文]. 合肥: 中國(guó)科學(xué)技術(shù)大學(xué), doi: 10.7666/d.y2029992.

    王馨悅, 荊濤, 張珅毅等. 2012. “嫦娥一號(hào)”衛(wèi)星太陽(yáng)高能粒子探測(cè)器的首次觀測(cè)結(jié)果. 地球物理學(xué)進(jìn)展, 2012, 27(6): 2289-2295, doi: 10.6038/j.issn.1004-2903.2012.06.003.

    (本文編輯何燕)

    Synchronization of energetic electron bursting and lunar orbiter surface charging to negative kilovolts

    WANG Xin-Yue1,2, ZHANG Ai-Bing1,2, JING Tao1,2, H. Reme3, KONG Ling-Gao1,2, ZHANG Shen-Yi1,2, LI Chun-Lai4

    1NationalSpaceScienceCenter,ChineseAcademyofSciences,Beijing100190,China2BeijingKeyLaboratoryofSpaceEnvironmentExploration,Beijing100190,China3CentreNationaldelaRecherchéScientifique,Toulouse31028,France4NationalAstronomicalObservatories,ChineseAcademyofSciences,Beijing100012,China

    Fifteen times of ambient plasma ions accelerated on lunar orbit during the bursts of 0.1~2 MeV energetic electrons (BEE) was observed when spacecraft Chang′E-1 and Chang′E-2 flight through the terrestrial magnetosheath or the transition region of the boundaries during 26 November 2007 to 5 February 2008 and 3 October 2010 to 28 February 2011. The 0.1~2 MeV BEE events were observed by High Energetic Particles Detector (HPD) on board Chang′E-1 in a 200 km lunar polar orbit and Chang′E-2 in a 100 km lunar polar orbit. And the ions energy increasing was observed by Solar Wind Ion Detector (SWID) on board Chang′E-1 and 2. These events were found under relatively steady solar wind conditions and weakly IMF values without strong environmental disturbances. The ions acceleration occurred after the energetic electrons bursting, and the accelerated ions energy was correlated positively with the fluxes of the energetic electrons with the spacecraft surface charging to negative kilovolts. We use the current balance equations to simulate the spacecraft and the lunar surface charging during the bursting of the energetic electrons. The observed data of the two lunar orbiters is used to derive the energetic electrons spectrum (>2 eV~2 MeV) with power-low distribution assumption during BEEs. The simulation results indicate that the energetic electron incident of BEE is the dominating cause of the spacecraft surface charging to negative kilovolts. The results also indicate that the lunar and the spacecraft surface will charge to negative kilovolts during the BEEs after the temporal integral of the energetic electrons fluxes reach up to >1011cm-2. The balance time for the spacecraft surface charging is about 10 times that the lunar surface charging. It is expected that the occurrence of the synchronization of energetic electrons bursting and surface charging to large negative voltage should be more frequent than the Chang′E-1 and 2 observations though the observed ions acceleration (25 times per year) are fewer than the BEEs (125 times per year).

    Chang′E-1; Chang′E-2; Lunar; Surface charging; Bursts of energetic electrons; Plasma

    10.6038/cjg20161001.

    國(guó)家自然科學(xué)基金(41204128),中國(guó)科學(xué)院戰(zhàn)略性先導(dǎo)科技專項(xiàng)(XDA04077100)資助.

    王馨悅,女,1977年生,博士,主要從事空間環(huán)境探測(cè)數(shù)據(jù)分析與探測(cè)器物理設(shè)計(jì)工作.E-mail: Orchard@nssc.ac.cn

    10.6038/cjg20161001

    P353,P354

    2015-11-02,2016-09-05收修定稿

    王馨悅, 張愛兵, 荊濤等. 2016. 高能電子爆發(fā)與繞月衛(wèi)星表面電位大幅下降的聯(lián)動(dòng)效應(yīng). 地球物理學(xué)報(bào),59(10):3533-3542,

    Wang X Y, Zhang A B, Jing T,et al. 2016. Synchronization of energetic electron bursting and lunar orbiter surface charging to negative kilovolts.ChineseJ.Geophys. (in Chinese),59(10):3533-3542,doi:10.6038/cjg20161001.

    猜你喜歡
    太陽(yáng)風(fēng)高能等離子體
    前方高能!戰(zhàn)機(jī)懟瞼
    連續(xù)磁活動(dòng)對(duì)等離子體層演化的影響
    基于低溫等離子體修飾的PET/PVC浮選分離
    搞笑秀
    意林彩版(2020年3期)2020-11-22 20:54:12
    多種觀測(cè)數(shù)據(jù)驅(qū)動(dòng)的三維行星際太陽(yáng)風(fēng)MHD模擬
    基于ACE飛船觀測(cè)的銀河宇宙線與太陽(yáng)風(fēng)變化的統(tǒng)計(jì)研究
    《高能少年團(tuán)》少年 未來(lái)可期
    等離子體種子處理技術(shù)介紹
    在太空聽聲音
    戈?duì)柕膲?mèng)想飛船啟航:美國(guó)發(fā)射首顆位于L1點(diǎn)的太陽(yáng)風(fēng)監(jiān)測(cè)衛(wèi)星
    太空探索(2015年5期)2015-07-12 12:52:26
    天堂动漫精品| 欧美成人免费av一区二区三区| 成人欧美大片| 成人av一区二区三区在线看| 国产午夜精品久久久久久| 中国美女看黄片| www国产在线视频色| 免费在线观看日本一区| 老熟妇仑乱视频hdxx| 午夜a级毛片| 日本免费a在线| 精品一区二区三区四区五区乱码| 色综合婷婷激情| 熟女电影av网| 在线免费观看的www视频| 精品久久久久久久末码| 中文字幕最新亚洲高清| 亚洲成人免费电影在线观看| 黄片小视频在线播放| 久久中文字幕人妻熟女| 亚洲一区高清亚洲精品| 欧美绝顶高潮抽搐喷水| 亚洲国产看品久久| 色精品久久人妻99蜜桃| 日韩欧美国产在线观看| 精品国产亚洲在线| 长腿黑丝高跟| 亚洲成人久久爱视频| 日本一本二区三区精品| 亚洲人成网站在线播放欧美日韩| 国产成年人精品一区二区| 1024视频免费在线观看| 国产又黄又爽又无遮挡在线| 亚洲中文字幕一区二区三区有码在线看 | 一二三四在线观看免费中文在| 丁香欧美五月| 夜夜爽天天搞| 亚洲天堂国产精品一区在线| 动漫黄色视频在线观看| 久久久水蜜桃国产精品网| 久久午夜亚洲精品久久| 黄色毛片三级朝国网站| 女人高潮潮喷娇喘18禁视频| 日本免费a在线| www国产在线视频色| 搞女人的毛片| av欧美777| 亚洲av熟女| 精品久久蜜臀av无| 精品福利观看| 岛国视频午夜一区免费看| 日韩高清综合在线| 怎么达到女性高潮| 国产三级中文精品| 久久精品国产亚洲av高清一级| 可以免费在线观看a视频的电影网站| 中文资源天堂在线| 99re在线观看精品视频| 狠狠狠狠99中文字幕| 成年女人毛片免费观看观看9| 特级一级黄色大片| 亚洲精品久久成人aⅴ小说| 99久久无色码亚洲精品果冻| 成人18禁在线播放| 中文亚洲av片在线观看爽| 丰满的人妻完整版| 一本一本综合久久| 成人国产综合亚洲| 91av网站免费观看| 这个男人来自地球电影免费观看| 中国美女看黄片| 99国产极品粉嫩在线观看| 99热这里只有精品一区 | 亚洲欧洲精品一区二区精品久久久| 亚洲精品粉嫩美女一区| 妹子高潮喷水视频| 亚洲国产中文字幕在线视频| 国产精品久久久久久人妻精品电影| 久热爱精品视频在线9| 美女扒开内裤让男人捅视频| 久久精品91无色码中文字幕| 欧美日本亚洲视频在线播放| 制服人妻中文乱码| 夜夜爽天天搞| 五月伊人婷婷丁香| 中文字幕av在线有码专区| 国产欧美日韩一区二区精品| 亚洲激情在线av| 成人特级黄色片久久久久久久| 中文字幕最新亚洲高清| 人成视频在线观看免费观看| 亚洲国产精品sss在线观看| 国产亚洲精品av在线| 欧美日韩一级在线毛片| 黄色视频不卡| 亚洲全国av大片| 在线观看舔阴道视频| 久久久精品国产亚洲av高清涩受| 亚洲五月天丁香| 国产黄a三级三级三级人| 亚洲美女视频黄频| 少妇粗大呻吟视频| 亚洲最大成人中文| 波多野结衣巨乳人妻| 欧美一级毛片孕妇| 母亲3免费完整高清在线观看| 国产成年人精品一区二区| 又粗又爽又猛毛片免费看| 国产乱人伦免费视频| 亚洲人成网站高清观看| 精品一区二区三区四区五区乱码| 午夜免费激情av| 91av网站免费观看| www.999成人在线观看| ponron亚洲| 日韩精品中文字幕看吧| 久久久久久久久免费视频了| 久久久久国产一级毛片高清牌| 亚洲精华国产精华精| 听说在线观看完整版免费高清| 国产欧美日韩精品亚洲av| 99久久精品热视频| 国产91精品成人一区二区三区| 男女午夜视频在线观看| 亚洲成人免费电影在线观看| 午夜日韩欧美国产| 亚洲avbb在线观看| 国产熟女午夜一区二区三区| 久久欧美精品欧美久久欧美| 免费在线观看视频国产中文字幕亚洲| 日本一本二区三区精品| 狂野欧美激情性xxxx| 可以在线观看的亚洲视频| 国产精品久久久久久人妻精品电影| 亚洲精品粉嫩美女一区| 免费电影在线观看免费观看| 人人妻,人人澡人人爽秒播| 小说图片视频综合网站| 亚洲一区中文字幕在线| 国产伦一二天堂av在线观看| 日本一本二区三区精品| 成人18禁高潮啪啪吃奶动态图| 少妇的丰满在线观看| bbb黄色大片| 亚洲欧美日韩高清专用| 欧美黑人巨大hd| 舔av片在线| 国产免费av片在线观看野外av| 啦啦啦免费观看视频1| 国产亚洲精品久久久久久毛片| 欧美 亚洲 国产 日韩一| 淫妇啪啪啪对白视频| 国产精品久久电影中文字幕| 国产麻豆成人av免费视频| 1024视频免费在线观看| 免费看美女性在线毛片视频| 亚洲熟妇中文字幕五十中出| 国产精品久久久av美女十八| 日韩欧美 国产精品| 亚洲熟女毛片儿| 日韩精品中文字幕看吧| 激情在线观看视频在线高清| netflix在线观看网站| 最近最新中文字幕大全电影3| 亚洲avbb在线观看| 婷婷丁香在线五月| 可以免费在线观看a视频的电影网站| 久久久精品国产亚洲av高清涩受| www.www免费av| 最近最新中文字幕大全免费视频| 极品教师在线免费播放| 床上黄色一级片| 久久香蕉国产精品| 久久久久久国产a免费观看| 久久久水蜜桃国产精品网| 国产高清视频在线观看网站| 美女高潮喷水抽搐中文字幕| 免费在线观看完整版高清| 国产精品久久久人人做人人爽| 少妇被粗大的猛进出69影院| 在线观看一区二区三区| 日本 av在线| 国产成人系列免费观看| 亚洲成人中文字幕在线播放| 女人高潮潮喷娇喘18禁视频| 国产男靠女视频免费网站| 神马国产精品三级电影在线观看 | 欧美日韩乱码在线| 成在线人永久免费视频| 国产高清有码在线观看视频 | 成人精品一区二区免费| 欧美精品亚洲一区二区| 亚洲自拍偷在线| 久久精品综合一区二区三区| 亚洲黑人精品在线| 久久久久性生活片| 国产aⅴ精品一区二区三区波| 久久精品亚洲精品国产色婷小说| ponron亚洲| 在线永久观看黄色视频| 少妇裸体淫交视频免费看高清 | 日日摸夜夜添夜夜添小说| 久久国产精品影院| 亚洲片人在线观看| 香蕉丝袜av| 非洲黑人性xxxx精品又粗又长| 日韩欧美在线二视频| 少妇熟女aⅴ在线视频| 天天躁夜夜躁狠狠躁躁| 99精品欧美一区二区三区四区| cao死你这个sao货| 成年女人毛片免费观看观看9| 啦啦啦韩国在线观看视频| 亚洲熟女毛片儿| 国产av一区二区精品久久| 国产伦人伦偷精品视频| 国产成人aa在线观看| 国产精品一区二区免费欧美| 两性午夜刺激爽爽歪歪视频在线观看 | 黄色毛片三级朝国网站| 成人精品一区二区免费| 麻豆成人午夜福利视频| 身体一侧抽搐| 欧美黄色淫秽网站| 天天躁夜夜躁狠狠躁躁| 久久人人精品亚洲av| 精品福利观看| 精品人妻1区二区| 日韩 欧美 亚洲 中文字幕| av在线天堂中文字幕| 在线观看日韩欧美| 日日爽夜夜爽网站| 色哟哟哟哟哟哟| 免费电影在线观看免费观看| 免费搜索国产男女视频| 欧美极品一区二区三区四区| 成人国产综合亚洲| 亚洲国产精品sss在线观看| 少妇裸体淫交视频免费看高清 | 一级毛片女人18水好多| 国产亚洲av高清不卡| 亚洲美女视频黄频| 色播亚洲综合网| 色综合欧美亚洲国产小说| 久久久久久久久中文| 精品一区二区三区视频在线观看免费| 色综合亚洲欧美另类图片| 两个人视频免费观看高清| 欧美日韩一级在线毛片| 午夜精品一区二区三区免费看| 不卡一级毛片| 国产私拍福利视频在线观看| 亚洲aⅴ乱码一区二区在线播放 | 国产男靠女视频免费网站| 国产一区二区激情短视频| 香蕉国产在线看| 中文字幕久久专区| 黄色 视频免费看| 欧美黑人精品巨大| 母亲3免费完整高清在线观看| 国产成人啪精品午夜网站| 亚洲精品一卡2卡三卡4卡5卡| 美女大奶头视频| 极品教师在线免费播放| 精品一区二区三区四区五区乱码| 一级毛片高清免费大全| 一区二区三区激情视频| 国产av一区二区精品久久| 久久亚洲真实| 久久午夜综合久久蜜桃| 亚洲专区国产一区二区| 色噜噜av男人的天堂激情| 成人午夜高清在线视频| 成人亚洲精品av一区二区| 国产乱人伦免费视频| 国产蜜桃级精品一区二区三区| а√天堂www在线а√下载| 校园春色视频在线观看| 欧美激情久久久久久爽电影| 在线视频色国产色| 俄罗斯特黄特色一大片| 婷婷丁香在线五月| 国产精品一区二区精品视频观看| 少妇被粗大的猛进出69影院| 嫁个100分男人电影在线观看| 一本精品99久久精品77| 三级男女做爰猛烈吃奶摸视频| 一进一出好大好爽视频| 丝袜人妻中文字幕| 国产精品永久免费网站| 亚洲中文字幕日韩| 美女黄网站色视频| 黄片大片在线免费观看| 高潮久久久久久久久久久不卡| 黄色 视频免费看| 午夜福利免费观看在线| 国产视频内射| 欧美乱码精品一区二区三区| 97人妻精品一区二区三区麻豆| 人成视频在线观看免费观看| 国产97色在线日韩免费| 久久久久国产精品人妻aⅴ院| 日韩精品免费视频一区二区三区| 欧美在线一区亚洲| 精品久久蜜臀av无| 国语自产精品视频在线第100页| 国产蜜桃级精品一区二区三区| 久久久久九九精品影院| 高潮久久久久久久久久久不卡| 神马国产精品三级电影在线观看 | 看黄色毛片网站| 欧美+亚洲+日韩+国产| 久久久久久大精品| 女人被狂操c到高潮| 欧美精品啪啪一区二区三区| 国产精品av视频在线免费观看| 欧美成狂野欧美在线观看| 免费在线观看日本一区| 国产三级中文精品| 国产伦人伦偷精品视频| 国产精品 欧美亚洲| 亚洲 欧美 日韩 在线 免费| 国产精品久久久久久精品电影| 午夜福利在线观看吧| 日本三级黄在线观看| 夜夜看夜夜爽夜夜摸| 国产一区二区激情短视频| aaaaa片日本免费| 欧美绝顶高潮抽搐喷水| 在线观看一区二区三区| 最近最新中文字幕大全免费视频| 国产精品 欧美亚洲| 色哟哟哟哟哟哟| 极品教师在线免费播放| 熟女少妇亚洲综合色aaa.| 国产亚洲精品久久久久久毛片| www.999成人在线观看| 日韩有码中文字幕| 国内精品一区二区在线观看| 亚洲,欧美精品.| 很黄的视频免费| 亚洲av五月六月丁香网| 亚洲精品色激情综合| 国产成人av教育| 亚洲专区字幕在线| 国产主播在线观看一区二区| 操出白浆在线播放| 亚洲av片天天在线观看| 五月伊人婷婷丁香| 俺也久久电影网| 欧美乱码精品一区二区三区| 日韩欧美一区二区三区在线观看| 欧美日韩福利视频一区二区| 亚洲午夜精品一区,二区,三区| 亚洲av成人av| 亚洲av电影不卡..在线观看| 精品国产亚洲在线| 亚洲熟妇熟女久久| 99riav亚洲国产免费| 欧美大码av| 女人高潮潮喷娇喘18禁视频| 俺也久久电影网| 看黄色毛片网站| 久久精品国产清高在天天线| 日本撒尿小便嘘嘘汇集6| 国产精品免费视频内射| 久久久久久国产a免费观看| 九色国产91popny在线| 男人舔女人的私密视频| 又大又爽又粗| 黑人操中国人逼视频| 亚洲精品色激情综合| 日韩欧美在线二视频| 久久久久久九九精品二区国产 | 好看av亚洲va欧美ⅴa在| 免费电影在线观看免费观看| 亚洲一区高清亚洲精品| 久久久久久免费高清国产稀缺| 国产亚洲欧美在线一区二区| 亚洲精品国产一区二区精华液| 亚洲中文字幕一区二区三区有码在线看 | 亚洲精品一区av在线观看| 亚洲五月婷婷丁香| 婷婷亚洲欧美| 九色成人免费人妻av| 色尼玛亚洲综合影院| 国产亚洲精品第一综合不卡| 十八禁人妻一区二区| 日韩国内少妇激情av| 波多野结衣高清作品| 国产精品久久电影中文字幕| 精品人妻1区二区| 亚洲色图 男人天堂 中文字幕| 国产精品一区二区三区四区免费观看 | 亚洲五月婷婷丁香| 999久久久精品免费观看国产| 好男人在线观看高清免费视频| 久久久久久人人人人人| 午夜老司机福利片| 美女大奶头视频| 亚洲avbb在线观看| 黄片小视频在线播放| 黑人欧美特级aaaaaa片| 哪里可以看免费的av片| 精品久久久久久久人妻蜜臀av| 嫩草影院精品99| 欧美绝顶高潮抽搐喷水| 亚洲五月天丁香| 全区人妻精品视频| 99在线视频只有这里精品首页| 日本免费一区二区三区高清不卡| 最新在线观看一区二区三区| 美女午夜性视频免费| 日本三级黄在线观看| 国产精品免费一区二区三区在线| 丁香六月欧美| 制服丝袜大香蕉在线| 色老头精品视频在线观看| 色播亚洲综合网| 亚洲人成网站在线播放欧美日韩| 国产av又大| 又粗又爽又猛毛片免费看| 夜夜夜夜夜久久久久| 日韩欧美精品v在线| 亚洲avbb在线观看| 88av欧美| 亚洲av成人一区二区三| 欧美zozozo另类| 精品少妇一区二区三区视频日本电影| 给我免费播放毛片高清在线观看| 午夜福利高清视频| 欧美中文日本在线观看视频| 日本 av在线| xxx96com| av有码第一页| 国产成+人综合+亚洲专区| 亚洲黑人精品在线| 妹子高潮喷水视频| 99久久国产精品久久久| 啦啦啦免费观看视频1| 国产片内射在线| 亚洲精品美女久久久久99蜜臀| 精品高清国产在线一区| 岛国在线免费视频观看| 小说图片视频综合网站| 日本一区二区免费在线视频| 欧美色欧美亚洲另类二区| 日韩成人在线观看一区二区三区| 熟女少妇亚洲综合色aaa.| 亚洲成人国产一区在线观看| 精品久久久久久久人妻蜜臀av| 少妇裸体淫交视频免费看高清 | 在线观看美女被高潮喷水网站 | 少妇裸体淫交视频免费看高清 | 哪里可以看免费的av片| 欧美一级a爱片免费观看看 | 婷婷亚洲欧美| 99国产精品一区二区蜜桃av| 在线观看舔阴道视频| 叶爱在线成人免费视频播放| av福利片在线观看| 可以在线观看的亚洲视频| 无遮挡黄片免费观看| 校园春色视频在线观看| 十八禁网站免费在线| 久久精品国产综合久久久| cao死你这个sao货| 老汉色∧v一级毛片| 国产视频内射| 在线永久观看黄色视频| 天堂av国产一区二区熟女人妻 | 亚洲中文日韩欧美视频| 色播亚洲综合网| 亚洲一区二区三区色噜噜| 97超级碰碰碰精品色视频在线观看| 在线观看一区二区三区| 免费观看精品视频网站| 国产黄色小视频在线观看| 国产亚洲av高清不卡| bbb黄色大片| 国产一区在线观看成人免费| 国产在线精品亚洲第一网站| 亚洲av电影在线进入| 在线十欧美十亚洲十日本专区| 日韩高清综合在线| 91字幕亚洲| 成熟少妇高潮喷水视频| 亚洲av成人av| 日本熟妇午夜| 国产av又大| 午夜影院日韩av| 亚洲欧洲精品一区二区精品久久久| 麻豆国产97在线/欧美 | or卡值多少钱| 曰老女人黄片| 好男人在线观看高清免费视频| 午夜影院日韩av| 伦理电影免费视频| 国产精品爽爽va在线观看网站| 亚洲午夜精品一区,二区,三区| 制服诱惑二区| 日韩欧美国产一区二区入口| 热99re8久久精品国产| 亚洲av成人不卡在线观看播放网| 亚洲一卡2卡3卡4卡5卡精品中文| 两个人的视频大全免费| 在线观看免费午夜福利视频| 欧美+亚洲+日韩+国产| 国产精品美女特级片免费视频播放器 | 国产av麻豆久久久久久久| 欧美黑人欧美精品刺激| 人妻丰满熟妇av一区二区三区| 一二三四社区在线视频社区8| e午夜精品久久久久久久| 亚洲欧美激情综合另类| 人人妻,人人澡人人爽秒播| 欧美3d第一页| av国产免费在线观看| 午夜福利成人在线免费观看| 午夜福利18| 桃色一区二区三区在线观看| 丰满人妻熟妇乱又伦精品不卡| 日韩 欧美 亚洲 中文字幕| 美女扒开内裤让男人捅视频| 久久精品aⅴ一区二区三区四区| 国产成人啪精品午夜网站| 免费av毛片视频| 精品第一国产精品| 国产精品一区二区三区四区久久| 中文字幕久久专区| 亚洲中文av在线| 成人精品一区二区免费| 最近视频中文字幕2019在线8| 一本综合久久免费| а√天堂www在线а√下载| 色老头精品视频在线观看| 欧美黑人巨大hd| 亚洲成人免费电影在线观看| 制服诱惑二区| 欧美+亚洲+日韩+国产| 97超级碰碰碰精品色视频在线观看| 午夜福利18| 久久久久国产一级毛片高清牌| 国产精品香港三级国产av潘金莲| 国产精品 欧美亚洲| 91在线观看av| 亚洲国产看品久久| 亚洲欧美精品综合久久99| 日本免费一区二区三区高清不卡| 在线观看日韩欧美| 亚洲精品一区av在线观看| 狂野欧美激情性xxxx| 99re在线观看精品视频| 精品日产1卡2卡| 9191精品国产免费久久| 国产三级黄色录像| 亚洲自偷自拍图片 自拍| 精品久久久久久久末码| 中文亚洲av片在线观看爽| 欧美成狂野欧美在线观看| 日韩欧美 国产精品| 老司机午夜十八禁免费视频| 搡老岳熟女国产| 国产av一区在线观看免费| 亚洲精品美女久久久久99蜜臀| 国产激情偷乱视频一区二区| 黄色视频,在线免费观看| 岛国在线观看网站| 俺也久久电影网| 精品电影一区二区在线| 欧美3d第一页| 男女做爰动态图高潮gif福利片| 久久久久久久久免费视频了| 欧美日韩中文字幕国产精品一区二区三区| 搡老熟女国产l中国老女人| 99国产精品一区二区蜜桃av| 精品福利观看| 两性夫妻黄色片| 日韩精品青青久久久久久| 成年女人毛片免费观看观看9| 午夜福利欧美成人| 成人18禁高潮啪啪吃奶动态图| 哪里可以看免费的av片| 久久久久久久久中文| 国产亚洲精品久久久久5区| 曰老女人黄片| 两性夫妻黄色片| 欧美又色又爽又黄视频| 成年女人毛片免费观看观看9| 国产蜜桃级精品一区二区三区| 国内精品久久久久精免费| 嫩草影院精品99| 亚洲av五月六月丁香网| 香蕉国产在线看| 亚洲真实伦在线观看| 成人三级做爰电影| 香蕉av资源在线| 黄色毛片三级朝国网站| 18禁美女被吸乳视频| 中文亚洲av片在线观看爽| 女警被强在线播放| 丁香欧美五月| 午夜久久久久精精品| 久久久水蜜桃国产精品网| 无限看片的www在线观看| 中国美女看黄片| 亚洲全国av大片| 老司机福利观看| 国产精品一区二区三区四区久久| 久久精品国产综合久久久| 欧美黑人巨大hd| 老熟妇仑乱视频hdxx| 老司机福利观看| 男女床上黄色一级片免费看| 91九色精品人成在线观看|