검색 전체 메뉴
PDF
맨 위로
OA 학술지
Variation of Magnetic Field (By , Bz) Polarity and Statistical Analysis of Solar Wind Parameters during the Magnetic Storm Period
  • 비영리 CC BY-NC
  • 비영리 CC BY-NC
ABSTRACT
Variation of Magnetic Field (By , Bz) Polarity and Statistical Analysis of Solar Wind Parameters during the Magnetic Storm Period
KEYWORD
disturbance storm time index , magnetic storms , solar wind parameters
  • 1. INTRODUCTION

    It is generally believed that the occurrence of a magnetic storm depends upon the solar conditions, particularly the southward interplanetary magnetic field (IMF) component. Some researchers have suggested a relationship between magnetic storms and solar wind parameters (Weigel 2010). Wu & Lepping (2006) investigated geomagnetic activity induced by interplanetary magnetic cloud (MC) during the past four solar cycles, 1965~1998, and found that the intensity of geomagnetic storms is more severe in a solar active period than in a solar quiet period.

    As suggested by some researchers, it is important to discuss the interplanetary drivers causing geomagnetic storms. Echer et al. (2008) also identified the interplanetary causes of intense geomagnetic storms and their solar dependence occurring during the solar cycle 23 (1996~2006). They concluded that the most important interplanetary structures leading to intense southward Bz are MC, which drove fast shocks causing 24% of the magnetic storms, and sheath fields, which also caused 24% of the storms. Borovsky & Denton (2006) showed twenty-one differences between coronal mass ejection (CME)-driven storms and co-rotating interaction region (CIR)-driven storms by table. According to their findings, CME-driven storms have denser plasma sheets, stronger ring currents, disturbance storm time (Dst) perturbation, and solar energetic particle events, while CIR-driven storms have longer duration and hotter plasma sheets.

    On the other hand, many studies have examined the role of solar wind parameters in the development of geomagnetic storms. Choi et al. (2008) investigated the statistical characteristics of solar wind dynamic pressure enhancements during geomagnetic storms. In their study, about 81% of storms studied indicate at least one event of dynamic pressure enhancements. Weigel (2010) quantified the influence of solar wind density using two statistical measures, and concluded that the solar wind density modifies the ability of a given value of the solar wind electric field to create a Dst perturbation. Lopez et al. (2004) also showed that high solar number density causes a change in the compression ratio of the bow shock for strong and southward IMF, which is typically associated with geomagnetic storms and an enhanced ring current.

    It is well known that magnetic storms are caused by enhanced magnetospheric convection during a prolonged period of southward IMF (Burton et al. 1975). On the other hand, Du et al. (2008) studied a geomagnetic storm that occurred on 21-22 January 2005, and found that magnetic storm is highly anomalous because the storm main phase developed during northward IMFs. They have interpreted this case as energy storage in the magnetotail and delayed release into the magnetosphere.

    After classifying MC polarity, Echer et al. (2005) presented a statistical study of MC parameters and geoeffectiveness.Considering that the IMF polarity of developing magnetic storms is important, it would be interesting to examine the influence of polarity change during magnetic storms. To identify the dependence of IMF polarity change, magnetic storms are divided into CME-driven storms, CIR-driven storms, and complicated types. Then, variation of magnetic field polarity (By, Bz) is analyzed during the main phase of magnetic storms. Echer et al. (2005) utilized only N-S and S-N polarity of MC, but in this research polarity change for the CME and CIR-driven storms is studied separately using the 12 types specified by IMF (By, Bz). The statistical relationship between solar wind parameters is also examined, including IMF Bz component and magnetic storms during the storm period by time-integrating. The time-integrated solar wind parameters are compared with the time-integrated Dst index during the storm period. The time lag in terms of hour from Bzmin to Dstmin during the main phases of CME and CIR-driven storms is also investigated.

    2. DATA AND METHOD

    Hourly average values of solar wind plasma and magnetic field data were obtained from National Aeronautics and Space Administration (NASA)'s OMNIWeb site (http://nssdc.gsfc.nasa.gov/omniweb) during 1997~2003. Magnetic field components are given in the GSM. To analyze the variation of IMF polarity and solar wind parameters, these data are used during the period of magnetic storms. Storms are classified according to their solar wind driver, including CME and CIRs. According to Borovsky & Denton(2006), CME-driven storms include CME sheaths, MC, and ejecta, while CIR-driven storms contain high-speed streams. Sheath fields present in the region between the interplanetary shock and the interplanetary CME (ICME). MC is characterized based on the following criteria: strong magnetic field, rotation of a large angle direction, low proton temperature, and low plasma beta (Burlaga et al. 1981). MC is classified as ejecta if it is hard to identify the rotation of magnetic field (Oh et al. 2007). CIR corresponds to an interaction region between low and high solar wind speed. It is also characterized by high proton temperature, and low number density.

    From some reference criteria, CME and CIR-driven storms were selected through an inspection of the data available OMNIWeb site. Through this process, a total of 156 storms were selected from the period of 1997~2003, that is, 109 (69.9%) CME-driven storms, 29 (18.6%) CIR-driven storms, and 18 (11.5%) complicated types. As mentioned above, complicated types were not included in this analysis.

    The Dst index data for this analysis are provided by the World Data Center for Geomagnetism, Kyoto (http://wdc.kugi.kyoto-u.ac.jp). Storms are divided into three categories in terms of the minimum value of the Dst index, Dstmin: severe (Dstmin < -200 nT), intense (-200 nT ≤ Dstmin < -100 nT), moderate (-100 nT ≤ Dstmin < -50 nT). A

    [Table 1.] List of magnetic storms used in this study.

    label

    List of magnetic storms used in this study.

    total of 156 magnetic storms, including 11 severe storms, 38 intense storms, and 107 moderate storms, were selected from the period of 1997~2003. Table 1 shows the list of storms used in this study during the main phase of magnetic storms. Assuming that the ring current provides the main contribution to the Dst index, Dst index is utilized for this analysis. The main phase of a magnetic storm is defined by the interval when the Dst index becomes less than -15 nT and reaches to the minimum value, while the late period is from the minimum to the moment when it recovers to -15 nT. As noted above, only magnetic storms with the Dst index exceeding -50 nT at the minimum value were used. Furthermore, only clearly recognizable storms with a single minimum and monotonically recovering storms in terms of the Dst index were chosen. Thus, the number of storms differs from the data presented in other papers because the criteria for storm selection are different (Choi et al. 2008).

    3. RESULTS

       3.1 Variation of Magnetic Field (By, Bz) Polarity during the Main Phase of Magnetic Storms

    Tsurutani et al. (1988) noted that magnetic storms are caused by long duration, intense southward IMFs. As mentioned in the introduction, Du et al. (2008) studied about a storm main phase during northward IMFs. Considering that the direction of IMF is important for magnetic storms, it would be worth investigating the interplanetary conditions of IMF By and Bz signs on the development of a storm. For this purpose, the direction change of IMF By and Bz components during the main phase of magnetic storm is examined.

    Fig. 1 shows temporal variations of (a) By, (b) Bz component of IMF, (c) solar wind temperature, (d) solar wind density, (e) solar wind speed, (f) beta, and (g) Dst index during 20-25 November 2003. The horizontal axis denotes universal time and the two vertical dashed lines represent the beginning and the maximum epoch of the storm. The dotted lines in Figs. 1a and b represent the zero values of IMF By and Bz, respectively. As shown in Figs. 1a and b, the IMF By component changes from a positive value to a negative value, and IMF Bz component remains negative during almost all of the main phase. In Fig. 1c, the solar wind temperature is somewhat high in the early stage of the main phase, but immediately becomes low. Values of solar wind density and speed in Figs. 1d and e are higher during the main phase than during the recovery phase. In contrast, the beta of plasma, symbolized by β, is the ratio of the plasma pressure to the magnetic pressure, which is lower during the main phase than during the recovery phase. Considering these conditions, Fig. 1 is an example of magnetic storms driven from MC. As can be seen from in Fig. 1, the Dst index decreased rapidly, and the minimum value of -422 nT at 2100 UT was recorded in November 2003.

    It is believed that the Bz component of IMF is more important than By on the development of magnetic storms. For this reason, Bz component was inserted in the first column of Table 2. Table 2 shows the classification of interplanetary conditions in terms of IMF Bz and By signs. F in Table 2 refers to the fluctuation of magnetic field, and the sign of '+ → -' or '- → +' stands for one polarity change. Magnetic storms are classified into 12 types according to IMF signs. No magnetic storm occurred during the positive IMF Bz period in this study. For this reason, those cases are not included in type classifications. As shown in Table 2, when IMF Bz remains steadily negative during the main phase of magnetic storms, they are classified as T1, T2, T3, and T4. A storm is included as one of these types if the percentage of positive values is below 10%. If the polarity of IMF Bz changes once, it belongs to T5~T8. Accordingly, if IMF Bz fluctuates, that type is one of T9~T12. For example, if the IMF condition is Bz < 0 and By > 0, it belongs to T1. When IMF Bz is fluctuating and By<0 during the main phase, type is T10. As can be seen from Table 2, almost half of CME-driven storms examined in this study, 44.2% (61 storms), occurred when IMF Bz remained steadily negative. Fig. 2 shows the

    histograms for 12 types during the main phase of CME-driven storms with the three storm categories. Since the number of storms in each category is not the same, the vertical axis is expressed in terms of a percentage value. For the CME-driven storms, it is found that the average number of storms for negative sign of IMF Bz (T1~T4) is high, at 56.4%, 53.0%, and 63.7% in each storm category,

    [Table 2.] Classification of interplanetary conditions in terms of IMF Bz and By signs.

    label

    Classification of interplanetary conditions in terms of IMF Bz and By signs.

    respectively. For the moderate storms, T2 shows the highest number of storms, while T3 has the highest percentage value for intense and severe storms. In particular, the sum of T3 and T4 for the severe storms is more than half in CME type storms. It is noted that IMF Bz seems to play a more important role in developing magnetic storms than By on the basis of Table 2 and Fig. 2. However, fluctuation of By component also seems to be important for big storms.

    On the other hand, the results for CIR-driven storms are different from those for the CME storm type. As shown in Fig. 3, the percentage of storms for T1~T4 is found to be low, at 29.2%, compared to that of T9~T12, which is 45.8% for moderate storms. When the condition is Bz < 0, the number of storms for CIR-driven type is approximately 52% lower than that of CME type for moderate storms. It was impossible to analyze storms this severe because there was no this category in CIR-driven type. For CIR-driven storms, a wide range of numbers was recorded for moderate storms, while there is no distribution in several types for intense storms.

       3.2 Correlations between the Time-Integrated Dst Index and Solar Wind Parameters during the Main Phase

    Echer et al. (2008) investigated the dependence of peak Dst on interplanetary parameters. This study also investigated the relationship between the magnetic storms and solar wind parameters. Taking into account the fact that the values of several solar wind parameters are temporary phenomena while a magnetic storm is continuous, it may be reasonable to compare time-integrated parameters rather than individual parameters. Thus, the difference between the present study and others is that the correlation used the time-integrated values is examined. In addition, I analyzed magnetic storms after dividing them into CME-driven storm and CIR-driven storms.

    Fig. 4 shows the correlations between the time-integrated Dst index and various solar wind parameters for the CME-driven storms: (a) time-integrated IMF Bz, (b) time-integrated temperature,(c) time-integrated solar wind density, and (d) time-integrated solar wind speed. 'r ' in each figure stands for correlation coefficient, and each filled circle represents one magnetic storm. As shown in Fig. 4a, the correlation coefficient is highest, at 0.77, between the time-integrated IMF Bz component and time-integrated Dst index during the main phase of magnetic storms. These indicate that IMF Bz is closely associated with the development of magnetic storms during

    the main phase for the CME-driven storms. On the other hand, no direct relation seems to exist between time-integrated temperature and time-integrated Dst index during the main phase. Fig. 5 shows the correlations for the CIR-driven storms. The vertical and horizontal axes of (a)~(d) indicate the same as Fig. 4. The highest correlation is found, with the correlation coefficient being 0.93, between time-integrated Dst index and time-integrated solar wind speed. In Fig. 5a, the correlation coefficient is lower than the coefficient of CME type.

    Table 3 shows the average values between the properties of CME-driven storms and CIR-driven storms. As can be seen from Table 3, for the CME-driven storms, |Dstmin| for the CME-driven storms is 34% higher than that of

    the CIR-driven storms, and |Bzmin| is also higher than for the CIR-driven type. Duration and time lag for the main phase is longer for the CIR-driven storms than for the

    [Table 3.] The average values between the properties of CME-driven storms and CIR-driven storms.

    label

    The average values between the properties of CME-driven storms and CIR-driven storms.

    CME-driven storms. Also, temperature for the CIR-driven storms is higher than that of the CME-driven storms. Average beta (β) is low, at 0.79, for CME-driven storms while beta is above 1.00 for the CIR-driven storms. These results correspond to the results obtained by Borovsky & Denton (2006). According to Borovsky & Denton (2006), plasma density is described as more superdense for CME-driven storms than for CME-driven storms. However, plasma density for the CME storms is similar to the CIR-driven type in the present work. This may be due to the selection criteria and periods used.

       3.3 Time Lag with Storm Size for CME-Driven Storms and CIR-Driven Storms

    Fig. 6 shows the relationship between Dstmin and Bzmin during the main phase for the CME-driven storms and CIR-driven storms. It would be meaningful to point out that the slopes of the linear regression lines for the two

    cases are different. The relationships between Dstmin and Bzmin during the main phase are expressed as follows:

    image
    image

    Eqs. (1) and (2) are for the CME-driven storms and CIR-driven storms. Compared with Figs. 6a and b, the slope in Fig. 6a is about two times higher than that of Fig. 6b. When Bzmin value is -20 nT, the corresponding Dstmin values are -155.4 nT and -115.2 nT for the CME and CIR storms, respectively. This means that IMF Bz for the CME-driven storms is more efficient than for the CIR-driven storms on the development of magnetic storms. The correlation coefficient for the CME-driven storms is significantly high, at 0.90, compared to for the CIR-driven storms.

    Fig. 7 shows the comparison of time lag for the CME and CIR-driven storms with storm size. For the CME-driven storms, time lag of 26% moderate storms is one hour, whereas time lag of 33% moderate storms is two hour for the CIR-driven storms. The pattern of intense storms is similar to that of moderate storms for the CME-driven storms, but the pattern of time lag for the severe storms is different. For CIR-driven storms, time lag of moderate storms is concentrated on two hours whereas that of intense storms is longer.

    4. SUMMARY AND DISCUSSION

    In this research, the manner in which the direction change of IMF By and Bz components (in GSM coordinate) during the main phase is related with the development of storms was examined. These were classified into 12 types according to the direction change of IMF By, Bz during the storm periods. Time-integrated solar wind parameters are compared to the time-integrated Dst index during the main phase for CME and CIR-driven storms, respectively. The differences between CME-driven storms and CIR-driven storms were also discussed.

    Several interesting results were summarized, and their implications briefly discussed. For the CME-driven storms, under steady conditions of Bz < 0, the number of storms represents more than half of storms. The percentage of storms was 56.4%, 53.0%, and 63.7% by storm category, respectively. One can note that the direction of Bz plays a more important role in the development of magnetic storms than By. However, for the CIR-driven storms, the number of moderate storms was only 29.2%, while the number of intense storms was over half (60.0%) under the Bz<0 condition. For CIR-driven storms, a wide range of numbers was recorded for the moderate storms, but there is no distribution in several types for the intense storms. Also, the fluctuations of IMF By and Bz components seem to be associated with the development of magnetic storms.

    The correlation coefficient is highest, at 0.77, between the time-integrated IMF Bz component and time-integrated Dst index for CME-driven storms, whereas the highest correlation is found, with a correlation coefficient of 0.93, between time-integrated Dst index and time-integrated solar wind speed for CIR-driven storms.

    The correlation coefficient between Dstmin and Bzmin is higher for the CME-driven storms than for the CIR-driven storms. Considering their slopes, IMF Bz of the CME-driven storms is more efficient on developing magnetic storms. Average time lag for the CIR-driven storms is longer than for CME-driven storms, and showed different aspects as shown in Fig. 7. As mentioned previously, Borovsky & Denton (2006) compared the properties between CME-driven storms and CIR-driven storms. In the present study, the results are approximately in accordance with Borovsky & Denton (2006), except for the solar wind density. This may be due to the period of data used and the criteria for storm selection.

참고문헌
  • 1. Borovsky JE, Denton MH 2006 Differences between CME-driven storms and CIR-driven storms [JGR] Vol.111 P.11447-11463 google cross ref
  • 2. Burlaga L, Sittler E, Mariani F, Schwenn R 1981 Magnetic loop behind an interplanetary shock: Voyager Helios and IMP 8 observations [JGR] Vol.86 P.6673-6684 google cross ref
  • 3. Burton RK, McPherron RL, Russell CT 1975 An empirical relationship between interplanetary conditions and Dst [JGR] Vol.80 P.4204-4214 google cross ref
  • 4. Choi C-R, Kim K-C, Lee D-Y, Kim J-H, Lee E 2008 Statistical characteristics of solar wind dynamic pressure enhancements during geomagnetic storms [JASS] Vol.25 P.113-128 google cross ref
  • 5. Du AM, Tsurutani BT, Sun W 2008 Anomalous geomagnetic storm of 21-22 January 2005: a storm main phase during northward IMFs [JGR] Vol.113 P.13284-13292 google cross ref
  • 6. Echer E, Alves MV, Gonzalez WD 2005 A statistical study of magnetic cloud parameters and geoeffectiveness [JASTP] Vol.67 P.839-852 google cross ref
  • 7. Echer E, Gonzalez WD, Tsurutani BT, Gonzalez ALC 2008 Interplanetary conditions causing geomagnetic storms (Dst ≤ -100 nT) during solar cycle 23 (1996-2006) [JGR] Vol.113 P.12744-12759 google cross ref
  • 8. Lopez RE, Wiltberger M, Hernandez S, Lyon JG 2004 Solar wind density control of energy transfer to the magnetosphere [GeoRL] Vol.31 P.8804-8807 google cross ref
  • 9. Oh SY, Yi Y, Kim YH 2007 Solar cycle variation of the interplanetary forward shock drivers observed at 1 AU [SoPh] Vol.245 P.391-410 google cross ref
  • 10. Tsurutani BT, Gonzalez WD, Tang F, Akasofu SI, Smith EJ 1988 Origin of interplanetary southward magnetic fields responsible for major magnetic storms near solar maximum (1978?1979) [JGR] Vol.93 P.8519-8531 google cross ref
  • 11. Weigel RS 2010 Solar wind density influence on geomagnetic storm intensity [JGR] Vol.115 P.15062-15072 google cross ref
  • 12. Wu C-C, Lepping RP 2006 Solar cycle effect on geomagnetic storms caused by interplanetary magnetic clouds [AnGeo] Vol.24 P.3383-3389 google cross ref
이미지 / 테이블
  • [ Table 1. ]  List of magnetic storms used in this study.
    List of magnetic storms used in this study.
  • [ Fig. 1. ]  Temporal variations of (a) By (b) Bz component of IMF (c) solar wind temperature (d) solar wind density (e) solar wind speed (f) beta and (g) Dst index during 20-25 November 2003. The two vertical dashed lines represent the beginning and the maximum epoch of the storm. The Dst index decreased rapidly and recorded the minimum value -422 nT at 2100 UT on November 2003. IMF: interplanetary magnetic field Dst: disturbance storm time.
    Temporal variations of (a) By (b) Bz component of IMF (c) solar wind temperature (d) solar wind density (e) solar wind speed (f) beta and (g) Dst index during 20-25 November 2003. The two vertical dashed lines represent the beginning and the maximum epoch of the storm. The Dst index decreased rapidly and recorded the minimum value -422 nT at 2100 UT on November 2003. IMF: interplanetary magnetic field Dst: disturbance storm time.
  • [ Table 2. ]  Classification of interplanetary conditions in terms of IMF Bz and By signs.
    Classification of interplanetary conditions in terms of IMF Bz and By signs.
  • [ Fig. 2. ]  Histograms for 12 types for CME-driven storms with the three storm categories. Since the number of storms for each category is not the same the vertical axis is expressed in terms of percentage values. CME: coronal mass ejection.
    Histograms for 12 types for CME-driven storms with the three storm categories. Since the number of storms for each category is not the same the vertical axis is expressed in terms of percentage values. CME: coronal mass ejection.
  • [ Fig. 3. ]  Same as Fig. 2 but for CIR-driven storms. CIR: co-rotating interaction region.
    Same as Fig. 2 but for CIR-driven storms. CIR: co-rotating interaction region.
  • [ Table 3. ]  The average values between the properties of CME-driven storms and CIR-driven storms.
    The average values between the properties of CME-driven storms and CIR-driven storms.
  • [ Fig. 4. ]  The correlations between the time-integrated Dst index and various solar wind parameters for the CME-driven storms: (a) time-integrated IMF Bz (b) time-integrated temperature (c) time-integrated solar wind density and (d) time-integrated solar wind speed. 'r ' in each figure stands for correlation coefficient. Each filled circle represents one magnetic storm. Dst: disturbance storm time CME: coronal mass ejection IMF: interplanetary magnetic field.
    The correlations between the time-integrated Dst index and various solar wind parameters for the CME-driven storms: (a) time-integrated IMF Bz (b) time-integrated temperature (c) time-integrated solar wind density and (d) time-integrated solar wind speed. 'r ' in each figure stands for correlation coefficient. Each filled circle represents one magnetic storm. Dst: disturbance storm time CME: coronal mass ejection IMF: interplanetary magnetic field.
  • [ Fig. 5. ]  Same as Fig. 4 but for CIR-driven storms. CIR: co-rotating interaction region.
    Same as Fig. 4 but for CIR-driven storms. CIR: co-rotating interaction region.
  • [ Fig. 6. ]  The relationship between Dstmin and Bzmin during the main phase for the CME and CIR-driven storms. CME: coronal mass ejection CIR: co-rotating interaction region.
    The relationship between Dstmin and Bzmin during the main phase for the CME and CIR-driven storms. CME: coronal mass ejection CIR: co-rotating interaction region.
  • [ Fig. 7. ]  Comparison of time lag with storm size for the CME and CIR-driven storms. CME: coronal mass ejection CIR: co-rotating interaction region.
    Comparison of time lag with storm size for the CME and CIR-driven storms. CME: coronal mass ejection CIR: co-rotating interaction region.
(우)06579 서울시 서초구 반포대로 201(반포동)
Tel. 02-537-6389 | Fax. 02-590-0571 | 문의 : oak2014@korea.kr
Copyright(c) National Library of Korea. All rights reserved.