A Novel Magnetosphere–Ionosphere Coupling Model for the Onset of Substrom Expansion Phase

Abstract

A novel magnetosphere–ionosphere (M-I) coupling model is proposed to simulate the brightening of the onset auroral arc of a magnetospheric substorm event. The new M-I coupling model is modified from the M-I coupling model proposed by the Alaska research team in 1988. We adjust the magnetospheric boundary conditions by including the Hall effects in the thin current sheet and allowing the spatial distributions of the reflection–transmission coefficient to vary with time. As a result, brightening and poleward drifting of multiple auroral arcs appear for the first time in an M-I coupling model. The new results indicate that the coupled Hall effects in the near-Earth plasma sheet and the E-region ionosphere play a vital role in triggering the onset of a magnetospheric substorm.

Introduction

Kan and Sun. (1985), Kan et al. (1988), and Kan and Sun. (1996) have proposed a highly simplified but an elegant simulation scheme to model the substorm-associated magnetosphere–ionosphere (M-I) coupling processes. Their results show a plasma flow pattern similar to the observed westward traveling surge during the substorm. The cross-polar-cap potential drop in the E-region ionosphere obtained in their model is slightly lower than the initial input cross-polar-cap potential drop. The nonuniform distributions of the conductivities and electric fields in the E-region ionosphere can result in the region-1 and region-2 field-aligned current distributions (e.g., Iijima and Potemra, 1976). Similar results have been obtained in the M-I coupling simulations studied by Miura and Sato. (1980); Miura and Sato. (1981). However, these early simulation studies (e.g., Miura and Sato, 1980, Miura and Sato, 1981); (Kan and Sun, 1985; Kan et al., 1988) fail to show onset aurora arc–associated upward field-aligned currents in the midnight region (e.g., Akasofu, 1964).

Kan and Sun. (1996) changed the conductivity enhancement scheme and modified the convection electric field by adding localized convection electric fields on the magnetotail boundary of the M-I coupling model. They successfully obtain upward field-aligned currents in the midnight region (Kan and Sun, 1996).

Based on satellite observations, the near-Earth plasma sheet shows a tail-like structure before the onset of a substorm. It changes to a dipole-like structure after the onset of the substorm (e.g., McPherron, 1972). Kaufmann. (1987) showed that the dipolarization of the near-Earth plasma sheet at the onset of a substorm is associated with the formation of current wedges. Dipolarization of the magnetic field occurs inside the current wedges, whereas the thinning processes continuously take place outside of the current wedges and tailward from the current wedges (Kaufmann, 1987). Ohtani et al. (1992) found several substorm onset events with explosive thinning of the plasma sheet before the onset of substorms.

It is believed that disruption of cross-tail currents in the near-Earth plasma sheet can trigger the dipolarization process. Theoretical models have been proposed to explain the current disruption in the near-Earth plasma sheet. These theoretical models include, but not limited to, the magnetic reconnection–associated resistive tearing-mode instability (e.g., Furth et al., 1963; Coppi et al., 1966; Schindler, 1974), the ion Weibel instability (e.g., Lui et al., 1993; Lui et al., 2008), the ballooning instability (e.g., Roux et al., 1991; Cheng and Lui, 1998), and the unloading instability triggered by the incident Alfvén waves from the ionosphere (Kan and Sun, 1996). Lyu and Chen. (2000) proposed another type of unloading instability due to the M-I coupling process and the nonuniform Hall effect in the explosively thinning region.

When a current sheet has a finite normal magnetic field component, and when the thickness of the current sheet is equal to or smaller than the gyroradius of the thermal ions, a sufficient amount of ions will become unmagnetized, whereas electrons are still magnetized. The unmagnetized ions will move across the thin current sheet along the electric field direction with a meandering trajectory. Magnetized electrons and unmagnetized ions will set up a Hall electric field along the drift direction of the magnetized electrons. As a result, the electric field will rotate left-handed with respect to the ambient magnetic field. If the thickness of the current sheet is nonuniform, it can lead to a nonuniform left-handed rotation of the electric field and results in localized field-aligned currents. The upward field-aligned current from the ionosphere to the thin current sheet in the pre-midnight region can increase the normal magnetic field component in the thin current sheet downward from the field-aligned current. Since the meandering motion of ions in the thin current sheet can enhance the cross-tail current intensity, increasing the current sheet thickness will reduce the number of unmagnetized ions. Decreasing the number of unmagnetized ions can reduce the cross-tail current intensity and further increase the current sheet thickness. This positive feedback process can result in current disruption and trigger the onset of a substorm (Lyu and Chen, 2000).

Unlike the other instabilities analyses, where waves of a given frequency are amplified with a well-determined growth rate, the unloading instability proposed by Lyu and Chen. (2000) does not have a well-defined growth rate or a well-defined frequency. Since the background magnetic field of the unloading instability is highly nonuniform and time dependent, we can only qualitatively show a positive feedback process to trigger the current disruption in the thin current sheet.

This study aimed to modify the M-I coupling model proposed by Kan et al. (1988) by including the Hall effects in the magnetotail as proposed by Lyu and Chen. (2000). The coupling between the nonuniform Hall effect in the high-latitude E-region ionosphere and the nonuniform Hall effect on the magnetospheric boundary will be examined.

Basic Equations of the Magnetosphere–Ionosphere Coupling Model

The new M-I coupling model is modified from the M-I coupling model proposed by Kan et al. (1988). For convenience, we shall call the M-I coupling model proposed by Kan et al. (1988) as the KZA′88 model. The high-latitude ionosphere in the KZA′88 model is a plane with uniform magnetic field B₀ perpendicular to the ionosphere. The nonuniform electric fields and conductivities yield nonuniform electrical currents in the high-latitude ionosphere. Since no charge accumulation could take place in the timescale of the KZA′88 model, the divergence of the height-integrated perpendicular current density in the E-region ionosphere implies a downward field-aligned current at the top of the E-region ionosphere. Likewise, the convergence of the height-integrated perpendicular current density can lead to an upward field-aligned current at the top of the E-region ionosphere. The upward field-aligned current at the top of the E-region ionosphere can be written in the following form (e.g., Kan et al., 1988).where is the unit vector parallel to , is the height-integrated Pedersen conductivity, is the height-integrated Hall conductivity, and is the vertical averaged electric field in the E-region ionosphere. The direction of is perpendicular to the background magnetic field in this study.

The coupling between the ionosphere and the magnetosphere in the KZA′88 model is achieved by the Alfvén waves that carry the electric field and field-aligned current and propagate back and forth between these two regions along the magnetic field lines. Based on the Walén relation (Walén, 1944), Kan and Sun. (1985) have shown that the Alfvén waves with group velocity parallel to the background magetic field will carry field-aligned currents:where , is the perturbed magnetic field in the Alfvén wave, is the effective height-integrated conductivity in the Alfvén wave, is the Alfvén speed, is the background plasma mass density, and is the perpendicular component of the perturbed electric field in the Alfvén wave. Likewise, Alfvén waves with group velocity antiparallel to the background magnetic field will carry field-aligned currents:

Since the background magnetic field is downward in the northern hemisphere high-latitude ionosphere, we have , and the group velocity of the incident Alfvén wave is parallel to the background magnetic field, whereas the group velocity of the reflected Alfvén wave is antiparallel to the background magnetic field. Thus, Eqs 2, 3 yieldwhere is the perpendicular component of the perturbed electric field in the incident Alfvén wave, is the upward field-aligned current carried by the incident Alfvén wave, is the perpendicular component of the perturbed electric field in the reflected Alfvén wave, and is the upward field-aligned current carried by the reflected Alfvén wave.

The reflected Alfvén waves from the ionosphere are determined in the following way in the KZA′88 model. When the sum of the carried by the incident Alfvén waves and the preexisting field-aligned currents above the E-region ionosphere do not match the required at the top of the E-region ionosphere, as shown in Eq. 1, Alfvén waves will be reflected from the ionosphere so that the upward field-aligned currents carried by the reflected Alfvén waves satisfy

According to Eq. 1, the preexisting field-aligned current satisfieswhere , , and are the preexisting , , and , respectively.

Substituting Eqs 1, 4, 5, 7 into Eq. 6 to eliminate , , , and , respectively, for the northern hemisphere ionosphere, it yields (Kan et al., 1988) the following equation:

The electric field in the E-region ionosphere is determined in the following way in the KZA′88 model.

Substituting Eq. 9 into Eq. 8 to eliminate for the northern-hemisphere ionosphere, it yields

Note that the KZA′88 model is modified from the M-I coupling model proposed by Kan and Sun. (1985). For convenience, we shall call the M-I coupling model proposed by Kan and Sun. (1985) as the KS′85 model. The matching conditions of the field-aligned currents and the perpendicular electric field used in the KS′85 model are given as follows:

Namely, both the preexisting field-aligned current and the preexisting ionospheric perpendicular electric field are ignored in the KS′85 model.

Since a current loop cannot last forever, especially when there is a finite Pederson conductivity in the E-region ionosphere, the current intensity in the preexisting current loop should decrease with time. In this study, we ignore the preexisting field-aligned current and adopt the matching condition of the field-aligned currents used in the KS′85 model. On the other hand, the preexisting electric field should be included in the M-I coupling model. When the preexisting electric field in the ionosphere is strong enough, the neutral particles might be pushed by the ions. Since the ionospheric conductivities are obtained in the neutral wind moving frame, we should determine the neutral wind flow velocity and modify the preexisting electric field to the neutral wind moving frame. However, modeling the enhanced neutral wind is beyond the scope of this study. We shall ignore the neutral wind effect and determine the electric field in the E-region ionosphere based on Eq. 9. Note that Kan and Sun (1996) also used Eqs 11, 9 instead of Eqs 6, 12 to match the field-aligned currents and to model the ionospheric electric field, respectively.

Substituting Eqs 1, 4, 5 into Eq. 11 to eliminate , , and , respectively, for the northern-hemisphere high-latitude ionosphere, it yields the following equation:Substituting Eq. 9 into Eq. 13 to eliminate for the northern-hemisphere ionosphere, it yields the following equation:Eq. 14 yieldswhere is a vector field. For simplicity, we choose in our simulation. Solving Eq. 15 for with , it yields

Note that due to the incompressibility nature of the Alfvén waves, the perpendicular electric fields in Alfvén waves are assumed to be curl-free electric fields in both the KZA′88 model and the KS′85 model. We will discuss in Section 4 on how the electric fields in the incident Alfvén waves are obtained from a reflection–transmission model. The electric fields in the incident Alfvén waves should not be curl-free electric fields everywhere due to the nonuniformity of the reflection–transmission coefficient . Since the electric fields in the incident Alfvén waves are not curl-free, we shall allow the electric fields obtained from Eq. 16 not to be curl-free either. Nonzero and indicate that the waves carrying field-aligned currents are not pure Alfvén waves. We recalled that Eqs 2, 3 are obtained based on the Walén relation of the Alfvén-mode wave (Kan and Sun, 1985). To justify the use of Eqs 2, 3 in the M-I coupling model, we need to show that Eqs 2, 3 can be obtained directly from Maxwell’s equations of the low-frequency waves with a field-aligned propagating speed equal to the Alfvén speed.

The time derivative of the low-frequency Ampere’s law yieldsSubstituting Faraday’s law into Eq. 17 to eliminate , it yieldsThe field-aligned component of Eq. 18 iswhere is a scalar differential operator. The first term on the right-hand side of Eq. 19, , is associated with the localized field-aligned potential jump. The formation of field-aligned potential jump is a kinetic process which is not included in this wave-associated M-I coupling model. Ignoring in Eq. 19, it yields

The field-aligned currents carried by parallel propagating waves with wave speed equal to the Alfvén speed should have a function form . The corresponding wave equation isComparing Eqs 20, 21, it yields, for parallel propagating waves, the following equation:

Likewise, the field-aligned currents carried by antiparallel propagating waves with a wave speed equal to the Alfvén speed should have a function form . The corresponding wave equation is given as follows:

Comparing Eqs 20, 23, it yields, for antiparallel propagating waves, the following equation:

Equations 22 and 24 are the same as Eqs 2, 3. Namely, Eqs 2, 3 are applicable to the field-aligned propagating fast-mode or slow-mode waves as long as the wave speed is equal to the Alfvén speed.

Modeling the Enhancement of the Hall Conductance by Upward Field-Aligned Currents

The enhancement of the conductance by the upward field-aligned current is included in the KZA′88 model, where the conductance is the height-integrated conductivity. For strong upward field-aligned current, the precipitating energetic electrons will bombard the neutral particles and increase the ionization rate. In the KZA′88 model, the Hall conductance is modeled bywhere is the initial Hall conductance at the given location, is the electron thermal flux in the loss cone, the coefficient is chosen to be , and the critical current density is chosen to be in the KZA′88 model. Based on Eq. 25, the maximum Hall conductance in the KZA′88 model is limited by the following equation:which is only slightly higher than the maximum of the initial Hall conductance in their simulation. Kan and Sun (1996) have removed the third condition in Eq. 25 in their new M-I coupling simulation model. As a result, their simulations show a maximum Hall conductance above 50 mho (Kan and Sun, 1996).

Based on the ground and satellite observations, the upward field-aligned current density can be as high as 5  (Kamide and Horwitz, 1978; Kamide, 1982; Podgorny et al., 2003; Wing and Johnson, 2015; Bunescu et al., 2019) or even 10  (e.g., Pitout et al., 2015). The enhancement of the Hall conductance can reach 40  or even 80  (Robinson et al., 1985). The ratio of the Hall conductance to the Pedersen conductance increases with increasing the Hall conductance (Robinson et al., 1985). For simplicity, we choose a constant ratio , which is the same as the one used in the KZA′88 model. To improve the conductance enhancement, we constructed a new model to determine . The field-aligned currents in our model do not include the preexisting field-aligned currents, but the preexisting conductance is included in our new model. The conductance enhancement is rewritten in the following form:where is the initial Hall conductance; is the preexisting Hall conductance; is the decaying parameter, which satisfies ; and the coefficient is chosen to be similar to the one used in Eq. 25. Note that the conductance enhancement shown in Eq. 25 was obtained based on the electron acceleration by a steady field-aligned potential drop (Knight, 1973; Fridman and Lemaire, 1980; Kan and Kamide, 1985). Therefore, to determine the field-aligned potential drop from the upward field-aligned current, only the electron thermal flux in the loss cone is removed from . However, the electron acceleration process in a steady field-aligned potential drop is different from the electron energization and acceleration by a propagating wave along the magnetic field line. In this study, we found that the relationship between and will lead to a similar relationship between and . To limit the conductance enhancement in a finite range of a given , we assume increases linearly with increasing Hall conductance. Based on the Hall conductance and the field-aligned current density reported in Robinson et al. (1985), we modeled the in the following way:where is the preexisting Hall conductance. In this study, we choose , , , and . As a result, for , the conductance will be enhanced only if .

Modeling the Wave Reflections at the Magnetospheric Boundary

According to the KZA′88 model, the incident perturbed electric field at the step satisfieswhere is the reflection–transmission coefficient (Kan and Sun, 1985; Kan et al., 1988).

Kan and Sun (1985) considered in the near-Earth plasma sheet, and in the low-latitude boundary layer and the distant plasma sheet. The boundary of the open-field-line region mapping onto the ionosphere () is located at a higher latitude in the nightside ionosphere than it is in the dayside ionosphere in the early simulation studies (Kan and Sun, 1985; Kan et al., 1988). Kan and Sun (1996) used a time-dependent distribution. The boundary of is nearly a circle at 80° magnetic latitude at the beginning of their simulation. The boundary of gradually expanded toward lower latitudes on the nightside. As a result, the boundary of the open-field-line region mapping onto the ionosphere is located at a lower latitude in the nightside ionosphere than it is in the dayside ionosphere (Kan and Sun, 1996).

In this study, we consider the magnetotail maps to ionosphere within an 8-h sector, from 20 MLT to 04 MLT in the nightside ionosphere, where MLT denotes the magnetic local time. Eq. 29 is applicable to the region outside the 8-h sector. Inside the 8-h sector, denotes regions map to the thick current sheet and denotes regions map to the thin current sheet. Namely, for MLT between 20 MLT and 04 MLT, we have

Since the thickness of the near-Earth plasma sheet varies during a substorm event, the distribution of R_m inside the 8-h sector is allowed to change with time in this study.

In addition to the reflection–transmission coefficient R_m, we introduce a new coefficient to mimic the proposed Hall effect in the thin current sheet region (Lyu and Chen, 2000). Since the Hall effect is a dynamo process, the Hall electric field and the Hall current should point in opposite directions. The Hall current is in the direction of . The Hall electric field should be in the direction of . Thus, we modify Eq. 28 to the following form:

Figure 1 shows two examples of and distributions used in this study. Panels (A) and (B) show the distributions of with and , respectively, where is the magnetic latitude of at 23 MLT. Panels (C) and (D) show the distributions of with and , respectively. In this study, we choose the equatorward boundary of the open-field-line region () at magnetic latitude and the poleward boundary of the dipole-field-line region at magnetic latitude similar to those used in previous simulation studies (Kan and Sun, 1985; Kan et al., 1988). The contour is located at magnetic latitude outside the 8-h sector and gradually shifts to magnetic latitude at 23 MLT inside the 8-h sector. We choose to be the magnetic latitude of at 23 MLT, but not at 00 MLT, because the cross-tail current is enhanced by the duskward motion of the unmagnetized ions, which are accelerated by the dawn-to-dusk electric field. Thus, we expect the thin current sheet to extend over the pre-midnight region.

FIGURE 1

Simulation Results

Table 1 lists the simulation parameters used in the 12 simulation cases to be present in this section. The magnetic latitude is 54° in Cases 1 and 2 but varies from 54° to 62° in Case 3. We choose 60° in Cases 4a∼6c. To show the importance of the Hall effect in the near-Earth plasma sheet, we set in Cases 2, 4c, 5c, and 6c. We also turn off the conductance enhancement scheme in Cases 4b, 5b, and 6b to show the effect of conductance enhancement in the E-region ionosphere. The initial distributions of the Hall conductance in Cases 1, 2, 3, 6a, 6b, and 6c are similar to the one shown in the KZA′88 model. The initial distribution of the Hall conductance is uniform in Cases 4a∼5c. The uniform conductance is in Cases 4a, 4b, and 4c but in Cases 5a, 5b, and 5c.

Figure 2 shows the simulation results obtained from Case 1 simulation. Case 1 is characterized by time-independent and distributions with . Panels (A) and (B) show color-level plots of obtained at the 10th step and the 15th step of the simulation, respectively. Panels (C) and (D) show color-level plots of the upward field-aligned current density obtained from the simulation at the 10th step and the 15th step, respectively. Strong field-aligned currents and enhanced Hall conductances appear in the midnight region.

FIGURE 2

Figure 3 shows color-level plots of (A) and (B) obtained from the 10th step of the simulation study of Case 2. Case 2 simulation has a time-independent distribution with . The Hall effect in the thin plasma sheet is turned off in Case 2. A weak upward field-aligned current can be found on the duskside of the polar cap boundary, which can slightly enhance the conductance in that region. The simulation has reached a steady state with very little change on and at the 8th step and beyond. The enhancement of on the duskside of the polar cap boundary is less than 10 mho, which is slightly above the maximum of the initial Hall conductance in the diffuse aurora region.

FIGURE 3

Figure 4 shows simulation results of Case 3, in the same format as those shown in Figure 2. Case 3 is characterized by time-dependent and distributions with changes from to and back to Since the distributions of and are time-dependent, the maximum of the Hall conductance and the maximum of the upward field-aligned current density obtained in Case 3 are lower than those shown in Case 1.

FIGURE 4

Figure 5 shows the time evolution of (A), (B), and (C) obtained from the simulations of Case 1 (the solid curves), Case 2 (the dashed curves), and Case 3 (the dotted curves). denotes the maximum upward field-aligned current density in the entire simulation domain obtained at the kth step of the simulation. The absolute value of denotes the maximum of the downward field-aligned current density in the entire simulation domain obtained at the kth step of the simulation. denotes the maximum Hall conductance in the entire simulation domain obtained at the kth step of the simulation. The time evolution of is similar to the time evolution of due to the choice of the conductance enhancement scheme discussed in Eqs 26, 27. The Hall effects in the magnetotail lead to a similar growth phase () in Case 1 and Case 3. The time-independent distributions of and used in the Case 1 simulation yield a second growth phase at in the Case 1 simulation. Ignoring the Hall effects in the magnetotail, the three dashed curves of Case 2 show very little enhancement on the field-aligned current density and the Hall conductance.

FIGURE 5

Figure 6 shows the time evolution of (A), (B), and (C) obtained from the simulations of Case 4a (the solid curves), Case 4b (the dashed curves), and Case 4c (the dotted curves). Simulation Cases 4a, 4b, and 4c are characterized by uniform high conductances at the beginning of the simulations. The initial Hall conductance is The conductance enhancement scheme is turned off in Case 4b. Without gradient on the ionospheric conductances, the field-aligned current density obtained in Case 4b is lower than it in Case 4a. The plasma sheet Hall effect is turned off in Case 4c. Without the Hall effect in the plasma sheet, very little enhancement on the field-aligned current density and the Hall conductance can be seen in the simulation of Case 4c.

FIGURE 6

Figure 7 shows the time evolution of (A), (B), and (C) obtained from the simulations of Case 5a (the solid curves), Case 5b (the dashed curves), and Case 5c (the dotted curves). Simulation Cases 5a, 5b, and 5c are characterized by uniform low conductances at the beginning of the simulations. The initial Hall conductance is . The conductance enhancement scheme is turned off in the simulation study of Case 5b. The plasma sheet Hall effect is turned off in the simulation study of Case 5c. Without gradient on the ionospheric conductances, the field-aligned current density obtained in Case 5b shows a minimal growth rate. Without considering the Hall effect in the plasma sheet, very little enhancement on the field-aligned current density and the Hall conductance can be found in the simulation study of Case 5c.

FIGURE 7

Let us compare the simulation results of Case 4a and Case 5a shown in Figures 6, 7. The curve obtained in Case 4a shows a short lag phase followed by a growth phase from the 4th step to the 10th step, but with a relatively low saturation level . The curve obtained in Case 5a shows a prolonged lag phase followed by a growth phase from the 8th step to the 16th step, but with a relatively high saturation level at . The curve obtained in the simulation study of Case 4a shows a second peak at the 14th step and a third peak at the 18th step. The highest upward field-aligned current density can be found at the 18th step with . The differences between the curves shown in Case 4a and Case 4b are less significant than the differences between the curves shown in Case 5a and Case 5b

Figure 8 shows the time evolution of (A)(B) and (C) obtained from Case 6a (the solid curves), Case 6b (the dashed curves), and Case 6c (the dotted curves) simulations. The initial distributions of the conductances of the three cases are nonuniform and are similar to the one used in Cases 1–3. The conductance enhancement scheme in the ionosphere is turned off in the simulation study of Case 6b. The plasma sheet Hall effect is turned off in the simulation study of Case 6c. Again, without the Hall effect in the plasma sheet, very little enhancement on the field-aligned current density and the Hall conductance can be found in the simulation study of Case 6c. But the differences between the curves obtained in Case 6a and Case 6b are relatively small due to the preexisting nonuniform conductance in the diffuse aurora region. The curves of Case 6a show a growth phase from the 6th step to the 12th step and with a saturation level at .

FIGURE 8

Summary and Discussion

In summary, we have constructed a new M-I coupling model modified from the M-I coupling model proposed by Kan et al. (1988). We adjust the magnetospheric boundary conditions by including the Hall effects in the thin current sheet. As a result, multiple brightening auroral arcs appear in the midnight region and gradually move poleward. This region is characterized by and .

A steady midnight auroral arc has also been found in the simulation study by Kan and Sun (1996). The formation of the midnight arc results from an enhanced localized convection electric field added in their simulation. A localized convection electric field has also been considered in the early simulation studies (e.g., Kan et al., 1988), but no midnight aurora arc can be found. These results indicate that only a particular type of convection electric field can result in a midnight aurora arc.

The enhancement of ionospheric conductances by the upward field-aligned currents can increase the nonuniformity of the ionospheric conductances. The nonuniform Hall effect in the plasma sheet can result in the nonuniform rotation of the electric field in the plasma sheet. The nonuniform rotation of the electric field in the plasma sheet obtained in this study and the localized convection electric field proposed by Kan and Sun (1996) can result in upward and downward field-aligned currents in the midnight region.

The simple M-I coupling model provides much helpful information for 3-dimensional global simulations of magnetospheric substorms. The simulation results shown in Figures 6–8 in the last section indicate that the high conductance in the ionosphere can speed up the growth phase of a substorm event but result in a relatively weak upward field-aligned current in the auroral arc. On the other hand, a substorm event with low preexisting conductance in the ionospheric boundary requires a longer time to complete the growth phase. Still, it can build up a much stronger upward field-aligned current in the auroral arc. The simulation results discussed in Figures 6–8 also indicate that including the conductance enhancement processes on the ionosphere boundary can increase the field-aligned current intensity in a global simulation of the magnetospheric substorm. The Hall effect in the thin current sheet should also be included in future simulation studies of magnetospheric substorms.

References

• 1

  AkasofuS.-I. (1964). The Development of the Auroral Substorm. Planet. Space Sci.12, 273–282. 10.1016/0032-0633(64)90151-5

  • CrossRef
  • Google Scholar

• 2

  BunescuC.VogtJ.MarghituO.BlagauA. (2019). Multiscale Estimation of the Field-Aligned Current Density. Ann. Geophys.37, 347–373. 10.5194/angeo-37-347-2019

  • CrossRef
  • Google Scholar

• 3

  ChengC. Z.LuiA. T. Y. (1998). Kinetic Ballooning Instability for Substorm Onset and Current Disruption Observed by AMPTE/CCE. Geophys. Res. Lett.25, 4091–4094. 10.1029/1998GL900093

  • CrossRef
  • Google Scholar

• 4

  CoppiB.LavalG.PellatR. (1966). Dynamics of the Geomagnetic Tail. Phys. Rev. Lett.16, 1207–1210. 10.1103/PhysRevLett.16.1207

  • CrossRef
  • Google Scholar

• 5

  FridmanM.LemaireJ. (1980). Relationship between Auroral Electrons Fluxes and Field Aligned Electric Potential Difference. J. Geophys. Res.85 (A2), 664–670. 10.1029/JA085iA02p00664

  • CrossRef
  • Google Scholar

• 6

  FurthH. P.KilleenJ.RosenbluthM. N. (1963). Finite-Resistivity Instabilities of a Sheet Pinch. Phys. Fluids6, 459. 10.1063/1.1706761

  • CrossRef
  • Google Scholar

• 7

  IijimaT.PotemraT. A. (1976). The Amplitude Distribution of Field-Aligned Currents at Northern High Latitudes Observed by Triad. J. Geophys. Res.81 (13), 2165–2174. 10.1029/JA081i013p02165

  • CrossRef
  • Google Scholar

• 8

  KamideY.HorwitzJ. L. (1978). Chatanika Radar Observations of Ionospheric and Field-Aligned Currents. J. Geophys. Res.83 (A3), 1063–1070. 10.1029/JA083iA03p01063

  • CrossRef
  • Google Scholar

• 9

  KamideY. (1982). The Relationship between Field-Aligned Currents and the Auroral Electrojets: A Review. Space Sci. Rev.31, 127–243. 10.1007/bf00215281

  • CrossRef
  • Google Scholar

• 10

  KanJ. R.KamideY. (1985). Electrodynamics of the Westward Traveling Surge. J. Geophys. Res.90 (A8), 7615–7619. 10.1029/JA090iA08p07615

  • CrossRef
  • Google Scholar

• 11

  KanJ. R.SunW. (1985). Simulation of the Westward Traveling Surge and Pi 2 Pulsations during Substorms. J. Geophys. Res.90 (A11), 10911–10922. 10.1029/JA090iA11p10911

  • CrossRef
  • Google Scholar

• 12

  KanJ. R.SunW. (1996). Substorm Expansion Phase Caused by an Intense Localized Convection Imposed on the Ionosphere. J. Geophys. Res.101 (A12), 27271–27281. 10.1029/96JA02426

  • CrossRef
  • Google Scholar

• 13

  KanJ. R.ZhuL.AkasofuS.-I. (1988). A Theory of Substorms: Onset and Subsidence. J. Geophys. Res.93 (A6), 5624–5640. 10.1029/JA093iA06p05624

  • CrossRef
  • Google Scholar

• 14

  KaufmannR. L. (1987). Substorm Currents: Growth Phase and Onset. J. Geophys. Res.92 (A7), 7471–7486. 10.1029/JA092iA07p07471

  • CrossRef
  • Google Scholar

• 15

  KnightS. (1973). Parallel Electric fields. Planet. Space Sci.21, 741–750. 10.1016/0032-0633(73)90093-7

  • CrossRef
  • Google Scholar

• 16

  LuiA. T. Y.YoonP. H.ChangC.-L. (1993). Quasi-linear Analysis of Ion Weibel Instability in the Earth's Neutral Sheet. J. Geophys. Res.98 (A1), 153–163. 10.1029/92JA02034

  • CrossRef
  • Google Scholar

• 17

  LuiA. T. Y.YoonP. H.MokC.RyuC.-M. (2008). Inverse cascade Feature in Current Disruption. J. Geophys. Res.113, a–n. 10.1029/2008JA013521

  • CrossRef
  • Google Scholar

• 18

  LyuL. H.ChenM. Q. (2000). A Kinetic M-I Coupling Model with Unloading Instability at Onset of Substorm. in Proc. 5th International Conference on Substorms, St. Petersburg, Russia, 315–318. 16-20 May 2000, ESA SP-443.

  • Google Scholar

• 19

  McPherronR. L. (1972). Substorm Related Changes in the Geomagnetic Tail: the Growth Phase. Planet. Space Sci.20, 1521–1539. 10.1016/0032-0633(72)90054-2

  • CrossRef
  • Google Scholar

• 20

  MiuraA.SatoT. (1981). Global Simulation of Auroral Arcs, in Physics of Auroral Arc Formation, ed. by AkasofuS.‐I.KanJ., pp. 321–332. Washington, DC: Geophysical Monograph Series, Volume 25. 10.1029/GM025p0321

  • CrossRef
  • Google Scholar

• 21

  MiuraA.SatoT. (1980). Numerical Simulation of Global Formation of Auroral Arcs. J. Geophys. Res.85 (A1), 73–91. 10.1029/JA085iA01p00073

  • CrossRef
  • Google Scholar

• 22

  OhtaniS.TakahashiK.ZanettiL. J.PotemraT. A.McEntireR. W.IijimaT. (1992). Initial Signatures of Magnetic Field and Energetic Particle Fluxes at Tail Reconfiguration: Explosive Growth Phase. J. Geophys. Res.97 (A12), 19311–19324. 10.1029/92JA01832

  • CrossRef
  • Google Scholar

• 23

  PitoutF.MarchaudonA.BlellyP. L.BaiX.FormeF.BuchertS. C.et al (2015). Swarm and ESR Observations of the Ionospheric Response to a Field‐aligned Current System in the High‐latitude Midnight Sector. Geophys. Res. Lett.42, 4270–4279. 10.1002/2015GL064231

  • CrossRef
  • Google Scholar

• 24

  PodgornyI. M.PodgornyA. I.MinamiS.RanaR. (2003). The Mechanism of Energy Release and Field-Aligned Current Generation during Substorms and Solar Flares. Adv. Polar Upper Atmos. Res.17, 77–83.

  • Google Scholar

• 25

  RobinsonR. M.VondrakR. R.PotemraT. A. (1985). Auroral Zone Conductivities within the Field-Aligned Current Sheets. J. Geophys. Res.90 (A10), 9688–9696. 10.1029/JA090iA10p09688

  • CrossRef
  • Google Scholar

• 26

  RouxA.PerrautS.RobertP.MoraneA.PedersenA.KorthA.et al (1991). Plasma Sheet Instability Related to the Westward Traveling Surge. J. Geophys. Res.96 (A10), 17697–17714. 10.1029/91JA01106

  • CrossRef
  • Google Scholar

• 27

  SchindlerK. (1974). A Theory of the Substorm Mechanism. J. Geophys. Res.79 (19), 2803–2810. 10.1029/JA079i019p02803

  • CrossRef
  • Google Scholar

• 28

  WalénC. (1944). On the Theory of Sunspots. Ark. Mat. Astron. Fys.30A, 15.

  • Google Scholar

• 29

  WingS.JohnsonJ. R. (2015). Theory and Observations of Upward Field-Aligned Currents at the Magnetopause Boundary Layer. Geophys. Res. Lett.42, 9149–9155. 10.1002/2015GL065464

  • CrossRef
  • Google Scholar