What caused the rapid recovery of the Carrington storm?

The geomagnetic storm during the Carrington event, which occurred on 2 September 1859, displayed extremely rapid recovery. The geomagnetic field increased by approximately 650 nT/h at Bombay, India, and by >300 nT/h in 1-h averaged data. Although the rapid recovery is considered due to a sudden increase in the magnetopause current, a sudden decrease of the ring current, or/and a sudden enhancement of the ionospheric currents, this study focuses on the ring current decay. The Carrington rapid recovery had a time constant (approximately 1 h) comparable to the storm development (i.e., decrease in the geomagnetic field), indicating that energy loss from the ring current region is predominantly controlled by E × B convection transport which is responsible for energy input during the storm main phase. This feature has led us to a hypothesis that the flow-out of dense ring current ions and injections of tenuous plasma sheet ions caused the rapid decay of the ring current and in turn the storm rapid recovery. This study examines whether the Carrington rapid recovery can be explained by the flow-out effect. We extend the empirical Burton’s model to a model that takes into consideration a sudden change in solar wind density which is correlated with plasma sheet density. We first apply the extended Burton’s model to previously observed four intense magnetic storms (Dst minimum < −200 nT) for which solar wind data are available. Using the best fit parameters found by forward modeling, the extended model estimates the recovery of the Carrington storm. The estimate indicates that a solar wind structure with a density bump by approximately 100 cm−3 (and southward interplanetary magnetic field (IMF) of 65 nT and solar wind speed of 1,500 km/s) can cause the rapid recovery under a continuous southward IMF condition. We conclude that the flow-out effect plays a significant role in producing the rapid recovery of the Carrington storm.


Background
A magnetic storm is characterized by a negative excursion of the geomagnetic field at low latitudes (e.g., Gonzalez et al. 1994).As a geomagnetic storm is a global phenomenon, its development and recovery are well represented by a decrease of the disturbance field (Dst) index, which is derived by averaging the geomagnetic field observed by four stations located at low latitudes (http://wdc.kugi.kyoto-u.ac.jp/dstdir/index.html).Dst decreases with increasing geomagnetic activity; it reaches −30 to −100 nT at minimum for small and moderate storms.In response to a powerful/large coronal mass ejection (CME) or high-speed solar wind stream, an intense storm with Dst decreasing down to < −100 nT can occur.
A magnetic storm recovers to the pre-storm level within several to 10 days.As far as the geomagnetic storms that have ever been observed, the recovery is more rapid for an intense storm than for a small/moderate storm (e.g., Yokoyama and Kamide 1997).Most of intense storms and some of moderate storms show a two-step recovery with a quick Dst increase followed by a slow increase to the pre-storm level.These characteristics indicate that the storm recovery is controlled by at least two or several physical processes that occur on different time scales.
On 1 September 1859, a huge solar flare was observed by Carrington (1859) and Hodgson (1859).About 17 h and 40 min later, a magnetogram installed at Bombay, India, exhibited an extremely large depression (approximately 1,600 nT in 2 h) (e.g., Tsurutani 2003;Siscoe et al. 2006).The extremely intense geomagnetic storm, known as the most intense storm ever recorded, called the Carrington storm, recovered by >1,300 nT in 2 h.The top panel of Figure 1 shows the time profile of the geomagnetic field at Bombay, India (after Siscoe et al. (2006), modified from Tsurutani et al. (2003)).Even the time profile of the 1-h averaged data shows that an approximately 850-nT decrease followed by an approximately 600-nT increase occurred on similar time scales.The rapid recovery preceded a several-hour-long slow recovery, and after a relatively small enhancement of another storm activity, a day-long slow recovery back to the pre-storm level.
It is well accepted that the global geomagnetic field disturbance predominantly results from changes in intensity of the current flowing in the Earth's inner magnetosphere, called the ring current, which is mainly controlled by the ion pressure and its spatial gradient (e.g., Ebihara and Ejiri 2003 and references therein).The ion pressure in the core ring current region (i.e., 3 < r < 8 R E , where r is the geocentric distance) is dominated by ions with energies of 1 to a few 100 s keV (e.g., Williams 1981).Loss of those energetic ions results in pressure decrease and in turn ring current decay and consequently causes recovery from the depression of the geomagnetic field.
A large number of observational and theoretical studies have investigated physical processes that can explain the storm rapid recovery.It remains unresolved what process dominates over the others, although the slow recovery is well recognized as a result of charge exchange collisions between ring current protons and low-energy neutrals called the geocorona.The major processes that could contribute most to the rapid recovery are (1) charge exchange loss of singly charged oxygen ions (O + ) (e.g., Hamilton et al. 1988;Daglis et al. 1999); (2) flow-out of ring current ions through the dayside magnetopause, under the conditions of continued/weakened convection where newly injected ions from the plasma sheet are more tenuous than those leaking out into the magnetosheath (e.g., Liemohn 1999Liemohn , 2001;;O'Brien et al. 2002;Keika et al. 2005); (3) pitch angle scattering of ring current ions in the stretched near-Earth magnetotail with small curvature and consequent precipitation into the atmosphere (e.g., Ebihara et al. 2011); and (4) pitch angle scattering of ring current ions through resonant interactions with electromagnetic ion cyclotron (EMIC) waves and consequent precipitation into the atmosphere (e.g., Cornwall et al. 1970;Kozyra et al. 1998a;Jordanova et al. 2001).The details of the processes and related previous studies are summarized in a review paper by Ebihara and Ejiri (2003).
The rapid recovery of the Carrington event (see Figure 1) remains an open question.Siscoe et al. (2006) and Li et al. (2006) calculated Carrington's recovery from the results of their simulation using a solar wind structure that could have caused the Carrington storm development; they made use of an empirical formula, Burton's model (Burton et al. 1975).The estimated recovery was much slower than the observed one.Aguado et al. (2010) proposed a hyperbolic model of storm recovery, and Cid et al. (2013) applied the model to intense storms including the Carrington storm.They showed high accuracy of the hyperbolic function to reproduce the rapid recovery.Aguado et al. (2010) and Cid et al. (2013) commented on the existence of diverse processes responsible for the two-step recovery.However, they did not mention about the relative significance of ion loss processes that could occur in the inner magnetosphere.Li et al. (2006) introduced a sudden enhancement of the solar wind density that starts at the time of the storm maximum.The density enhancement, which enhances the magnetopause current, produces recovery from the storm development.They suggested that a density sudden increase up to 1,800 cm −3 can reproduce the observed rapid recovery.Such an extreme solar wind density enhancement is, however, unrealistic because it is 1 order magnitude higher than the observed limit at 1 AU distance.
In this paper, we examine whether the Carrington rapid recovery can be explained by the flow-out effect associated with changes in the plasma sheet density (i.e., leakage of denser ions and injections of more tenuous ions), by extending Burton's model and applying it to previously observed intense storms.There have been a large number of modeling studies that demonstrated a significant role of plasma sheet density in variations of the ring current intensity and the storm magnitude.A dense plasma sheet during a period of enhanced convection results in higher ion pressure in the inner magnetosphere and in turn stronger ring current (e.g., Chen et al. 1994;Kozyra et al. 1998b;Ebihara and Ejiri 2000;Ebihara et al. 2005).
A storm recovery due to the flow-out effect can be explained mostly by the rate of plasma sheet density change and the drift speed of energetic ions.A density change determines how much the magnetic field magnitude changes, and the ion drift speed determines the time scale of the magnetic field change.The ion drift is mostly controlled by E × B drift for both ion injections and flow-out.The Carrington event consists of the rapid main phase and rapid recovery phase, both of which have comparable time constants (approximately 1 h).This feature indicates that ion transport during the rapid recovery is controlled by a process that dominates ion transport during the main phase, namely, the E × B drift.We thus hypothesize that the ion flow-out effect is a plausible process that causes the rapid recovery.
The fact that the Carrington event recovered on a similar time scale to its development concerns us about possible enhancements of strong ground-induced currents (GIC) (e.g., Pulkkinen et al. 2012) even for the recovery phase.GIC could cause severe impact on our modern technological infrastructure such as pipelines and power grids as occurred during the 13 March 1989 storm (Bolduc 2002).Rapid recovery is thus a space weather issue that needs to be resolved like large geomagnetic sudden commencements (Araki 2014).
In the ' Analysis and results' section, we extend Burton's model to take into consideration the effects of solar wind density changes on storm development and decay.We then apply the extended model to four intense magnetic storms that occurred in years 2001, 2003, and 2004.The 1-h averaged symmetric disturbance field in the horizontal direction (SYM-H) during the storms is presented in Figure 1.For each event, by forward modeling, we find a pair of the model parameters that best fits the observed SYM-H temporal profile.We then apply the model to the Carrington event, using solar wind drivers estimated/used in several previous papers and the best fit parameters derived for the four intense storms.We present the results and discuss about preferred conditions for the flow-out process and other significant processes in the 'Results and discussion' section.We summarize our conclusions in the 'Conclusions' section.

Extended Burton's model
In order to examine the flow-out effects on the storm rapid recovery, we extend Burton's empirical model by adding a new term that represents a change in the solar wind density.Temporal variations of the Dst index based on the original Burton's equation (Burton et al. 1975) are expressed as: where Dst 0 is the corrected Dst for which the contributions from the current other than the ring current are subtracted.In this study, we subtracted three types of currents: the magnetopause current, the tail current, and the inductive current on/under the ground.It is well accepted that the contribution from the magnetopause current is proportional to the square root of the solar wind dynamic pressure (e.g., Siscoe et al. 1968).The relationship is expressed as: where b and c are constants.In this study, we use b = 7.26 nT/nPa 1/2 and c = 11.0 nT, which were derived by O'Brien and McPherron (2000).The contributions from the ground-induced current and the tail current have been controversial issues (e.g., Langel and Estes 1985;Alexeev et al. 1996;Turner et al. 2000).In this study, we assume 25% contributions from both currents, taking numbers close to the estimates by Langel and Estes (1985) and Turner et al. (2000).
The energy input term, Q, of Equation 1 is defined as: where a and E c are constants, and V and B s are the solar wind speed and southward component of the interplanetary magnetic field (IMF), respectively.E c was estimated as 0.49 mV/m by O' Brien and McPherron (2000).τ is the recovery time constant.The Dst recovery rate represented by the loss term increases with increasing Dst, indicating that the loss process(es) considered here is(are) more effective as a storm develops and in turn the plasma pressure in the inner magnetosphere increases.Charge exchange loss can be categorized in this group.Previous studies showed that a constant value of τ cannot reproduce the rapid recovery and two-step recovery of intense storms (e.g., Hamilton et al. 1988;Gonzalez et al. 1989;Liemohn et al. 1999;Monreal MacMahon and Llop-Romero 2008).In this study, we choose τ of 1 day to represent charge exchange collisions of <50 keV protons or approximately 100 keV O+ ions with the geocorona at L = 3 to 4 (Ebihara and Ejiri 2003).This study focuses most on evaluating how quickly the flow-out effect can recover the Carrington storm.For simplicity, we take into consideration the most general, best accepted loss process, that is, charge exchange of energetic (1 to a few 100 s keV) ions.
The loss term of the original Burton's model assumes the rate of Dst recovery proportional to the plasma pressure in the inner magnetosphere.However, a change in the plasma sheet density, which is the key phenomenon that makes the flow-out process effective, is independent of plasma pressure already stored in the inner magnetosphere.The rate of plasma sheet density variations and the drift speed of energetic ions determine the rate of storm recovery.Thus, in order to describe the flow-out effect, we need to add a new term that changes as a function of plasma sheet density variations and the global convection electric field.Since the plasma sheet density and the convection electric field are well correlated with the solar wind density and electric field, respectively (Terasawa et al. 1997;Borovsky et al. 1998;Ebihara and Ejiri 2000), we propose an extended Burton's model as: where τ is the time constant; the input term, Q, is the same as the original Burton's model: and the new term, Ψ, is expressed as: where N SW is the solar wind density, B s is the IMF southward component, and E c is a constant.The time variation of N SW is defined as: where Δt r is the time required for the plasma sheet density to respond to changes in the solar wind density.Δt d is chosen as the typical time required for ring current ions to drift from the near-Earth magnetotail to the dayside magnetopause.When energetic ions are injected into the inner magnetosphere (i.e., crossed the outer edge of the inner magnetosphere on the night side) at t = t 0 , the injected ions drift to the dayside magnetopause for a time period of Δt and eventually flow out through the dayside magnetopause at t = t 0 + Δt d .What affects the ring current intensity is the difference in density between the outflowing ring current ions and ions newly injected from the plasma sheet at t = t 0 + Δt d .Thus, for a Dst change at time t, the density of the outflowing ring current ions and the injected ions are correlated with the solar wind density at t = t − Δt r − Δt d and t = t − Δt r .In this study, we chose Δt d = 4 h.

Application to observed intense storms
We surveyed intense storm events during which the solar wind and IMF were monitored.We found four events, the first three of which are listed in the top tenth intense storms recorded since year 1957.The bottom four panels of Figure 1 present the time profile of the 1-h averaged SYM-H index for the four storms.We used the SYM-H index instead of Dst index because it has the same time resolution as the OMNI data.The minimum of 1-h averaged SYM-H is −379.

Search of the best fit parameters
We performed forward modeling to find a set of the parameters (α and β) of Equations 5 and 6 that best matches our extended model to the 1-h averaged SYM-H profile of each storm.Our forward modeling takes the following steps.
A. The SYM-H time profile based on our extended model (Equation 4) is calculated for different parameter pairs [α i , β i ]. α i is chosen from −5.0 to −0.1 with 0.1 increment; β i is chosen from −0.1 to 0.0 with 0.01 increment.B. [α j , β j ] pairs are found that provide the SYM-H minimum within 5 nT of the observed SYM-H minimum.C. The pair [α 0 , β 0 ] that minimizes the mean square (σ 2 ) of the differences between the modeled SYM-H and the observed SYM-H is selected from the [α j , β j ] pairs.σ 2 is calculated as:

2004-
where t start and t end are the start and end times of a storm interval of interest, respectively.In this study, we defined t start as 00 UT of the storm start date and t end as 24 UT of the storm end date which we chose as 5 days after the start date.We conducted the same analysis for different values of Δt r .We chose Δt r from 0 to 5 h with 1-h increment for the November 2003, April 2001, and April 2003 events.We chose Δt r from 10 to 20 h with 1-h increment for the November 2004 event because it is apparent that a >10-h time lag is needed for the solar wind density decrease to account for the storm rapid recovery.The parameter pair [α 0 , β 0 ] and Δt r that returns the smallest σ 2 are regarded as the best fit.

Results and discussion
Figure 2 presents IMF Z component of the magnetic field (Bz) in GSM coordinates, solar wind speed, solar wind density, and observed SYM-H, all of which are averaged over 1 h, for the November 2004 storm.The bottom panel also shows SYM-H derived from the original Burton's model (blue line) and SYM-H derived from our extended model (red line) with the best fit parameters.We chose the constant a of Equation 3 so that the original Burton's model is consistent with the observed SYM-H at its minimum.Data and model results for the November 2003 storm, the April 2001 storm, and the March 2001 storm are shown in Figures 3, 4, and 5, respectively, in the same format as Figure 2. All the storms were caused mostly by an increase of solar wind speed and a southward excursion of IMF (and in turn an increase of the dawn-dusk electric field).The pair of the best fit parameters α 0 , β 0 , Δt r is listed in Table 1, along with the SYM-H minimum, the IMF Bz minimum, and the solar wind density maximum.For all events, the SYM-H profile derived from the original Burton's model does not show rapid recovery similar to the observed SYM-H.Our extended model, however, displays rapid recovery; the rapidness of the early recovery is comparable to that of the observed SYM-H.The extended model also reproduces two-step recovery.For the March 2001 event, the original Burton's model displays two-step development rather than two-step recovery; our extended model shows rapid development followed by rapid recovery.
All events except for the November 2004 event display an increase in solar wind density almost simultaneously with an increase in the solar wind speed and a decrease in IMF Bz.The enhanced density reduced down to the pre-event level in 12 to 24 h, that is, when or before the IMF Bz negative excursion recovers to the pre-event level.These features indicate enhanced density in the magnetic cloud generated by a CME event.Such a density bump plays the key role in the extended model producing rapid recovery.
For the November 2004 event, a density bump appears in front of the southward IMF region (probably the sheath region of the interplanetary shock).IMF was northward during the period of the density enhancement.Since plasma entry across the flank-side magnetopause needs longer time than transport along the filed lines reconnected at the dayside magnetopause, the dense solar wind plasma could have spent a longer time to enter into the plasma sheet than the southward IMF cases.

Modeling of the Carrington recovery
We examine what SYM-H time profile the extended model produces with solar wind inputs that could cause the Carrington event.Figure 6 shows previously modeled solar wind inputs (speed, IMF Bz, and density) and the SYM-H time profiles estimated by the original Burton's model and our extended model.We used the best fit parameter pair obtained for the November 2004 storm.The solar wind speed and IMF Bz are the same as those used by Li et al. (2006) except for the duration of southward IMF.IMF Bz is 2 nT before the event, followed by a decrease down to −65 nT and recovery to 2 nT.The duration of southward IMF is 2 h, 1 h longer than that used by Li et al. (2006).The solar wind speed is 450 km/s before the event, increases to 1,500 km/s, and recovers to 450 km/s.We used Δt r = 1 h and Δt d = 4 h.We chose the density profile such that the density enhancement appears in the early part of the southward IMF interval; this is a reasonable choice because all of the four intense storms examined in this study were accompanied by a density enhancement during the part/entire interval of strong southward IMF (see Figure 2).The density is 5 cm −3 except for a 1-h bump of 100 cm −3 .
Our model derived a SYM-H time profile showing a rapid recovery in the time scale comparable to the storm development.The original Burton's model did not reproduce such a rapid recovery.Figure 7 presents the SYM-H profiles derived from the original Burton's model and our model with the best fit parameter set obtained for each storm event examined in this study.The parameter pair is denoted at the top of each panel.Although the storm maximum level (i.e., SYM-H minimum) differs among cases, it is noted that our extended model reproduced two-step recovery for all cases.

Preferred conditions for the flow-out effects
The flow-out process is effective in decaying the ring current when (1) plasma sheet ions injecting into the inner magnetosphere are more tenuous than ring current ions and (2) drift trajectories encounter the dayside magnetopause.These conditions require (i) a decrease in plasma sheet density and (ii) strong and prolonged E × B convection.Since plasma sheet density is correlated with solar wind density (Terasawa et al. 1997;Borovsky et al. 1998;Ebihara and Ejiri 2000), a decrease in solar wind density is needed.For the continuous E × B drift, the dawn-dusk convection electric field needs to remain strong, and thus, solar wind speed needs to keep fast and/or the IMF southward component remains large.In summary, the flow-out effect is significant when solar wind density decreases during a period of fast solar wind and/or southward IMF.
For three of the four storm events examined in this study (the November 2003, April 2001, and March 2001 storm), an increase in solar wind density occurs during the earlier part of a period of strong southward IMF, and the enhanced density decreases during the later part.IMF was southward even after the density peak.For the November 2004 event, solar wind density sharply increases before IMF turns southward; we speculate that plasma sheet density becomes high before the convection electric field becomes large.We believe that solar wind structures with a density enhancement only in the earlier part of or prior to a southward IMF period are not unusual.
Our estimates of the Carrington recovery based on the extended Burton's model suggest that a solar wind density decrease by approximately 100 cm −3 is high enough to produce the observed rapid recovery.We believe that the number is more reasonable than a density of approximately 2,000 cm −3 estimated by Li et al. (2006), who proposed that such a huge increase in solar wind density drastically increases the magnetopause current and in turn recovers the Dst/SYM-H index quickly.

Response time of plasma sheet density
For three of the four storm events examined in this study (the November 2003, April 2001, and March 2001 storm), the best fit Δt r was ≤2 h.The quick response is supported by the density bump that occurred during a period of southward IMF.We believe that the dense solar wind that entered into the plasma sheet through  (2003), Siscoe et al. (2006), andLi et al. (2006) except for the duration of southward IMF.A longer duration was used for this study.IMF, interplanetary magnetic field; Bz, the Z component of the magnetic field; GSM, geocentric solar magnetospheric; SYM-H, symmetric disturbance field in the horizontal direction.
the reconnected field lines convected at a speed comparable to the magnetosheath solar wind.For the November 2004 storm, on the other hand, we needed to adopt the response time of 12 h to best fit our model to the observed SYM-H.This is probably because the dense solar wind took a long time to enter into the plasma sheet under the condition of northward IMF.We assume plasma transport across the flank-side magnetopause due to instabilities such as the Kelvin-Helmholtz instability caused by flow shear at the magnetopause boundary.A larger number of intense storm events are necessary to identify the typical response time and/or its relationship with other solar wind parameters such as solar wind speed and IMF directions.

Effects of ionospheric outflow
Outflow of ionospheric ions into the magnetotail is enhanced after an interplanetary shock hit the magnetosphere and when energetic ions/electrons precipitate into the ionosphere particularly during substorms (e.g., Moore et al. 1999).Such outflowing ionospheric ions increase plasma density.This effect was not taken into consideration in our model.This study is the first step to incorporate the flow-out effect associated with solar wind density variations into Burton's model.It is well accepted that plasma sheet density is statistically well correlated with solar wind density (e.g., Terasawa et al. 1997).Effects of ionospheric outflow will be addressed in our future models.

Ring current decay proportional to ring current energy content
O' Brien and McPherron (2000) found a correlation between τ and VB s and suggested that VB s controls the position of the ring current.They assumed that the position is controlled by the convection electric field, which is correlated with the solar wind electric field.As the convection electric field increases, ring current ions are transported closer to the Earth, and the charge exchange loss rate increases because of higher density of the geocorona.We chose a constant value of τ for our model in this study to focus on evaluating whether the flow-out effect due to solar wind density decrease contributes to ring current rapid decay.We used a typical loss time scale of the most general, best accepted ring current ion loss process, charge exchange with the geocorona at L = 3 to 4. Time-varying τ will be in our future studies.Wang et al. (2003) and Xie et al. (2008) proposed τ as a function of P d for northward IMF and attributed the correlation to the size of the magnetosphere.The new loss term of our model depends on ΔN SW and in turn ΔP d .We believe that the flow-out effect expressed by the new loss term can be distinguished from the effect of a change in the magnetospheric size.The change does not affect the densities of both injected and outflowing ions but the energy of outflowing ions, because the magnetic field strength adjacent to the dayside magnetopause changes.

Conclusions
We examined intense geomagnetic storms in order to address an open question: why did the Carrington storm recover very rapidly (>800 nT in a few hours)?We tested a hypothesis that the rapid recovery is due to the flow-out of energetic ions that generate the westward ring current, by extending Burton's model.We added a new term representing variations of solar wind density which is correlated with plasma sheet density.We then applied the extended model to previously observed intense magnetic storms.Plasma sheet density variations make the flow-out process effective under the condition where the convection electric field remains strong.When the density drops, energetic ions contributing to the ring current flow out through the dayside magnetopause and are replaced by low-density plasma transported from a tenuous plasma sheet.We compared the time profile of the observed SYM-H (SYM-Hx obs ) with that derived from our extended model (SYM-H ext ).We performed forward modeling to find the best fit of SYM-H ext to SYM-H obs .The best fit SYM-H ext displayed more rapid recovery than SYM-H based on the original Burton's model.Furthermore, our model reproduced a two-step recovery, a rapid initial recovery followed by a slow one.We also applied the best fit parameters to the Carrington event, utilizing solar wind data used by Li et al. (2006).Our model produced rapid recovery similar to the Carrington event recovery when the solar wind density increases up to approximately 100 cm −3 and suddenly drops under a condition of continuous southward IMF of 65 nT and solar wind speed of 1,500 km/s.Such a solar wind density enhancement, as well as the IMF strength and the solar wind speed, is not beyond the observed limit.We conclude that the flow-out effect plays a significant role in the rapid recovery of the Carrington storm.

Figure 1
Figure1Geomagnetic field variations for the Carrington event and the intense storms examined.The top panel shows geomagnetic field variations at Bombay, India, during the Carrington event (afterSiscoe et al. (2006), modified fromTsurutani et al. (2003)).The bottom four panels present 1-h averaged SYM-H index for the intense storms examined in this study.SYM-H, symmetric disturbance field in the horizontal direction.
3 nT for the November 2004 storm, −461.3 nT for the November 2003 storm, −246.3 nT for the April 2001 storm, and −412.0 nT for the March 2001 storm.We used 1-h averaged OMNI solar wind data for the parameters required for the extended Burton's model, specifically solar wind density, speed, and IMF north-south component.

Figure 2
Figure2Solar wind speed, solar wind density, IMF Bz, and 1-h averaged SYM-H for the November 2004 storm.The bottom panel also includes 1-h SYM-H derived from the original Burton's model (blue line) and 1-h SYM-H derived from our extended Burton's model (red line).IMF, interplanetary magnetic field; Bz, the Z component of the magnetic field; GSM, geocentric solar magnetospheric; SYM-H, symmetric disturbance field in the horizontal direction.

Figure 3
Figure 3 Same as Figure 2 except for the November 2003 storm.IMF, interplanetary magnetic field; Bz, the Z component of the magnetic field; GSM, geocentric solar magnetospheric; SYM-H, symmetric disturbance field in the horizontal direction.

Figure 4
Figure 4 Same as Figure 2 except for the April 2001 storm.IMF, interplanetary magnetic field; Bz, the Z component of the magnetic field; GSM, geocentric solar magnetospheric; SYM-H, symmetric disturbance field in the horizontal direction.

Figure 5
Figure 5 Same as Figure 2 except for the March 2001 storm.IMF, interplanetary magnetic field; Bz, the Z component of the magnetic field; GSM, geocentric solar magnetospheric; SYM-H, symmetric disturbance field in the horizontal direction.

Table 1 Figure 6
Figure6Solar wind conditions and the SYM-H time profiles estimated from the original Burton's model and our model.We used the parameter pair that best fits the observed SYM-H for the November 2004 storm.The solar wind inputs are almost the same as those estimated byTsurutani (2003),Siscoe et al. (2006), andLi et al. (2006) except for the duration of southward IMF.A longer duration was used for this study.IMF, interplanetary magnetic field; Bz, the Z component of the magnetic field; GSM, geocentric solar magnetospheric; SYM-H, symmetric disturbance field in the horizontal direction.

Figure 7
Figure 7 SYM-H time profiles of the Carrington storm modeled by the original Burton's model and our model.We used each pair of the best fit parameters derived for the intense storms examined.