A&A 376, 288-291 (2001)
P. L. Israelevich1 - T. I. Gombosi2 - A. I. Ershkovich1 - K. C. Hansen2 - C. P. T. Groth2 - D. L. DeZeeuw2 - K. G. Powell3
1 - Department of Geophysics and Planetary Sciences, Raymond and Beverly Sackler Faculty of Exact Sciences,
Tel Aviv University, Ramat Aviv, 69978 Tel Aviv, Israel
2 - Space Physics Research Laboratory, Department of Atmospheric, Oceanic and Space Sciences, University
of Michigan, Ann Arbor, MI 48109, USA
3 - Department of Aerospace Engineering, University of Michigan, Ann Arbor, MI 48109, USA
Received 18 April 2001 / Accepted 15 June 2001
The existence of the radial component of the electric current flowing toward the Sun is revealed in numerical simulation. The total strength of the radial current is A. The only way to fulfil the electric current continuity is to close the radial electric current by means of field- aligned currents at the polar region of the Sun. Thus, the surface density of the closure current flowing along the solar surface can be estimated as 4 A/m, and the magnetic field produced by this current is T, i.e. several percent of the intrinsic magnetic field of the Sun. This seems to mean that any treatment of the solar magnetic field generation should take into account the heliospheric current circuit as well as the currents flowing inside the Sun.
Key words: solar wind - Sun: magnetic fields
Solar wind outflow from the Sun causes the magnetic field lines originating from opposite hemispheres to extend radially outward and brings them together near the equatorial plane (e.g. Pneuman & Kopp 1971). These field lines are of opposite magnetic polarities, and thus the interplanetary magnetic field must change sign abruptly within the narrow layer near the equatorial plane. This implies the presence of a thin sheet of a very high current density. The current circulates around the dipole axis in the same direction as the original current generating the dipole field. It is this heliospheric current sheet that separates fields and plasma flows from different hemispheres. Rotation of the Sun twists the interplanetary magnetic field lines in such a way that they take the shape of Archimedean spiral (Parker 1958). Another effect of the rotation is warping of the current sheet plane into a structure that resembles a ballerinaís twirling skirt. Because of this warping, the Earth crosses the current sheet at least twice during each solar rotation (sometimes more than twice, if the current sheet is wavy enough). These current sheet crossings are observed at the Earth as crossings of interplanetary sector boundaries, when the polarity of the magnetic field changes its sign (as observed from the Earth) (Wilcox & Ness 1965).
Another consequence of solar rotation and magnetic field line
twisting was predicted by Alfvén (1981). The spiral form of
the magnetic field
lines means that there is a significant radial component of the
electric current in
the sheet, along with the azimuthal component. This point is
|Figure 1: a) Interplanetary magnetic field lines and electric current lines in the solar wind near the heliospheric current sheet surface. Solid lines show the magnetic field lines slightly above the surface, and dashed lines correspond to the field lines slightly below the current sheet surface. Empty lines show the direction of the electric current in the sheet. b) Heliospheric current circuit in the meridional plane. Thick lines show the closure of the radial component of the heliospheric electric current.|
|Open with DEXTER|
Since the publication of Alfvén's conceptual model of the heliospheric current system very little work has been done to follow up on this idea. Recently, however, several efforts are under way to model the global three-dimensional heliosphere. These efforts are motivated by the space weather and Living with a Star initiatives and they, for the first time, attempt the understand the global current system connecting the sun and the heliosphere. The purpose of this note is to show that the three dimensional structure of the heliospheric current system obtained from a self-consistent, first-principles based numerical model of the solar wind outflow with realistic intrinsic solar magnetic field is consistent with Alfvén's conceptual model.
The work published by Groth et al. (2000) included both a model of the solar corona and solar wind during quiet solar conditions as well as a simulation of a space weather event and its interaction with the Earth. Here we concentrate on the steady state solar wind like solution. The solar wind model is generated by solving the governing equations of ideal MHD using the BATS-R-US code. A detailed description of the code can be found in the papers by Powell et al. (1999). The code employs a high resolution finite volume scheme which is based on an approximate Riemann solver for MHD. In this approach, the single fluid MHD equations are solved for in a fully three-dimensional (no dimensional splitting) tightly coupled manner, rather than in separate steps. BATS-R-US makes use of a block based adaptive grid which allows the model to resolve structures of widely varying length scales.
The computational domain used in the calculation was a rectangular box defined by , , and . The adapted computational grid for the solar wind solution consisted of over cells with 8 refinement levels. The minimum cell size at the solar surface was 1/16 , while the cells at 1 AU where approximately 10 in size.
Boundary conditions at the solar surface were dependent on local flow conditions: plasma could freely leave the reservoir, but no "backflow" was allowed. In addition, no deviation was allowed from the intrinsic magnetic field at the surface. Simple zero-gradient boundary conditions are posed at the outer boundaries. The intrinsic solar magnetic field at the solar surface was assumed to be a combination of a dipole and octupole both tilted from the rotation axis. In order to obtain a nominal solar wind-like solution, Groth et al. (2000) included a volumetric heat source beyond, but near the solar surface. The role of the volumetric heating function is to mimic the combined effects of energy absorption above the transition region, heat conduction (which is not included in ideal MHD) and radiative losses. A complete description of the solar wind model as well as the simulation of the space weather event and its interaction with Earth may be found in Groth et al. (2000).
Figure 2a shows isolines of
(a product of the azimuthal component of
magnetic field by the radial distance) in the meridional plane XZ. The thick
corresponds to the current sheet.
|Figure 2: Isolines of the product of the azimuthal component of the magnetic field by the radial distance from the Sun. a) meridional plane and b) equatorial plane.|
|Open with DEXTER|
Figure 2b shows the isolines of the same quantity, in the equatorial plane XY. The intersection of the current sheet with the equatorial plane is given by the condition , and is shown by the thick line. Note that the isolines of do not necessarily coincide with the magnetic field lines. The quantity equals zero at the surface of the current sheet, and grows with the distance from the surface. Therefore -value in the equatorial plane (Fig. 2b) can be considered as a measure of the current sheet distance from the equatorial plane. In other words, it corresponds to the altitude of warps of the current sheet. Thus, the current sheet warping clearly exhibits a spiral structure (the ballerina skirt model).
The distributions of the electric current density in the plane
to the solar magnetic dipole axis (
are given in
|Figure 3: Isolines of products of the components of the electric current density a) and b) by the radial distance from the Sun shown in the plane perpendicular to the solar magnetic dipole axis.|
|Open with DEXTER|
Continuity of the electric current requires the radial component closure
through the solar atmosphere in a way sketched in Fig. 1b. We are able to
follow this closure, and Fig. 4 gives a polar view of isolines of the radial
current on the Sun's surface.
|Figure 4: Polar view of the distribution of the electric current density radial component on the solar surface.|
|Open with DEXTER|
In contrast to simulation result (Fig. 4), the positive radial component of the current, according to Alfvén, is limited to a single point (the pole) (see Fig. 1). Such a picture is an idealization which is not confirmed by our numerical simulation. In fact, the results of calculations (presented in Fig. 4) show that there is a continuous transition from positive to negative value of the radial electric current. All the region within the ring of downward (negative) field-aligned current is filled by the positive (upward) current. Unfortunately, more detailed comparison with the Alfvén prediction is impossible since Alfvén (1981) did not discuss this question.
The numerical simulation has been performed for the main solar magnetic field configuration typical for Solar Minimum conditions. During the next Solar Minimum, the solar magnetic dipole polarity is reversed. Obviously, the radial component of the electric current also reverses its sign (Alfvén 1981). Near the Solar Maximum, the dipole component of the main magnetic field almost vanishes. Domination of higher harmonics makes the main field irregular, and, therefore, the regular structure of the heliospheric current system (Fig. 1) disappers. However, radial electric currents should be present in the solar wind current sheets if the plane of the sheet is not parallel to the Sun's rotation axis. The heliospheric magnetic field behavior during the solar cycle can be described as a magnetic flux transport from the equatorial to polar regions of the Sun. During the 11-year solar cycle, the solar wind transfers the magnetic flux equal to the main solar magnetic flux through a hemisphere (Israelevich et al. 2000).
The three-dimensional heliospheric current system revealed in our simulation may play an important role in the processes on the Sun. The total strength of the radial current is A. Thus the surface density of the closure current flowing along the solar surface can be estimated as A/m, and the magnetic field produced by this current is T, i.e. several percent of the intrinsic magnetic field of the Sun. This seems to mean that any consideration of the solar magnetic field generation should take into account the heliospheric current circuit as well as the currents flowing inside the Sun. Such a conlusion corroborates our recent result (Israelevich et al. 2000) that there is a feedback between the solar wind flow and the main solar magnetic field.
This work was supported by the U.S.-Israel Binational Science Foundation (BSF).