程國(guó)勝，苑順周，趙蕾

(南京信息工程大學(xué)數(shù)理學(xué)院，江蘇南京210044)

磁層頂日下點(diǎn)距離R0與磁暴Dst指數(shù)的相關(guān)性

程國(guó)勝，苑順周，趙蕾

(南京信息工程大學(xué)數(shù)理學(xué)院，江蘇南京210044)

利用2004—2006年ACE、WIND衛(wèi)星觀測(cè)的太陽(yáng)風(fēng)數(shù)據(jù)和相應(yīng)時(shí)期反映磁暴大小的Dst指數(shù)，針對(duì)200個(gè)不同級(jí)別的磁暴事件，分析了磁層頂日下點(diǎn)距離R0與磁暴Dst指數(shù)的線性相關(guān)性。分析顯示，在極端太陽(yáng)風(fēng)條件下，Dst指數(shù)時(shí)間序列比借助于Chao Model計(jì)算出的磁層頂日下點(diǎn)距離R0的時(shí)間序列延遲了約3 h。經(jīng)修正時(shí)間延遲后，對(duì)磁層頂日下點(diǎn)距離R0與磁暴Dst指數(shù)的線性相關(guān)性進(jìn)行分析。結(jié)果表明:發(fā)生超級(jí)磁暴時(shí)，二者線性相關(guān)系數(shù)的均值為0.77;大磁暴時(shí)，線性相關(guān)系數(shù)的均值為0.74;中等磁暴時(shí)，線性相關(guān)系數(shù)的均值為0.47。此外，用最小二乘數(shù)據(jù)擬合得到了不同級(jí)別磁暴時(shí)的R0-Dst關(guān)系圖。

磁層頂日下點(diǎn)距離;磁暴;線性相關(guān)系數(shù);最小二乘擬合方法

0 引言

磁層頂是磁層與太陽(yáng)風(fēng)等離子體的邊界，決定著磁層的形狀，所以磁層頂?shù)奈恢脜?shù)在空間天氣研究中有著至關(guān)重要的作用。Chapman and Ferraro(1931)根據(jù)太陽(yáng)風(fēng)動(dòng)壓與磁層磁壓的平衡，首次提出磁層頂邊界的存在;隨后，F(xiàn)erraro(1952)描述了磁層頂?shù)男螤?。Aubry et al.(1970)注意到，行星際磁場(chǎng)的強(qiáng)度大小和方向也會(huì)影響磁層頂?shù)降厍虻木嚯x。隨著科技發(fā)展，人們采集的衛(wèi)星數(shù)據(jù)越來(lái)越多，越來(lái)越精確，一些空間物理學(xué)家們借助于已有數(shù)據(jù)分析，提出了多個(gè)磁層頂模型(Holzer and Slavin，1978;Sibeck et al.，1991;Roelof and Sibeck，1993;Shue et al.，1998;Boardsen et al.，2000;Kalegaev and Lyutov.，2000;Chao et al.，2002)，模擬了在不同太陽(yáng)風(fēng)條件下的磁層頂位置。通常太陽(yáng)風(fēng)條件下，磁層頂日下點(diǎn)距離R0大約有10RE～11RE(地球半徑);太陽(yáng)風(fēng)條件極弱時(shí)，R0可達(dá)14RE;極端太陽(yáng)風(fēng)條件時(shí)，磁層頂日下點(diǎn)會(huì)被壓縮到地球同步軌道以內(nèi)，即R0小于6.6RE(Kuznetsov and Suvorova，1998)。Plaschke et al.(2009)給出了磁層頂運(yùn)動(dòng)的統(tǒng)計(jì)分析研究。Yang et al.(2002)發(fā)現(xiàn)PR96模型(Petrinec and Russell，1996)、Shue98模型(Shue et al.，1998)和C02模型(Chao et al.，2002)適合應(yīng)用于模擬極端太陽(yáng)風(fēng)條件下磁層頂位形，其中C02模型能夠更好的預(yù)測(cè)極端太陽(yáng)風(fēng)條件下的磁層頂穿越事件。

磁層中的環(huán)電流是赤道面附近圍繞地球的一個(gè)環(huán)形電流帶，其主體部分是西向電流，磁暴期間總強(qiáng)度可達(dá)幾百萬(wàn)安培，因而引起地球磁場(chǎng)水平分量大幅度減小。然而，磁層頂日下點(diǎn)距離與磁層頂電流有關(guān)，根據(jù)等離子體物理的單粒子理論可知，磁層頂電流是由太陽(yáng)風(fēng)等離子體產(chǎn)生的一種抗磁電流，亦會(huì)引起地球磁場(chǎng)水平分量減小。擾動(dòng)一般限于高緯極區(qū)，但在行星際磁場(chǎng)具有長(zhǎng)期的南向分量且具有較大的幅度(大于10～15 nT)時(shí)，磁層連續(xù)受到壓力，磁場(chǎng)擾動(dòng)到達(dá)赤道區(qū)域。在中低緯度的測(cè)站使用的每一小時(shí)量測(cè)地磁擾動(dòng)強(qiáng)度的一種分級(jí)指標(biāo)，稱之為Dst指數(shù)，主要是量測(cè)地磁水平分量強(qiáng)度變化的一種物理量。Dst指數(shù)是研究磁暴期間磁場(chǎng)擾動(dòng)特征的一個(gè)關(guān)鍵參數(shù)，對(duì)分析地磁場(chǎng)活動(dòng)有重要的作用(劉振興和濮祖蔭，1997)，通常Dst＜－200 nT為超級(jí)磁暴，－200 nT≤Dst＜－100 nT為大磁暴，－100 nT≤Dst＜－50 nT為中等磁暴，－50 nT≤Dst＜－30 nT為小磁暴(Gonzalez et al.，1999)。同時(shí)，太陽(yáng)風(fēng)動(dòng)壓和行星際磁場(chǎng)是預(yù)測(cè)磁暴模型的主要輸入值(Wei et al.，2007)，也是日側(cè)磁層頂位置的主要影響因素，因此磁層頂日下點(diǎn)距離與磁暴指數(shù)之間必然存在某種關(guān)系。Shue et al.(2000)比較了低緯磁層頂位形模型后發(fā)現(xiàn)，眾多模型只考慮磁層頂位置變化受外在因素的影響，而沒(méi)有考慮磁層頂位置是否與磁層磁暴活動(dòng)有關(guān)聯(lián)的問(wèn)題。

針對(duì)Shue提出的上述問(wèn)題，迄今還未見(jiàn)到相關(guān)的研究。本文利用衛(wèi)星觀測(cè)到的數(shù)據(jù)，對(duì)磁層頂日下點(diǎn)距離與磁暴指數(shù)之間的關(guān)系進(jìn)行了統(tǒng)計(jì)分析。通過(guò)比較PR96、Shue98和C02模型所得的R0，發(fā)現(xiàn)C02模型更適合用來(lái)計(jì)算磁層頂日下點(diǎn)距離R0。利用2004年ACE、WIND衛(wèi)星提供的太陽(yáng)風(fēng)流和行星際磁場(chǎng)數(shù)據(jù)，基于C02模型，分析了不同級(jí)別磁暴發(fā)生時(shí)R0-Dst的相關(guān)性;然后，利用2005—2006年相應(yīng)數(shù)據(jù)，對(duì)所得結(jié)論進(jìn)行驗(yàn)證，發(fā)現(xiàn)在發(fā)生不同級(jí)別磁暴時(shí)，R0-Dst之間存在著線性相關(guān)性，磁暴越強(qiáng)，線性相關(guān)性越大。最后，按磁暴級(jí)別，對(duì)2004—2006年200個(gè)磁暴事件進(jìn)行統(tǒng)一處理，得到在超級(jí)、大、中等磁暴發(fā)生時(shí)線性相關(guān)系數(shù)的均值分別為0.77、0.74、0.47，并用最小二乘擬合得到不同級(jí)別磁暴時(shí)的R0-Dst關(guān)系式。

1 數(shù)據(jù)

所用的太陽(yáng)風(fēng)參數(shù)和行星際磁場(chǎng)的數(shù)據(jù)來(lái)源于ACE、WIND衛(wèi)星，分辨率為1 h。在GSM(geocentric solar magnetic system)坐標(biāo)系下，行星際磁場(chǎng)(IMF)的南向分量Bz、太陽(yáng)風(fēng)速度Vsw和密度Nsw數(shù)據(jù)主要來(lái)自ACE衛(wèi)星，對(duì)ACE衛(wèi)星缺失的數(shù)據(jù)通過(guò)WIND衛(wèi)星數(shù)據(jù)進(jìn)行補(bǔ)充。實(shí)時(shí)的地磁Dst指數(shù)數(shù)據(jù)來(lái)自于WDC(World Data Center for geomagnetism，Kyoto)，分辨率為1 h。以上資料均取自NASA的Space Physics Data Facility(SPDF)網(wǎng)站(http://cdaweb.gsfc.nasa.gov/cnofs/)。

2004年太陽(yáng)活動(dòng)極為頻繁，存在不同級(jí)別的磁暴，數(shù)據(jù)比較完整。根據(jù)磁暴級(jí)別定義進(jìn)行統(tǒng)計(jì)可知，2004年的82個(gè)磁暴，有小磁暴54個(gè)、中等磁暴22個(gè)、大磁暴5個(gè)和超級(jí)磁暴1個(gè)。利用2005—2006年的118個(gè)磁暴(小磁暴76個(gè)、中等磁暴31個(gè)、大磁暴10個(gè)和超級(jí)磁暴1個(gè))發(fā)生時(shí)的數(shù)據(jù)對(duì)2004年的結(jié)論進(jìn)行驗(yàn)證。

2 日側(cè)磁層頂模型比較

在眾多的磁層頂位置模型中，大都體現(xiàn)了磁層頂?shù)奈锢磉^(guò)程，并把這些物理過(guò)程轉(zhuǎn)化成為數(shù)學(xué)模型(Russell，2003)。例如，磁層頂日下點(diǎn)距離對(duì)于太陽(yáng)風(fēng)壓力的指數(shù)依賴性體現(xiàn)了磁偶極場(chǎng)的性質(zhì);磁層頂日下點(diǎn)距離對(duì)IMF Bz的非線性依賴性體現(xiàn)了磁層頂侵蝕的非線性飽和。大部分磁層頂位置模型應(yīng)用的是通常太陽(yáng)風(fēng)條件下磁層頂穿越事件的數(shù)據(jù)(Dmitriev et al.，2005)，而模型用在極端太陽(yáng)風(fēng)條件下時(shí)需要進(jìn)一步考慮。PR96模型、Shue98模型和C02模型是擬合磁層頂位置較為合適的模型(Yang et al.，2002)。下面分別介紹這3種模型的日側(cè)部分。

2.1 日側(cè)PR96模型

PR96模型的向陽(yáng)一側(cè)假設(shè)磁層頂關(guān)于XGSE(geocentric solar ecliptic system)軸旋轉(zhuǎn)對(duì)稱，且地心與橢圓交點(diǎn)重合，使用的是XGSE≥－10RE空間范圍內(nèi)低緯磁層頂穿越數(shù)據(jù)，給出向陽(yáng)側(cè)磁層頂位形模型如下:

其中:IMF北向時(shí)，m1=0;IMF南向時(shí)，m1=0.16。這里r代表磁層頂?shù)膹较蚓嚯x，θ是太陽(yáng)的天頂角;當(dāng)θ=0時(shí)，r=r0代表日下點(diǎn)距離。

2.2 日側(cè)Shue98模型

在處理上游太陽(yáng)風(fēng)數(shù)據(jù)時(shí)，采用固定太陽(yáng)風(fēng)傳輸時(shí)間，結(jié)合ISEE和IMP 8衛(wèi)星的數(shù)據(jù)，Shue98模型采用方程

描述了磁層頂?shù)奈恢煤托螤睢Ｆ渲?地球作為原點(diǎn);r代表磁層頂?shù)膹较蚓嚯x;θ是太陽(yáng)的天頂角;參數(shù)r0和α分別代表磁層頂日下點(diǎn)距離和磁層頂張角。

2.3 日側(cè)C02模型

C02模型采用了Shue98模型方程，且用分段函數(shù)來(lái)描述不同太陽(yáng)風(fēng)條件下磁層頂位形。假設(shè)在通常和極端太陽(yáng)風(fēng)條件下模型具有相同的表達(dá)式，r0的具體表達(dá)式為:

在通常太陽(yáng)風(fēng)條件下(r0≥7.0RE)，推導(dǎo)出的系數(shù)為:a1=11.646，a2=0.216，a3=0.122，a4=6.215;在極端太陽(yáng)風(fēng)條件下(r0≤6.4RE)，推導(dǎo)出的系數(shù)為:a1=11.646，a2=0.169，a3=0.158，a4=6.800;如果6.7RE＜r0＜7RE，則ln(r0)－ln(7)=c1(Dpc2);如果6.4RE＜r0＜6.7RE，則r0－6.4=c3(Dpc4)。由于r0和在r=7R處以及r=6.4R處連續(xù)，可分別求出c1和c2以及c3和c4。

2.4 模型的比較

下面將2004年4月的前40 h數(shù)據(jù)應(yīng)用于PR96、Shue98和C02模型，分別計(jì)算磁層頂日下點(diǎn)距離，然后進(jìn)行比較(圖1)。

圖1 PR96、Shue98和C02模型所計(jì)算的磁層頂日下點(diǎn)距離R0的比較(虛線:PR96;實(shí)線:Shue98;星形線:C02)Fig.1 Comparison of R0calculated by PR96，Shue98 and C02 models，respectively(The dash line:PR96;the solid line:Shue98;the asteroid line:C02)

Yang et al.(2002)對(duì)3種模型進(jìn)行比較后發(fā)現(xiàn)，C02模型能夠更好地預(yù)測(cè)極端太陽(yáng)風(fēng)條件下的磁層頂穿越事件。圖1顯示，3種模型在通常太陽(yáng)風(fēng)條件下區(qū)別不大，而在極端太陽(yáng)風(fēng)條件下區(qū)別明顯，其中C02模型能夠更好地描述磁層頂日下點(diǎn)位置。基于上述原因，本文采用C02模型來(lái)計(jì)算磁層頂日下點(diǎn)距離R0。

3 R0-Dst的相關(guān)性

下面分別考慮2004年的82個(gè)磁暴(超級(jí)磁暴1個(gè)、大磁暴5個(gè)、中等磁暴22個(gè)和小磁暴54個(gè))的Dst指數(shù)與磁層頂日下點(diǎn)距離R0的相關(guān)性。在用計(jì)算R-D之間的相關(guān)系數(shù)后，得出不同級(jí)別磁暴情況下的線性相關(guān)系數(shù)均值。

3.1 超級(jí)磁暴

對(duì)2004年太陽(yáng)風(fēng)數(shù)據(jù)流和地磁指數(shù)數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析，只找到2004年11月7—18日發(fā)生的一個(gè)超級(jí)磁暴。通過(guò)對(duì)此次磁暴Dst指數(shù)時(shí)間序列和R0時(shí)間序列比較發(fā)現(xiàn)，磁暴Dst指數(shù)時(shí)間序列比磁層頂日下點(diǎn)距離的時(shí)間序列延遲了約3 h。修正時(shí)間延遲后，磁暴Dst指數(shù)極小值與磁層頂日下點(diǎn)距離R0的極小值有明顯的對(duì)應(yīng)。圖2給出了修正時(shí)間延遲后磁層頂日下點(diǎn)距離R0與磁暴Dst指數(shù)的對(duì)比。

圖2 修正時(shí)間延遲后的磁暴Dst指數(shù)與磁層頂日下點(diǎn)距離R0的對(duì)比Fig.2 Contrast between Dstindex and magnetopause standard off distance R0after the lag time corrected

比較圖2a與2b發(fā)現(xiàn)，在這次超級(jí)磁暴中隨著時(shí)間的變化，磁暴Dst指數(shù)和磁層頂日下點(diǎn)距離R0的變化曲線趨勢(shì)一致。進(jìn)一步考慮Dst指數(shù)和R0之間的相關(guān)性，得到二者線性相關(guān)系數(shù)的均值是0.762。圖3給出了對(duì)此次超級(jí)磁暴266個(gè)R0-Dst數(shù)據(jù)集用最小二乘擬合方法得到的結(jié)果(直線為最小二乘擬合得到的結(jié)果;方程式y(tǒng)=0.018x+11是擬合關(guān)系式)。

圖3 超級(jí)磁暴期間Dst指數(shù)和磁層頂日下點(diǎn)距離R0的線性擬合Fig.3 The linear fitting between Dstindex and magnetopause standard off distance R0during the super storm

3.2 大磁暴

統(tǒng)計(jì)分析2004年太陽(yáng)風(fēng)數(shù)據(jù)流和地磁指數(shù)數(shù)據(jù)，一共有5個(gè)大磁暴。大磁暴發(fā)生時(shí)R0-Dst之間的線性相關(guān)系數(shù)均值為0.724(表1)。

表1 5個(gè)大磁暴發(fā)生時(shí)R0-Dst之間的線性相關(guān)系數(shù)Table 1 The linear correlation coefficients of R0-Dstduring 5 big storms，respectively

下面研究這5個(gè)大磁暴。對(duì)每個(gè)磁暴發(fā)生時(shí)Dst指數(shù)的時(shí)間序列和R0的時(shí)間序列進(jìn)行了對(duì)比分析，發(fā)現(xiàn)把Dst指數(shù)時(shí)間序列延遲R0時(shí)間序列3 h的情況下，它們之間相關(guān)系數(shù)最高;而延遲2 h或4 h的情況下，相關(guān)系數(shù)會(huì)明顯下降。因此得到結(jié)果如下:大磁暴發(fā)生時(shí)，磁暴Dst指數(shù)的時(shí)間序列比R0的時(shí)間序列延遲了約3 h。修正時(shí)間延遲后，Dst指數(shù)與R0變化曲線趨勢(shì)一致。進(jìn)一步對(duì)所有大磁暴事件的303個(gè)R0-Dst數(shù)據(jù)集，用最小二乘擬合方法進(jìn)行擬合。圖4表示的是大磁暴事件Dst指數(shù)和R0的線性擬合結(jié)果(擬合得到的函數(shù)關(guān)系式是y=0.028x+12)。

3.3 中等磁暴

圖4 大磁暴期間Dst指數(shù)和磁層頂日下點(diǎn)距離R0的線性擬合Fig.4 The linear fitting between Dstindex and magnetopause standard off distance R0during big storms

現(xiàn)在考慮2004年22個(gè)典型的中等磁暴，如表2所示。可見(jiàn)，中等磁暴發(fā)生時(shí)，R0-Dst之間的線性相關(guān)系數(shù)均值為0.453 4。分析這22個(gè)中等磁暴，對(duì)每個(gè)磁暴發(fā)生時(shí)Dst指數(shù)的時(shí)間序列和R0的時(shí)間序列進(jìn)行了對(duì)比分析，發(fā)現(xiàn)把Dst指數(shù)時(shí)間序列延遲R0時(shí)間序列3 h的情況下，它們之間相關(guān)系數(shù)最高;而延遲2 h或4 h的情況下，相關(guān)系數(shù)會(huì)明顯下降。因此得到結(jié)果如下:中等磁暴發(fā)生時(shí)，磁暴Dst指數(shù)時(shí)間序列比R0的時(shí)間序列也延遲了約3 h。修正時(shí)間延遲后，磁暴Dst指數(shù)和磁層頂日下點(diǎn)距離R0的變化曲線具有相當(dāng)程度的一致性。

進(jìn)一步利用所有中等磁暴事件的587個(gè)R0-Dst數(shù)據(jù)集，用最小二乘線性擬合方法進(jìn)行擬合。圖5表示的是中等磁暴事件Dst指數(shù)和R0之間的線性擬合結(jié)果(擬合得到的函數(shù)關(guān)系式是y=0.027x+11)。

表2 22個(gè)中等磁暴發(fā)生時(shí)R0-Dst之間的線性相關(guān)系數(shù)Table 2 The linear correlation coefficients of R0-Dstduring 22 middle storms，respectively

對(duì)2004年的54個(gè)小磁暴進(jìn)行研究。結(jié)果顯示，小磁暴發(fā)生時(shí)R0-Dst之間的線性相關(guān)系數(shù)均值僅為0.24左右，因此小磁暴事件發(fā)生時(shí)，Dst指數(shù)和R0之間沒(méi)有明顯的線性相關(guān)性。

圖5 中等磁暴期間Dst指數(shù)和磁層頂日下點(diǎn)距離R0的線性擬合Fig.5 The linear fitting between Dstindex and magnetopause standard off distance R0during middle storms

4 驗(yàn)證與討論

根據(jù)上節(jié)分析發(fā)現(xiàn)，在磁暴發(fā)生時(shí)R0-Dst之間存在著線性相關(guān)性，磁暴級(jí)別越強(qiáng)，線性相關(guān)性越大。那么，這個(gè)結(jié)果是否具有普適性呢?針對(duì)這個(gè)問(wèn)題，用2005—2006年118個(gè)不同級(jí)別磁暴數(shù)據(jù)，對(duì)上述結(jié)論進(jìn)行驗(yàn)證，處理方法與上節(jié)類(lèi)似，結(jié)果見(jiàn)表3。表3顯示，2005—2006年R0-Dst之間平均相關(guān)系數(shù)與2004年的結(jié)果相近，誤差值都在0.05以內(nèi);在超級(jí)磁暴、大磁暴、中等磁暴和小磁暴發(fā)生的情況下，誤差率分別為1.3%、5.5%、9.1%和8.2%。由此得出結(jié)論:在發(fā)生不同級(jí)別磁暴時(shí)，R0-Dst之間確實(shí)存在著線性相關(guān)性，且級(jí)別越強(qiáng)，相關(guān)性越大。

最后，將2004—2006年200個(gè)不同級(jí)別磁暴事件統(tǒng)一進(jìn)行處理。對(duì)表3中2004年和2005/2006年二者的平均相關(guān)系數(shù)取均值，得到如下結(jié)果:發(fā)生超級(jí)磁暴時(shí)，R0-Dst線性相關(guān)系數(shù)均值為0.77;發(fā)生大磁暴時(shí)，R0-Dst線性相關(guān)系數(shù)均值為0.74;發(fā)生中等磁暴時(shí)，R0-Dst線性相關(guān)系數(shù)均值為0.47。

表3 2004與2005/2006年R0-Dst之間平均相關(guān)系數(shù)的比較Table 3 Comparison of the average correlation coefficients in 2004 and 2005/2006

5 結(jié)論

本文利用C02模型得到極端太陽(yáng)風(fēng)條件下的磁層頂日下點(diǎn)距離R0，對(duì)R0與磁暴Dst指數(shù)進(jìn)行相關(guān)性分析后發(fā)現(xiàn)，兩者存在約3 h的延遲，也就是說(shuō)，從太陽(yáng)風(fēng)激波突然壓縮磁層頂?shù)接绊懘艑与娏黧w系(主要包括磁層頂電流和環(huán)電流)產(chǎn)生抗磁效應(yīng)所需要的時(shí)間大概為3 h。修正時(shí)間延遲后R0-Dst之間存在著線性相關(guān)性，磁暴級(jí)別越強(qiáng)，線性相關(guān)性越大，且得到了不同級(jí)別磁暴時(shí)的R0-Dst擬合關(guān)系。

本文結(jié)果將為研究磁暴活動(dòng)提供了一個(gè)新途徑，有助于認(rèn)識(shí)太陽(yáng)風(fēng)—磁層耦合過(guò)程。R0-Dst的線性相關(guān)性研究，為研究磁層頂位置和磁暴活動(dòng)提供了一個(gè)新思路，使得通過(guò)磁層磁暴活動(dòng)強(qiáng)度大小來(lái)反演日側(cè)磁層頂位置的遠(yuǎn)近成為可能。對(duì)此問(wèn)題，我們將進(jìn)行后續(xù)研究。另外，伴隨著Cluster、THEMIS等衛(wèi)星的升空，空間物理學(xué)家已經(jīng)對(duì)磁層頂位置進(jìn)行動(dòng)態(tài)研究(Zhang et al.，2009);以前大部分模型把磁層頂?shù)奈恢每醋鲀H僅依賴于上游太陽(yáng)風(fēng)和行星際磁場(chǎng)，事實(shí)上，隨著上游太陽(yáng)風(fēng)環(huán)境的改變，導(dǎo)致了磁層頂?shù)奈恢迷谄胶鈶B(tài)附近來(lái)回震蕩。考慮磁層頂運(yùn)動(dòng)速度時(shí)，磁層頂位置變化與磁暴的關(guān)系有待進(jìn)一步研究。

致謝:對(duì)提供幫助的中國(guó)科學(xué)院空間科學(xué)與應(yīng)用研究中心王赤研究員，臺(tái)灣國(guó)立中央大學(xué)太空科學(xué)研究所Jih-Kwan Chao(趙寄昆)教授，及提出建設(shè)性修改意見(jiàn)的審稿人，一并表示衷心感謝!

劉振興，濮祖蔭.1997.我國(guó)磁層物理研究的進(jìn)展和展望［J］.地球物理學(xué)報(bào)，40(增刊):7-19.

Aubry M B，Russell C T，Kivelson M G.1970.Inward motion of the magnetopause before substorm［J］.J Geophys Res，75(34):7018-7031.

Boardsen S A，Eastman T E，Sotirelis T，et al.2000.An empirical model of the high-latitude magnetopause［J］.J Geophys Res，105(A10):193-219.

Chao J K，Wu D J，Lin C H，et al.2002.Models for the size and shape of the earth's magnetopause and bow shock［C］//Lyu Ling-Hsiao.Cospar Colloquia series Vol.12:Space weather study using multipoint techniques.Pergamon:Elsevier Science Ltd.:127-134.

Chapman S，F(xiàn)erraro V C A.1931.A new theory of magnetic storm:1.The initial phase［J］.J Geophys Res，36(77/79):171-186.

Dmitriev A，Chao J K，Thomsen M，et al.2005.Geosynchronous magnetopause crossings on 29—31 October 2003［J］.J Geophys Res，110，A08209，doi:10.1029/2004JA010582.

Ferraro V C A.1952.On the theory of the first phase of a geomagnetic storm:A new illustrative calculation based on idealized(plane not cylindrical)model field distribution［J］.J Geophys Res，57(1):15-49.

Gonzalez W D，Tsurutani B T，Cluade A L.1999.Interplanetary origin of geomagnetic storms［J］.Space Science Reviews，88(3):529-562.

Holzer R E，Slavin J A.1978.Magnetic flux transfer associated with expansions and contractions of the dayside magnetopause［J］.J Geophys Res，83(A8):3831-3839.

Kalegaev V，Lyutov Y G.2000.The solar wind control of the magnetopause［J］.Adv Space Res，25(3/4):1489-1492.

Kuznetsov S N，Suvorova A V.1998.An empirical model of the magnetopause for board ranges of solar wind pressure and Bz IMF［C］//Polar cap boundary phenomena.Netherlands:Academic Publishers:120-131.

Petrinec S M，Russell C T.1996.Near-earth magnetotail shape and size as determined from the magnetopause flaring angle［J］.J Geophys Res，101(A1):137-152.

Plaschke F，Glassmeier K H，Auster H U，et al.2009.Statistical study of the magnetopause motion:First results from THEMIS［J］.J Geophys Res，114，A00C10，doi:10.1029/2008JA013423.

Roelof E C，Sibeck D G.1993.Magnetopause shape as a bivariate function of interplanetary magnetic field Bzand solar wind dynamic pressure［J］.J Geophys Res，98(A12):421-450.

Russell C T.2003.The structure of the magnetopause［J］.Planetary and Space Science，51(1):731-744.

Shue J H，Song P，Russell C T，et al.1998.Magnetopause location under extreme solar wind conditions［J］.J Geophys Res，103(A8):17691-17700.

Shue J H，Russell C T，Song P.2000.Shape of the low-latitude magnetopause:Comparison of models［J］.Adv Space Res，25(7):1471-1484.

Sibeck D G，Lopez R E，Roelof E C.1991.Solar wind control of the magnetopause shape，location，and motion［J］.J Geophys Res，96(A4):5489-5495.

Wei H L，Zhu D Q，Billings S A，et al.2007.Forecasting the geomagnetic activity of the Dst index using multiscale radial basis function networks［J］.Adv Space Res，40(12):1863-1870.

Yang Y H，Chao J K，Lin C H，et al.2002.Comparison of three magnetopause prediction models under extreme solar-wind conditions［J］.JGeophysRes，107(A1)，1008，doi:10.1029/2001JA000079.

Zhang H，Zong Q G，Sibeck D G，et al.2009.Dynamic motion of the bow shock and magnetopause observed by THEMIS spacecraft［J］.J Geophys Res，114，A00C12，doi:10.1029/2008JA013488.

Correlation between magnetopause standard off distance R0and geomagnetic storm Dstindex

CHENG Guo-sheng，YUAN Shun-zhou，ZHAO Lei

(School of Mathematics＆Physics，NUIST，Nanjing 210044，China)

With the solar wind data observed by ACE/WIND spacecraft and the geomagnetic storm Dstindex data during 2004—2006，the linear correlation between magnetopause standard off distance R0and Dstindex，in which the grade of geomagnetic storm is characterized，is investigated in terms of 200 different graded storms.It is a fact that the time series of Dstindex lag about 3 h behind the time series of the magnetopause standard off distance R0calculated by Chao Model under the extreme solar wind.After the time lag corrected，the linear correlation between magnetopause standard off distance R0and geomagnetic storm Dstindex is revealed.The results show that the average linear correlation coefficients are 0.77，0.74 and 0.47，respectively，corresponding to the super，the big and the middle storms.Under different graded storms，the charts of R0-Dstfitted by the least square method are shown.

magnetopause standard off distance;geomagnetic storm;linear correlation coefficient;least square method

P353

A

1674-7097(2011)04-0504-06

2010-10-15;改回日期:2010-12-09

國(guó)家自然科學(xué)基金資助項(xiàng)目(60973157)

程國(guó)勝(1963—)，男，安徽肥西人，博士，教授，研究方向?yàn)榭臻g天氣學(xué)、智能計(jì)算，chenggs@nuist.edu.cn.

程國(guó)勝，苑順周，趙蕾.2011.磁層頂日下點(diǎn)距離R0與磁暴Dst指數(shù)的相關(guān)性［J］.大氣科學(xué)學(xué)報(bào)，34(4):504-509.

Cheng Guo-sheng，Yuan Shun-zhou，Zhao Lei.2011.Correlation between magnetopause standard off distance R0and geomagnetic storm Dstindex［J］.Trans Atmos Sci，34(4):504-509.

(責(zé)任編輯:倪東鴻)