The substorm current wedge and field‐aligned currents in MHD simulations of magnetotail reconnection

We present results of a three-dimensional MHD simulation of magnetotail reconnection, which, in contrast to earlier ones, starts from an equilibrium configuration that contains a region 1 type field-aligned current system at the plasma sheet/lobe boundary. This current system is found from equilibrium theory as the consequence of boundary conditions, which require that field lines at large distances from the Earth become aligned with the tail axis, consistent with observations and the expected consequences of the interaction with the solar wind. As in earlier simulations, the dynamic evolution develops from a slow diffusion to the fast growth of a three-dimensional tearing mode, initiated by the sudden occurrence or increase of resistivity, leading again to the formation and subsequent tailward ejection of a plasmoid, associated with fast plasma flows. In addition, we find changes of the electric current system which are consistent with the substorm current wedge picture inferred from observations, i.e., a deviation of cross-tail current from a near-Earth tail region toward the Earth on the dawn side and away on the dusk side, associated with field-aligned currents of the same signature. While the total cross-tail current is reduced at the location of the near-Earth reconnection site, the local current density at the near-Earth X line actually increases, consistent with the strong plasma sheet thinning in that region. The source regions of the field-aligned currents, indicated by strong peaks of ▽·j∥, are identified as regions of strong plasma flow vorticity near the boundaries of the fast earthward flow earthward of the reconnection site. This flow shear increases preexisting magnetic shear resulting in an increased twist of magnetic flux tubes connected with these regions and thus a generation or increase of the field-aligned currents. We note that the current diversion in our simulation is not caused by a locally enhanced resistivity (although in reality this might be the case) but by the localization of reconnection in a three-dimensional tail configuration.

[1]  Harold P. Furth,et al.  Finite‐Resistivity Instabilities of a Sheet Pinch , 1963 .

[2]  Robert L. McPherron,et al.  Satellite studies of magnetospheric substorms on August 15, 1968. IX - Phenomenological model for substorms. , 1973 .

[3]  T. Potemra,et al.  The amplitude distribution of field-aligned currents at northern high latitudes observed by TRIAD. Interim report , 1975 .

[4]  T. Potemra,et al.  Field‐aligned currents in the dayside cusp observed by Triad , 1976 .

[5]  N. A. Saflekos,et al.  Large-scale characteristics of Birkeland currents , 1979 .

[6]  E. W. Hones,et al.  Three‐dimensional computer modeling of dynamic reconnection in the geomagnetic tail , 1981 .

[7]  R. Walker,et al.  Neutral sheet current interruption and field-aligned current generation by three dimensional driven reconnection , 1983 .

[8]  L. Zanetti,et al.  Large-scale Birkeland currents in the dayside polar region during strongly northward IMF: A new Birkeland current system , 1984 .

[9]  E. W. Hones,et al.  Field‐aligned plasma flow in MHD simulations of magnetotail reconnection and the formation of boundary layers , 1986 .

[10]  T. Ogino A three-dimensional MHD simulation of the interaction of the solar wind with the earth's magnetosphere - The generation of field-aligned currents , 1986 .

[11]  J. Lyon,et al.  The solar wind-magnetosphere-ionosphere current-voltage relationship , 1987 .

[12]  J. Birn Magnetotail equilibrium theory: The general three-dimensional solution , 1987 .

[13]  J. Birn Three-dimensional equilibria for the extended magnetotail and the generation of field-aligned current sheets , 1989 .