Storm-time Substorms Research Articles

MAGNETOSPHERIC SUBSTORMS AND RELATIVISTIC ELECTRONS

The most recent findings on the dynamics of the outer radiation belt (ORB) and the physics of magnetospheric substorms are examined. Specifically, we investigate the relationship between storm time substorms and the energetic electron population that forms the ORB. Traditionally, storm time substorms have been considered as the primary source of energetic electrons, which are further accelerated during storms to contribute to the formation of the ORB. However, several observations have demonstrated that large magnetospheric substorms can generate high-energy electrons even in the absence of magnetic storms. Substorms introduce dispersionless injections of energetic electrons deep into the magnetosphere from the geosynchronous orbit during storm times. The injected electrons undergo additional acceleration via the betatron mechanism during the storm recovery phase, thus increasing the ORB population. To gain a better understanding of this process, it is crucial to study plasma sheet turbulence, substorm onset processes, and the brightening of auroral arcs. By analyzing the aforementioned findings, this study aims to highlight the need for reanalyzing of the role of auroral processes in the formation of the ORB.

Concurrent Empirical Magnetic Reconstruction of Storm and Substorm Spatial Scales Using Data Mining and Virtual Spacecraft

Data mining (DM) has ushered in a new era of empirical magnetic reconstructions of the magnetosphere via application of the k-nearest neighbors (kNN) method. In this approach, the combined magnetosphere storm-substorm state is characterized by the Sym-H and AL indices, their time derivatives, and the solar wind electric field vBzIMF. However, using the DM reconstructions to account for the substorm contributions to the ring current as well as describing storm-time substorms remains a problem. The inner region r ≤ 12RE, where the ring current develops, has a much higher density of data than the tail region 12RE ≤ r ≤ 22RE, where substorms operate. This results in two models inconsistent in their scales dictated by the corresponding data densities. The inner model reconstructs storm time dynamics, including the formation of the westward and eastward ring current and pressure distributions. The outer model captures substorm features, including the thinning and rapid dipolarization of the tail sheet during the growth and expansion phases, respectively. However, the substorm model is insufficient to reconstruct the eastward ring current while the storm model cannot fully reproduce substorm effects because it overfits in the tail region. This issue is addressed by constructing a hybrid model which is fit using virtual magnetic field observations generated by sampling the other two models. The resulting merged resolution model concurrently captures the spatial scales associated with both storms in the inner region and substorms in the near-tail region. Hence it is particularly useful for investigation of the storm-substorm relationship, including storm-time substorms and the impact of individual substorm injections to the buildup of the storm-time ring current.

Spatially Quasi-periodic Finger-like Auroras during Substorms

Abstract North–south-aligned discrete auroras are common phenomena during substorms. A quasi-periodic form of such auroras distributed azimuthally in the nightside resembles a finger-like morphology. They were noticed in sawtooth-time substorms, but have not yet been paid much attention during other geomagnetic times. In this study, we investigate the longitudinal periodicity of those auroras for different types of substorms using IMAGE data from 2001 to 2003. Using an automated detection algorithm, we identified multiple north–south-aligned thin auroral structures in auroral images and found that their average center and length in latitude, and longitudinal wavelength (∼25°) are essentially the same. However, the occurrent rates for finger-like auroras are substantially different, 83%, 64%, and 18% for sawtooth-time substorms, storm-time substorms and isolated substorms, respectively. Therefore, statistical results indicate that finger-like auroras occur in favor of geomagnetically more active times. Combining our data analysis and previous simulations, we suggest that interchange instability is a possible mechanism of spatially quasi-periodic finger-like auroras when a substorm injection happens in a very wide local time. It is particularly important to point out that the sawtooth events not only have a temporal periodicity as defined in particle injections, but also have a spatial periodicity as observed in auroral morphology for the majority of them.

Penetration electric fields observed at middle and low latitudes during the 22 June 2015 geomagnetic storm

Prompt penetration electric fields were detected during the intense geomagnetic storm on 22 June, 2015 by the HF Doppler sounders at middle and low latitudes; Prague (Czech Republic), Iitate, Oarai, Onna (Japan), and Zhongli (Republic of China). The storm was initiated by the storm sudden commencement (SC) at 1833 UT, immediately followed by the main phase with the minimum SYM-H of − 139 nT at 2017 UT. The convection electric field was found to be up to 5.1 mV/m at low latitudes, which is enough to cause ionospheric disturbances such as the enhancement of the total electron content, equatorial ionization anomaly and so on. The electric fields of the SC and main phase are westward in the night (0330 MLT), while eastward in the evening (1930 MLT). The main phase electric field intensified the eastward equatorial electrojets (EEJ) at Huancayo in South America on the dayside (1330 MLT) and the westward EEJ on the nightside at Guam in the western Pacific (0330 MLT) and Tirunelveli in India (00 MLT). The eastward direction of the evening-time electric field agrees with a feature of the evening anomaly of the penetration electric field, while the intensity is suppressed because of the sunlit condition in the evening in the summer solstice. The sudden northward turning of the interplanetary magnetic field (IMF) caused a transition from the main phase to the recovery phase, when the EEJ turned into the counterelectrojet (CEJ). It was found that the HF Doppler sounders detected an impulsive eastward overshielding electric field at 2005 UT on the nightside during the northward IMF-driven CEJ with the duration of 17 min and intensity of 8.0 mV/m at Iitate (0500 MLT) and 11.8 mV/m at Prague (2130 MLT). The impulsive overshielding electric field was accompanied by the anti-sunward plasma flows at middle latitude (50°–60° GML) equatorward of the dusk sunward convection as observed with SuperDARN, which signifies intensification of the Region-2 field-aligned currents at the flow reversal. Substorm signatures are identified with intensified westward auroral electrojet and magnetic dipolarization as detected by the geosynchronous orbit satellite, ETS-VIII. Consequently, the stormtime substorm may provide strong overshielding electric fields to the middle latitude ionosphere on the nightside.

Energy Budget of High-speed Plasma Flows in the Terrestrial Magnetotail

Abstract Significant advances have been presented for the theoretical model and quantitative investigation of the energy input from the solar wind and its subsequent release into the ring current, the Joule heating, and particle precipitation energy flux that are closely related to high-speed flows in the plasma sheet of the terrestrial magnetotail during both isolated and storm-time substorms. Here we further determine the plasma flow vorticity/shear and heating observed in the terrestrial magnetotail by the Magnetospheric Multiscale (MMS) mission by deducing that other energy may be dissipated along with the above-mentioned three energy releases. The energy transported by the flow is also estimated here. Three-dimensional observations have shown new detailed information that can be compared with quasi-2D simulations. In this investigation, we calculate vorticity/shear directly from simultaneous observation of four MMS satellites. Our results generally indicate that plasma heating or temperature enhancements are related to both the flow vorticity/shear and current density, but more strongly with flow vorticity/shear. This research suggests that the ubiquitous energy input from the solar wind can be used to estimate plasma heating or temperature enhancements in the absence of any process or phenomenon relating to explosive energy release in planetary magnetospheres.

Swarm field-aligned currents during a severe magnetic storm of September 2017

Abstract. Swarm satellite observations are used to characterize the extreme behavior of large- and small-scale field-aligned currents (FACs) during the severe magnetic storm of September 2017. Evolutions of the current intensities and the equatorward displacement of FACs are analyzed while the satellites cross the pre-midnight, pre-noon, dusk and dawn sectors in both hemispheres. The equatorward boundaries of FACs mainly follow the dynamics of the ring current as monitored in terms of the SYM-H index. The minimum latitude of the FAC boundaries is limited to 50∘ magnetic latitude (MLat). The FAC densities are very variable and may increase dramatically, especially in the nightside ionosphere during the storm-time substorms. At the peak of substorms, the average FAC densities reach >3 µA m−2. The dawn–dusk asymmetry is manifested in the enhanced dusk-side R2 FACs in both hemispheres. In the 1 Hz data filamentary high-density structures are always observed. In the pre-noon sector, the bipolar structures (7.5 km width FACs of opposite polarities adjacent to each other) dominate, while at the other local times the upward and downward FACs tend to be latitudinally separated. The most intense small-scale FACs, up to ∼80 µA m−2, are observed just in the post-midnight sector. Simultaneous magnetic and plasma perturbations indicate that this structure is likely a current system of a mesoscale auroral arc.

Large‐Scale Ducting of Pc1 Pulsations Observed by Swarm Satellites and Multiple Ground Networks

AbstractLow Earth orbit satellites frequently encounter Pc1 pulsations, but most have been observed with limited latitudinal extent or short lifetime. In this study we analyze two large‐scale Pc1 pulsations (both latitudinally wide and long‐lasting) generated by ionospheric ducting effect using Swarm and ground magnetometers on 25 June and 3 September 2015. Swarm observed the 25 June pulsations on both dayside and nightside during the storm time substorm (a strong geomagnetic storm on 23 June with Dst = − 204 nT). We found the Pc1 pulsations were pervasive in both magnetic local time sectors of dayside and nightside for at least 2 hr. Another large Pc1 pulsation on 3 September was observed during a nonstorm substorm period. We conclude that (1) ionospheric ducting can transmit Pc1 waves to a wide range of L shells, (2) geomagnetic storm is not a prerequisite for such large‐scale ducting, and (3) wave intensity can abruptly decrease across sharp gradients in the ionospheric plasma density.

Geomagnetically Induced Currents Caused by Interplanetary Shocks With Different Impact Angles and Speeds

AbstractThe occurrence of geomagnetically induced currents (GICs) poses serious threats to modern technological infrastructure. Large GICs result from sharp variations of the geomagnetic field (dB/dt) caused by changes of large‐scale magnetospheric and ionospheric currents. Intense dB/dt perturbations are known to occur often in high‐latitude regions as a result of storm time substorms. Magnetospheric compressions usually caused by interplanetary shocks increase the magnetopause current leading to dB/dt perturbations more evident in midlatitude to low‐latitude regions, while they increase the equatorial electrojet current leading to dB/dt perturbations in dayside equatorial regions. We investigate the effects of shock impact angles and speeds on the subsequent dB/dt perturbations with a database of 547 shocks observed at the L1 point. By adopting the threshold of dB/dt = 100 nT/min, identified as a risk factor to power systems, we find that dB/dt generally surpasses this threshold when following impacts of high‐speed and nearly frontal shocks in dayside high‐latitude locations. The same trend occurs at lower latitudes and for all nightside events but with fewer high‐risk events. Particularly, we found nine events in equatorial locations with dB/dt > 100 nT/min. All events were caused by high‐speed and nearly frontal shock impacts and were observed by stations located around noon local time. These high‐risk perturbations were caused by sudden strong and symmetric magnetospheric compressions, more effectively intensifying the equatorial electrojet current, leading to sharp dB/dt perturbations. We suggest that these results may provide insights for GIC forecasting aiming at preventing degradation of power systems due to GICs.

Stormtime substorm onsets: occurrence and flow channel triggering

Bright auroral emissions during geomagnetic storms provide a good opportunity for testing the proposal that substorm onset is frequently triggered by plasma sheet flow bursts that are manifested in the ionosphere as auroral streamers. We have used the broad coverage of the ionospheric mapping of the plasma sheet offered by the high-resolution THEMIS all-sky-imagers (ASIs) and chose the main phases of 9 coronal mass ejection (CME) related and 9 high-speed stream (HSS)-related geomagnetic storms, and identified substorm auroral onsets defined as brightening followed by poleward expansion. We found a detectable streamer heading to near the substorm onset location for all 60 onsets that we identified and were observed well by the ASIs. This indicates that substorm onsets are very often triggered by the intrusion of plasma with lower entropy than the surrounding plasma to the onset region, with the caveat that the ASIs do not give a direct measure of the intruding plasma. The majority of the triggering streamers are “tilted streamers,” which extend eastward as their eastern tip tilts equatorward to near the substorm onset location. Fourteen of the 60 cases were identified as “Harang streamers,” where the streamer discernibly turns toward the west poleward of reaching to near the onset latitude, indicating flow around the Harang reversal. Using the ASI observations, we observed substantially less substorm onsets for CME storms than for HSS storms, a result in disagreement with a recent finding of approximately equal substorm occurrences. We suggest that this difference is a result of strong non-substorm streamers that give substorm-like signatures in ground magnetic field observations but are not substorms based on their auroral signature. Our results from CME storms with steady, strong southward IMF are not consistent with the ~ 2–4 h repetition of substorms that has been suggested for moderate to strong southward IMF conditions. Instead, our results indicate substantially lower substorm occurrence during such steady driving conditions. Our results also show the much more frequent occurrence of substorms during HSS period, which is likely due to the highly fluctuating IMF.

Simultaneous Remote Observations of Intense Reconnection Effects by DMSP and MMS Spacecraft During a Storm Time Substorm.

During a magnetic storm on 23 June 2015, several very intense substorms took place, with signatures observed by multiple spacecraft including DMSP and Magnetospheric Multiscale (MMS). At the time of interest, DMSP F18 crossed inbound through a poleward expanding auroral bulge boundary at 23.5 h magnetic local time (MLT), while MMS was located duskward of 22 h MLT during an inward crossing of the expanding plasma sheet boundary. The two spacecraft observed a consistent set of signatures as they simultaneously crossed the reconnection separatrix layer during this very intense reconnection event. These include (1) energy dispersion of the energetic ions and electrons traveling earthward, accompanied with high electron energies in the vicinity of the separatrix; (2) energy dispersion of polar rain electrons, with a high‐energy cutoff; and (3) intense inward convection of the magnetic field lines at the MMS location. The high temporal resolution measurements by MMS provide unprecedented observations of the outermost electron boundary layer. We discuss the relevance of the energy dispersion of the electrons, and their pitch angle distribution, to the spatial and temporal evolution of the boundary layer. The results indicate that the underlying magnetotail magnetic reconnection process was an intrinsically impulsive and the active X‐line was located relatively close to the Earth, approximately at 16–18 RE.

Contribution of storm time substorms to the prompt electric field disturbances in the equatorial ionosphere

AbstractThis study tries to bring out the fact that storm time substorms can compete and at times significantly contribute to the geomagnetically disturbed time prompt penetration electric field effects on low and equatorial latitudes. Observations of unusual equatorial plasma drift data from Jicamarca Unattended Long‐term Investigations of the Ionosphere and Atmosphere during two space weather events show that substorms can induce both eastward and westward penetration electric fields under steady southward interplanetary magnetic field (IMF Bz) conditions. During the first event on 2 January 2005, the enhancement of the daytime eastward electric field over Jicamarca due to substorm is found to be comparable with the Sq and interplanetary electric field (IEFy) generated electric fields combined. During the second event on 19 August 2006, the substorm is seen to weaken the daytime eastward field thereby inducing a westward field in spite of the absence of northward turning of IMF Bz (overshielding). The westward electric field perturbation in the absence of any overshielding events is observationally sparse and contrary to the earlier results. Further, the substorm‐induced field is found to be strong enough to compete or almost nullify the effects of storm time IEFy fields. This study also shows quantitatively that at times substorm contribution to the disturbed time prompt electric fields can be significant and thus should be taken into consideration in evaluating penetration events over low latitudes.

Contributions of substorm injections to SYM‐H depressions in the main phase of storms

AbstractSubstorm injections bring energetic particles to the inner magnetosphere. But the role of the injected population in building up the storm time ring current is not well understood. By surveying Los Alamos National Laboratory geosynchronous data during 34 storm main phases, we show evidence that at least some substorm injections can contribute to substorm‐time scale SYM‐H/Dst depressions in the main phase of storms. For event studies, we analyze two typical events in which the main‐phase SYM‐H index exhibited stepwise depressions that are correlated with particle flux enhancement due to injections and with AL index. A statistical study is performed based on 95 storm time injection events. The flux increases of the injected population (50–400 keV) are found proportional to the sharp SYM‐H depressions during the injection interval. By identifying dispersionless and dispersive injection signals, we estimate the azimuthal extent of the substorm injection. Statistical results show that the injection regions of these storm time substorms are characterized with an azimuthal extent larger than 06:00 magnetic local time. These results suggest that at least some substorm injections may mimic the large‐scale enhanced convection and contribute to sharp decreases of Dst in the storm main phase.

Geomagnetically induced currents around the world during the 17 March 2015 storm

AbstractGeomagnetically induced currents (GICs) represent a significant space weather issue for power grid and pipeline infrastructure, particularly during severe geomagnetic storms. In this study, magnetometer data collected from around the world are analyzed to investigate the GICs caused by the 2015 St. Patrick's Day storm. While significant GIC activity in the high‐latitude regions due to storm time substorm activity is shown for this event, enhanced GIC activity was also measured at two equatorial stations in the American and Southeast Asian sectors. This equatorial GIC activity is closely examined, and it is shown that it is present both during the arrival of the interplanetary shock at the storm sudden commencement (SSC) in Southeast Asia and during the main phase of the storm ∼10 h later in South America. The SSC caused magnetic field variations at the equator in Southeast Asia that were twice the magnitude of those observed only a few degrees to the north, strongly indicating that the equatorial electrojet (EEJ) played a significant role. The large equatorial magnetic field variations measured in South America are also examined, and the coincident solar wind data are used to investigate the causes of the sudden changes in the EEJ ∼10 h into the storm. From this analysis it is concluded that sudden magnetopause current increases due to increases in the solar wind dynamic pressure, and the sudden changes in the resultant magnetospheric and ionospheric current systems, are the primary drivers of equatorial GICs.

The proton and electron radiation belts at geosynchronous orbit: Statistics and behavior during high‐speed stream‐driven storms

AbstractThe outer proton radiation belt (OPRB) and outer electron radiation belt (OERB) at geosynchronous orbit are investigated using a reanalysis of the LANL CPA (Charged Particle Analyzer) 8‐satellite 2‐solar cycle energetic particle data set from 1976 to 1995. Statistics of the OPRB and the OERB are calculated, including local time and solar cycle trends. The number density of the OPRB is about 10 times higher than the OERB, but the 1 MeV proton flux is about 1000 times less than the 1 MeV electron flux because the proton energy spectrum is softer than the electron spectrum. Using a collection of 94 high‐speed stream‐driven storms in 1976–1995, the storm time evolutions of the OPRB and OERB are studied via superposed epoch analysis. The evolution of the OERB shows the familiar sequence (1) prestorm decay of density and flux, (2) early‐storm dropout of density and flux, (3) sudden recovery of density, and (4) steady storm time heating to high fluxes. The evolution of the OPRB shows a sudden enhancement of density and flux early in the storm. The absence of a proton dropout when there is an electron dropout is noted. The sudden recovery of the density of the OERB and the sudden density enhancement of the OPRB are both associated with the occurrence of a substorm during the early stage of the storm when the superdense plasma sheet produces a “strong stretching phase” of the storm. These storm time substorms are seen to inject electrons to 1 MeV and protons to beyond 1 MeV into geosynchronous orbit, directly producing a suddenly enhanced radiation belt population.

Transient, small-scale field-aligned currents in the plasma sheet boundary layer during storm time substorms.

We report on field‐aligned current observations by the four Magnetospheric Multiscale (MMS) spacecraft near the plasma sheet boundary layer (PSBL) during two major substorms on 23 June 2015. Small‐scale field‐aligned currents were found embedded in fluctuating PSBL flux tubes near the separatrix region. We resolve, for the first time, short‐lived earthward (downward) intense field‐aligned current sheets with thicknesses of a few tens of kilometers, which are well below the ion scale, on flux tubes moving equatorward/earthward during outward plasma sheet expansion. They coincide with upward field‐aligned electron beams with energies of a few hundred eV. These electrons are most likely due to acceleration associated with a reconnection jet or high‐energy ion beam‐produced disturbances. The observations highlight coupling of multiscale processes in PSBL as a consequence of magnetotail reconnection.

Oxygen impacts on dipolarization fronts and reconnection rate

AbstractSpacecraft observations near a magnetotail X line show that oxygen (O+) ions are minor species during nonstorm substorms, but they can become major species during some of the storm time substorms. Dipolarization fronts (DFs), which are characterized by a sharp increase northward magnetic field in the magnetotail, are commonly observed during magnetospheric substorms. In this study, we investigated the O+ effects on DFs and the reconnection rate during magnetotail reconnection. We used a 2.5‐D implicit particle‐in‐cell simulation in a 2‐D Harris current sheet in the presence of H+ and O+ ions. Simulation runs with equal number densities of O+ and H+ (O+ run) and with pure H+ ion species (H+ run) were performed. Comparing the two different runs, we found that both the reconnection rate and the DF speed in the O+ run were much less than those in the H+ run. By studying the force balance and plasma composition at the DF, we found that the outflow magnetic flux and DF propagation were encumbered by the current sheet O+ inertia, which reduced the DF speed and delayed the reconnection rate in the O+ run. We also found an ambipolar electric field in the O+ run due to the different inflow and outflow speeds of O+ and electrons in the O+ diffusion region. As a result, this ambipolar electric field induced O+ drag on the convective magnetic field in the O+ diffusion region. The small reconnection rate determined in the O+ run can be attributed to the current sheet inertia and the O+ drag on the convective magnetic flux.

Substorm related CNA near equatorward boundary of the auroral oval in relation to interplanetary conditions

Cosmic noise absorption (CNA) at high latitudes is a typical manifestation of enhanced precipitation of energetic charged particles during the course of a magnetospheric substorm. Present analysis demonstrates the energetic particles precipitate to the high latitude ionosphere during substorms, affecting upper and lower regions of the ionosphere simultaneously. Previous studies have reported that intense and short-lived CNA events associated with substorms are mostly observed in the midnight sector of the auroral oval. In the current study, we have examined such type of CNA events predominantly occurring during 0000–0600UT (2300–0500MLT) at an Indian Antarctic station Maitri (corrected geomagnetic (CGM) coordinates 62.59°S, 53.59°E), which is located at the equatorward edge of the auroral oval. Absorption events related to isolated substorm and storm-time substorms exhibit distinct features in terms of their intensity and extent in latitude and longitude. Our study suggests that the maximum intensity of CNAs depends on the interplanetary conditions, such as, the solar wind speed, southward component of IMF Bz, and duskward component of IEF Ey. Moreover, the role of duskward component of IEF Ey is more noteworthy than other interplanetary parameters.

First results from imaging riometer installed at Indian Antarctic station Maitri

Cosmic noise absorption (CNA) measurred by imaging riometer, is an excellent tool to passively study the high latitude D-region ionospheric conditions and dynamics. An imaging riometer has been installed at Indian Antarctic station Maitri (geographic 70.75°S, 11.75°E; corrected geomagnetic 63.11°S, 53.59°E) in February 2010. This is the first paper using the imaging riometer data from Maitri. The present paper introduces the details of this facility, including its instrumentation, related CNA theory and its applications. Sidereal shift of around 2 hours in the diurnal pattern validates the data obtained from the newly installed instrument. Moreover, the strength of cosmic noise signal on quiet days also varies with months. This is apparently due to solar ionization of D-region ionosphere causing enhanced electron density where collision frequency is already high. The main objective of installing the imaging riometer at Maitri is to study magneotspheric–ionospheric coupling during substorm processes. In the current study, we present two typical examples of disturbed time CNA associated with storm-time and non-storm time substorm. Results reveal that CNA is more pronounced during storm-time substorm as compared to non-storm time substorm. The level of CNA strongly depends upon the strengthening of convectional electric field and the duration of southward turning of interplanetary magnetic field before the substorm onset.

Magnetotail current contribution to the Dst Index Using the MT Index and the WINDMI model

In this paper, we have improved the capabilities of a low dimensional nonlinear dynamical model called WINDMI to determine the state of the global magnetosphere by employing the magnetotail (MT) index as a measurement constraint during large geomagnetic storms. The MT index is derived from particle precipitation measurements made by the NOAA/POES satellites. This index indicates the location of the nightside ion isotropic boundary, which is then used as a proxy for the strength of the magnetotail current in the magnetosphere. In Asikainen et al. (2010), the contribution of the tail current to the Dst index is estimated from an empirical relationship based on the MT index. Here the WINDMI model is used as a substitute to arrive at the tail current and ring current contribution to the Dst index, for comparison purposes. We run the WINDMI model on 7 large geomagnetic storms, while optimizing the model state variables against the Dst index, the MT index, and the AL index simultaneously. Our results show that the contribution from the geotail current produced by the WINDMI model and the MT index are strongly correlated, except during some periods when storm time substorms are observed. The inclusion of the MT index as an optimization constraint in turn increases our confidence that the ring current contribution to the Dst index calculated by the WINDMI model is correct during large geomagnetic storms.

Substorms under northward interplanetary magnetic field: Statistical study

It was recently noted that substorms can occur even under prolonged northward interplanetary magnetic field (IMF) conditions. Based on the substorm list obtained from the IMAGE spacecraft, we perform a statistical study on the features of substorms during northward IMF interval. The strength of the substorm is represented by the AL index decrease and the total intensity of the auroral bulge. Four main features have been found as follows: (1) Most substorms occur soon after a southward IMF, and intense substorms are more likely to occur for short duration of northward IMF period (2 ~ 5 h), whereas no intense substorms occur after prolonged northward IMF condition. (2) There is a positive correlation between the strength of the substorm and the two solar wind parameters (the IMF |By| and the solar wind dynamic pressure Pd). (3) The average strength of the substorms during the storm period is much larger than that of the substorms during the period without storm. Meanwhile, nearly all strong storm time substorms occurred either during the intense storm period, or during the late main phase or the early recovery phase of the storm. (4) About half of substorms, either an increase or a decrease in the solar wind dynamic pressure, are found within 30 min preceding each onset time. Such features indicate that the energy stored in the magnetotail during a previous southward IMF period is the main energy source for substorms under northward IMF condition, especially for intense substorms, and both the IMF |By| and the solar wind dynamic pressure play an important role in the energy accumulation during the northward IMF period.