• <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
    亚洲欧美日韩高清在线视频| 国产xxxxx性猛交| 午夜久久久在线观看| 久久久久久久午夜电影 | 人妻丰满熟妇av一区二区三区| 天天添夜夜摸| 亚洲国产精品合色在线| av福利片在线| 久久中文字幕一级| 少妇粗大呻吟视频| 亚洲中文av在线| 激情视频va一区二区三区| 在线十欧美十亚洲十日本专区| 国产1区2区3区精品| 国产精品亚洲av一区麻豆| 丝袜美腿诱惑在线| 在线视频色国产色| 国产精品免费视频内射| 久久午夜综合久久蜜桃| 桃红色精品国产亚洲av| 亚洲精品美女久久久久99蜜臀| 大型黄色视频在线免费观看| 亚洲av五月六月丁香网| 久久久久久久久免费视频了| 亚洲人成电影观看| 成人免费观看视频高清| 日韩精品青青久久久久久| 男女高潮啪啪啪动态图| 夜夜爽天天搞| 国产又色又爽无遮挡免费看| 久久人妻av系列| 中国美女看黄片| 国产成年人精品一区二区 | 久久久水蜜桃国产精品网| 国产av在哪里看| 亚洲人成网站在线播放欧美日韩| 国产精品一区二区免费欧美| 99精品久久久久人妻精品| 亚洲av电影在线进入| av在线播放免费不卡| 亚洲中文字幕日韩| 99久久99久久久精品蜜桃| 99香蕉大伊视频| 精品福利永久在线观看| 88av欧美| 91成人精品电影| 国产伦一二天堂av在线观看| 91大片在线观看| 亚洲免费av在线视频| 亚洲五月天丁香| 亚洲精品美女久久av网站| 亚洲少妇的诱惑av| 国产精品久久久久成人av| 成人亚洲精品av一区二区 | 两性午夜刺激爽爽歪歪视频在线观看 | 日本一区二区免费在线视频| 免费高清视频大片| 色哟哟哟哟哟哟| 国产精品98久久久久久宅男小说| 亚洲情色 制服丝袜| 黄色片一级片一级黄色片| 精品福利永久在线观看| 国产精品爽爽va在线观看网站 | 999久久久国产精品视频| 美女国产高潮福利片在线看| 日日摸夜夜添夜夜添小说| 无人区码免费观看不卡| cao死你这个sao货| 亚洲片人在线观看| 亚洲三区欧美一区| 热re99久久国产66热| 99久久综合精品五月天人人| 一级作爱视频免费观看| 在线看a的网站| 女性生殖器流出的白浆| 精品人妻在线不人妻| 波多野结衣高清无吗| 精品少妇一区二区三区视频日本电影| 亚洲中文日韩欧美视频| 色综合欧美亚洲国产小说| 久99久视频精品免费| 国产黄a三级三级三级人| 精品无人区乱码1区二区| 美女国产高潮福利片在线看| 黄片小视频在线播放| 精品免费久久久久久久清纯| 老汉色∧v一级毛片| 欧美性长视频在线观看| 美国免费a级毛片| av电影中文网址| 国产欧美日韩综合在线一区二区| 露出奶头的视频| 另类亚洲欧美激情| 国产精品乱码一区二三区的特点 | 99riav亚洲国产免费| 国产黄a三级三级三级人| 黄色女人牲交| 丁香欧美五月| 欧美日韩亚洲高清精品| 国产精品一区二区精品视频观看| 国产亚洲精品综合一区在线观看 | 99国产精品一区二区蜜桃av| 免费在线观看完整版高清| 日本免费一区二区三区高清不卡 | 91精品国产国语对白视频| av片东京热男人的天堂| a在线观看视频网站| 午夜福利在线免费观看网站| 日本精品一区二区三区蜜桃| 亚洲国产看品久久| 97碰自拍视频| 午夜视频精品福利| 欧美日韩黄片免| 成人特级黄色片久久久久久久| 又黄又爽又免费观看的视频| svipshipincom国产片| 久久久国产一区二区| 婷婷丁香在线五月| 久久久久国产一级毛片高清牌| 免费在线观看黄色视频的| 一级作爱视频免费观看| 69精品国产乱码久久久| 真人一进一出gif抽搐免费| 成年人黄色毛片网站| 免费看十八禁软件| 中国美女看黄片| 日韩大码丰满熟妇| 日韩高清综合在线| 18美女黄网站色大片免费观看| 久久亚洲真实| www.www免费av| 国产成人精品久久二区二区91| 亚洲狠狠婷婷综合久久图片| 久久精品影院6| 亚洲av美国av| 久久精品aⅴ一区二区三区四区| 精品国产一区二区三区四区第35| 男女床上黄色一级片免费看| 久久天堂一区二区三区四区| 欧美黑人精品巨大| 日韩成人在线观看一区二区三区| 精品电影一区二区在线| 国产精品98久久久久久宅男小说| 岛国在线观看网站| 在线十欧美十亚洲十日本专区| 中文字幕精品免费在线观看视频| 亚洲av成人一区二区三| 久久人人爽av亚洲精品天堂| 国产精品一区二区在线不卡| 高清在线国产一区| 亚洲欧美精品综合久久99| 国产极品粉嫩免费观看在线| 又大又爽又粗| 亚洲av第一区精品v没综合| 午夜福利,免费看| 少妇被粗大的猛进出69影院| 欧美+亚洲+日韩+国产| 夜夜看夜夜爽夜夜摸 | 日本黄色视频三级网站网址| 99国产精品99久久久久| 国产成人av激情在线播放| 久久国产精品人妻蜜桃| 正在播放国产对白刺激| 制服人妻中文乱码| 国产片内射在线| 精品久久久久久久久久免费视频 | 超碰成人久久| 色婷婷久久久亚洲欧美| av电影中文网址| 久久精品亚洲精品国产色婷小说| 久久精品亚洲精品国产色婷小说| 757午夜福利合集在线观看| 不卡av一区二区三区| 男人舔女人下体高潮全视频| 村上凉子中文字幕在线| 国产精品香港三级国产av潘金莲| 精品一品国产午夜福利视频| 国产精品自产拍在线观看55亚洲| 欧美日韩乱码在线| 久久精品影院6| 男女床上黄色一级片免费看| 99在线人妻在线中文字幕| 国产三级黄色录像| 亚洲人成77777在线视频| 日韩欧美免费精品| 黑丝袜美女国产一区| 大型av网站在线播放| 一个人观看的视频www高清免费观看 | 97人妻天天添夜夜摸| 色播在线永久视频| 欧美性长视频在线观看| 欧美日韩黄片免| 12—13女人毛片做爰片一| 亚洲精品国产区一区二| 男人舔女人的私密视频| 无遮挡黄片免费观看| 国产成人系列免费观看| 国产高清激情床上av| 黄片小视频在线播放| 99国产精品免费福利视频| 欧美激情极品国产一区二区三区| 在线天堂中文资源库| 美女高潮到喷水免费观看| 精品国产乱码久久久久久男人| 久久99一区二区三区| 免费女性裸体啪啪无遮挡网站| 99国产精品一区二区蜜桃av| 母亲3免费完整高清在线观看| 日韩视频一区二区在线观看| 天堂中文最新版在线下载| 99国产综合亚洲精品| 日韩av在线大香蕉| a级毛片在线看网站| 亚洲一区二区三区欧美精品| 国产精品成人在线| 精品一区二区三区四区五区乱码| 操美女的视频在线观看| 精品午夜福利视频在线观看一区| 亚洲欧洲精品一区二区精品久久久| 黄片播放在线免费| 亚洲一区二区三区不卡视频| 亚洲色图av天堂| 国产99白浆流出| 乱人伦中国视频| 国产一区在线观看成人免费| 88av欧美| 久久中文字幕人妻熟女| 久久精品91无色码中文字幕| 亚洲精华国产精华精| 欧美人与性动交α欧美精品济南到| 国产成人一区二区三区免费视频网站| 97超级碰碰碰精品色视频在线观看| 亚洲avbb在线观看| 久久国产精品人妻蜜桃| 日日摸夜夜添夜夜添小说| 香蕉久久夜色| 亚洲三区欧美一区| 日韩欧美国产一区二区入口| av在线播放免费不卡| 日本黄色日本黄色录像| 精品国产乱码久久久久久男人| 又大又爽又粗| 热99re8久久精品国产| 成人亚洲精品av一区二区 | 老熟妇仑乱视频hdxx| 亚洲国产欧美日韩在线播放| 欧美成人午夜精品| 变态另类成人亚洲欧美熟女 | 香蕉久久夜色| 一二三四在线观看免费中文在| 伊人久久大香线蕉亚洲五| 亚洲片人在线观看| 亚洲精品国产精品久久久不卡| 51午夜福利影视在线观看| 一区在线观看完整版| 国产区一区二久久| 欧美日韩亚洲国产一区二区在线观看| 国产精品一区二区免费欧美| 日韩免费高清中文字幕av| 欧美日韩黄片免| 男人舔女人下体高潮全视频| 日韩欧美一区二区三区在线观看| 国产一区二区三区综合在线观看| 日日干狠狠操夜夜爽| 高清毛片免费观看视频网站 | 麻豆一二三区av精品| 99国产精品一区二区三区| 国产亚洲精品久久久久久毛片| 日韩免费av在线播放| 欧美黑人精品巨大| 亚洲熟妇中文字幕五十中出 | 国产精品 欧美亚洲| 免费高清视频大片| 岛国视频午夜一区免费看| 欧美日韩中文字幕国产精品一区二区三区 | 国产黄色免费在线视频| 久久精品91无色码中文字幕| 国产熟女午夜一区二区三区| 欧美中文日本在线观看视频| 日本三级黄在线观看| 50天的宝宝边吃奶边哭怎么回事| 精品卡一卡二卡四卡免费| 久久香蕉精品热| 国产亚洲精品第一综合不卡| 国产一区二区激情短视频| 电影成人av| 一二三四社区在线视频社区8| 午夜视频精品福利| 夜夜看夜夜爽夜夜摸 | 啦啦啦在线免费观看视频4| 中文字幕色久视频| 男人操女人黄网站| 久久久久久久久久久久大奶| 嫩草影院精品99| 国产精品二区激情视频| 国产av一区二区精品久久| 亚洲国产毛片av蜜桃av| 神马国产精品三级电影在线观看 | 国产一区二区三区在线臀色熟女 | 国产精品久久久久久人妻精品电影| 老汉色av国产亚洲站长工具| 啦啦啦免费观看视频1| 涩涩av久久男人的天堂| 免费在线观看完整版高清| 黄色视频不卡| 一级片免费观看大全| 精品国产一区二区久久| 国产精品一区二区在线不卡| 好看av亚洲va欧美ⅴa在| 亚洲 欧美 日韩 在线 免费| av在线播放免费不卡| 国产av精品麻豆| 国产成人一区二区三区免费视频网站| av网站在线播放免费| 在线观看免费午夜福利视频| 午夜激情av网站| 18禁国产床啪视频网站| 亚洲伊人色综图| 欧美日韩av久久| 国产伦一二天堂av在线观看| 一进一出抽搐gif免费好疼 | 亚洲成国产人片在线观看| 在线观看免费高清a一片| 欧美成人午夜精品| 一a级毛片在线观看| 99国产综合亚洲精品| 成年人黄色毛片网站| 成人亚洲精品av一区二区 | 精品人妻在线不人妻| 久久久久久大精品| 青草久久国产| 欧美一级毛片孕妇| 又黄又粗又硬又大视频| 国产99白浆流出| 精品国产一区二区三区四区第35| 亚洲人成77777在线视频| 女警被强在线播放| 两个人看的免费小视频| 国产免费男女视频| 男人舔女人的私密视频| www.999成人在线观看| 日本a在线网址| 黄色片一级片一级黄色片| 丰满迷人的少妇在线观看| 欧美人与性动交α欧美软件| 免费看a级黄色片| 欧美激情极品国产一区二区三区| 亚洲精华国产精华精| 亚洲成a人片在线一区二区| 一区二区三区国产精品乱码| 18禁观看日本| a级片在线免费高清观看视频| 性色av乱码一区二区三区2| 波多野结衣一区麻豆| 99久久综合精品五月天人人| 脱女人内裤的视频| 激情视频va一区二区三区| 久久中文看片网| 丝袜人妻中文字幕| 一级片'在线观看视频| 色老头精品视频在线观看| 亚洲第一av免费看| 大陆偷拍与自拍| bbb黄色大片| 国产成人av激情在线播放| 亚洲中文字幕日韩| 亚洲第一青青草原| 宅男免费午夜| 国产高清激情床上av| 欧美黑人欧美精品刺激| 免费人成视频x8x8入口观看| 国内毛片毛片毛片毛片毛片| 成人免费观看视频高清| 欧美黑人欧美精品刺激| 亚洲精品美女久久av网站| 99在线人妻在线中文字幕| 精品国产乱码久久久久久男人| 久久精品成人免费网站| 后天国语完整版免费观看| 热99re8久久精品国产| 精品日产1卡2卡| 久久影院123| 中文字幕高清在线视频| 又黄又爽又免费观看的视频| 精品少妇一区二区三区视频日本电影| 亚洲伊人色综图| 嫩草影院精品99| 精品午夜福利视频在线观看一区| 国产精品99久久99久久久不卡| 两个人免费观看高清视频| 夜夜爽天天搞| 淫妇啪啪啪对白视频| 18禁观看日本| 午夜福利免费观看在线| 精品国产一区二区久久| 一本大道久久a久久精品| 巨乳人妻的诱惑在线观看| 大型av网站在线播放| 国产亚洲欧美精品永久| 又黄又爽又免费观看的视频| 欧美午夜高清在线| 亚洲人成电影观看| 精品一区二区三区av网在线观看| 每晚都被弄得嗷嗷叫到高潮| 亚洲狠狠婷婷综合久久图片| 免费不卡黄色视频| 在线观看舔阴道视频| 国产真人三级小视频在线观看| 一个人观看的视频www高清免费观看 | 亚洲人成电影免费在线| 黑人操中国人逼视频| 亚洲片人在线观看| 国产精品成人在线| 在线国产一区二区在线| 日韩 欧美 亚洲 中文字幕| 亚洲人成电影观看| 国产精品野战在线观看 | 久久久精品国产亚洲av高清涩受| 欧美日韩国产mv在线观看视频| 丁香欧美五月| 日韩av在线大香蕉| 制服人妻中文乱码| 亚洲全国av大片| 久久亚洲真实| 国产精品电影一区二区三区| 日本五十路高清| 中国美女看黄片| 男男h啪啪无遮挡| 国产野战对白在线观看| 久久精品成人免费网站| 真人做人爱边吃奶动态| 国产在线观看jvid| 中文欧美无线码| 老司机深夜福利视频在线观看| 国产一卡二卡三卡精品| 麻豆国产av国片精品| 91九色精品人成在线观看| avwww免费| 别揉我奶头~嗯~啊~动态视频| 一本综合久久免费| 亚洲第一青青草原| 午夜老司机福利片| 国产乱人伦免费视频| 成人免费观看视频高清| 亚洲 国产 在线| 看免费av毛片| 岛国在线观看网站| 午夜福利,免费看| 亚洲人成伊人成综合网2020| 91精品国产国语对白视频| 男女午夜视频在线观看| 丝袜美腿诱惑在线| av天堂久久9| 亚洲 欧美 日韩 在线 免费| 琪琪午夜伦伦电影理论片6080| 99精品久久久久人妻精品| 国产精品久久久久成人av| 村上凉子中文字幕在线| 99riav亚洲国产免费| 天堂影院成人在线观看| 精品无人区乱码1区二区| 欧美最黄视频在线播放免费 | 免费人成视频x8x8入口观看| 欧美大码av| 欧美亚洲日本最大视频资源| av在线天堂中文字幕 | 日韩人妻精品一区2区三区| 欧美久久黑人一区二区| 热re99久久精品国产66热6| 午夜福利在线免费观看网站| 后天国语完整版免费观看| 国产在线精品亚洲第一网站| 久久久久亚洲av毛片大全| 国产午夜精品久久久久久| 少妇的丰满在线观看| 国产伦一二天堂av在线观看| 精品人妻1区二区| 99国产精品一区二区蜜桃av| www.精华液| 成人av一区二区三区在线看| 国产成人影院久久av| 两个人看的免费小视频| 亚洲成人精品中文字幕电影 | a级毛片在线看网站| 久久热在线av| 日本免费一区二区三区高清不卡 | 国内久久婷婷六月综合欲色啪| 亚洲专区字幕在线| 最好的美女福利视频网| 亚洲精品在线观看二区| 一a级毛片在线观看| 久久久久久久久久久久大奶| 国产精品永久免费网站| 99国产极品粉嫩在线观看| 黄频高清免费视频| 欧美乱妇无乱码| 欧美精品啪啪一区二区三区| 亚洲第一av免费看| 如日韩欧美国产精品一区二区三区| 亚洲 欧美一区二区三区| 91成人精品电影| 亚洲狠狠婷婷综合久久图片| 国产欧美日韩综合在线一区二区| 一边摸一边抽搐一进一出视频| 99久久人妻综合| 国产成人影院久久av| 免费在线观看日本一区| 日韩精品免费视频一区二区三区| 久久久久九九精品影院| 午夜福利影视在线免费观看| 天天影视国产精品| 欧美日韩福利视频一区二区| 国产亚洲精品综合一区在线观看 | 最近最新中文字幕大全免费视频| 精品午夜福利视频在线观看一区| 亚洲在线自拍视频| 欧美日韩亚洲高清精品| 后天国语完整版免费观看| 97碰自拍视频| xxx96com| 人人澡人人妻人| 欧美日韩视频精品一区| 亚洲欧美精品综合久久99| 免费一级毛片在线播放高清视频 | av国产精品久久久久影院| 日韩免费av在线播放| 国内毛片毛片毛片毛片毛片| 中国美女看黄片| 国产色视频综合| 日韩大码丰满熟妇| 在线观看免费日韩欧美大片| 大陆偷拍与自拍| 亚洲欧美激情综合另类| 搡老熟女国产l中国老女人| 涩涩av久久男人的天堂| 国产精品爽爽va在线观看网站 | 人人澡人人妻人| 成人亚洲精品一区在线观看| 亚洲一区二区三区不卡视频| 国产伦人伦偷精品视频| 亚洲熟女毛片儿| 视频区图区小说| 嫩草影视91久久| 欧美丝袜亚洲另类 | 涩涩av久久男人的天堂| 黄色视频不卡| 一边摸一边抽搐一进一小说| 久久热在线av| 无遮挡黄片免费观看| 人成视频在线观看免费观看| 精品一区二区三区视频在线观看免费 | 午夜免费成人在线视频| 在线观看一区二区三区| 国产精品九九99| x7x7x7水蜜桃| 成人亚洲精品av一区二区 | 新久久久久国产一级毛片| 制服诱惑二区| 午夜视频精品福利| 欧美成人免费av一区二区三区| 变态另类成人亚洲欧美熟女 | 国产区一区二久久| 丰满饥渴人妻一区二区三| 色婷婷av一区二区三区视频| 国产一区在线观看成人免费| 高潮久久久久久久久久久不卡| 精品国产乱码久久久久久男人| 欧美国产精品va在线观看不卡| 最近最新中文字幕大全免费视频| 老司机在亚洲福利影院| 69av精品久久久久久| 在线观看66精品国产| 国产黄色免费在线视频| 成人特级黄色片久久久久久久| 久久久久久亚洲精品国产蜜桃av| 国产精品久久久人人做人人爽| 看片在线看免费视频| 久久中文字幕人妻熟女| 一区二区三区激情视频| 久久久久久亚洲精品国产蜜桃av| 午夜福利欧美成人| 一级毛片女人18水好多| 一个人免费在线观看的高清视频| 免费一级毛片在线播放高清视频 | 电影成人av| 亚洲专区中文字幕在线| 国产精品美女特级片免费视频播放器 | 老熟妇乱子伦视频在线观看| 高清黄色对白视频在线免费看| 变态另类成人亚洲欧美熟女 | av网站在线播放免费| 亚洲av成人av| 99精品在免费线老司机午夜| av网站免费在线观看视频| 好男人电影高清在线观看| 女人精品久久久久毛片| 国产欧美日韩精品亚洲av| 日韩精品中文字幕看吧| av网站免费在线观看视频| 久久久久亚洲av毛片大全| www.熟女人妻精品国产| 国产精品一区二区三区四区久久 | 欧美日韩福利视频一区二区| 亚洲国产欧美日韩在线播放| 日韩高清综合在线| 真人做人爱边吃奶动态| 亚洲少妇的诱惑av| 国产精品二区激情视频| 国产黄色免费在线视频| 日韩欧美在线二视频| 亚洲专区国产一区二区| 亚洲av日韩精品久久久久久密| 久久热在线av| 亚洲精品一区av在线观看| 午夜精品久久久久久毛片777|