Statistical analysis of the mid-latitude trough position during different categories of magnetic storms and different storm intensities

Abstract

The ionospheric mid-latitude trough minimum position as a function of geomagnetic storm time is identified and investigated statistically in terms of the category and the intensity of storms. The data of ion density derived from DMSP and DEMETER satellites were used to extract the trough position. The variations of mid-latitude trough in 41 moderate magnetic storms and 88 intense magnetic storms in the 23rd solar cycle were studied. The results show that the trough moves toward the equator as Dst index decreases and toward the pole as Dst index increases. Compared with the ICME, MC and CIR storms, in sheath storms the trough shifts to lower latitude at the end of the main phase, although the average storm intensity is weak. During the storm recovery phase, the rapid recovery of the trough position can be seen at the start of the recovery phase for moderate CIR storms. We also calculated the correlation between the minimum latitude of the trough position and the storm magnitude as well as other related main phase parameters during all storms. We found that the minimum latitude of the trough position exhibits a strong correlation with the storm magnitude during magnetic storms. However, the correlation coefficients between the trough position and other related main phase parameters are very low.

Introduction

The mid-latitude trough is a typical feature of the nightside F region and lies at the interface between the mid-latitude and auroral ionospheres. Throughout the current paper, ‘trough’ refers specifically to the mid-latitude trough. The trough has been studied for several decades (Muldrew 1965; Horvath and Essex 2003; Nilsson et al. 2005; Pryse et al. 2006; Middleton et al. 2008; He et al. 2011; Lee et al. 2011; Ishida et al. 2014). Most studies have displayed the quiet-time features of the trough (e.g., Rodger et al. 1992; Scali 1992; Rodger 2008). Some studies investigate the dynamics of the trough during substorms (e.g., Rodger et al. 1986; Rodger and Dudeney 1987; Zou et al. 2011). Only a few studies focus on the trough dynamics during magnetic storms. These studies mainly studied the relationship between trough position and Kp index and Dst index for all phases of storms (e.g., Deminov et al. 1995a, b, 1996; Karpachev et al. 1996). However, the trough dynamics have not been investigated in the catalog of different magnetic storms.

In the near-Earth space, magnetic storms are predominantly driven by two different types of solar disturbances: coronal mass ejections (CMEs) and corotating interaction regions (CIRs). A CIR is formed at the interface between the high-speed solar wind and preceding low-speed solar wind. The high-speed solar wind originates from the regions of open magnetic field in the solar corona (Krieger et al. 1973). When a single coronal hole persists for many solar rotation periods, CIRs may recur with a ~27-day period (or more frequently if more than one coronal hole is presented on the sun) (Tsurutani et al. 2006). CMEs with vast explosions of plasma and embedded magnetic flux erupt from magnetic fields associated with sunspots (Webb and Howard 2012). Interplanetary coronal mass ejections (ICMEs) are generally believed to be the counterparts of CMEs at the sun (Richardson and Cane 2010). ICMEs occur predominantly during solar maximum and generally cause the largest excursions in the Dst and Kp indices (Richardson and Cane 2012). ICME storms are further classified as either magnetic clouds (MC) storms or sheath storms by the geoeffective structure. Magnetic clouds are a subset of CME ejecta. Sheath regions are ahead of the CME ejecta. Magnetic clouds are identified in the solar wind by a strong magnetic field with a large rotation of the magnetic field direction (Burlaga et al. 1981; Klein and Burlaga 1982). In the sheath regions, the solar wind plasma is compressed and heated.

Several studies have contrasted the effects of ICMEs and CIRs (Wilken et al. 1999; Dmitriev et al. 2005). Huttunen et al. (2002) found that sheath regions drove stronger auroral activity, while magnetic clouds drove stronger ring current enhancement. Borovsky and Denton (2006) showed 21 differences between ICME and CIR storms in a number of regions, including the solar wind, the magnetosphere and the auroral zone. Denton et al. (2006) compared average plasma sheet parameters for the two types of events. They found that CIR storms had a more important modulation in the plasma sheet temperature, while ICME storms had a more important modulation in the plasma sheet density. Katus et al. (2015b) showed that the responses of the plasmapause location were also very different.

The aim of the present study was to systematically examine the general features of the trough position by undertaking a superposed epoch analysis. We focus on the differences seen in different categories of magnetic storms and different storm intensities. The paper begins with reviewing satellite data, identification of storms and method of superposed epoch analysis in “Data and Methodology” section. In “Results” section, we present the results for different categories of magnetic storms and different storm intensities. Fourth and fifth sections present the “Discussion” and “Conclusions” of the paper.

Data and methodology

Observations

We examined nocturnal latitudinal profiles of the total ion density in 1996–2005 measured by the Special Sensor-Ions, Electrons, Scintillation (SSIES) instrument package onboard the Defense Meteorological Satellite Program’s (DMSP) F12, F14 and F15 satellites. The DMSP satellites are in Sun-synchronous circular polar orbits at about 850 km altitude. The orbits of these DMSP satellites are on the 2100 LT meridian when in the nightside sector. In addition, we also examined nocturnal latitudinal profiles of the total ion density in 2005–2009 measured by the Plasma Analysis Instrument (IAP) onboard the DEMETER satellite. The French microsatellite DEMETER was launched in June 2004, and it was operated until December 2010. The spacecraft had a nearly Sun-synchronous circular orbit (LT ≈ 10.30 and ≈ 22.30) at an altitude of about 700 km.

Selection of trough-like phenomena

The trough minimum position can be found by analyzing latitudinal profiles of ion density between 40° and 70° magnetic latitude during each storm. The latitudinal range between 40° and 70° magnetic latitude has been demonstrated to be suitable by Yang et al. (2015). Deminov et al. (1995a) demonstrated that the trough position primarily depends on the actual magnetic activity level. So, the variation in the trough position during each storm is inspected visually. If the variation in the trough position is not consistent with the variation in the Dst index, corresponding trough signatures are removed. When the DMSP and DEMETER satellites are placed on the 2100 and 2230 LT meridians in the nightside sector, then the local time of the trough position is in the pre-midnight region. Therefore, the local time distribution of the trough position is small and the local time effect of the trough position can be neglected.

Identification of storms

Zhang et al. (2007a, b) classified 88 intense storms (Dstmin ≤ −100 nT) in the 23rd solar cycle (1996–2005) according to their solar wind driver. These 88 intense storms include 11 CIR storms and 77 ICME storms, of which 33 storms are driven by MC and 19 storms by sheath.

In order to further examine the response of trough position to storm intensity, we applied the following selection criteria. The moderate ICME storms (−100 nT ≤ Dstmin ≤ −50 nT) studied in this study are selected on the basis of an initial ICME list (1996–2009) of Richardson and Cane (2010). The moderate CIR storms (−100 nT ≤ Dstmin ≤ −50 nT) studied in this study are selected on the basis of an initial CIR list (1995–2004) of Jian et al. (2006). The CIR list (including 2005–2009) is later updated at http://www.srl.caltech.edu/ACE/ASC/DATA/level3/index.html.

The Dst index is then examined during each storm. Selection of moderate storms began with a broad search for all storms with the peak Dst between −50 and −100 nT. According to the selection criteria provided by Denton et al. (2006), storms are rejected when either (1) one or more events occur within the 72 h after the initial event, or (2) clear minimum storm Dst is not found during each event occurred. These criteria provide a list of 16 CIR and 25 ICME moderate storms. It should be noted that these moderate ICME storms contain MC storms due to the nature of the initial ICME list. Therefore, the storms studied in “Superposed epoch analysis results for CIR storms and ICME storms” and “Superposed epoch analysis results for moderate storms and intense storms” sections contain intense MC storms.

Figure 1 shows the distribution of ICME and CIR storms as well as the corresponding minimum Dst index for each event. It is clear that the majority of intense storms are ICME storms, and the moderate storms show a significant proportion of CIR storms. During solar maximum, ICME storms are typically most numerous; in contrast, CIR storms are sparse.

Fig. 1

Distribution of ICME storms and CIR storms used in this study

Superposed epoch analysis

In this study, a superposed epoch analysis similar in method to other recent studies (Katus et al. 2013, 2015a, b) is employed. This method uses multiple markers in each storm to normalize the time line. First, the time of the storm peak Dst is defined as the end of the main phase. The end of the recovery phase is then defined as the maximum Dst in the following 96 h. The start of the main phase is selected as the maximum Dst in the 24 h prior to the storm peak. In this study, the identification of a storm sudden commencement is not required. The start of the initial phase is defined as the largest rate of Dst increase in the 8 h prior to the start of the main phase. Finally, the 6-h period before the start of the initial phase is concatenated to present information concerning the initial state of the magnetosphere. The epoch markers are then used to calculate the duration of each storm phase. These values are then utilized to normalize the epoch time line. Each storm phase is essentially stretched or compressed to the average duration of each phase using linear interpolation. The superposition used the 1-h time-resolution bins in the normalized time period. The process is repeated for trough positions.

Correlation coefficient and P value

The correlation coefficients and P values are calculated to determine the relationship between the minimum latitude of the trough position and the storm magnitude as well as other related main phase parameters during all storms. It is known that the confidence in a relationship is formally determined not just by the correlation coefficient but also by the number of pairs in the data. The number of pairs determines the ‘significance’ of the statistical analyses. The P value is widely used to determine statistical significance in statistical hypothesis testing, specifically in null hypothesis significance testing. The P value is defined as the probability of getting a correlation as large as the observed value by random chance, when the null hypothesis is true. A small P value (typically <0.05) is taken as strong evidence against the null hypothesis, then the null hypothesis can be rejected, and it can be concluded that the correlation R is significant. P values are usually found on a reference table by first calculating a Chi-square value.

Results

Storm example

The storm event on May 8, 2005, was associated with the passage of a CIR. Figure 2 describes time series of the Dst index and the Kp index, as well as the trough position derived from different satellites for this storm. The red vertical lines show the start or end of each phase. The red circles represent variations of trough position with longitude between 0 and 200. The black circles represent variations of trough position with longitude between 200 and 360. There are little black circles in Fig. 2d, which is because DMSP ion density data over this longitude range don’t cover the trough region in the Southern Hemisphere. It is clear that the trough position depends on the current magnetic activity level. The trough moves toward the equator as Dst index decreases and toward the pole as Dst index increases. The behavior of the trough position mirrors that of the Dst index. However, it should be pointed out that the trough position reaches its minimum value in a few to several hours prior to the minimum in Dst index, and the changes in Kp index are ahead of variations in trough position (see Fig. 2e). The results are consistent with the findings of Deminov et al. (1995a). During the storm recovery phase, the trough position recovers to pre-event values in about 20 h (see Fig. 2c). In general, the recovery of plasmapause has a larger timescale (Rodger. 2008). Therefore, the relationship between the trough and plasmapause breaks down during the storm recovery phase. It can be seen from Fig. 2c that the variations in the Dst index during the recovery phase of storm can result in large variations in the trough position. In addition, it is clear that the longitudinal effect is weak and can be neglected during the magnetic storm. This result is consistent with the finding of Karpachev et al. (1996). So, in the later study, the longitudinal effect of trough is neglected.

Fig. 2

May 8, 2005, CIR storm. a Dst index, b Kp index, c Northern Hemisphere trough positions derived from DMSP-F15 satellite data, d Southern Hemisphere trough positions derived from DMSP-F14 satellite data, e Northern Hemisphere trough positions derived from DEMETER satellite data. The red vertical lines show the start or end of each phase

Variations in the trough position of hemispheres

To reveal the differences between the Northern and Southern Hemispheres in trough’s response to magnetic storm, Fig. 3 shows the result for all of the MC storms because of their largest number of events. Figure 3a displays the Dst index data density as a function of epoch time. The color of the plot presents the density of data points in each bin. Each bin is 3 h across the x axis and 5 nT across the y axis. The black and white lines indicate the mean and the median, respectively. Figure 3b shows the trough positions as a function of epoch time in the Northern Hemisphere. The pink, green and black circles indicate the data derived from DMSP-F12, DMSP-F14 and DMSP-F15, respectively. The black and blue lines indicate the mean and the median, respectively. For clarity, Fig. 3c only shows the average trough positions and the upper and lower quartiles in both hemispheres along the normalized epoch time line. The black and blue lines indicate the results of the Northern Hemisphere and Southern Hemisphere, respectively. The red vertical lines show the start or end of each phase. As shown in Fig. 3a, the average Dst index increases with time during the initial phase, then drops to approximately −130 nT at the storm peak and then slowly recovers to initial value. Figure 3b shows that the average trough position peaks at the peak of Dst index. It is striking that there is almost the same temporal evolution of the trough position in the Northern Hemisphere and Southern Hemisphere (see Fig. 3c). In quiet conditions, the difference between the trough positions of the Northern Hemisphere and that of the Southern Hemisphere is mainly due to the longitudinal difference. During magnetic storms, the longitudinal effect of the trough position can be neglected. Therefore, the trough position follows a similar variation throughout the storm sequence regardless of hemispheres. The result is also correct for the other three categories of magnetic storms. In the later text, we only present the results of the Northern Hemisphere.

Fig. 3

Superposed epoch plots of the Dst index (a), Northern Hemisphere trough positions (b) and the average trough positions in both hemispheres (c) for 33 MC storms. The color of the plot in a presents the density of data points in 5 nT by 3-h bins. The black and white lines in a indicate the mean and the median, respectively. The pink, green and black circles in b indicate the data derived from DMSP-F12, DMSP-F14 and DMSP-F15, respectively. The black and blue lines in b indicate the mean and the median, respectively. The black and blue lines in c indicate the results of the Northern Hemisphere and Southern Hemisphere, respectively. The short vertical lines in each figure indicate the range bounded by the upper and lower quartiles. The red vertical lines show the start or end of each phase

Superposed epoch analysis results for CIR storms and ICME storms

Figure 4 shows the average Dst index and trough position for CIR (a) and ICME (b) storms with the peak Dst between −50 and −150 nT, as a function of epoch time. The numbers of CIR and ICME storms are 27 and 60, respectively. It is clear that ICME storms cause larger excursions in the Dst index and trough position. During the storm main phase, changes in the trough position follow Dst index variations irrespective of storm drivers. During the storm recovery phase, the rapid recovery of the trough position can be found at the start of the CIR storm recovery phase, and the rate of recovery slows markedly with time. However, during the ICME storm recovery phase, changes in the trough position are consistent with those in the Dst index.

Fig. 4

Superposed epoch analysis of Dst index and trough position for CIR (a) and ICME (b) storms using the normalized time line. The solid blue lines indicate the average Dst index. The solid green lines indicate the average trough position. The short vertical lines indicate the range bounded by the upper and lower quartiles. The red vertical lines show the start or end of each phase

Superposed epoch analysis results for moderate storms and intense storms

Figure 5 shows the statistic results in the same format as Fig. 4, but for moderate (a) and intense (b) storms. The numbers of moderate and intense storms are 41 and 46, respectively. It is evident that intense storms have longer durations, both for the main phase and for the recovery phase. The trough position is closer to the equator at the end of the intense storm main phase. During the storm recovery phase, the trough position recovers rapidly at the start of the moderate storm recovery phase, but this feature disappears at the start of the intense storm recovery phase. In addition, this feature is also pronounced during the CIR storm recovery phase (see Fig. 4a). Therefore, this feature is associated with moderate CIR storms. This finding is verified through the analysis of the trough position for moderate CIR storms.

Fig. 5

Superposed epoch analysis of Dst index and trough position for moderate (a) and intense (b) storms using the normalized time line. The solid blue lines indicate the average Dst index. The solid green lines indicate the average trough position. The short vertical lines indicate the range bounded by the upper and lower quartiles. The red vertical lines show the start or end of each phase

Superposed epoch analysis results for different solar wind drivers

Figure 6 shows the results of the superposed epoch analysis for the four categories of magnetic storms with the peak Dst between −100 and −150 nT. Figure 2c presents the fact that the variations in the Dst index during the storm recovery phase can result in large variations in the trough position. Therefore, the storms with large variations in the Dst index during the recovery phase are removed. Table 1 shows several parameters for the four categories of magnetic storms. It is evident that these plots are quite similar to each other. The trough moves toward the equator as Dst index decreases and toward the pole as Dst index increases. For CIR, ICME and MC storms, the minimum latitude of the trough position decreases as the storm magnitude increases. The minimum latitude of the trough position for sheath storms is closer to the equator, although the average storm intensity is weak. In addition, it is striking that the average duration of the main phase is shortest for sheath storms. It seems that the minimum latitude of the trough position correlates with the duration of the main phase, or the rate of Dst decreases during the main phase, or the magnitude of Dst decreases during the main phase.

Fig. 6

Superposed epoch analysis of Dst index and trough position for the four categories of magnetic storms using the normalized time line. The solid blue lines indicate the average Dst index. The solid green lines indicate the average trough position. The short vertical lines indicate the range bounded by the upper and lower quartiles. The red vertical lines show the start or end of each phase

Table 1 Parameters characterizing the four categories of magnetic storms

Figure 7 shows the relationship between the minimum latitude of the trough position and the storm magnitude. The data in Fig. 7a, b are derived from DMSP and DEMETER satellites, respectively. The correlation coefficients and P values between the minimum latitude of the trough position and the Dst index at the start of main phase (Dstmax, second row), the magnitude of Dst decrease during the main phase (ΔDst, third row), the duration of main phase (ΔT, forth row), the rate of Dst decrease during the main phase (ΔDst/ΔT, fifth row) and the Dst peak value (Dstmin, sixth row) are summarized in Table 2. It is evident that there is a strong (correlation coefficients >0.8) and significant (P value <0.0001) correlation between the minimum latitude of the trough position and the Dstmin. Because ΔDst correlates with Dstmin, the correlation with P value <0.001 between the minimum latitude of the trough position and ΔDst are also statistically significant. However, the correlation between the minimum latitude of the trough position and Dstmax and ΔT is weak (correlation coefficient <0.3) and not significant (P value >0.1). So, the minimum latitude of the trough position does not correlate with other parameters associated with the main phase. The minimum latitude of the trough position is largely dependent on Dstmin for all storms. As the Dstmin decreases, the trough moves to lower latitude.

Fig. 7

Relationship between the minimum latitude of the trough position and the storm magnitude. The data in a and b are derived from DMSP and DEMETER satellites, respectively

Table 2 Correlation coefficients (cc) and P values for different satellite data

Discussion

In this paper, we specify the statistical characteristics of the trough position for different categories of magnetic storms and different storm intensities. During magnetic storms, the longitudinal effect in the trough position is weak and can be neglected. The result is consistent with the finding of Karpachev et al. (1996). In addition, the temporal evolution of the trough position in the Northern and Southern Hemispheres during magnetic storms is symmetric.

During the main phase of magnetic storms, the average trough position moves equatorward with decreasing Dst index. During this phase, an increase in the intensity of the ring current after the rotation of the vertical component of interplanetary magnetic field (IMF) Bz from north to south is accompanied by the displacement of the inner boundary of plasma sheet ions toward the Earth (Hamilton et al. 1988). Since this boundary is connected with the trough, the trough position follows the variation and moves equatorward (Deminov et al. 1995a). Compared with the ICME, MC and CIR storms, in sheath storms the trough shifts to lower latitude at the end of the main phase. This phenomenon is associated with the solar wind properties of sheath storms. Previous studies have addressed the similarities and differences of the dominant solar wind parameters for the four categories of magnetic storms with the peak Dst between −100 and −150 nT (Katus et al. 2015b). During the main phase of sheath storms, the greater inner magnetospheric convection associated with larger southward IMF Bz and solar wind velocity exists. The greater inner magnetospheric convection can result in the larger displacement of the inner boundary of plasma sheet and the trough position. In addition, it has also been demonstrated in our study that the minimum latitude of the trough position is largely dependent on the magnitude of the storms. The minimum latitude of the trough position does not correlate with other parameters associated with the main phase.

During the recovery phase of magnetic storms, the average trough position moves poleward with increasing Dst index. The timescale for recovery of the trough is of the order of hours (see Fig. 2c). It is striking that the rapid recovery of the trough position occurs at the start of the recovery phase for moderate CIR storms. For strong CIR storms and ICME storms, changes in the trough position follow Dst index variations. To understand the cause of the differences, this study examines the dominant solar wind parameters for moderate and intense CIR storms. Figure 8 shows superposed epoch solar wind parameters (solar wind electric field, solar wind dynamic pressure and IMF Bz) for moderate (left panels) and intense (right panels) CIR storms. These hourly averaged data are taken from the high-time-resolution (1 min) OMNI database (http://omniweb.gsfc.nasa.gov/). The red vertical lines show the start or end of each phase. It is striking that these parameters continue to vary at the start of the recovery phase for moderate storms, but remain constant for intense storms. This phenomenon is not found at the start of the ICME storm recovery phase. These differences of solar wind parameters may lead to the rapid recovery of the trough position for moderate CIR storms.

Fig. 8

Superposed epoch plots of solar wind parameters (solar wind electric field, solar wind dynamic pressure and IMF Bz) for moderate (left panels) and intense (right panels) CIR storms. The black and white lines indicate the mean and the median, respectively. The red vertical lines show the start or end of each phase

Although theoretically CIRs are essentially sheaths, their solar wind properties are very different. These differences may lead to different behaviors of the trough. Therefore, future work should be carried out to examine the relationship between the various solar wind parameters and the trough.

Conclusions

This study shows evidence that the trough dynamics investigated are different for different categories of magnetic storms and different storm intensities. The properties of trough position during magnetic storms may be summarized as follows.

• (1) The temporal evolution of the trough position in the Northern and Southern Hemispheres during magnetic storms is symmetric.

• (2) During the main phase of magnetic storms, the average trough position moves equatorward with decreasing Dst index. Compared with the ICME, MC and CIR storms, in sheath storms the trough shifts to lower latitude at the end of the main phase, although the average storm intensity is weak.

• (3) During the recovery phase of magnetic storms, the average trough position moves poleward with increasing Dst index. The rapid recovery of the trough position can be seen at the start of the recovery phase for moderate CIR storms.

• (4) As the magnitude of magnetic storm increases, the trough moves to lower latitude.

As a conclusion, during the sheath and CIR storms, the trough’s behavior is different from that during the ICME and MC storms. During the ICME and MC storms, changes in the trough position are consistent with those in the Dst index. Compared with the ICME, MC and CIR storms, in sheath storms the trough shifts to lower latitude at the end of the main phase, although the average storm intensity is weak. The rapid recovery of the trough position appears at the start of the recovery phase in the moderate CIR storms but not in ICME and MC storms. This information is likely to assist in validating and testing ionospheric models where the trough structure is included.