黃子萱

摘要：本文選擇2008年至2012年，行星際磁場北向期間15個太陽風(fēng)動壓突增事件，研究地面磁場擾動信號（初始脈沖（PI）和主要脈沖（MI））的變化規(guī)律和空間分布。統(tǒng)計結(jié)果顯示地磁場H分量脈沖信號極性分布規(guī)律在向陽側(cè)（6：00-18：00MLT）與Araki模型預(yù)言的一致。

關(guān)鍵詞：初始脈沖;主要脈沖;空間分布;Araki模型

1 緒論

空間中太陽風(fēng)擾動，會引起地球磁層和電離層的大規(guī)模擾動，當(dāng)太陽風(fēng)動壓突增作用于地球磁層時，磁層頂首先被壓縮，隨后在向陽側(cè)磁層頂?shù)叵騻鞑サ目炷４帕黧w波被激發(fā)，當(dāng)磁流體波傳播到地面時，中低緯度區(qū)域地磁臺站一般會觀測到地磁場水平分量（H分量）響應(yīng)增強，即觸發(fā)地磁場 SI（Sudden Impulse），而在高緯地區(qū)的地磁臺站則觀測到地磁場H分量呈雙極性變化，即地磁場H分量先突然降低（增加）并快速恢復(fù)（平均持續(xù)2～3min），然后突然增加（降低）并持續(xù)一段時間，分布在不同磁地方時的地磁臺站觀測到的地磁場極性變化不同^[1]。通常，地磁場H分量的響應(yīng)可以分為兩個部分，其中前者稱為初始脈沖響應(yīng)（The Preliminary Impulse，PI），后者則稱為主脈沖響應(yīng)（The Main Impulse，MI）。針對這種地磁響應(yīng)現(xiàn)象，Araki提出了著名的磁暴急始（SC）模型來試圖解釋這種全球的現(xiàn)象，稱為Araki SC模型^[2]，解釋了PI和MI在不同磁地方時的分布變化規(guī)律。本文的目就是通過對2008-2012年間所有同時在磁層有THEMIS衛(wèi)星觀測的太陽風(fēng)動壓突增事件中地磁響應(yīng)的特征統(tǒng)計分析來檢驗Araki物理模型。

1.1?太陽風(fēng)動壓脈沖簡介

太陽上的活動過程，如耀斑爆發(fā)、日冕物質(zhì)拋射、冕洞（Coronal Hole）等不停地擾動著平靜的背景太陽風(fēng)，在行星際空間形成了很多擾動結(jié)構(gòu)，如動壓脈沖（Dynamic Pressure Pulse），行星際激波（Shock），磁云（Magnetic Cloud，MC）行星際日冕物質(zhì)拋射（ICME），磁流體力學(xué)波動（MHD Wave），以及各種間斷面結(jié)構(gòu)^[3]。其中最受人們關(guān)注的就是太陽風(fēng)動壓脈沖結(jié)構(gòu)（Dynamic Pressure Pulse），即太陽風(fēng)動壓會在很短時間內(nèi)（小于1min） 突然大幅的增強，即發(fā)生我們所謂的太陽動壓突增事件，突增的太陽風(fēng)動壓會引起地球磁層-電離層大范圍的擾動，該擾動的影響是遍布全球的，它能影響包括磁層以及電離層在內(nèi)幾乎所有的等離子區(qū)域和電流體系。

1.2 Araki 物理模型簡介

由行星際激波和各間斷面結(jié)構(gòu)所造成的太陽風(fēng)動壓突增，會引起全球的地磁急始脈沖事件（geomagnetic sudden commencement，SC），Araki[1994]曾對這種現(xiàn)象提出了最簡單并且最全面的解釋，他通過總結(jié)大量的觀測結(jié)果與對別人的理論研究，提出太陽風(fēng)動壓突增所引起的SC事件中在高緯度地區(qū)磁場的H分量即水平分量，由兩個連續(xù)的、方向相反的脈沖組成。其中如果PI出現(xiàn)的是上升趨勢，則記為正，簡稱PPI（the positive preliminary impulse），如果PI出現(xiàn)的是下降的趨勢，則記為負，簡稱PRI（the reverse preliminary impulse）。在極區(qū)較低緯度，晨側(cè)地磁響應(yīng)特征為先正后負，但是在昏側(cè)則為先負后正。在更高緯度，地磁響應(yīng)特征又不一樣，在晨側(cè)的時候則是先負后正，在昏側(cè)的時候是先正再負。中低緯度，地磁響應(yīng)呈現(xiàn)一個類似斜坡的增長，但是仍疊加兩個雙極脈沖結(jié)構(gòu)，只是隨緯度的降低，脈沖結(jié)構(gòu)的振幅也降低。在赤道地區(qū)，兩個脈沖結(jié)構(gòu)只出現(xiàn)在日側(cè)，夜側(cè)則很不明顯。

Araki將地磁突增事件中的地磁變化的特征響應(yīng)分解為兩部分^[2]：

公式中DL表示地磁場對太陽風(fēng)動壓的突增時的一個類似于直流信號的階躍響應(yīng)，DP表示的是兩個脈沖信號響應(yīng)，這兩個脈沖信號就指的是前文所提到的初始脈沖和主要脈沖響應(yīng)：

角的值表示的是圖的刻度值（nT），較細的實線和虛線分別代表的是DP和DL，粗線代表的是總的疊加場，水平虛線表示H=0。

2 數(shù)據(jù)與處理

2.1 太陽風(fēng)動壓脈沖結(jié)構(gòu)事件選擇

2.2?地磁場數(shù)據(jù)處理

研究使用了來自加拿大的地磁臺站（CARISMA）^[5]、THEMIS地磁臺站^[6]、阿拉斯加和國際極光地磁效應(yīng)監(jiān)測（IMAGE）的地磁臺站陣列、格林蘭地磁臺站陣列、Magnetometer Array for Cusp and Cleft Studies（MACCS）^[7]的地磁數(shù)據(jù)。本文選擇的90個地磁臺站全球分布情況見圖2。

利用程序下載上述各地磁臺站陣列中各臺站的地磁數(shù)據(jù)，時間分辨率為0.5至1min。同時為了更加清晰的辨別地磁響應(yīng)中PI和MI的極性，我們對所有地磁數(shù)據(jù)都進行了高通濾波，截止周期為1200s。

由于受環(huán)境或其他條件的影響，并不是所有地磁臺站的地磁場H、D分量對于太陽風(fēng)動壓突增的響應(yīng)都有明顯的結(jié)構(gòu)，為了避免人為錯判所帶來的誤差，我們對于波動較大，或者結(jié)構(gòu)特征不夠明顯的地磁響應(yīng)，均予以排除。在這里列舉了7個典型的地磁場H分量響應(yīng)作為例子說明，如圖3所示，臺站PINA、PTRS的地磁響應(yīng)從12：00UT開始，表現(xiàn)為突然的下降，然后再上升并持續(xù)一段時間，根據(jù)前文所描述，第一個突然下降的脈沖即為初始脈沖，第二個上升的脈沖即為主要脈沖，而臺站NAQ、SCO的響應(yīng)則是剛剛相反，表現(xiàn)為先上升再下降，在制作PI和MI的極性響應(yīng)記錄表時，我們將上升脈沖信號記為“+”，下降脈沖信號則記為“-”，而對于類似與KAR、FYKN、GJOA臺站這樣響應(yīng)波動較大或地磁響應(yīng)結(jié)構(gòu)不明顯，無法通過肉眼直接辨別。

3結(jié)果與討論

我們選擇發(fā)生在2012年5月21日19點37分的太陽風(fēng)動壓突增事件進行分析。此次事件中向日面部分臺站的地磁響應(yīng)列于圖4中，我們可以看到隨著臺站所在地磁緯度增高，地磁場H分量的PI/MI的極性逐漸發(fā)生反轉(zhuǎn)。晨側(cè)，緯度較低的地區(qū)，地磁場H分量的響應(yīng)特征為先正后負，而緯度較高的地區(qū)，地磁場H分量的響應(yīng)特征為先負后正?；鑲?cè)，地磁場H分量的PI和MI的極性變化與晨側(cè)剛好相反，即在緯度較低的區(qū)域，地磁場H分量的響應(yīng)特征為先負后正，緯度較高的地區(qū)，H分量的響應(yīng)特征為先正后負。為了對PI/MI的分布規(guī)律，有更直觀的了解，我們將本次事件中，地磁場H分量的PI極性展示在地磁坐標(biāo)系中（圖5），并將臺站GAKO、KIAN、GILL等10個臺站在地磁坐標(biāo)系中的具體位置標(biāo)注在其中。其中X軸和Y軸的分別代表磁當(dāng)?shù)貢r和地磁緯度，圖中黑點表示極性響應(yīng)為負，加號表示極性響應(yīng)為正。

圖5?左圖是2012年5月21日19：37的太陽風(fēng)動壓突增事件地磁場H分量的PI極性在地磁坐標(biāo)系中的分布。

右圖為模擬結(jié)果，引自Lam et al.[2001]^[4]。

為便于比較，將Lam et al.[2001]給出的地磁場H分量的PI極性的模擬計算結(jié)果列在圖5中右圖。可以看出，該事件與模擬計算結(jié)果在向日側(cè)（06：00-12：00 MLT、12：00-18：00 MLT）呈現(xiàn)相同的分布規(guī)律，而在夜側(cè)（00：00-06：00 MLT和18：00-24：00 MLT）兩者并不一致。此外，該事件晨側(cè)的PI極性發(fā)生反轉(zhuǎn)的分界線的位置～67°N（AACGM），低于模型所給的～75°N，因此可能位于閉合磁力線區(qū)。

4 結(jié)論

本文首先給出典型事例2012年5月21日19：37地磁場H分量的PI正負極性在地磁坐標(biāo)系中的分布規(guī)律，并與Lam et al.[2001]的模擬結(jié)果作了對比。發(fā)現(xiàn)在向陽側(cè)（06：00-12：00 MLT、12：00-18：00 MLT），極性分布和模型的結(jié)果一致，但是晨側(cè)正負極性的分界線（～67°N）低于模擬值（75°N），因此可能位于閉合磁力線區(qū)。

參考文獻

[1]姚麗，左平兵，劉振興，陳化然，地球磁層對太陽風(fēng)動壓脈沖結(jié)構(gòu)響應(yīng)的研究進展天文學(xué)進展，2012.28（4）

[2]Araki，T.（1994），A physical model of the geomagnetic sudden commencement，in Solar Wind Sources of Magnetospheric Ultra-Low-Frequency Waves，Geophys. Monogr. Ser.，vol. 81，edited by M. J Engebretson，K. Takahashi，and M. Scholer，pp. 183–?200，AGU，Washington，D.C

[3]李暉，行星際擾動對地球空間環(huán)境的影響研究，博士論文，空間科學(xué)與應(yīng)用研究中心，2011

[4]Lam，M. M.，and A. S. Rodger（2001），A case study test of Arakis physical model of geomagnetic sudden commencement，J .Geophys. Res.，106，13，135.

[5]Mann，I. R.，et al.（2008），The Upgraded?CARISMA?Magnetometer Array In the THEMIS?Era，Space?Sci. Rev.，141，413–451，doi：10.1007/s 11214-008-9457-6.

[6]Russell，C. T.，et al.（2008），THEMIS ground-based magnetometers，Space Sci. Rev.，141，389–412，doi：10.1007/s11214-008-9337-0.

[7]Hughes，W. J.，and M. J. Engebretson（1997），Spatial，MACCS：Magnetometer array for cusp and cleft studies，in Satellite-Ground Based Coordination Sourcebook，edited by M. Lockwood，M. N. Wild，and H. J. Opgenoorth，Rep. ESA SP-1198，p. 119.

（作者單位：北京市第十二中學(xué)）