?？宋? 郭建鵬*, 何林峽, 林海博, 雷桓, 陳艷, 孔令高

1 北京師范大學(xué)天文與天體物理前沿科學(xué)研究所, 北京 1022062 北京師范大學(xué)天文系行星與空間物理研究團(tuán)組, 北京 1008753 南京大學(xué)(蘇州校區(qū)) 深空探測(cè)科學(xué)與技術(shù)研究院, 江蘇蘇州 2151634 中國(guó)科學(xué)院國(guó)家空間科學(xué)中心, 北京 100190

0 引言

火星引力場(chǎng)大約只有地球的2/5,并且火星弓激波的空間尺度小,距離火星近,因此火星的氫(H)、氧(O)逃逸層可以延伸至弓激波外.O逃逸層達(dá)到1～2個(gè)火星半徑;H的相對(duì)原子質(zhì)量較O更小,H逃逸層可至十幾個(gè)火星半徑.水手6號(hào)(Mariner-6)和水手7號(hào)(Mariner-7)衛(wèi)星上的紫外光譜儀在Lyman-α波段(121.6 nm附近)的觀測(cè)結(jié)果,首次證實(shí)了火星周圍存在延展的氫逃逸層(“氫冕”; Anderson and Hord, 1971;Anderson, 1974;Chaufray et al., 2008).逃逸層中的中性原子可以直接和太陽(yáng)風(fēng)相互作用,通過(guò)光致電離、電荷交換和電子碰撞形成新生離子,隨即被太陽(yáng)風(fēng)拾起,稱為“拾起”離子.拾起離子在洛倫茲力作用下圍繞行星際磁場(chǎng)做回旋運(yùn)動(dòng),同時(shí)受太陽(yáng)風(fēng)對(duì)流電場(chǎng)加速,運(yùn)動(dòng)軌跡呈現(xiàn)擺線.拾起離子是火星大氣逃逸的一個(gè)重要通道.火星向陽(yáng)側(cè)產(chǎn)生的拾起O+,由于其回旋半徑大于火星半徑,故會(huì)在太陽(yáng)風(fēng)對(duì)流電場(chǎng)作用下加速逃離火星.Jakosky等(2018)根據(jù)火星大氣和揮發(fā)物演化任務(wù)(Mars Atmosphere and Volatile EvolutioN, MAVEN)最新觀測(cè)結(jié)果得到,當(dāng)前拾起O+的損失率為5×1024s-1,約占總逃逸的8%左右.

火星弓激波上游的拾起離子在朝向火星運(yùn)動(dòng)過(guò)程中,可能在弓激波處被反射.Ip等(1992)使用Phobos-2攜帶的ASPERA質(zhì)譜儀的數(shù)據(jù),分析了火星附近的拾起H+的密度分布,并結(jié)合太陽(yáng)風(fēng)速度的變化,認(rèn)為被弓激波反射后的拾起H+會(huì)導(dǎo)致太陽(yáng)風(fēng)減速.隨后,Dubinin等(2006)利用MEX(Mars Express)衛(wèi)星,在火星弓激波上游觀測(cè)到拾起H+和能量大于10 keV的反射H+,他們認(rèn)為當(dāng)拾起H+以較小法向速度入射到弓激波時(shí),被弓激波的靜電勢(shì)反射回到太陽(yáng)風(fēng)中.而Yamauchi等(2011)認(rèn)為靠近火星弓激波的拾起H+可以進(jìn)行多重加速,多重加速的離子速度與被弓激波多次反射后的速度一致,進(jìn)而推測(cè)拾起H+是發(fā)生了多次鏡面反射后逃離弓激波.

O+的相對(duì)原子質(zhì)量比H+大,回旋半徑可達(dá)數(shù)個(gè)火星半徑,可以克服靜電勢(shì)穿越火星弓激波.Yamauchi等(2015)分析MEX數(shù)據(jù)時(shí)發(fā)現(xiàn)火星弓激波上游存在O+束(>1 keV),且O+的能量隨距離弓激波高度的增加而增大,他們認(rèn)為O+被磁鞘和太陽(yáng)風(fēng)的對(duì)流電場(chǎng)加速,其能量隨高度的增加而增大,逃離火星.隨后,基于MAVEN探測(cè)器數(shù)據(jù),Masunaga等(2016)在火星磁鞘中觀測(cè)到部分環(huán)狀分布的O+,同時(shí)在弓激波上游觀測(cè)到更高能的反射O+束(>10 keV),他們推斷弓激波上游的拾起O+進(jìn)入磁鞘后,在強(qiáng)磁場(chǎng)區(qū)域發(fā)生反射,在磁鞘電磁場(chǎng)的作用下,繞新的引導(dǎo)中心運(yùn)動(dòng),被對(duì)流電場(chǎng)加速返回到上游太陽(yáng)風(fēng)中.但該觀測(cè)只在短時(shí)間內(nèi)(128 s)在火星磁鞘中看到反射O+,難以判定發(fā)生反射的位置.隨后,Masunaga等(2017)對(duì)O+的反射率進(jìn)行估算,發(fā)現(xiàn)拾起O+在火星弓激波下游被反射的反射比約為14.1%±6.7%,極端太陽(yáng)風(fēng)條件下反射比增加到38.9%±10.2%.

關(guān)于火星拾起O+反射現(xiàn)象的解釋,存在兩種觀點(diǎn).一種觀點(diǎn)認(rèn)為,高能O+束可能是朝向火星方向沉降的太陽(yáng)風(fēng)拾起O+與逃逸層中O+碰撞產(chǎn)生的.在不顯著改變沉降離子總能量的情況下,碰撞可以使其速度方向發(fā)生偏轉(zhuǎn),故沉降的O+可能回旋到太陽(yáng)風(fēng)中.理論計(jì)算上,沉降的太陽(yáng)風(fēng)拾起O+在火星逃逸面(～200 km)之上發(fā)生碰撞的概率很小,故而產(chǎn)生的反射O+相空間密度應(yīng)小于沉降O+密度(Yagi et al., 2012).但觀測(cè)顯示反射O+和沉降O+密度相當(dāng)(Masunaga et al., 2016).另一種觀點(diǎn)認(rèn)為,能量較高的太陽(yáng)風(fēng)沉降O+和火星逃逸層中性粒子之間先發(fā)生電荷交換(Chen et al., 2001),成為高能中性原子后,再次被光致電離或者與激波下游的熱電子發(fā)生電子碰撞后被電離(Crider et al., 2000),兩次電離過(guò)程就可能使得高能拾起O+做回旋運(yùn)動(dòng)返回到太陽(yáng)風(fēng)中.這種情況下發(fā)生碰撞的平均自由程(≈1800RM,RM為火星半徑)遠(yuǎn)大于火星磁鞘的厚度(～1RM),電荷交換過(guò)程可以忽略,通過(guò)這種機(jī)制產(chǎn)生的反射O+難以被觀測(cè)到(昋若寓等, 2022),所以目前O+反射機(jī)制尚不明確.

相較于缺少磁力計(jì)的MEX,MAVEN攜帶高精度的磁場(chǎng)和離子探測(cè)儀,有助于我們進(jìn)一步認(rèn)識(shí)火星空間拾起離子的動(dòng)力學(xué)過(guò)程.若反射O+持續(xù)進(jìn)入探測(cè)器,將為我們追蹤反射發(fā)生的空間位置提供重要機(jī)遇.為此,我們分析2014—2021期間的 MAVEN觀測(cè)數(shù)據(jù),找到了一例這樣的觀測(cè)事件.具體而言,在2016年9月25日04∶18∶30—04∶24∶54 UT期間,MAVEN運(yùn)行在-E半球(電流電場(chǎng)E指向火星)的磁鞘區(qū),軌道高度逐漸降低,持續(xù)觀測(cè)到背離火星運(yùn)動(dòng)的O+(103～104eV)束流,即正在反射的拾起O+.本文將對(duì)該事件中反射O+束的來(lái)源進(jìn)行分析,并就拾起O+可能發(fā)生反射的位置提出一種合理的推斷.

1 MAVEN儀器數(shù)據(jù)

MAVEN衛(wèi)星在2013年11月成功發(fā)射,并于2014年9月進(jìn)入火星軌道.MAVEN衛(wèi)星為橢圓軌道,近火點(diǎn)約150 km,遠(yuǎn)火點(diǎn)約6400 km,軌道周期為4.5 h,軌道覆蓋面較廣,故可以對(duì)太陽(yáng)風(fēng)和火星磁鞘進(jìn)行全面采樣.在本文工作中,我們使用了超熱和熱離子成分分析儀(Supra thermal and Thermal Ion Composition, STATIC;McFadden et al., 2015)、磁力計(jì)(Magnetometer, MAG;Connerney et al., 2015)、太陽(yáng)風(fēng)離子分析儀(Solar Wind Ion Analyzer, SWIA;Halekas et al., 2015)、太陽(yáng)風(fēng)電子分析儀(Solar Wind Electron Analyzer, SWEA;Mitchell et al., 2016)四個(gè)儀器的原始數(shù)據(jù).

STATIC采用超環(huán)面靜電分析器,可以測(cè)量不同能量和角度分布中的離子組成.探測(cè)視場(chǎng)為360°×90°,角分辨率為22.5°×22.5°,探測(cè)能量范圍為0.1 eV～30 keV,能量分辨率為15%.STATIC 具有不同類型的數(shù)據(jù)產(chǎn)品:“d0”數(shù)據(jù)的時(shí)間分辨率為 32 s或128 s,含有32個(gè)能量段和8個(gè)質(zhì)量段;“c6”數(shù)據(jù)的時(shí)間分辨率為4 s,有32個(gè)能量段和64個(gè)質(zhì)量段.在本項(xiàng)工作中,“d0”數(shù)據(jù)用于繪制拾起O+的速度分布函數(shù)(Velocity Distribution Functions,VDFs)圖以及O+能譜圖,而“c6”數(shù)據(jù)主要給出了質(zhì)量通量圖,可以幫助我們識(shí)別火星感應(yīng)磁層邊界的位置,這是由于MAVEN進(jìn)入感應(yīng)磁層后,源自火星電離層的行星重離子通量升高.

MAG由兩個(gè)獨(dú)立的三軸磁通門磁力計(jì)傳感器組成,用于測(cè)量在火星-太陽(yáng)-軌道(Mars-Sun-Orbital, MSO)坐標(biāo)系下的磁場(chǎng)三分量,測(cè)量范圍在0.1～60000 nT,分辨率高達(dá)0.008 nT.MAG可以幫助識(shí)別MAVEN穿越弓激波和感應(yīng)磁層邊界的時(shí)刻.

SWIA采用帶靜電偏轉(zhuǎn)板的超環(huán)面靜電分析器,用于測(cè)量太陽(yáng)風(fēng)和磁鞘離子(主要是質(zhì)子)的數(shù)據(jù),幫助研究太陽(yáng)風(fēng)與火星高層大氣的相互作用.SWIA的時(shí)間分辨率為4 s,視場(chǎng)為360°×90°,角分辨率為22.5°×22.5°,能量范圍為25 eV～25 keV.在本文中,我們使用SWIA的“onboardsvymom”和“onboardsvyspec”數(shù)據(jù),在磁鞘中太陽(yáng)風(fēng)離子通量會(huì)顯著增強(qiáng),幫助我們識(shí)別MAVEN經(jīng)過(guò)火星弓激波的時(shí)間.

SWEA采用帶靜電偏轉(zhuǎn)板的半球形靜電分析器,用于測(cè)量5 eV～4.6 keV的電子能量和角度分布.SWEA的能量分辨率為17%,在1.6 keV能量處,探測(cè)視場(chǎng)為360°×120°,角度分辨率是22.5°×20°,另外,需要注意的是,SWEA部分視場(chǎng)會(huì)被衛(wèi)星本體遮擋.

為了研究拾起離子的運(yùn)動(dòng)狀態(tài)隨時(shí)間和方位的變化,我們引入了VDF.假設(shè)入射微分粒子通量在儀器中的能量響應(yīng)和角響應(yīng)內(nèi)是均勻的,采樣時(shí)間τ(s)內(nèi)檢測(cè)粒子數(shù)為C(vpqr),則

(1)

(2)

其中數(shù)值系數(shù)為(1×10-4)由m2到cm2變化而來(lái),由(1)和(2)式可得速度分布函數(shù)

(3)

再結(jié)合MAVEN觀測(cè)數(shù)據(jù),STATIC-“d0”提供了32個(gè)能量段、4個(gè)仰俯角、16個(gè)方位角的計(jì)數(shù),計(jì)算可得所有角度的離子相空間密度,進(jìn)而投影至垂直于磁場(chǎng)方向繪制出VDF圖,VDF圖中縱軸代表電場(chǎng)方向速度,橫軸代表垂直于磁場(chǎng)和電場(chǎng)方向的速度.

2 結(jié)果

圖1展示了2016年9月25日03∶50—04∶40 UT期間MAVEN的觀測(cè)數(shù)據(jù).圖中紅色豎直虛線(04∶05∶30 UT)標(biāo)識(shí)出MAVEN穿越火星弓激波的時(shí)刻,MAVEN從火星弓激波上游穿越至下游,觀測(cè)到磁場(chǎng)強(qiáng)度突然增強(qiáng),如圖1h磁場(chǎng)強(qiáng)度從5 nT增加到40 nT,且有劇烈擾動(dòng),太陽(yáng)風(fēng)速度減慢(圖1i),離子能譜(圖1a)和電子能譜(圖1b)顯著變寬,STATIC所有質(zhì)量能譜(圖1c)也顯著增寬,如圖1藍(lán)色豎直虛線(04∶29 UT)標(biāo)識(shí)出火星感應(yīng)磁層邊界.在感應(yīng)磁層邊界,可以看到主要離子組份從太陽(yáng)風(fēng)質(zhì)子(圖1e中的1 amu)轉(zhuǎn)變?yōu)樾行侵仉x子(O+和O2+對(duì)應(yīng)圖1e中的16 amu和32 amu;Brain et al., 2010),并且磁場(chǎng)強(qiáng)度擾動(dòng)變小,磁場(chǎng)增強(qiáng)(Nagy et al.,2004; Franco et al.,2022).MAVEN衛(wèi)星進(jìn)入感應(yīng)磁層后,由于主要粒子成分發(fā)生改變,這時(shí)計(jì)算的電場(chǎng)存在誤差,但不影響對(duì)本文事件的研究.

圖1 MAVEN在2016年9月25日03∶50—04∶40 UT的觀測(cè)結(jié)果,包含火星磁鞘內(nèi)氧離子(O+)反射事件從上至下依次為:(a)SWIA離子能譜、(b)SWEA電子能譜、(c)STATIC離子能譜、(d)STATIC中相對(duì)原子質(zhì)量為12～20的離子能譜、(e)STATIC質(zhì)譜、(f)—(g)STATIC中A3D1和A2D0方向的O+能譜、(h)—(j)MSO坐標(biāo)系下的磁場(chǎng)、太陽(yáng)風(fēng)速度、對(duì)流電場(chǎng)、(k)MAVEN衛(wèi)星高度和(l)MAVEN衛(wèi)星與火星之間的角度(即衛(wèi)星的太陽(yáng)天頂角). 紅色和藍(lán)色豎虛線分別表示弓激波位置和感應(yīng)磁層邊界,MAVEN在太陽(yáng)風(fēng)、磁鞘和磁層的探測(cè)時(shí)間段分別用粉色、青色和紫色條帶表示.

通過(guò)識(shí)別出弓激波和感應(yīng)磁層邊界,我們確定了磁鞘的位置.與太陽(yáng)風(fēng)中相比,磁鞘中離子能譜增寬,超熱電子顯著,由于磁場(chǎng)的壓縮導(dǎo)致對(duì)流電場(chǎng)的增強(qiáng)(Fang et al., 2008).在磁鞘中,04∶18∶30—04∶24∶54 UT期間STATIC儀器可以持續(xù)觀測(cè)到反射拾起O+事件,圖1f和圖1g選取了STATIC中兩個(gè)方向,我們發(fā)現(xiàn)在A2D0這個(gè)探測(cè)方向出現(xiàn)能量高于Emax(理論上拾起O+被對(duì)流電場(chǎng)加速可獲得的最大能量;白色實(shí)線),白色箭頭所指的拾起O+,我們認(rèn)為這部分是磁鞘中反射的拾起O+.

通過(guò)分析此期間MAVEN在MSO坐標(biāo)系下的運(yùn)行軌道,可以確定反射O+事件發(fā)生的空間區(qū)域.從圖2a可以看到反射事件發(fā)生在感應(yīng)磁層邊界的上方.圖2c和2d顯示MAVEN的空間位置位于火星日側(cè)南半球下方,且其運(yùn)行軌道朝向火星感應(yīng)磁層內(nèi)部移動(dòng).從圖1-l可知衛(wèi)星的太陽(yáng)天頂角大概位于60°左右.

圖2 MAVEN衛(wèi)星在2016年9月25日00∶07∶24—04∶40∶48 UT的運(yùn)行軌道圖,以火星半徑(RM)為單位此軌道#3874依次在(a)柱坐標(biāo)系平面、(c)YMSO-ZMSO平面和(d)XMSO-ZMSO平面的投影圖,分別用粉色、青色和紫色條帶表示其在太陽(yáng)風(fēng)、磁鞘和磁層的探測(cè)時(shí)間段. 紅色虛線和藍(lán)色虛線分別表示弓激波和感應(yīng)磁層邊界,紅色正方形和三角形分別表示磁鞘內(nèi)出現(xiàn)高能O+反射的開(kāi)始和結(jié)束時(shí)刻.

圖3(a—d)顯示此期間局地對(duì)流電場(chǎng)、磁場(chǎng)以及速度的方向,局地對(duì)流電場(chǎng)指向火星(-E半球).為了詳細(xì)研究O+的反射,圖3(e—h)是選取了04∶18∶30—04∶24∶54 UT期間的四個(gè)時(shí)間段內(nèi)的VDF.從VDF中我們能看到兩個(gè)離子群.第一個(gè)離子群是出現(xiàn)在紅色環(huán)上的新生拾起離子,這部分離子的相空間密度高(10-12s3·m-6)、能量較低,主要源自觀測(cè)點(diǎn)附近局地新生的拾起離子.第二個(gè)離子群是沿著反向火星和局地對(duì)流電場(chǎng)運(yùn)動(dòng)的O+束,可判斷為在磁鞘中發(fā)生反射的拾起O+,這部分反射O+的能量范圍較大.能量范圍1～10 keV的反射拾起O+的相空間密度較高(10-13s3·m-6);能量大于10 keV的反射拾起O+的相空間密度較低(10-14s3·m-6).此外,在圖3g中還出現(xiàn)一個(gè)朝向火星運(yùn)動(dòng)的高能低相空間密度的離子束,可判斷為從太陽(yáng)風(fēng)向火星沉降的O+束.

圖3 2016年9月25日04∶18∶30—04∶24∶54 UT期間,STATIC觀測(cè)得到的O+速度分布函數(shù)圖(VDFs)以及此時(shí)MAVEN在MSO坐標(biāo)系中的位置(a)—(d)紅線和藍(lán)線分別表示在XMSO-ZMSO平面內(nèi)弓激波和感應(yīng)磁層邊界的投影,該位置上的綠色、橘色和紫色線分別表示電場(chǎng)、質(zhì)子流速度和磁場(chǎng)的方向.(e)—(h)投影在垂直于局地磁場(chǎng)的平面上的O+ VDF圖,VE表示沿電場(chǎng)方向速度,VE×B表示垂直于磁場(chǎng)和電場(chǎng)方向速度,灰色圓從內(nèi)到外分別對(duì)應(yīng)100 eV、1 keV和10 keV,紅色虛線圓表示局地產(chǎn)生的拾起O+的環(huán)分布,由垂直于局地磁場(chǎng)的速度計(jì)算,空白部分為超出儀器視場(chǎng).

根據(jù)STATIC的64個(gè)探測(cè)方向的O+能譜(圖4),我們找出與VDF中反射拾起O+對(duì)應(yīng)的離子群,即A2D0、A2D1中出現(xiàn)的能量高于Emax的離子,對(duì)應(yīng)VDF中低相空間密度高能的反射O+;A4D1、A4D0、A3D1出現(xiàn)的離子可能對(duì)應(yīng)VDF中較高相空間密度低能的反射O+.

3 討論

基于MAVEN探測(cè)器數(shù)據(jù),Masunaga等(2016)在磁鞘中觀測(cè)到O+呈部分環(huán)狀分布,并且在火星弓激波上游觀測(cè)到高能O+束,他們給出解釋認(rèn)為太陽(yáng)風(fēng)中的O+注入磁鞘后,在磁鞘電磁場(chǎng)作用下,繞新的引導(dǎo)中心回旋,呈部分環(huán)狀分布,其中一部分O+被反射回到弓激波上游形成高能O+束,并被對(duì)流電場(chǎng)加速,逃離火星.然而需要注意的是,他們的觀測(cè)只在短時(shí)間內(nèi)看到了反射O+,故不能明確拾起O+在火星磁鞘內(nèi)發(fā)生反射的空間位置.本文中的事件則在火星磁鞘中持續(xù)觀測(cè)到能量超過(guò)Emax的反射O+束,對(duì)研究火星拾起O+的反射提供了新的觀測(cè)證據(jù).

對(duì)于拾起O+可能發(fā)生反射的空間位置,我們做如下討論:在04∶18∶30—04∶24∶54 UT期間,MAVEN衛(wèi)星位于-E半球,反射O+束的運(yùn)動(dòng)方向幾乎與對(duì)流電場(chǎng)方向相反,所以該O+束不會(huì)是受對(duì)流電場(chǎng)加速形成.從軌道運(yùn)行圖看,MAVEN衛(wèi)星的太陽(yáng)天頂角位于60°左右,排除了該O+是沿著磁尾逃逸的情況.Masunaga等(2016)的觀點(diǎn)認(rèn)為拾起O+可能是在磁鞘強(qiáng)磁場(chǎng)區(qū)域發(fā)生反射,但他們并未看到反射的具體位置.在本事件時(shí)間段中,MAVEN朝向火星感應(yīng)磁層邊界運(yùn)行,且從04∶19 UT開(kāi)始,磁場(chǎng)強(qiáng)度逐漸增強(qiáng)超過(guò)20 nT,磁場(chǎng)y分量的方向發(fā)生了改變,并持續(xù)觀測(cè)到反射O+束,由此我們推斷拾起O+是在火星感應(yīng)磁層邊界上方強(qiáng)磁場(chǎng)區(qū)域發(fā)生了反射.

4 結(jié)論

基于MAVEN觀測(cè),我們分析了一例在火星磁鞘中持續(xù)觀測(cè)到的反射O+事件,進(jìn)而推斷了拾起O+發(fā)生反射的可能位置.反射拾起O+大致可分為兩部分:一部分源自上游太陽(yáng)風(fēng)拾起的O+,具有一定的初速度,注入磁鞘后被反射,對(duì)應(yīng)于低相空間密度和較高能量的部分;另一部分源自火星磁鞘中局地新生的拾起O+,反射后形成高相空間密度和較低能量的部分.我們推斷拾起O+發(fā)生反射的位置在感應(yīng)磁層邊界附近的強(qiáng)磁場(chǎng)區(qū)域.本文的觀測(cè)事件補(bǔ)充了拾起O+在火星磁鞘中反射的觀測(cè)證據(jù),給出了拾起O+可能的反射位置,相關(guān)的反射機(jī)制還有待進(jìn)一步深入研究.未來(lái),結(jié)合中國(guó)的“天問(wèn)一號(hào)”對(duì)火星空間的探測(cè),將有助于揭示火星空間拾起離子動(dòng)力學(xué)過(guò)程及其物理機(jī)制.

致謝本工作使用的MAVEN 數(shù)據(jù)由美國(guó)宇航局(https:∥pds-ppi.igpp.ucla.edu/)提供.